Hi, I'm hoping you can help me with the vertical datums feature, which will be a huge help. I have a project referenced in WGS84 Lat/Lng and meters HAE, and would like to export it ideally in a state plane with NAVD88 feet.
I see that when I use the 'Convert' button in the reference pane, and pick the WGS84 + EGM96 geoid height, the z values all change to reasonable geoid heights in meters. The value of "chunk.crs.wkt" changes to be a COMPD_CS with VERT_CS/VERT_DATUM components that make sense. However, I can't seem to build my own PhotoScan.CoordinateSystem object with a vertical datum component. Even trying to create a new CoordinateSystem with the WKT that was chosen automatically fails (with a non-pythonish error in the binding):
>>> PhotoScan.CoordinateSystem(PhotoScan.app.document.chunk.crs.wkt)
Traceback (most recent call last):
File "<console>", line 1, in <module>
SystemError: NULL result without error in PyObject_Call
What I would ideally like to do is build a CoordinateSystem with something like the following, to get state plane data in US Survey Foot units in all three axes, and use it with chunk.exportDem.
COMPD_CS["California zone 4 NAVD88 ftUS",
PROJCS["NAD83(NSRS2007) / California zone 4 (ftUS)",
GEOGCS["NAD83(NSRS2007)",
DATUM["NAD83_National_Spatial_Reference_System_2007",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6759"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4759"]],
UNIT["US survey foot",0.3048006096012192,
AUTHORITY["EPSG","9003"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",37.25],
PARAMETER["standard_parallel_2",36],
PARAMETER["latitude_of_origin",35.33333333333334],
PARAMETER["central_meridian",-119],
PARAMETER["false_easting",6561666.667],
PARAMETER["false_northing",1640416.667],
AUTHORITY["EPSG","3496"],
AXIS["X",EAST],
AXIS["Y",NORTH]],
VERT_CS["EGM96 geoid height",
VERT_DATUM["EGM96 geoid",2005,
AUTHORITY["EPSG","5171"]],
UNIT["US survey foot",0.3048006096012192,
AUTHORITY["EPSG","9003"]],
AUTHORITY["EPSG","5773"]]
]
Thanks,
John