proj4 = '+proj=longlat +ellps=aust_SA +nadgrids=SAD69_003.gsb +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)