proj4 = '+proj=longlat +ellps=intl +towgs84=-96.062,-82.428,-121.753,-4.801,-0.345,1.376,1.496 +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)