proj4 = '+proj=lcc +lat_0=44 +lon_0=-70 +lat_1=50 +lat_2=46 +x_0=800000 +y_0=0 +ellps=clrk66 +nadgrids=NTv2_0.gsb +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)