from mapscript import mapObj,layerObj wkt = '''PROJCS["Lambert_Conformal_Conic",GEOGCS["Unknown",DATUM["NTF_Paris_Meridian",SPHEROID["Clarke_1880_modified_for_IGN",6378249.2,293.4660213]],PRIMEM["non-Greenwich",2.337229166667],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["standard_parallel_1",43.199291275544],PARAMETER["standard_parallel_2",44.996093814511],PARAMETER["latitude_of_origin",44.1],PARAMETER["central_meridian",0],PARAMETER["false_easting",600000],PARAMETER["false_northing",200000],UNIT["Meter",1.0]]''' m = mapObj('') m.setWKTProjection(wkt) lyr = layerObj(m) lyr.setWKTProjection(wkt)