harmonica.filters.derivative_easting_kernel

harmonica.filters.derivative_easting_kernel(fft_grid, order=1)

Filter for easting derivative in frequency domain

Return a xarray.DataArray with the values of the frequency domain filter for computing the easting derivative. The filter is built upon the frequency coordinates of the passed fft_grid and is defined as follows:

\[g(\mathbf{k}) = (i k_e)^n\]

where \(\mathbf{k}\) is the wavenumber vector (\(\mathbf{k} = 2\pi \mathbf{f}\) where \(\mathbf{f}\) is the frequency vector), \(k_e\) is the easting wavenumber component of \(\mathbf{k}\), \(i\) is the imaginary unit and \(n\) is the order of the derivative.

Parameters:

fft_gridxarray.DataArray

Array with the Fourier transform of the original grid. Its dimensions should be in the following order: freq_northing, freq_easting. Use xrft.xrft.fft and xrft.xrft.ifft functions to compute the Fourier Transform and its inverse, respectively.

orderint

The order of the derivative. Default to 1.

Returns:

da_filterxarray.DataArray

Array with the kernel for the easting derivative filter in frequency domain.

See also

harmonica.derivative_easting

References

[Blakely1995]