harmonica.tilt_angle

harmonica.tilt_angle(grid)

Calculates the tilt angle of a potential field grid

Compute the tilt of a regular gridded potential field \(M\). The horizontal derivatives are calculated through finite-differences while the upward derivative is calculated using FFT.

Parameters:

gridxarray.DataArray

A two dimensional xarray.DataArray whose coordinates are evenly spaced (regular grid). Its dimensions should be in the following order: northing, easting. Its coordinates should be defined in the same units.

Returns:

tilt_gridxarray.DataArray

A xarray.DataArray with the calculated tilt in radians.

Notes

The tilt is calculated as:

\[\text{tilt}(f) = \tan^{-1} \left( \frac{ \frac{\partial M}{\partial z} }{ \sqrt{ \left( \frac{\partial M}{\partial x} \right)^2 + \left( \frac{\partial M}{\partial y} \right)^2 } } \right)\]

where \(M\) is the regularly gridded potential field.

References

[Blakely1995] [MillerSingh1994]

Examples using harmonica.tilt_angle

Tilt of a regular grid

Tilt of a regular grid