I have a specific problem re: projected coordinate systems. I am working with point data from IceSat2 and a LiDAR-derived DEM of Daviess county in Indiana. I am trying to compare height values for points (in shapefile format).
The DEM coordinate system comes from the data portal with reference system:
EPSG Code 2966, NAD_1983_StatePlane_Indiana_West_FIPS_1302_Feet Units of Measure: us_survey_feet Horizontal Datum: NAD83 Spheroid: GRS 1980 Vertial: GEOID NAVD1988 (elevation units of us_survey_feet)
The point data was originally taken from ICESAT2 data files which are said to have a WGS84 reference
This data set (ATL03) contains height above the WGS 84 ellipsoid (ITRF2014 reference frame)
When I write out the point shapefile, I use EPSG:4326. I then re-projected the DEM from EPSG:2966 to EPSG:4326.
When I finally sample the raster with the point shapefile, the heights are similar but offset by some ~34 meters...
I understand it's likely something simple with a disconnect between the spatial reference system and the projected coordinate system or something, but am at a loss. Please let me know a way forward!