from mapscript import mapObj,layerObj wkt = '''PROJCS["GRS80_127",GEOGCS["GCS_ITRF_2000",DATUM["D_ITRF_2000 (EPSG ID 9028080)",SPHEROID["GRS_1980 (EPSG ID 9021980)",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["degree",0.017453292519943295],AXIS["Longitude",EAST],AXIS["Latitude",NORTH]],PROJECTION["Transverse_Mercator"],PARAMETER["central_meridian",127.0],PARAMETER["latitude_of_origin",38.0],PARAMETER["scale_factor",1.0],PARAMETER["false_easting",200000.0],PARAMETER["false_northing",600000.0],UNIT["m",1.0],AXIS["x",EAST],AXIS["y",NORTH]]''' m = mapObj('') m.setWKTProjection(wkt) lyr = layerObj(m) lyr.setWKTProjection(wkt)