proj4 = '+proj=longlat +ellps=bessel +towgs84=565.4171,50.3319,465.5524,1.9342,-1.6677,9.1019,4.0725 +no_defs +type=crs' m = Map(256,256,proj4) lyr = Layer('Name',proj4)