cosmic_ray_reject — ASTROSTACK cosmic op

데이터 종류: image2dimage2d

호출: import astrostack; astrostack.cosmic_ray_reject(frame, sigma=5.0, f_lim=2.0, replace_box=5, iters=1)(또는 opsastrostack.get("cosmic_ray_reject"))

사용법

단일 프레임의 우주선 제거(라플라시안 첨예도).

> 아래 상세 설명은 원문입니다 —— 요약과 제목은 번역되어 있습니다.

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 未満または偶数の場合。

참고(샘플 데이터·문헌)

• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.

• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.

• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.

실행 가능한 예제(이 연산자를 실제로 호출하는 검증된 샘플)

astro_stackingpy -3.11 examples/astro_stacking.py

타입이 이어지는 다음 연산자(image2d 를 입력으로 받는 것)

frame_quality · noise_sigma · star_detect · psf_fit · aperture_photometry · frame_align

같은 카테고리(cosmic)

cosmic_ray_reject_stack


*Provenance: astrostack.py — ASTROSTACK 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

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