Files

15 lines
384 B
C++

#include "cgeoid.h"
#include "GeographicLib/Geoid.hpp"
extern "C"
double HeightAboveEllipsoid(double lat, double lon, double h) {
try {
// Declare static so that g is only constructed once
static const GeographicLib::Geoid g("egm2008-1");
return h + GeographicLib::Geoid::GEOIDTOELLIPSOID * g(lat, lon);
}
catch (...) {
return GeographicLib::Math::NaN();
}
}