mswkt = '''VERT_CS["NGVD29 height (ftUS)",VERT_DATUM["National Geodetic Vertical Datum 1929",2005,AUTHORITY["EPSG","5102"]],UNIT["US survey foot",0.304800609601219,AUTHORITY["EPSG","9003"]],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","5702"]]''' m = mapObj('') m.setWKTProjection(mswkt) lyr = layerObj(m) lyr.setWKTProjection(mswkt)