proj4 = '+proj=utm +zone=29 +ellps=intl +towgs84=-83.11,-97.38,-117.22,0.0276,-0.2167,0.2147,0.1218 +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)