corrected_geoid_undulation

  • Calculates the topographically corrected geoidal undulation at a given latitude and longitude using an iterative approach

Calling Sequence

from geoid_toolkit.corrected_geoid_undulation import corrected_geoid_undulation
N = corrected_geoid_undulation(lat, lon, 'WGS84', clm, slm, tclm, tslm, lmax, R, GM)

Source code

geoid_toolkit.corrected_geoid_undulation(lat, lon, refell, clm, slm, tclm, tslm, lmax, R, GM, density, GAUSS=0, EPS=1e-08)[source]

Calculates the topographically corrected geoidal undulation using the iterative approach described in Barthelmes [1], Moazezi and Zomorrodian [8]

Parameters:
lat: float

latitude in degrees

lon: float

longitude in degrees

refell: str

Reference ellipsoid name

  • 'CLK66': Clarke 1866

  • 'GRS67': Geodetic Reference System 1967

  • 'GRS80': Geodetic Reference System 1980

  • 'HGH80': Hughes 1980 Ellipsoid

  • 'WGS72': World Geodetic System 1972

  • 'WGS84': World Geodetic System 1984

  • 'ATS77': Quasi-earth centred ellipsoid for ATS77

  • 'NAD27': North American Datum 1927

  • 'NAD83': North American Datum 1983

  • 'INTER': International

  • 'KRASS': Krassovsky (USSR)

  • 'MAIRY': Modified Airy (Ireland 1965/1975)

  • 'TOPEX': TOPEX/POSEIDON ellipsoid

  • 'EGM96': EGM 1996 gravity model

clm: float

cosine spherical harmonics for a gravity model

slm: float

sine spherical harmonics for a gravity model

tclm: float

cosine spherical harmonics for a topographic model

tslm: float

sine spherical harmonics for a topographic model

lmax: int

maximum spherical harmonic degree

R: float

average radius used in gravity model

GM: float

geocentric gravitational constant used in gravity model

GAUSS: float, default 0

Gaussian Smoothing Radius in km

EPS: float, default 1e-8

level of precision for calculating geoid height

Returns:
N: float

geoidal undulation for a given ellipsoid in meters