proj4 = '+proj=lcc +lat_0=31.1666666666667 +lon_0=-100 +lat_1=27.4166666666667 +lat_2=34.9166666666667 +x_0=1000000 +y_0=1000000 +ellps=GRS80 +nadgrids=AB_CSRS.DAC +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)