invert4geom.DatasetAccessorInvert4Geom.regional_trend

invert4geom.DatasetAccessorInvert4Geom.regional_trend#

DatasetAccessorInvert4Geom.regional_trend(trend, regional_shift=0, mask_column=None, reverse_regional_residual=False)[source]#

Calculate the gravity misfit as the difference between dataset variables gravity_anomaly and forward_gravity. Then separate the misfit into regional and residual components where residual = misfit - regional.

Approximate the regional field by fitting a polynomial trend to the gravity misfit using verde.Trend.

The resulting regional field can be shifted with regional_shift, and the calculated residual field can be multiplied by the values in mask_column.

The following new variables are added to the dataset:

misfit, reg, res, starting_forward_gravity, starting_misfit, starting_reg, starting_res

Parameters:
  • trend (int) – order of the polynomial trend to fit to the data

  • regional_shift (float) – shift to add to the regional field, by default 0

  • mask_column (str | None) – Name of optional dataset variable with values to multiply the calculated residual gravity field by, should have values of 1 or 0, by default None.

  • reverse_regional_residual (bool) – if True, reverse the regional and residual fields after calculation, by default False

Return type:

None