ct_dual_energy_reconstruct — TOMOGRAPHY 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.