Multi-Material Decomposition for Single Energy CT Using Material Sparsity Constraint

Prior image-based generative adversarial learning for multi-material decomposition in photon counting computed tomography

BACKGROUND

Photon counting detector computed tomography (PCD-CT) is a novel promising technique providing higher spatial resolution, lower radiation dose and greater energy spectrum differentiation, which create more possibilities to improve image quality. Multi-material decomposition is an attractive application for PCD-CT to identify complicated materials and provide accurate quantitative analysis. However, limited by the finite photon counting rate in each energy window of photon counting detector, the noise problem hinders the decomposition of high-quality basis material images.

METHODS

To address this issue, an end-to-end multi-material decomposition network based on prior images is proposed in this paper. First, the reconstructed images corresponding to the full spectrum with less noise are introduced as prior information to improve the overall signal-to-noise ratio of the data. Then, a generative adversarial network is designed to mine the relationship between reconstructed images and basis material images based on the information interaction of material decomposition. Furthermore, a weighted edge loss is introduced to adapt to the structural differences of different basis material images.

RESULTS

To verify the performance of the proposed method, simulation and real studies are carried out. In simulation study of structured fibro-glandular tissue model, the results show that the proposed method decreased the root mean square error by 67 % and 26 % on adipose, 66 % and 28 % on fibroglandular, 52 % and 8 % on calcification, compared to butterfly network and dual interactive Wasserstein generative adversarial network.

CONCLUSION

Experimentally, the proposed method shows certain advantages over other methods on noise suppression effect, detail retention ability and decomposition accuracy.

Image Domain Multi-Material Decomposition Noise Suppression Through Basis Transformation and Selective Filtering

Spectral CT can provide material characterization ability to offer more precise material information for diagnosis purposes. However, the material decomposition process generally leads to amplification of noise which significantly limits the utility of the material basis images. To mitigate such problem, an image domain noise suppression method was proposed in this work. The method performs basis transformation of the material basis images based on a singular value decomposition. The noise variances of the original spectral CT images were incorporated in the matrix to be decomposed to ensure that the transformed basis images are statistically uncorrelated. Due to the difference in noise amplitudes in the transformed basis images, a selective filtering method was proposed with the low-noise transformed basis image as guidance. The method was evaluated using both numerical simulation and real clinical dual-energy CT data. Results demonstrated that compared with existing methods, the proposed method performs better in preserving the spatial resolution and the soft tissue contrast while suppressing the image noise. The proposed method is also computationally efficient and can realize real-time noise suppression for clinical spectral CT images.

Quantitative multi-energy micro-CT: A simulation and phantom study for simultaneous imaging of four different contrast materials using an energy integrating detector

Emerging from the development of single-energy Computed Tomography (CT) and Dual-Energy Computed Tomography, Multi-Energy Computed Tomography (MECT) is a promising tool allowing advanced material and tissue decomposition and thereby enabling the use of multiple contrast materials in preclinical research. The scope of this work was to evaluate whether a usual preclinical micro-CT system is applicable for the decomposition of different materials using MECT together with a matrix-inversion method and how different changes of the measurement-environment affect the results. A matrix-inversion based algorithm to differentiate up to five materials (iodine, iron, barium, gadolinium, residual material) by applying four different acceleration voltages/energy levels was established. We carried out simulations using different ratios and concentrations (given in fractions of volume units, VU) of the four different materials (plus residual material) at different noise-levels for 30 keV, 40 keV, 50 keV, 60 keV, 80 keV and 100 keV (monochromatic). Our simulation results were then confirmed by using region of interest-based measurements in a phantom-study at corresponding acceleration voltages. Therefore, different mixtures of contrast materials were scanned using a micro-CT. Voxel wise evaluation of the phantom imaging data was conducted to confirm its usability for future imaging applications and to estimate the influence of varying noise-levels, scattering, artifacts and concentrations. The analysis of our simulations showed the smallest deviation of 0.01 (0.003-0.15) VU between given and calculated concentrations of the different contrast materials when using an energy-combination of 30 keV, 40 keV, 50 keV and 100 keV for MECT. Subsequent MECT phantom measurements, however, revealed a combination of acceleration voltages of 30 kV, 40 kV, 60 kV and 100 kV as most effective for performing material decomposition with a deviation of 0.28 (0-1.07) mg/ml. The feasibility of our voxelwise analyses using the proposed algorithm was then confirmed by the generation of phantom parameter-maps that matched the known contrast material concentrations. The results were mostly influenced by the noise-level and the concentrations used in the phantoms. MECT using a standard micro-CT combined with a matrix inversion method is feasible at four different imaging energies and allows the differentiation of mixtures of up to four contrast materials plus an additional residual material.

Imaging-proteomic analysis for prediction of neoadjuvant chemotherapy responses in patients with breast cancer

BACKGROUND

Optimizing patient selection for neoadjuvant chemotherapy in patients with breast cancer remains an unmet clinical need. Quantitative features from medical imaging were reported to be predictive of treatment responses. However, the biologic meaning of these latent features is poorly understood, preventing the clinical use of such noninvasive imaging markers. The study aimed to develop a deep learning signature (DLS) from pretreatment magnetic resonance imaging (MRI) for predicting responses to neoadjuvant chemotherapy in patients with breast cancer and to further investigate the biologic meaning of the DLS by identifying its underlying pathways using paired MRI and proteomic sequencing data.

METHODS

MRI-based DLS was constructed (radiogenomic training dataset, n = 105) and validated (radiogenomic validation dataset, n = 26) for the prediction of pathologic complete response (pCR) to neoadjuvant chemotherapy. Proteomic sequencing revealed biological functions facilitating pCR (n = 139). Their associations with DLS were uncovered by radiogenomic analysis.

RESULTS

The DLS achieved a prediction accuracy of 0.923 with an AUC of 0.958, higher than the performance of the model trained by transfer learning. Cellular membrane formation, endocytosis, insulin-like growth factor binding, protein localization to membranes, and cytoskeleton-dependent trafficking were differentially regulated in patients showing pCR. Oncogenic signaling pathways, features correlated with human phenotypes, and features correlated with general biological processes were significantly correlated with DLS in both training and validation dataset (p.adj < 0.05).

CONCLUSIONS

Our study offers a biologically interpretable DLS for the prediction of pCR to neoadjuvant chemotherapy in patients with breast cancer, which may guide personalized medication.

A generalizable new figure of merit for dose optimization in dual energy cone beam CT scanning protocols

Objective. This study proposes and evaluates a new figure of merit (FOMn) for dose optimization of Dual-energy cone-beam CT (DE-CBCT) scanning protocols based on size-dependent modeling of radiation dose and multi-scale image quality.Approach. FOMn was defined using Z-score normalization and was proportional to the dose efficiency providing better multi-scale image quality, including comprehensive contrast-to-noise ratio (CCNR) and electron density (CED) for CatPhan604 inserts of various materials. Acrylic annuluses were combined with CatPhan604 to create four phantom sizes (diameters of the long axis are 200 mm, 270 mm, 350 mm, and 380 mm, respectively). DE-CBCT was decomposed using image-domain iterative methods based on Varian kV-CBCT images acquired using 25 protocols (100 kVp and 140 kVp combined with 5 tube currents).Main results. The accuracy of CED was approximately 1% for all protocols, but degraded monotonically with the increased phantom sizes. Combinations of lower voltage + higher current and higher voltage + lower current were optimal protocols balancing CCNR and dose. The most dose-efficient protocols for CED and CCNR were inconsistent, underlining the necessity of including multi-scale image quality in the evaluation and optimization of DE-CBCT. Pediatric and adult anthropomorphic phantom tests confirmed dose-efficiency of FOMn-recommended protocols.Significance. FOMn is a comprehensive metric that collectively evaluates radiation dose and multi-scale image quality for DE-CBCT. The models and data can also serve as lookup tables, suggesting personalized dose-efficient protocols for specific clinical imaging purposes.

Statistical iterative spectral CT imaging method based on blind separation of polychromatic projections

Spectral computed tomography (CT) can provide narrow-energy-width reconstructed images, thereby suppressing beam hardening artifacts and providing rich attenuation information for component characterization. We propose a statistical iterative spectral CT imaging method based on blind separation of polychromatic projections to improve the accuracy of narrow-energy-width image decomposition. For direct inversion in blind scenarios, we introduce the system matrix into the X-ray multispectral forward model to reduce indirect errors. A constrained optimization problem with edge-preserving regularization is established and decomposed into two sub-problems to be alternately solved. Experiments indicate that the novel algorithm obtains more accurate narrow-energy-width images than the state-of-the-art method.