Elevation transforms and HAE vs MSL #1987
Replies: 2 comments 1 reply
-
Do me a favor - this is too much to digest and parse as-is. Please give us one example of one |
Beta Was this translation helpful? Give feedback.
-
Here is a more concise example:
Gives the output (with input LL 40.5745329750, -111.7853343571): [osgEarth] !^!^!^! mapPoint elev(m): 1597.71 Basically, I would expect these elevation values to differ. The goal is to get MSL and HAE from the same elevation pool sample by using GeoPoint::transform with different SRS/vdatums. |
Beta Was this translation helpful? Give feedback.
-
I'm trying to figure out how best to consistently get HAE and MSL. I wrote some test code from some previous forum discussions and I'm not getting the results I anticipated. Can you guide me on how to get MSL and HAE elevations and, if possible, how to get the layer from which the elevation value was retrieved from?
Note that GeoPoint::transform doesn't seem to transform the elevation from one vdatum to another. I'm also confused on why the last two VerticalDatum::transform calls always come back with the same results.
Any insight on what I'm doing wrong would be greatly appreciated!
Here is the sample code, the outputs, and the .earth file I am using:
Inputs:
const std::vector< std::pair<double, double> > LAT_LON_LIST_UTAH
{
{40.5745329750, -111.7853343571 }, //LCC
{40.5726146663, -111.7765861811 }, //Little Cottonwood Road fork
{40.5888667984, -111.6388564531 }, //Alta
{40.5817361188, -111.6569948744 }, //Snowbird
{40.3761752391, -105.5236120331 }, //Estes Park
{45.4108516473, -75.6875127888 }, //Ottowa
{38.9315317322, -119.9790876095} //South Lake Tahoe
};
Test code:
Output:
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 1614.74
[osgEarth] !^!^!^! mapPointEGM96(m): 1614.74 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 1614.74 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 1614.74 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 1614.74 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 1631.77
[osgEarth] !^!^!^! EGM84: 1629.87
[osgEarth] !^!^!^! EGM08: 1614.74
[osgEarth] !^!^!^! GEOID: 1614.74
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 1649.66
[osgEarth] !^!^!^! mapPointEGM96(m): 1649.66 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 1649.66 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 1649.66 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 1649.66 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 1666.68
[osgEarth] !^!^!^! EGM84: 1664.78
[osgEarth] !^!^!^! EGM08: 1649.66
[osgEarth] !^!^!^! GEOID: 1649.66
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 2609.62
[osgEarth] !^!^!^! mapPointEGM96(m): 2609.62 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 2609.62 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 2609.62 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 2609.62 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 2626.31
[osgEarth] !^!^!^! EGM84: 2624.47
[osgEarth] !^!^!^! EGM08: 2609.62
[osgEarth] !^!^!^! GEOID: 2609.62
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 2484.32
[osgEarth] !^!^!^! mapPointEGM96(m): 2484.32 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 2484.32 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 2484.32 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 2484.32 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 2501.07
[osgEarth] !^!^!^! EGM84: 2499.21
[osgEarth] !^!^!^! EGM08: 2484.32
[osgEarth] !^!^!^! GEOID: 2484.32
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 2304.33
[osgEarth] !^!^!^! mapPointEGM96(m): 2304.33 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 2304.33 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 2304.33 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 2304.33 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 2317.47
[osgEarth] !^!^!^! EGM84: 2316.54
[osgEarth] !^!^!^! EGM08: 2304.33
[osgEarth] !^!^!^! GEOID: 2304.33
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 66.3461
[osgEarth] !^!^!^! mapPointEGM96(m): 66.3461 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 66.3461 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 66.3461 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 66.3461 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 99.9251
[osgEarth] !^!^!^! EGM84: 99.1827
[osgEarth] !^!^!^! EGM08: 66.3461
[osgEarth] !^!^!^! GEOID: 66.3461
[osgEarth] geodetic (wgs 84)
[osgEarth] !^!^!^! mapPoint elev(m): 1905.18
[osgEarth] !^!^!^! mapPointEGM96(m): 1905.18 (WGS84:EGM96-wgs 84)
[osgEarth] !^!^!^! mapPointEGM84(m): 1905.18 (WGS84:EGM84-wgs 84)
[osgEarth] !^!^!^! mapPointHAE (m): 1905.18 (WGS84:geodetic-wgs 84)
[osgEarth] !^!^!^! mapPointEGM96(m): 1905.18 (WGS84:EGM2008-wgs 84)
[osgEarth] !^!^!^! EGM96: 1928.94
[osgEarth] !^!^!^! EGM84: 1928.02
[osgEarth] !^!^!^! EGM08: 1905.18
[osgEarth] !^!^!^! GEOID: 1905.18
Earth file:
Beta Was this translation helpful? Give feedback.
All reactions