stopping op• 資料種類:image2d → table
• 呼叫: import fullseye as fs; fs.ledger.iv_richardson_lucy_until(image=None, psf_sigma=1.5, psf=None, noise='poisson', sigma=None, gain=1.0, tau=1.0, max_iter=2000, plateau_window=25, plateau_rtol=0.0001, keep_history=True, callback=None) -> 'dict'(要直接呼叫實作,import deconvstop; deconvstop.iv_richardson_lucy_until(image=None, psf_sigma=1.5, psf=None, noise='poisson', sigma=None, gain=1.0, tau=1.0, max_iter=2000, plateau_window=25, plateau_rtol=0.0001, keep_history=True, callback=None) -> 'dict';從台帳取用則 opsconverge.get("iv_richardson_lucy_until"))
依偏差原理停止的 Richardson-Lucy 反卷積（卜瓦松：Bertero 2010，τ = 1；高斯：Morozov）；拒絕模型值域之外的資料。
> 以下的詳細說明為原文 —— 摘要與標題已翻譯。
image: the blurred, noisy 2-D image. psf: a 2-D PSF (normalised here), or a Gaussian of psf_sigma px. The blur is periodic (FFT), so the RL step x <- x * H^T(y / H x) keeps sum(H x) = sum(y) exactly (Lemma 3.2 of Benvenuto & Piana).
• noise="poisson": image * gain must be photon counts (non-negative integers within 1e-6); stop at the first k with the Poisson discrepancy D_k <= tau (Benvenuto & Piana 2012 eq. (30); tau = 1 as proposed by Bertero et al. 2010).
• noise="gaussian": sigma (the noise standard deviation, in image units) is required; stop at the first k with ||H x_k - y||^2 / (N sigma^2) <= tau (Morozov, eq. (18)). The image must be non-negative (add the known offset first).
The iteration starts from the flat image of the data's mean (k = 0). Returns image (x_k in the input units), stopped, status, k_stop, discrepancy_history (D_0, D_1, ...), mass_rel_error_max (Lemma 3.2), tau and noise.
Not stopped -> image = None and status:
• "outside_model_range" — the discrepancy levelled off above tau (relative decrease below plateau_rtol per iteration over plateau_window iterations): y is not in the range of the model, so the criterion can never be met (Benvenuto & Piana: the rule is a regularisation iff y is in H(C)). Typical causes: the assumed PSF is wider than the real blur, a background is missing, the noise level is wrong.
• "max_iter" — still decreasing when max_iter ran out.
The last iterate is under last_iterate (not an estimate). callback(k, x_k, D_k), if given, is called after every iteration (x_k in the input units) — for plotting the run.
• convergent_computation 族使用指南
• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。
• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。
• 演算法的正典(作者、年份)與用途見上面的族使用指南。
• converge_tour — py -3.11 examples/converge_tour.py
table 作為輸入)—
stopping)—
*Provenance: deconvstop.py — CONVERGE 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.