proj4 = '+proj=longlat +ellps=intl +towgs84=-81.1,-89.4,-115.8,0.485,0.024,0.413,-0.54 +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)