harmonica.upward_continuation#

harmonica.upward_continuation(grid, height_displacement, *, pad=True, pad_kwargs=None)[source]#

Calculate the upward continuation of a potential field grid.

Compute the upward continuation of regular gridded data using frequency domain calculations through Fast Fourier Transform.

Note

Any non-dimensional coordinates of the grid will be dropped since upward continuation may have made them no longer correct.

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.
height_displacementfloat: The height displacement of upward continuation. For upward continuation, the height displacement should be positive. Its units are the same units of the grid coordinates.
padbool, optional: If True, will add padding to the grid before taking the Fourier Transform and applying the filter and remove it after the inverse Fourier Transform. Adding padding usually helps reduce edge effects from signal truncation. Default is True.
pad_kwargsdict or None, optional: Any additional keyword arguments that should be passed to the xarray.DataArray.pad function in the form of a dictionary. If none are given, the default padding of 25% the dimensions of the grid will be added using the “edge” method.

Returns:

upward continuationxarray.DataArray: A xarray.DataArray after upward continuation of the passed grid.

Examples using `harmonica.upward_continuation`#