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.