Variance Estimation

cwitools.analysis.estimate_variance(inputfits, zwindow=10, rescale=True, sigmaclip=4, zmask=(0, 0), fmin=0.9, fmax=10)

Estimates the 3D variance cube of an input cube.

Parameters:

• cubePath (str) – Path to the input cube.
• zWindow (int) – Size of z-axis bins to use for 2D variance estimation. Default: 10.
• zmask (int tuple) – Wavelength layers to exclude while estimating variance.
• rescale (bool) – Set to TRUE to perform layer-by-layer rescaling of 2D variance.
• fMin (float) – The minimum rescaling factor (Default 0.9)
• fMax (float) – The maximum rescaling factor (Default: 10)
• fileExt (str) – The extension to use for the output cube (Default .var.fits)

Returns:

Estimated variance cube

Return type:

NumPy ndarray