Defining vertical datum when exporting point clouds

[ivan-zd]

Hi,

When exporting the point cloud using exportPoints() function I can define only the output CRS but not the vertical datum as well. Is there any possibility to define the vertical datum (using EPSG ideally) so that it's also written to the header of the LAS/LAZ file?

[Paulo]

Hi Ivan-ZD,

you just need to define the CRS that includes your geoid. Supposing your chunk crs is already defined with geoid then you coud just do:

Code: [Select]

crs = chunk.crs
chunk.exportPoints(...,crs=crs,...)
Hope this gets you along the way....

[ivan-zd]

Hi Paulo,

Thank you for your reply. I will try once with geoid file but I'm afraid it won't be good approach for me on the long run as I was hoping for rather generic approach where I could just define EPSG code and that it's then written to the header of the LAZ file. The reason behind is that I need to have LAZ file with this vertical datum info when I pass it over to Cesium ION for tiling because otherwise Cesium assume that heights are relative to ellipsoid height and it's not fitting the globe.

I see PDAL can help with setting correctly vertical datum info but I was hoping I could do that directly in Agisoft to avoid additional steps.

Anyway, I'll look more in details there, thanks again for your answer.


Best regards,
Ivan

[Alexey Pasumansky]

Hello Ivan,

If you need to export the dense cloud in the compound coordinate system that uses height above geoid, the corresponding coordinate system definition should be passed to the export task (can be loaded from external PRJ file or as a string variable) and also the corresponding geoid file should be initialized in the code using addGeoid function before any operations for the related coordinate systems are applied. Instead of geoid file, however, you can use fixed offset defined via crs.geoid_height method.