proj4 = '+proj=utm +zone=29 +ellps=intl +towgs84=-157.89,-17.16,-78.41,2.118,2.697,-1.434,-5.38 +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)