mascado.distortions.analysis module

mascado.distortions.analysis.analyze_contributions_over_order(vf, minorder=0, maxorder=None, gridsize=100, units='arcseconds', info=' ')[source]

Get contributions of the terms of a specific degree.

The contribution is measured as RMS of the vector field over the normalized coordinates \([-1,1]\times[-1,1]\). The model is evaluated on a regular grid to calculate the RMS value. The grid has to be well sampled for the result to be correct. For higher orders, a larger gridsize is needed.

Parameters:
  • vf (mascado.distortions.polynomials.PolyVectorField) – Vector field with internal parameters.
  • minorder (int)
  • maxorder (int or None) – If None, the order of the vector field is used.
  • gridsize (int) – Number of points along each axis in normalized coordinates.
  • units (str) – Units of vector field. Used for console output.
  • info (str or False) – Indent of console output. Set to False to suppress output.
Returns:

DataFrame with two columns: order and modelrms.

Return type:

pandas.DataFrame

mascado.distortions.analysis.analyze_residuals_over_order(positions, distortions, maxorder, minorder=0, poly=<class 'mascado.distortions.polynomials.Legendre'>, distortionunits='arcseconds', maxcondition=100.0, info=' ')[source]

Get RMS of residuals after nth order fit.

Parameters:
  • positions ((N, 2)-shaped array) – Normalized positions.
  • distortions ((N, 2)-shaped array) – Distortions.
  • maxorder (int) – End of iteration.
  • minorder (int) – Start of iteration.
  • poly (mascado.distortions.polynomials.Polynomial subclass) – Class of polynomials to use (not an instance). By default, Legendre polynomials are used.
  • distortionunits (str) – Units of distortions. Used for console output.
  • info (str or False) – Indent of console output. Set to False to suppress output.
Returns:

DataFrame with two columns: order and resrms.

Return type:

pandas.DataFrame