proj4 = '+proj=utm +zone=18 +ellps=intl +towgs84=221.899,274.136,-397.554,1.361573e-05,-2.174431e-06,-1.36241e-05,-2.199943 +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)