cosmic op• 資料種類:image2d → image2d
• 呼叫: import fullseye as fs; fs.ledger.cosmic_ray_reject(frame, sigma=5.0, f_lim=2.0, replace_box=5, iters=1)(要直接呼叫實作,import astrostack; astrostack.cosmic_ray_reject(frame, sigma=5.0, f_lim=2.0, replace_box=5, iters=1);從台帳取用則 opsastrostack.get("cosmic_ray_reject"))
• 經台帳呼叫的回傳值: fullseye.ledger.cosmic_ray_reject(...) 只回傳**宣告的 out 型別 image2d 的值**(底層函式還會回傳輔助資訊)。需要被捨棄的部分時,請用 fullseye.ledger.cosmic_ray_reject.raw(...),或直接呼叫 astrostack.cosmic_ray_reject。
• 底層回傳: (cleaned, mask) -> image2d
單影格的宇宙線剔除(拉普拉斯銳度)。
> 以下的詳細說明為原文 —— 摘要與標題已翻譯。
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 未満または偶数の場合。
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• 演算法的正典(作者・年份)與用途見上面的族使用指南。
• astro_stacking — py -3.11 examples/astro_stacking.py
image2d 作為輸入)frame_quality · noise_sigma · star_detect · psf_fit · aperture_photometry · frame_align
cosmic)*Provenance: astrostack.py — ASTROSTACK 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.