harmonica.reduction_to_pole

harmonica.reduction_to_pole

harmonica.reduction_to_pole(grid, inclination, declination, magnetization_inclination=None, magnetization_declination=None)[source]

Calculate the reduction to the pole of a magnetic field grid

Compute the reduction to the pole of regular gridded magnetic data using frequency domain calculations through Fast Fourier Transform.

Parameters
  • grid (xarray.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.

  • inclination (float in degrees) – The inclination of the inducing Geomagnetic field.

  • declination (float in degrees) – The declination of the inducing Geomagnetic field.

  • magnetization_inclination (float in degrees or None) – The inclination of the total magnetization of the anomaly source. If None, the magnetization_inclination will be set equal to the inclination, neglecting remanent magnetization and self demagnetization. Default None.

  • magnetization_declination (float in degrees) – The declination of the total magnetization of the anomaly source. If None, the magnetization_declination will be set equal to the declination, neglecting remanent magnetization and self demagnetization. Default None.

Returns

reduced_to_pole_grid (xarray.DataArray) – A xarray.DataArray after reduction to the pole of the passed grid.

References

[Blakely1995]

Examples using harmonica.reduction_to_pole