proj4 = '+proj=utm +zone=29 +ellps=intl +towgs84=-181.5,-90.3,-187.2,0.144,0.492,-0.394,17.57 +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)