harmonica.derivative_northing

harmonica.derivative_northing

harmonica.derivative_northing(grid, order=1, method='finite-diff')[source]

Calculate the derivative of a regular grid in the northing direction

Compute the spatial derivative in the northing direction of regular gridded data. It can compute using accurate central differences using xarray.differentiate or through frequency domain calculations through Fast Fourier Transform.

Important

Choosing the finite differences option produces more accurate results without border effects.

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.

  • order (int) – The order of the derivative. Default to 1.

  • method (str (optional)) – Method that will be used for computing the easting derivative. It can be either "finite-diff", for computing using xarray.differentiate, or "fft", for using FFT-based filters. Default "finite-diff".

Returns

derivative (xarray.DataArray) – A xarray.DataArray with the northing derivatives of the passed grid. Its units are the same units of the grid per units of its coordinates to the power of the passed order.

References

[Blakely1995]