proj4 = '+proj=longlat +ellps=intl +towgs84=-178.4,-83.2,-221.3,0.54,-0.532,-0.126,21.2 +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)