This module provides comprehensive functions for extracting 1D spectra from IFU data cubes using various extraction methods.

Key functionalities include:

• Simple aperture extraction with user-defined or automatic masks
• Optimal extraction following Horne (1986) algorithm
• Multiple mask generation methods (circular, Gaussian fit, position-based)
• Bad pixel handling and quality flag propagation
• Variance and error propagation throughout extraction
• Adaptive radius estimation for point sources
• Weighted extraction with running median smoothing
• Cosmic ray rejection and outlier clipping

Extraction methods supported:

• MASK: Use provided external mask
• MAX: Automatic center detection from maximum flux
• POSITION: User-specified center position
• OPTIMAL: Horne (1986) optimal extraction with PSF modeling
• FIT: 2D Gaussian fitting for PSF determination

The optimal extraction implements a modified version of Horne (1986):

• PSF model constructed from data
• Variance taken from input cube (not updated iteratively)
• Slice-by-slice outlier clipping
• Running weighted average smoothing instead of polynomial fitting
• Quality flags incorporated into rejection scheme

Note

All extraction methods properly propagate errors and handle bad pixels

Mask convention: 1 = good pixel, 0 = bad pixel (opposite of CPL binary masks)