kernels.ccdclip_mask_1d

kernels.ccdclip_mask_1d(
    values,
    *,
    mask=None,
    sigma=(3.0, 3.0),
    maxiters=5,
    ddof=0,
    nkeep=1,
    maxrej=None,
    cenfunc='median',
    clip_cen=None,
    revert_on_nkeep=True,
    rdnoise=0.0,
    gain=1.0,
    snoise=0.0,
    scale_ref=1.0,
    zero_ref=0.0,
    validate=True,
)

CCD noise-model clipping (1-D mask-only variant).

Parameters

values : (ndarray, shape(N))

One-dimensional value vector. Accepted dtypes are uint8, uint16, int16, int32, float32, and float64. Integer inputs are promoted to the package’s floating workspace when validate is True.

mask : ndarray of bool = None

Input mask; True means already masked. Must have shape (N,).

sigma : float or tuple of float = (3.0, 3.0)

User-supplied clipping multiplier, not the measured data spread itself. Values are clipped when their residual from the clipping center exceeds sigma * spread. A scalar applies the same multiplier to both tails; a 2-tuple (sigma_lower, sigma_upper) clips asymmetrically. Thresholds must be finite and non-negative.

maxiters : int = 5

Maximum number of clipping iterations per pixel.

ddof : int = 0

Delta degrees of freedom for the per-pixel spread estimator.

nkeep : int = 1

Minimum number of unmasked values to preserve at each pixel. An iteration that would drop below this count is reverted in full. The default, 1, prevents clipping from rejecting every finite sample in a per-output stack. Set 0 only when all samples may be rejected.

maxrej : int or None = None

Maximum number of values that may be rejected at each pixel. None means no limit.

cenfunc : ('median', 'lmedian', 'mean') = "median"

Center estimator used each iteration. lmedian accepts the alias lmed and uses the lower of the two middle samples for even valid counts.

clip_cen : ('median', 'lmedian', 'mean') = "median"

Center used when computing the spread for iterative rejection. None uses the same center as cenfunc; lmedian also accepts lmed.

revert_on_nkeep : bool = True

If True, an iteration that would leave fewer than nkeep usable samples at a pixel is reverted in full for that pixel. Input masks, threshold masks, and non-finite samples remain excluded. The default is True. If False, clipping is strict and may leave fewer than nkeep samples.

rdnoise : float = 0.0

Read noise in electrons.

snoise : float = 0.0

Sky-noise coefficient. Adds a Poisson-like contribution scaled by the local count level.

scale_ref : float = 1.0

Reference scale factor. Typically the inverse of the exposure-time normalization applied during reduction.

zero_ref : float = 0.0

Reference zero-point offset. Accounts for a known DC bias that shifts the effective noise level.

validate : bool = True

If True, check dimensionality and normalize dtype/contiguity before entering the Rust kernel. If False, callers must provide inputs that satisfy the compiled kernel assumptions.

Returns

mask_rej : ndarray of bool, shape (N,)

True where a value was rejected by this kernel.

Notes

The kernel evaluates CCD rejection on gain-corrected scratch values (electrons). The per-pixel noise threshold is approximately::

sqrt((1 + snoise) * abs(noise_center + zero_ref) * scale_ref + rdnoise**2)