cosmic_ray_reject — ASTROSTACK cosmic op

Datenarten: image2dimage2d

Aufruf: import astrostack; astrostack.cosmic_ray_reject(frame, sigma=5.0, f_lim=2.0, replace_box=5, iters=1) (oder opsastrostack.get("cosmic_ray_reject"))

Verwendung

Entfernung kosmischer Strahlen aus einem Einzelbild (Laplace-Schärfe).

> Die ausführliche Beschreibung unten ist der Originaltext — Zusammenfassung und Überschriften sind übersetzt.

van Dokkum, *Cosmic-Ray Rejection by Laplacian Edge Detection*,

PASP 113, 1420 (2001) の考え方 —— 宇宙線は星より鋭い。星は PSF で

ぼけているので必ず数画素にまたがるが、宇宙線は光学系を通っていないので

1〜数画素で立ち上がる。そこで

1. ラプラシアン `L` の正の成分を雑音で規格化した有意度

`S = L / (2 sigma_noise)` を作り、

2. 微細構造像 `F = median3 - median7(median3) と比べて L / F` が

`f_lim` を超えるものだけを宇宙線とする。

2 番目の条件が無いと星の中心が必ず宇宙線に見える(星も局所的には

尖っている)。`f_lim` はその境目で、原論文の推奨は 2.0。

`iters` を増やすと、除去 → 再測定を繰り返す(大きなヒットの裾が残るとき)。

置換は `replace_box` の窓での非汚染画素の中央値

Returns `(cleaned, mask)`:

• `cleaned —— (H, W)` float64、宇宙線画素を置換した像。

• `mask —— (H, W) bool、True` = 宇宙線と判定した画素。

Raises `ValueError`: 2-D でない / 非有限を含む / *sigma* が非正 /

*f_lim* が非正 / *replace_box* が 3 未満または偶数の場合。

Referenzen (Beispieldaten, Literatur)

• Katalog der Beispieldaten (Download-URLs / Lizenzen) — 2-D nutzt skimage.data (BSD/Public Domain) plus synthetische Bilder, 3-D nennt Download-URLs echter Datenquellen (Stanford, PDS, …).

• Herkunft und Literatur der Operatoren — die Quellen der Forschung/Verfahren, auf denen diese Operatorfamilie beruht.

• Der kanonische Algorithmus (Autor, Jahr) und seine Anwendungen stehen im Familienleitfaden oben.

Ausführbare Beispiele (verifizierte Samples, die diesen Operator wirklich aufrufen)

astro_stackingpy -3.11 examples/astro_stacking.py

Typkompatible Folge-Operatoren (nehmen image2d als Eingabe)

frame_quality · noise_sigma · star_detect · psf_fit · aperture_photometry · frame_align

Gleiche Kategorie (cosmic)

cosmic_ray_reject_stack


*Provenance: astrostack.py — ASTROSTACK Operator-Registry. Diese Notiz wird von tools/opdocs.py md erzeugt (nicht von Hand bearbeiten).*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.