cwitools.subtraction.psf_sub_1d

cwitools.subtraction.psf_sub_1d(fits_in, pos, fit_rad=2, sub_rad=15, slice_rad=4, wl_window=150, wmasks=[], slice_axis=2)

Models and subtracts a point-source on a slice-by-slice basis.

Args:

fits_in (astrop FITS object): Input data cube/FITS. pos (float tuple): (x,y) position of the source to subtract. fit_rad (float): Inner radius, in arcsec, to use for fitting PSF within a slice. sub_rad (float): Outer radius, in arcsec, over which to subtract PSF within a slice. slice_rad (int): Number of slices to subtract over, as distance from the

slice indicated by the `pos’ argument.

wl_window (float): Size of white-light window (in Angstrom) to use.

This is the window used to form a white-light image centered on each wavelength layer. Default: 150A.

wmasks (list): List of wavelength tuples to exclude when making

white-light images. Use to exclude nebular emission or sky lines.

slice_axis (int): Which axis represents the slices of the image.

For KCWI default data cubes, slice_axis = 2. For PCWI data cubes, slice_axis = 1.

Returns:

numpy.ndarray: PSF-subtracted data cube numpy.ndarray: PSF model data