mswkt = '''VERT_CS["Chatham Island 1959 height",VERT_DATUM["Waitangi (Chatham Island) 1959",2005,AUTHORITY["EPSG","5169"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Gravity-related height",UP],AUTHORITY["EPSG","5771"]]''' m = mapObj('') m.setWKTProjection(mswkt) lyr = layerObj(m) lyr.setWKTProjection(mswkt)