gsolve.reductions.terrain_corrections.tcorr_harmonica_topography

gsolve.reductions.terrain_corrections.tcorr_harmonica_topography(point, topography, topography_density, parallel=False, disable_checks=False)

Calculate topographic terrain correction using harmonica.

This is called by calculate_terrain_correction() to compute terrain corrections for the topography portions of a DEM. You should not need to call this function directly, however it is exposed for advanced use cases.

Parameters:

pointtuple of floats (x, y, z)

The (x, y, z) coordinates of the observation point.

topographyxarray.DataArray

The topographic surface elevations in meters. The DataArray should have no values below ‘sea level’

topography_densityxarray.DataArray

The density of the topographic surface in kg/m^3, with the same coordinates as topography. Must be compatible with topography.

parallelbool, optional

Whether to use parallel processing, by default False.

disable_checksbool, optional

Whether to disable input checks, by default False

Returns:

float

The calculated topographic terrain correction at point

See also

harmonica.prism_layer.gravity

the underkying harmonica method used to compute the terrain correction.