Atmosphere correction (pygeoid.reduction.atmosphere)¶
Calculate atmospheric correction for the gravity anomalies.
-
pygeoid.reduction.atmosphere.ussa76_density(alt_arr: Unit("km") = <Quantity 0. km>) -> Unit("kg / m3")[source]¶ Return atmospheric density from USSA76 model.
- Refer to the following document (2 doc codes) for details of this model:
- NOAA-S/T 76-1562 NASA-TM-X-74335
All assumptions are the same as those in the source documents.
Derived from: https://github.com/mattljc/atmosphere-py
Parameters: alt_arr (Quantity) – Altitude above sea level.
-
pygeoid.reduction.atmosphere.iag_atm_corr_sph(density_function: Callable[[astropy.units.quantity.Quantity], astropy.units.quantity.Quantity], height: Unit("m"), height_max: Unit("m"), samples=10000.0) -> Unit("mGal")[source]¶ Return atmospheric correction to the gravity anomalies by IAG approach.
This function numerically integrates samples from density function by trapezoidal rule. The spherical layering of the atmosphere is considered.
IAG approach:
g_atm = G*M(r) / r**2
inf /- M(r) = 4*pi*| rho(r) * r**2 dr
- / r
Parameters: - density_function (callable) – The density_funtion is called for all height samples to calculate density of the atmosphere.
- height (Quantity) – Height above sea level.
- height_max (Quantity) – Maximum height of the atmosphere layer above sea level.
- samples (float) – Number of samples for integration. Default is 1e4.
- Quantity – Atmospheric correction.
-
pygeoid.reduction.atmosphere.grs80_atm_corr_interp(height: Unit("m"), kind: str = 'linear') -> Unit("mGal")[source]¶ Return GRS 80 atmospheric correction, in mGal.
Interpolated from the table data [1]_.
Note: If height < 0 m or height > 40000 m, then correction is extrapolated
Parameters: - height (Quantity) – Height above sea level.
- kind (str or int, optional) – Specifies the kind of interpolation as a string (‘linear’, ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, ‘cubic’ where ‘zero’, ‘slinear’, ‘quadratic’ and ‘cubic’ refer to a spline interpolation of zeroth, first, second or third order) or as an integer specifying the order of the spline interpolator to use. Default is ‘linear’.
Returns: Atmospheric correction.
Return type: Quantity
References
[1] Moritz, H. (1980). Geodetic reference system 1980. Bulletin Géodésique, 54(3), 395-405
-
pygeoid.reduction.atmosphere.wenzel_atm_corr(height: Unit("m")) -> Unit("mGal")[source]¶ Return atmospheric correction by Wenzel, in mGal.
Parameters: height (Quantity) – Height above sea level. Returns: Atmospheric correction. Return type: Quantity References
[1] Wenzel, H., 1985, Hochauflosende Kugelfunktionsmodelle fur des Gravitationspotential der Erde [1]: Wissenschaftliche arbeiten der Fachrichtung Vermessungswesen der Universitat Hannover, 137