ct_scan op• 数据种类:projstack × projstack × ctgeometry → table
• 调用: import fullseye as fs; fs.ledger.ct_dual_energy_reconstruct(proj_low, proj_high, geometry, mu_low=(1.0, 0.4), mu_high=(0.6, 0.3), method='fdk', device='auto', **options)(要直接调用实现,import ctscan; ctscan.ct_dual_energy_reconstruct(proj_low, proj_high, geometry, mu_low=(1.0, 0.4), mu_high=(0.6, 0.3), method='fdk', device='auto', **options);从台账取用则 opstomography.get("ct_dual_energy_reconstruct"))
• GPU:无(仅 CPU 参考实现)
由同一几何的两次扫描得到两个基材料体数据。
> 以下的详细说明为原文 —— 摘要与标题已翻译。
At each detector sample the two line integrals are two linear equations in the two material path lengths, L = M t (Alvarez & Macovski 1976). The system is solved in the projection domain, where beam hardening has not yet been smeared across the volume, and each material is reconstructed on its own. *mu_low* and *mu_high* are effective attenuation coefficients; for a polychromatic beam, call dualxray.basis_material_decompose per view first and reconstruct its output.
The noise of the two path lengths grows with the condition number of M, which is returned.
:param proj_low: (view, V, U) line integrals of the low-energy scan. :param proj_high: the high-energy scan, same shape. :param geometry: dict from ct_scan_geometry. :param mu_low: attenuation of the two basis materials at the low energy. :param mu_high: attenuation of the two basis materials at the high energy. :param method: "fdk", "sirt", "os_sart" or "tv". :param device: "auto", "cpu" or "cuda". :param options: passed to the reconstruction. :returns: dict with basis_0 and basis_1 (path-length densities, volumes),
and condition.
:raises ValueError: on mismatched shapes or a singular matrix. Empty input
(proj_low with 0 elements) raises ValueError.
• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。
• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。
• 算法的正典(作者、年份)与用途见上面的族使用指南。
• 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
table 作为输入)—
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 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.