ms_clean#

xrayvision.clean.ms_clean(dirty_map, dirty_beam, pixel_size, scales=None, clean_beam_width=<Quantity 4. arcsec>, gain=0.1, thres=0.01, niter=5000)[source]#

Clean the map using a multiscale clean algorithm.

Parameters:
  • dirty_map (Quantity) – The 2D dirty map to be cleaned

  • dirty_beam (Quantity) – The 2D dirty beam should have the same dimensions as dirty_map

  • pixel_size (Annotated[Quantity, Unit("arcsec / pix")]) – The pixel size in arcsec

  • scales (array-like, optional, optional) – The scales to use eg [1, 2, 4, 8]

  • clean_beam_width (Quantity) – The width of the gaussian to convolve the model with. If set to 0.0 the gaussian convolution is disabled

  • gain (float) – The gain per loop or loop gain

  • thres (float) – Terminates clean when residuals.max() <= thres`

  • niter (int) – Maximum number of iterations to perform

Returns:

Cleaned image

Return type:

Quantity | ndarray[Any, dtype[float64]]

Notes

This is an implementation of the multiscale clean algorithm as outlined in [R1] adapted for x-ray Fourier observations.

It is based on the on the implementation in the CASA software which can be found here ~L956.

References

[R1]

Cornwell, T. J., “Multiscale CLEAN Deconvolution of Radio Synthesis Images”, IEEE Journal of Selected Topics in Signal Processing, vol 2, p793-801, Paper #noqa