When trying to tile a LiDAR dataset that has a compound CS, a geoid-height Z adjustment is applied even without the flag being set. This dataset (from UND in Grand Forks, ND) is adjusted by -26.672m.
$ tilepoints -f test_aaa.laz & tail -f ~/viz-points/pdgpoints.log
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: File: /mnt/ceph/home/nesbitt/lidar/test_aaa.laz
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: Merge: False
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: Intensity > RGB: False
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: Intens. scalar: 1.0x
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: Translate Z: +0.0
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: From geoid: None
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: Archive input: False
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: LLVS CSV input: False
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: LLV color scale: None
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: Given name: test_aaa
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.log_init_stats: File extension: .laz
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.run: Creating dir /mnt/ceph/home/nesbitt/lidar/rewrite
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.run: Creating dir /mnt/ceph/home/nesbitt/lidar/archive
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.run: Creating dir /mnt/ceph/home/nesbitt/lidar/3dtiles
2025-03-28T00:32:58:INFO:pdgpoints.pipeline.run: Rewriting file with new OGC WKT... (step 1 of 4)
2025-03-28T00:32:58:INFO:pdgpoints.lastools_iface.log_subprocess_output: subprocess output: "WARNING: no projection information. ignoring '-set_ogc_wkt' for '/mnt/ceph/home/nesbitt/lidar/test_aaa.laz'"
2025-03-28T00:34:45:INFO:pdgpoints.lastools_iface.las2las_ogc_wkt: Finished las2las (107 sec / 1.8 min)
2025-03-28T00:34:45:INFO:pdgpoints.pipeline.run: Doing lasinfo dump... (step 2 of 4)
2025-03-28T00:34:45:INFO:pdgpoints.utils.get_epsgs_from_wkt: Found compound coordinate system (COMPD_CS): 2 entries
2025-03-28T00:34:45:INFO:pdgpoints.utils.get_epsgs_from_wkt: Found horizontal EPSG: 32614 (WGS 84 / UTM zone 14N)
2025-03-28T00:34:45:INFO:pdgpoints.utils.get_epsgs_from_wkt: Found vertical EPSG: 3855 (EGM2008 height)
2025-03-28T00:34:45:INFO:pdgpoints.lastools_iface.lasinfo: Compound CRS: WGS 84 / UTM zone 14N
2025-03-28T00:34:45:INFO:pdgpoints.lastools_iface.lasinfo: Compound VRS: EGM2008 height
2025-03-28T00:34:45:INFO:pdgpoints.lastools_iface.lasinfo: Writing WKT to /mnt/ceph/home/nesbitt/lidar/test_aaa-wkt.laz-wkt.txt
2025-03-28T00:34:45:INFO:pdgpoints.lastools_iface.lasinfo: Finished lasinfo (0 sec / 0.0 min)
2025-03-28T00:34:45:INFO:pdgpoints.pipeline.run: self.from_geoid="None", las_vrs="3855"
2025-03-28T00:34:45:INFO:pdgpoints.pipeline.run: Getting mean lat/lon from las file... (step 3 of 4)
2025-03-28T00:34:45:INFO:pdgpoints.lastools_iface.lasmean: Writing abridged XY file to /mnt/ceph/home/nesbitt/lidar/test_aaa-wkt.laz-xy.txt
2025-03-28T00:35:54:INFO:pdgpoints.lastools_iface.lasmean: X mean: 642672.533 Y mean: 5309301.630 (WGS 84 / UTM zone 14N)
2025-03-28T00:35:54:INFO:pdgpoints.lastools_iface.lasmean: Finished las2las (68 sec / 1.1 min)
2025-03-28T00:35:54:INFO:pdgpoints.pipeline.run: Resolving geoid/tidal model... (step 3 of 4)
2025-03-28T00:35:54:INFO:pdgpoints.pipeline.run: Looking up ellipsoid height of EGM2008 at (47.921, -97.090)... (step 3 of 4)
2025-03-28T00:35:55:INFO:pdgpoints.pipeline.run: Manual Z transformation: 0.000
2025-03-28T00:35:55:INFO:pdgpoints.pipeline.run: Geoid height adjustment: -26.672
2025-03-28T00:35:55:INFO:pdgpoints.pipeline.run: Translating Z values by -26.672
2025-03-28T00:35:55:INFO:pdgpoints.pipeline.run: Starting las2las rewrite... (step 4 of 4)
2025-03-28T00:35:55:INFO:pdgpoints.lastools_iface.las2las: Rewriting LAS to avoid VLR size errors (e.g. PDAL reading QTModeler files)
2025-03-28T00:37:45:INFO:pdgpoints.lastools_iface.las2las: Finished las2las (110 sec / 1.8 min)
2025-03-28T00:37:45:INFO:pdgpoints.pipeline.run: Starting tiling process... (step 5 of 4)
2025-03-28T00:37:45:INFO:pdgpoints.py3dtiles_iface.tile: File: /mnt/ceph/home/nesbitt/lidar/rewrite/test_aaa.laz
2025-03-28T00:37:45:INFO:pdgpoints.py3dtiles_iface.tile: Creating tile directory
2025-03-28T00:37:45:INFO:pdgpoints.py3dtiles_iface.tile: CRS to convert from: EPSG:32614
2025-03-28T00:37:45:INFO:pdgpoints.py3dtiles_iface.tile: CRS to convert to: EPSG:4978
2025-03-28T00:40:30:INFO:pdgpoints.py3dtiles_iface.tile: Finished tiling (165 sec / 2.8 min)
2025-03-28T00:40:30:INFO:pdgpoints.pipeline.run: Cleaning up processing artifacts.
2025-03-28T00:40:30:INFO:pdgpoints.pipeline.run: Finished processing test_aaa.laz (451 sec / 7.5 min)
When trying to tile a LiDAR dataset that has a compound CS, a geoid-height Z adjustment is applied even without the flag being set. This dataset (from UND in Grand Forks, ND) is adjusted by -26.672m.