make_mask
- acstools.satdet.make_mask(filename, ext, trail_coords, sublen=75, subwidth=200, order=3, sigma=4, pad=10, plot=False, verbose=False)
Create DQ mask for an image for a given satellite trail. This mask can be added to existing DQ data using
acstools.utils_findsat_mrt.update_dq()
.Note
Unlike
detsat()
, multiprocessing is not available for this function.- Parameters:
filename (str) – FITS image filename.
ext (int, str, or tuple) – Extension for science data, as accepted by
astropy.io.fits
.trail_coords (ndarray) – One of the trails returned by
detsat()
. This must be in the format of[[x0, y0], [x1, y1]]
.sublen (int, optional) – Length of strip to use as the fitting window for the trail.
subwidth (int, optional) – Width of box to fit trail on.
order (int, optional) – The order of the spline interpolation for image rotation. See
skimage.transform.rotate()
.sigma (float, optional) – Sigma of the satellite trail for detection. If points are a given sigma above the background in the subregion then it is marked as a satellite. This may need to be lowered for resolved trails.
pad (int, optional) – Amount of extra padding in pixels to give the satellite mask.
plot (bool, optional) – Plot the result.
verbose (bool, optional) – Print extra information to the terminal, mostly for debugging.
- Returns:
mask – Boolean array marking the satellite trail with
True
.- Return type:
ndarray
- Raises:
ImportError – Missing scipy or skimage>=0.11 packages.
IndexError – Invalid subarray indices.
ValueError – Image has no positive values, trail subarray too small, or trail profile not found.