CT Image Reconstruction Using Polychromatic Model and Optimization
Find Innovative SolutionsGenerate Solutions
Solution Overview
Problem
Existing CT image reconstruction methods, particularly for dual-energy CT, rely on approximate line integral models that ignore X-ray energy information, leading to quantitatively inaccurate results and significant beam-hardening artifacts due to the mismatch between physical models.
Innovation Solution
Implementing a realistic polychromatic physical model combined with an analytical algorithm and a single-variable optimization method to solve the non-linear polychromatic X-ray integral model in the projection domain, allowing for efficient and accurate decomposition of sinograms into two physical basis components.
Engineering Contradictions & Design Principles
Engineering Contradiction Analysis
1Ease of manufacture
If an approximate line integral model is used for dual-energy CT reconstruction, then the computation is simpler, but the measurement precision and manufacturing precision deteriorate due to quantitative inaccuracy and beam-hardening artifacts
Solution Approach 1:
The patent transforms the non-linear polychromatic X-ray integral model into a linear model by changing parameters through logarithmic transformation and material decomposition. This allows the use of efficient linear reconstruction algorithms while maintaining quantitative accuracy by properly accounting for polychromatic effects through the transformation: I(E) = I0(E) * exp(-∫μ(E,r)dl), where the linearized form enables accurate material decomposition without beam-hardening artifacts
Solution Approach 2:
The patent introduces basis material decomposition as an intermediary step between raw projection data and final images. By decomposing the attenuation coefficient into basis materials (μ(E,r) = Σwi(r)μi(E)), the system mediates between the complex polychromatic physics and the simplified linear reconstruction, enabling accurate quantitative results while maintaining computational efficiency
2Measurement precision
If a realistic polychromatic physical model with single-variable optimization is used, then the measurement precision and manufacturing precision improve, but the device complexity increases
Solution Approach 1:
The patent extracts the non-linear optimization problem into a separate material decomposition step, isolating the complexity from the main reconstruction pipeline. By separating the linear reconstruction (FBP) from the non-linear material decomposition, the system achieves high quantitative accuracy while keeping the overall device complexity manageable through modular architecture
Solution Approach 2:
The patent segments the reconstruction process into distinct stages: (1) linearized reconstruction using FBP, (2) material decomposition through basis materials, and (3) iterative optimization for refinement. This segmentation allows each component to be optimized independently, reducing overall system complexity while maintaining high measurement precision
3Ease of operation
If existing image-domain reconstruction methods are used, then the ease of operation improves, but the measurement precision deteriorates due to substantial approximations in energy spectra
Solution Approach 1:
The patent inverts the traditional approach by performing material decomposition in the projection domain rather than the image domain. This inversion allows the decomposition to occur before any approximations are made, preserving energy spectrum accuracy while maintaining operational simplicity through the linearized model: ln[I0(E)/I(E)] = ∫μ(E,r)dl, which can be directly solved for basis material concentrations
Applied Scientific Principles
This section explains which scientific principles are used to turn an abstract innovation direction into a practical engineering solution.
Function Achieved in This Case
This approach significantly reduces computation costs and increases convergence speed, resulting in a more practical and efficient image reconstruction method that accurately decomposes components and minimizes beam-hardening artifacts, producing high-quality CT images.
Implementation Method 1
The attenuation coefficient can be represented by a linear combination of photoelectric and Compton scattering components
Implementation Method 2
The attenuation coefficient can be represented by a linear combination of photoelectric and Compton scattering components
Implementation Method 3
X-ray detectors are currently operated in a current-integrating mode
Data Source
AI summary
Systems and methods for reconstructing images for computed tomography are provided. Image reconstruction can be based on a realistic polychromatic physical model, and can include use of both an analytical algorithm and a single-variable optimization method. The optimization method can be used to solve the non-linear polychromatic X-ray integral model in the projection domain, resulting in an accurate decomposition for sinograms of two physical basis components.