harmonica.total_gradient_amplitude

harmonica.total_gradient_amplitude(grid)

Calculates the total gradient amplitude of a potential field grid

Compute the total gradient amplitude of a regular gridded potential field M. The horizontal derivatives are calculated though 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:

total_gradient_amplitude_gridxarray.DataArray

A xarray.DataArray after calculating the total gradient amplitude of the passed grid.

Notes

The total gradient amplitude is calculated as:

\[A(x, y) = \sqrt{ \left( \frac{\partial M}{\partial x} \right)^2 + \left( \frac{\partial M}{\partial y} \right)^2 + \left( \frac{\partial M}{\partial z} \right)^2 }\]

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

References

[Blakely1995]

Examples using harmonica.total_gradient_amplitude

Total gradient amplitude of a regular grid

Total gradient amplitude of a regular grid