proj4 = '+proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +ellps=intl +nadgrids=nzgd2kgrid0005.gsb +units=m +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)