from mapscript import mapObj,layerObj wkt = '''PROJCS["Lambert Maroc Zone I (Nord)",GEOGCS["NTF (Paris meridian)",DATUM["NTF (Paris meridian)",SPHEROID["Clarke 1880 (IGN)",6378249.20000000,293.466021293627]],PRIMEM["Paris",0.0],UNIT["Decimal Degree",0.01745329251994330]],PROJECTION["Lambert Conformal Conic"],PARAMETER["Standard_Parallel_1",31.7278664],PARAMETER["Standard_Parallel_2",34.8717272],PARAMETER["Latitude_Of_Origin",33.3000000],PARAMETER["Central_Meridian",-3.062770833],PARAMETER["False_Easting",500000.00],PARAMETER["False_Northing",300000.00],UNIT["Meter",1.000000000000]]''' m = mapObj('') m.setWKTProjection(wkt) lyr = layerObj(m) lyr.setWKTProjection(wkt)