ct_scan op• 데이터 종류: projstack × ctgeometry → voxel
• 호출: import fullseye as fs; fs.ledger.os_sart_reconstruct(projections, geometry, n_iter=10, n_subsets=16, relaxation=0.8, initial=None, nonnegative=True, device='auto', on_iter=None)(구현을 직접 호출하려면 import ctscan; ctscan.os_sart_reconstruct(projections, geometry, n_iter=10, n_subsets=16, relaxation=0.8, initial=None, nonnegative=True, device='auto', on_iter=None), 원장에서 가져오려면 opstomography.get("os_sart_reconstruct"))
• GPU: 없음(CPU 참조 구현만)
OS-SART: SIRT 갱신을 순서가 있는 뷰 부분집합마다 적용한다.
> 아래 상세 설명은 원문입니다 —— 요약과 제목은 번역되어 있습니다.
One iteration passes over all subsets, so it costs about one SIRT iteration and moves about *n_subsets* times as far early on. With inconsistent data it settles into a limit cycle rather than a point.
:param projections: (view, V, U) line integrals. :param geometry: dict from ct_scan_geometry. :param n_iter: passes over all subsets, 1 .. 5000. :param n_subsets: subsets, 1 .. n_views. :param relaxation: step lam, (0, 2). :param initial: starting volume; None gives zeros. :param nonnegative: clip to >= 0 after each subset. :param device: "auto", "cpu" or "cuda". :param on_iter: optional on_iter(iteration, volume), called after each full
pass; returning a boolean True (bool or numpy.bool_) stops
there, any other return value is ignored.
:returns: (Z, Y, X) float volume. :raises ValueError: on a shape that disagrees with the geometry, a relaxation
outside (0, 2), or more subsets than views. Empty input (projections
with 0 elements) raises ValueError.
The default 10 passes is not converged. On consistent data (24^3 Shepp-Logan, cone, 20 views, 16 subsets) the relative residual is 0.0135 at 10 passes and 0.0044 at 40, nRMS to the truth 0.0278 -> 0.0176 (37 %); on the measured walnut (fan, 50 views, 10 subsets) the RMSE to the reference is 0.0048 /mm at 10 passes and 0.0036 at 40 (2026-11-09 audit). Pick the count from a curve (on_iter, see sirt_reconstruct).
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.
• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.
• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.
• ct_scan_modes — py -3.11 examples/ct_scan_modes.py
• poc_ct_casting_sparse_views — py -3.11 examples/poc_ct_casting_sparse_views.py
• poc_real_walnut_cbct — py -3.11 examples/poc_real_walnut_cbct.py
voxel 를 입력으로 받는 것)radon_volume · ct_forward_project
ct_scan)ct_scan_geometry · shepp_logan_3d · ellipsoid_phantom · ellipsoid_projections · ct_forward_project · ct_adjoint_project · fdk_reconstruct · sirt_reconstruct
*Provenance: ctscan.py — TOMOGRAPHY 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.