oemAddErrorPolyFit

Workspace.oemAddErrorPolyFit(self, jac_targets: JacobianTargets = self.jac_targets, oem: OptimalEstimationData = self.oem, measurement_sensor: ArrayOfSensorObsel = self.measurement_sensor, t: Vector, sensor_elem: Index, polyorder: Index = 0, matrix: BlockMatrix, inverse: BlockMatrix = []) → None

Set a measurement error to polynomial fit.

This is a generic error that is simply added to measurement_vec as if

\[y = y_0 + \epsilon(p_0,\; p_1,\; \cdots,\; p_n),\]

where \(y\) represents measurement_vec and \(y_0\) is the measurement vector without any errors)

Order 0 means constant: \(y = y_0 + a\)

Order 1 means linear: \(y = y_0 + a + b t\)

and so on. The derivatives that are added to the model_state_vec are those with regards to a, b, etc..

Note

The rule for the sensor_elem GIN is a bit complex. Generally, methods such as measurement_sensorAddSimple() will simply add a single unique frequency grid to all the different SensorObsel that they add to the measurement_sensor. The GIN sensor_elem is 0 for the first unique frequency grid, 1 for the second, and so on. See ArrayOfSensorObsel member methods in python for help identifying and manipulating how many unique frequency grids are available in measurement_sensor.

This method wraps jac_targetsAddErrorPolyFit() together with adding the covariance matrices, to oem.covmat_diagonal_blocks for assembly by oemFinalizeDiagonal().

The input covariance matrices must fit the size of the later computed model state represented by jac_targets. The inverse block is optional.

Warning

Automatic size constraints are not checked for this method. Group-invariant checks on read-only inputs still apply.

Author: Richard Larsson

Parameters:

  • jac_targets (~pyarts3.arts.JacobianTargets, optional) – A list of targets for the Jacobian Matrix calculations. Defaults to self.jac_targets. [INOUT]

  • oem (~pyarts3.arts.OptimalEstimationData, optional) – Numerical problem and results for oemCalc() and oemCalcReduced(). Defaults to self.oem. [INOUT]

  • measurement_sensor (~pyarts3.arts.ArrayOfSensorObsel, optional) – A list of sensor elements that fully describe one or more observing sensor(s). Defaults to self.measurement_sensor. [IN]

  • t (Vector) – The grid of \(y\). As \(t\) above. [IN]

  • sensor_elem (Index) – The sensor element whose frequency grid to use. [IN]

  • polyorder (~pyarts3.arts.Index, optional) – The order of the polynomial fit. Maximum \(n\) above. Defaults to 0 [IN]

  • matrix (BlockMatrix) – The covariance diagonal block matrix. [IN]

  • inverse (~pyarts3.arts.BlockMatrix, optional) – The inverse covariance diagonal block matrix. Defaults to pyarts3.arts.BlockMatrix() [IN]