topic Reprojecting LAS datasets using Extract LAS Tool not transforming Z values in LIDAR Resources Questions

Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

James_Whitacre — Wed, 23 Sep 2020 15:32:32 GMT

I am trying to reproject LAS datasets derived from Pix4D photogrammetry output of a UAV flight that were created with the WGS_1984_UTM_Zone_17N / VCS:WGS_1984 (Z unit is meters above ellipsoid) spatial reference so that it can align with other third-party LAS dataset (traditional lidar; not photogrammetry) in the NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet / VCS:NAVD88_height_ (Z unit is Ft_US) spatial reference.

The Extract LAS (3D Analyst)—ArcGIS Pro | Documentation, it explains that the LAS data can be reprojected using the Output Coordinate System environment and it is claimed:

'If the spatial reference contains a vertical coordinate system, the LAS files can also be reprojected to another height reference system if the vertical datum transformation grids are installed and a transformation is available from the source reference and the target coordinate system.'

I have the Coordinate System for ArcGIS grids installed and the appropriate transformation is set in the Map/Scene Properties and everything is projected on-the-fly brilliantly in the map and scenes when viewed. The highlighted transformation below is what is selected in the Map/Scene Properties > Transformation for the above coordinate systems.

However, I cannot set the Geographic Transformation environment setting to match the above compound transformation, so I left it blank hoping the software would choose the appropriate XY and Z transformations. However, the Geographic Transformations (Environment setting)—ArcGIS Pro | Documentation says that 'compound transformations' are not valid in the environment setting (which makes me think that this particular transformation is not possible in this tool).

So, when I run the Extract LAS tool, the vertical location of the points are not transformed/reprojected (while the horizontally they are), and instead appear to be located in the original Z coordinate system, which in this case is below the expected vertical location.

So is the claim that the Z coordinate system for LAS data can be reprojected false? Am I missing something? Any guidance would be appreciated. This will likely be mission critical for me. Thanks.

Re: Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

amann — Fri, 13 Nov 2020 15:58:18 GMT

This can be tricky within the Esri software environment. Are you able to use LAStools?

This will do the job without the headache. Depending on your preferences, the CLI can be easier to use for setting VCS if they do not show up in the GUI.

Re: Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

MarcSwartz — Sat, 19 Dec 2020 18:07:55 GMT

Please do provide a CLI example using (LASTools) las2las.exe going from WGS_1984_UTM_Zone_17N / VCS:WGS_1984 to NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet / VCS:NAVD88_height_ (Z unit is Ft_US). Extra credit if you include Geoid 18

Re: Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

amann — Mon, 21 Dec 2020 02:10:28 GMT

I took a shot at your request and was able to produce decent results using the following command:

D:\LAStools\bin>las2las -i ..\PAUTMv1\PA_UTMtoNAD83.las

-vertical_navd88_geoid12b

-set_ogc_wkt -target_sp83 PA_S

-target_feet -target_elevation_feet -target_precision 0.001

-odix _PA_S_foot -olas

Run LASInfo to see WKT OGC CS details.

I will take a look later this week to work in NAD 83 FIPS instead of NAD 83 and Geoid18 instead of Geoid12b. This should get you started.

As far as I can tell, you can achieve the same results (except for baking the Geoid in the header) using the "Extract LAS" tool in ArcGIS Pro and set your CS/transformation in Environments.

Re: Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

MarcSwartz — Mon, 21 Dec 2020 09:50:21 GMT

Great job. Now others can benefit from having an example available.

Re: Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

JimZack — Fri, 13 Sep 2024 22:25:04 GMT

I know it's probably four years too late, but here's a solution/workaround:

Problem: Reprojecting the Vertical Coordinate System of Point Cloud Data in an LAS Dataset (esri.com)

Re: Reprojecting LAS datasets using Extract LAS Tool not transforming Z values

tikola — Wed, 04 Jun 2025 08:29:01 GMT

And even easier in 3.5: Project LAS (Data Management)—ArcGIS Pro | Documentation

In nutshell that SLPK trick has worked well a while but it has bit uneccessary steps as we just use some Esri geoprocessing functions just as they come from factory. Next step uses simply previous step output as an input. I thought that this might be the case of model builder to queue functions in the row, but even better as 3.5 contains factory made function for this task.

So this has been doable in most recent ArcGIS Pro version and the latest 3.5 has a single tool for this task. My case is typically a drone based LAS in ellipsoidal heights that needs to be converted to local correct orthometric height system. Just yesterday did one in old way and learned after that there is this function now.