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 패밀리 가이드
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.
• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.
• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.
• converge_tour — py -3.11 examples/converge_tour.py
table 를 입력으로 받는 것)—
stopping)—
*Provenance: deconvstop.py — CONVERGE 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.