ct_scan op• データ種: projstack × ctgeometry → voxel
• 呼び出し: import fullseye as fs; fs.ledger.tv_reconstruct(projections, geometry, tv_weight=0.0003, n_iter=200, initial=None, nonnegative=True, device='auto', on_iter=None) (実装を直接呼ぶなら import ctscan; ctscan.tv_reconstruct(projections, geometry, tv_weight=0.0003, n_iter=200, initial=None, nonnegative=True, device='auto', on_iter=None)、台帳から引くなら opstomography.get("tv_reconstruct"))
• GPU: なし(CPU 参照実装のみ)
Chambolle-Pock 主双対法による全変動(TV)再構成。
> 以下の詳細説明は原文のままです —— 要約と見出しは訳出済み。
Solves min 0.5 ||A x - p||^2 + w ||grad x||_1 with x >= 0, in the form Sidky, Jorgensen & Pan 2012 give for CT. The prior says that the object is piecewise constant, which is what lets it fill data that was never measured: few views and limited angle. On an object that is not piecewise constant it paints the same flat patches, so a low residual here is not evidence of truth.
*tv_weight* is in units of ||A||^2 times the volume scale, so that one value behaves alike across grid sizes.
:param projections: (view, V, U) line integrals. :param geometry: dict from ct_scan_geometry. :param tv_weight: strength of the prior, >= 0. :param n_iter: primal-dual iterations, 1 .. 5000. :param initial: starting volume; None gives zeros. :param nonnegative: enforce x >= 0. :param device: "auto", "cpu" or "cuda". :param on_iter: optional on_iter(iteration, volume); returning a boolean
True (bool or numpy.bool_) stops after that iteration, any other
return value is ignored (same rule as sirt_reconstruct).
:returns: (Z, Y, X) float volume. :raises ValueError: on a shape that disagrees with the geometry or a negative
weight. Empty input (projections with 0 elements) raises ValueError.
• サンプルデータ カタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
• 演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。
• アルゴリズムの正典(著者・年)と用途は上記ファミリ使い方ガイドに記載。
• 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 operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.