proj4 = '+proj=longlat +ellps=intl +towgs84=-116.641,-56.931,-110.559,0.893,0.921,-0.917,-3.52 +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)