from mapscript import mapObj,layerObj wkt = '''PROJCS["UTM Zone 31, Northern Hemisphere",GEOGCS["Undefined",DATUM["Undefined",SPHEROID["Clarke 1866",6378206.4,294.9786982139006]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",3],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]''' m = mapObj('') m.setWKTProjection(wkt) lyr = layerObj(m) lyr.setWKTProjection(wkt)