TIGRE: a MATLAB-GPU toolbox for CBCT image reconstruction

Sparse-view CBCT reconstruction using meta-learned neural attenuation field and hash-encoding regularization

Cone beam computed tomography (CBCT) is an emerging medical imaging technique to visualize the internal anatomical structures of patients. During a CBCT scan, several projection images of different angles or views are collectively utilized to reconstruct a tomographic image. However, reducing the number of projections in a CBCT scan while preserving the quality of a reconstructed image is challenging due to the nature of an ill-posed inverse problem. Recently, a neural attenuation field (NAF) method was proposed by adopting a neural radiance field algorithm as a new way for CBCT reconstruction, demonstrating fast and promising results using only 50 views. However, decreasing the number of projections is still preferable to reduce potential radiation exposure, and a faster reconstruction time is required considering a typical scan time. In this work, we propose a fast and accurate sparse-view CBCT reconstruction (FACT) method to provide better reconstruction quality and faster optimization speed in the minimal number of view acquisitions (< 50 views). In the FACT method, we meta-trained a neural network and a hash-encoder using a few scans (= 15), and a new regularization technique is utilized to reconstruct the details of an anatomical structure. In conclusion, we have shown that the FACT method produced better, and faster reconstruction results over the other conventional algorithms based on CBCT scans of different body parts (chest, head, and abdomen) and CT vendors (Siemens, Phillips, and GE).

Shape-matching-based fracture reduction aid concept exemplified on the proximal humerus-a pilot study

PURPOSE

Optimizing fracture reduction quality is key to achieve successful osteosynthesis, especially for epimetaphyseal regions such as the proximal humerus (PH), but can be challenging, partly due to the lack of a clear endpoint. We aimed to develop the prototype for a novel intraoperative C-arm-based aid to facilitate true anatomical reduction of fractures of the PH.

METHODS

The proposed method designates the reduced endpoint position of fragments by superimposing the outer boundary of the premorbid bone shape on intraoperative C-arm images, taking the mirrored intact contralateral PH from the preoperative CT scan as a surrogate. The accuracy of the algorithm was tested on 60 synthetic C-arm images created from the preoperative CT images of 20 complex PH fracture cases (Dataset A) and on 12 real C-arm images of a prefractured human anatomical specimen (Dataset B). The predicted outer boundary shape was compared with the known exact solution by (1) a calculated matching error and (2) two experienced shoulder trauma surgeons.

RESULTS

A prediction accuracy of 88% (with 73% 'good') was achieved according to the calculation method and an 87% accuracy (68% 'good') by surgeon assessment in Dataset A. Accuracy was 100% by both assessments for Dataset B.

CONCLUSION

By seamlessly integrating into the standard perioperative workflow and imaging, the intuitive shape-matching-based aid, once developed as a medical device, has the potential to optimize the accuracy of the reduction of PH fractures while reducing the number of X-rays and surgery time. Further studies are required to demonstrate the applicability and efficacy of this method in optimizing fracture reduction quality.

X-ray CT metal artifact reduction using neural attenuation field prior

BACKGROUND

The presence of metal objects in computed tomography (CT) imaging introduces severe artifacts that degrade image quality and hinder accurate diagnosis. While several deep learning-based metal artifact reduction (MAR) methods have been proposed, they often exhibit poor performance on unseen data and require large datasets to train neural networks.

PURPOSE

In this work, we propose a sinogram inpainting method for metal artifact reduction that leverages a neural attenuation field (NAF) as a prior. This new method, dubbed NAFMAR, operates in a self-supervised manner by optimizing a model-based neural field, thus eliminating the need for large training datasets.

METHODS

NAF is optimized to generate prior images, which are then used to inpaint metal traces in the original sinogram. To address the corruption of x-ray projections caused by metal objects, a 3D forward projection of the original corrupted image is performed to identify metal traces. Consequently, NAF is optimized using a metal trace-masked ray sampling strategy that selectively utilizes uncorrupted rays to supervise the network. Moreover, a metal-aware loss function is proposed to prioritize metal-associated regions during optimization, thereby enhancing the network to learn more informed representations of anatomical features. After optimization, the NAF images are rendered to generate NAF prior images, which serve as priors to correct original projections through interpolation. Experiments are conducted to compare NAFMAR with other prior-based inpainting MAR methods.

RESULTS

The proposed method provides an accurate prior without requiring extensive datasets. Images corrected using NAFMAR showed sharp features and preserved anatomical structures. Our comprehensive evaluation, involving simulated dental CT and clinical pelvic CT images, demonstrated the effectiveness of NAF prior compared to other prior information, including the linear interpolation and data-driven convolutional neural networks (CNNs). NAFMAR outperformed all compared baselines in terms of structural similarity index measure (SSIM) values, and its peak signal-to-noise ratio (PSNR) value was comparable to that of the dual-domain CNN method.

CONCLUSIONS

NAFMAR presents an effective, high-fidelity solution for metal artifact reduction in 3D tomographic imaging without the need for large datasets.

DLPVI: Deep learning framework integrating projection, view-by-view backprojection, and image domains for high- and ultra-sparse-view CBCT reconstruction

This study proposes a deep learning framework, DLPVI, which integrates projection, view-by-view backprojection (VVBP), and image domains to improve the quality of high-sparse-view and ultra-sparse-view cone-beam computed tomography (CBCT) images. The DLPVI comprises a projection domain sub-framework, a VVBP domain sub-framework, and a Transformer-based image domain model. First, full-view projections were restored from sparse-view projections via the projection domain sub-framework, then filtered and view-by-view backprojected to generate VVBP raw data. Next, the VVBP raw data was processed by the VVBP domain sub-framework to suppress residual noise and artifacts, and produce CBCT axial images. Finally, the axial images were further refined using the image domain model. The DLPVI was trained, validated, and tested on CBCT data from 163, 30, and 30 real patients respectively. Quantitative metrics including root-mean-square error (RMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and feature similarity (FSIM) were calculated to evaluate the method performance. The DLPVI was compared with 15 state-of-the-art (SOTA) methods, including 2 projection domain models, 10 image domain models, and 3 projection-image dual-domain frameworks, on 1/8 high-sparse-view and 1/16 ultra-sparse-view reconstruction tasks. Statistical analysis was conducted using the Kruskal-Wallis test, followed by the post-hoc Dunn's test. Experimental results demonstrated that the DLPVI outperformed all 15 SOTA methods for both tasks, with statistically significant improvements (p < 0.05 in Kruskal-Wallis test and p < 0.05/15 in Dunn's test). The proposed DLPVI effectively improves the quality of high- and ultra-sparse-view CBCT images.

Strategies for generating synthetic computed tomography-like imaging from radiographs: A scoping review

BACKGROUND

Advancements in tomographic medical imaging have revolutionized diagnostics and treatment monitoring by offering detailed 3D visualization of internal structures. Despite the significant value of computed tomography (CT), challenges such as high radiation dosage and cost barriers limit its accessibility, especially in low- and middle-income countries. Recognizing the potential of radiographic imaging in reconstructing CT images, this scoping review aims to explore the emerging field of synthesizing 3D CT-like images from 2D radiographs by examining the current methodologies.

METHODS

A scoping review was carried out following PRISMA-SR guidelines. Eligibility criteria for the articles included full-text articles published up to September 9, 2024, studying methodologies for the synthesis of 3D CT images from 2D biplanar or four-projection x-ray images. Eligible articles were sourced from PubMed MEDLINE, Embase, and arXiv.

RESULTS

76 studies were included. The majority (50.8 %, n = 30) were published between 2010 and 2020 (38.2 %, n = 29) and from 2020 onwards (36.8 %, n = 28), with European (40.8 %, n = 31), North American (26.3 %, n = 20), and Asian (32.9 %, n = 25) institutions being primary contributors. Anatomical regions varied, with 17.1 % (n = 13) of studies not using clinical data. Further, studies focused on the chest (25 %, n = 19), spine and vertebrae (17.1 %, n = 13), coronary arteries (10.5 %, n = 8), and cranial structures (10.5 %, n = 8), among other anatomical regions. Convolutional neural networks (CNN) (19.7 %, n = 15), generative adversarial networks (21.1 %, n = 16) and statistical shape models (15.8 %, n = 12) emerged as the most applied methodologies. A limited number of studies included explored the use of conditional diffusion models, iterative reconstruction algorithms, statistical shape models, and digital tomosynthesis.

CONCLUSION

This scoping review summarizes current strategies and challenges in synthetic imaging generation. The development of 3D CT-like imaging from 2D radiographs could reduce radiation risk while simultaneously addressing financial and logistical obstacles that impede global access to CT imaging. Despite initial promising results, the field encounters challenges with varied methodologies and frequent lack of proper validation, requiring further research to define synthetic imaging's clinical role.

RNAF: Regularization neural attenuation fields for sparse-view CBCT reconstruction

Cone beam computed tomography (CBCT) is increasingly used in clinical settings, with the radiation dose incurred during X-ray acquisition emerging as a critical concern. Traditional algorithms for reconstructing high-quality CBCT images typically necessitate hundreds of X-ray projections, prompting a shift towards sparse-view CBCT reconstruction as a means to minimize radiation exposure. A novel approach, leveraging the Neural Attenuation Field (NAF) based on neural radiation field algorithms, has recently gained traction. This method offers rapid and promising CBCT reconstruction outcomes using a mere 50 views. Nonetheless, NAF tends to overlook the inherent structural properties of projected images, which can lead to shortcomings in accurately capturing the structural essence of the object being imaged. To address these limitations, we introduce an enhanced method: Regularization Neural Attenuation Fields (RNAF). Our approach includes two key innovations. First, we implement a hash coding regularization technique designed to retain low-frequency details within the reconstructed images, thereby preserving essential structural information. Second, we incorporate a Local Patch Global (LPG) sampling strategy. This method focuses on extracting local geometric details from the projection image, ensuring that the intensity variations in randomly sampled X-rays closely mimic those in the actual projection image. Comparative analyses across various body parts (Chest, Jaw, Foot, Abdomen, Knee) reveal that RNAF substantially outperforms existing algorithms. Specifically, its reconstruction quality exceeds that of previous NeRF-based, optimization-based, and analysis algorithms by margins of at least 2.09 dB, 3.09 dB, and 13.84 dB respectively. This significant enhancement in performance underscores the potential of RNAF as a groundbreaking solution in the realm of CBCT imaging, offering a path towards achieving high-quality reconstructions with reduced radiation exposure. Our implementation is publically available at https://github.com/springXIACJ/FRNAF.

Cone-beam computed tomography (CBCT) image-quality improvement using a denoising diffusion probabilistic model conditioned by pseudo-CBCT of pelvic regions

Cone-beam computed tomography (CBCT) is widely used in radiotherapy to image patient configuration before treatment but its image quality is lower than planning CT due to scattering, motion, and reconstruction methods. This reduces the accuracy of Hounsfield units (HU) and limits its use in adaptive radiation therapy (ART). However, synthetic CT (sCT) generation using deep learning methods for CBCT intensity correction faces challenges due to deformation. To address these issues, we propose enhancing CBCT quality using a conditional denoising diffusion probability model (CDDPM), which is trained on pseudo-CBCT created by adding pseudo-scatter to planning CT. The CDDPM transforms CBCT into high-quality sCT, improving HU accuracy while preserving anatomical configuration. The performance evaluation of the proposed sCT showed a reduction in mean absolute error (MAE) from 81.19 HU for CBCT to 24.89 HU for the sCT. Peak signal-to-noise ratio (PSNR) improved from 31.20 dB for CBCT to 33.81 dB for the sCT. The Dice and Jaccard coefficients between CBCT and sCT for the colon, prostate, and bladder ranged from 0.69 to 0.91. When compared to other deep learning models, the proposed sCT outperformed them in terms of accuracy and anatomical preservation. The dosimetry analysis for prostate cancer revealed a dose error of over 10% with CBCT but nearly 0% with the sCT. Gamma pass rates for the proposed sCT exceeded 90% for all dose criteria, indicating high agreement with CT-based dose distributions. These results show that the proposed sCT improves image quality, dosimetry accuracy, and treatment planning, advancing ART for pelvic cancer.

Multi-stage deep learning artifact reduction for parallel-beam computed tomography

Computed tomography (CT) using synchrotron radiation is a powerful technique that, compared with laboratory CT techniques, boosts high spatial and temporal resolution while also providing access to a range of contrast-formation mechanisms. The acquired projection data are typically processed by a computational pipeline composed of multiple stages. Artifacts introduced during data acquisition can propagate through the pipeline and degrade image quality in the reconstructed images. Recently, deep learning has shown significant promise in enhancing image quality for images representing scientific data. This success has driven increasing adoption of deep learning techniques in CT imaging. Various approaches have been proposed to incorporate deep learning into computational pipelines, but each has limitations in addressing artifacts effectively and efficiently in synchrotron CT, either in properly addressing the specific artifacts or in computational efficiency. Recognizing these challenges, we introduce a novel method that incorporates separate deep learning models at each stage of the tomography pipeline - projection, sinogram and reconstruction - to address specific artifacts locally in a data-driven way. Our approach includes bypass connections that feed both the outputs from previous stages and raw data to subsequent stages, minimizing the risk of error propagation. Extensive evaluations on both simulated and real-world datasets illustrate that our approach effectively reduces artifacts and outperforms comparison methods.

Deep Learning-Based Multi-View Projection Synthesis Approach for Improving the Quality of Sparse-View CBCT in Image-Guided Radiotherapy

While radiation hazards induced by cone-beam computed tomography (CBCT) in image-guided radiotherapy (IGRT) can be reduced by sparse-view sampling, the image quality is inevitably degraded. We propose a deep learning-based multi-view projection synthesis (DLMPS) approach to improve the quality of sparse-view low-dose CBCT images. In the proposed DLMPS approach, linear interpolation was first applied to sparse-view projections and the projections were rearranged into sinograms; these sinograms were processed with a sinogram restoration model and then rearranged back into projections. The sinogram restoration model was modified from the 2D U-Net by incorporating dynamic convolutional layers and residual learning techniques. The DLMPS approach was trained, validated, and tested on CBCT data from 163, 30, and 30 real patients respectively. Sparse-view projection datasets with 1/4 and 1/8 of the original sampling rate were simulated, and the corresponding full-view projection datasets were restored via the DLMPS approach. Tomographic images were reconstructed using the Feldkamp-Davis-Kress algorithm. Quantitative metrics including root-mean-square error (RMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and feature similarity (FSIM) were calculated in both the projection and image domains to evaluate the performance of the DLMPS approach. The DLMPS approach was compared with 11 state-of-the-art (SOTA) models, including CNN and Transformer architectures. For 1/4 sparse-view reconstruction task, the proposed DLMPS approach achieved averaged RMSE, PSNR, SSIM, and FSIM values of 0.0271, 45.93 dB, 0.9817, and 0.9587 in the projection domain, and 0.000885, 37.63 dB, 0.9074, and 0.9885 in the image domain, respectively. For 1/8 sparse-view reconstruction task, the DLMPS approach achieved averaged RMSE, PSNR, SSIM, and FSIM values of 0.0304, 44.85 dB, 0.9785, and 0.9524 in the projection domain, and 0.001057, 36.05 dB, 0.8786, and 0.9774 in the image domain, respectively. The DLMPS approach outperformed all the 11 SOTA models in both the projection and image domains for 1/4 and 1/8 sparse-view reconstruction tasks. The proposed DLMPS approach effectively improves the quality of sparse-view CBCT images in IGRT by accurately synthesizing missing projections, exhibiting potential in substantially reducing imaging dose to patients with minimal loss of image quality.

Synthesis and Characterization of Gold-Nanoparticle-Loaded Block Copolymer Vectors for Biomedical Applications: A Multivariate Analysis

Gold nanoparticles (GNPs) encapsulated in amphiphilic block copolymers are a promising system for numerous biomedical applications, although critical information on the effects of various preparation variables on the structure and properties of this unique type of nanomaterial is currently missing from the literature. In this research, we synthesized GNPs functionalized with thiol-terminated polycaprolactone (PCL-GNPs) before encapsulating them into poly(ε-caprolactone)-b-poly(ethylene glycol) (PCL-b-PEG) micellar nanoparticles via nanoprecipitation to yield GNP-loaded polymeric nanoparticles (GNP-PNPs). We explored the role of different manufacturing variables (water volume, PCL-b-PEG to PCL-GNP ratio, and PEG block length) on the sizes, morphologies, GNP occupancies, colloidal gold concentrations, and time stability of GNP-PNPs. Despite our motivation to increase colloidal gold concentrations for K-edge CT imaging applications, there was only moderate variation in the concentration of colloidal gold (cAu = ∼100-150 μg/mL) over the range of investigated experimental variables; however, postformulation exposure to compressed air flow provided samples with increased gold concentrations and CT contrast above the visual threshold in imaging phantoms. The range of formulation variables also had only a weak effect on mean effective hydrodynamic diameters (dh,eff = ∼150 nm). Statistical analysis of TEM images revealed that the mean number of GNPs within GNP-PNP vectors smaller than 50 nm, Zave,d<50 nm, is generally higher for preparations involving PCL-b-PEG with the shorter of the two different PEG block lengths. Preparations with the shorter PEG block length copolymer were also found to produce GNP-PNP colloids with greater time stability in dh,eff and cAu. Consistent with our previous study using MB-MDA-231 cells, we found increased gold uptake in MCF-7 cells with increasing Zave,d<50 nm. This study provides a roadmap for optimizing important figures of merit for existing biomedical applications, including CT imaging and radiotherapy sensitization, and for developing new diagnostic and therapeutic strategies using GNP-PNPs.

Single-View Fluoroscopic X-Ray Pose Estimation: A Comparison of Alternative Loss Functions and Volumetric Scene Representations

Many tasks performed in image-guided procedures can be cast as pose estimation problems, where specific projections are chosen to reach a target in 3D space. In this study, we construct a framework for fluoroscopic pose estimation and compare alternative loss functions and volumetric scene representations. We first develop a differentiable projection (DiffProj) algorithm for the efficient computation of Digitally Reconstructed Radiographs (DRRs) from either Cone-Beam Computerized Tomography (CBCT) or neural scene representations. We introduce two innovative neural scene representations, Neural Tuned Tomography (NeTT) and masked Neural Radiance Fields (mNeRF). Pose estimation is then performed within the framework by iterative gradient descent using loss functions that quantify the image discrepancy of the synthesized DRR with respect to the ground-truth, target fluoroscopic X-ray image. We compared alternative loss functions and volumetric scene representations for pose estimation using a dataset consisting of 50 cranial tomographic X-ray sequences. We find that Mutual Information significantly outperforms alternative loss functions for pose estimation, avoiding entrapment in local optima. The alternative discrete (CBCT) and neural (NeTT and mNeRF) volumetric scene representations yield comparable performance (3D angle errors, mean ≤ 3.2° and 90% quantile ≤ 3.4°); however, the neural scene representations incur a considerable computational expense to train.

Fast computational approach with prior dimension reduction for three-dimensional chemical component analysis using CT data of spectral imaging

Spectral image (SI) measurement techniques, such as X-ray absorption fine structure (XAFS) imaging and scanning transmission electron microscopy (STEM) with energy-dispersive X-ray spectroscopy (EDS) or electron energy loss spectroscopy (EELS), are useful for identifying chemical structures in composite materials. Machine-learning techniques have been developed for automatic analysis of SI data and their usefulness has been proven. Recently, an extended measurement technique combining SI with a computed tomography (CT) technique (CT-SI), such as CT-XAFS and STEM-EDS/EELS tomography, was developed to identify the three-dimensional (3D) structures of chemical components. CT-SI analysis can be conducted by combining CT reconstruction algorithms and chemical component analysis based on machine-learning techniques. However, this analysis incurs high-computational costs owing to the size of the CT-SI datasets. To address this problem, this study proposed a fast computational approach for 3D chemical component analysis in an unsupervised learning setting. The primary idea for reducing the computational cost involved compressing the CT-SI data prior to CT computation and performing 3D reconstruction and chemical component analysis on the compressed data. The proposed approach significantly reduced the computational cost without losing information about the 3D structure and chemical components. We experimentally evaluated the proposed approach using synthetic and real CT-XAFS data, which demonstrated that our approach achieved a significantly faster computational speed than the conventional approach while maintaining analysis performance. As the proposed procedure can be implemented with any CT algorithm, it is expected to accelerate 3D analyses with sparse regularized CT algorithms in noisy and sparse CT-SI datasets.

Machine learning-driven image synthesis and analysis applications for inertial confinement fusion (invited)

Recent fusion breakeven [Abu-Shawareb et al., Phys. Rev. Lett. 132, 065102 (2024)] in the National Ignition Facility (NIF) motivates an integrated approach to data analysis from multiple diagnostics. Deep neural networks provide a seamless framework for multi-modal data fusion, automated data analysis, optimization, and uncertainty quantification [Wang et al., arXiv:2401.08390 (2024)]. Here, we summarize different neural network methods for x-ray and neutron imaging data from NIF. To compensate for the small experimental datasets, both model based physics-informed synthetic data generation and deep neural network methods, such as generative adversarial networks, have been successfully implemented to allow a variety of automated workflows in x-ray and neutron image processing. We highlight results in noise emulation, contour analysis for low-mode analysis and asymmetry, denoising, and super-resolution. Further advances in the integrated multi-modal imaging, in sync with experimental validation and uncertainty quantification, will help with the ongoing experimental optimization in NIF, as well as the maturation of alternate inertial confinement fusion (ICF) platforms such as double-shells.

Source-detector trajectory optimization for FOV extension in dental CBCT imaging

In dental imaging, Cone Beam Computed Tomography (CBCT) is a widely used imaging modality for diagnosis and treatment planning. Small dental scanning units are the most popular due to their cost-effectiveness. However, these small systems have the limitation of a small field of view (FOV) as the source and detector move at a limited angle in a circular path. This often limits the FOV size. In this study, we addressed this issue by modifying the source-detector trajectory of the small dental device. The main goal of this study was to extend the FOV algorithmically by acquiring projection data with optimal projection angulation and isocenter location rather than upgrading any physical parts of the device. A novel algorithm to implement a Volume of Interest (VOI) guided trajectory is developed in this study based on the small dental imaging device's geometry. In addition, this algorithm is fused with a previously developed off-axis scanning method which uses an elliptical trajectory, to compensate for the existing constraints and to further extend the FOV. A comparison with standard circular trajectory is performed. The FOV of such a standard trajectory is a circle of 11 cm diameter in the axial plane. The proposed novel trajectory extends the FOV significantly and a maximum FOV of 19.5 cm is achieved with the Structural Similarity Index Measure (SSIM) score ranging between (≈98-99%) in different VOIs. The study results indicate that the proposed source-detector trajectory can extend dental imaging FOV and increase imaging performance, which ultimately results in more precise diagnosis and enhanced patient outcomes.

Evaluation of a compact cone beam CT concept with high image fidelity for point-of-care brain imaging

Cone beam computed tomography (CBCT) has potential advantages for developing portable, cost-effective point-of-care CT systems for intracranial imaging, such as early stroke diagnosis, hemorrhage detection, and intraoperative navigation. However, large volume imaging with flat panel detector based CBCT significantly increases the scattered radiation fluence which reduces its image quality and utility. To address these issues, a compact CBCT concept with enhanced image quality was investigated for intracranial imaging. The new system features a novel antiscatter collimator and data correction method to address the challenges in imaging large volumes with CBCT. A benchtop CBCT prototype was constructed. Imaging studies with anthropomorphic phantoms showed that soft tissue visualization, Hounsfield Unit (HU) accuracy, contrast, and spatial resolution increased significantly with the proposed CBCT concept, and they were comparable to the values measured in the gold standard multidetector-row CT (MDCT) images. Contrast-to-noise ratio (CNR) in CBCT images was within 12-31% of the CNR in MDCT images. These findings indicate that a compact CBCT system integrated with effective scatter suppression techniques may have increased utility in the context of brain imaging, and the proposed approach may enable the development of point-of-care CT systems for head imaging based on flat panel detector based CBCT technology.

A simulation framework for preclinical proton irradiation workflow

Objective.The integration of proton beamlines with x-ray imaging/irradiation platforms has opened up possibilities for image-guided Bragg peak irradiations in small animals. Such irradiations allow selective targeting of normal tissue substructures and tumours. However, their small size and location pose challenges in designing experiments. This work presents a simulation framework useful for optimizing beamlines, imaging protocols, and design of animal experiments. The usage of the framework is demonstrated, mainly focusing on the imaging part.Approach.The fastCAT toolkit was modified with Monte Carlo (MC)-calculated primary and scatter data of a small animal imager for the simulation of micro-CT scans. The simulated CT of a mini-calibration phantom from fastCAT was validated against a full MC TOPAS CT simulation. A realistic beam model of a preclinical proton facility was obtained from beam transport simulations to create irradiation plans in matRad. Simulated CT images of a digital mouse phantom were generated using single-energy CT (SECT) and dual-energy CT (DECT) protocols and their accuracy in proton stopping power ratio (SPR) estimation and their impact on calculated proton dose distributions in a mouse were evaluated.Main results.The CT numbers from fastCAT agree within 11 HU with TOPAS except for materials at the centre of the phantom. Discrepancies for central inserts are caused by beam hardening issues. The root mean square deviation in the SPR for the best SECT (90 kV/Cu) and DECT (50 kV/Al-90 kV/Al) protocols are 3.7% and 1.0%, respectively. Dose distributions calculated for SECT and DECT datasets revealed range shifts <0.1 mm, gamma pass rates (3%/0.1 mm) greater than 99%, and no substantial dosimetric differences for all structures. The outcomes suggest that SECT is sufficient for proton treatment planning in animals.Significance.The framework is a useful tool for the development of an optimized experimental configuration without using animals and beam time.

Feasibility study of photon-counting CT for material identification based on YSO/SiPM detector: A proof of concept

BACKGROUND

Current photon-counting computed tomography (CT) systems utilize semiconductor detectors, such as cadmium telluride (CdTe), cadmium zinc telluride (CZT), and silicon (Si), which convert x-ray photons directly into charge pulses. An alternative approach is indirect detection, which involves Yttrium Orthosilicate (YSO) scintillators coupled with silicon photomultipliers (SiPMs). This presents an attractive and cost-effective option due to its low cost, high detection efficiency, low dark count rate, and high sensor gain.

OBJECTIVE

This study aims to establish a comprehensive quantitative imaging framework for three-energy-bin proof-of-concept photon-counting CT based on YSO/SiPM detectors developed in our group using multi-voltage threshold (MVT) digitizers and assess the feasibility of this spectral CT for material identification.

METHODS

We developed a proof-of-concept YSO/SiPM-based benchtop spectral CT system and established a pipeline for three-energy-bin photon-counting CT projection-domain processing. The empirical A-table method was employed for basis material decomposition, and the quantitative imaging performance of the spectral CT system was assessed. This evaluation included the synthesis errors of virtual monoenergetic images, electron density images, effective atomic number images, and linear attenuation coefficient curves. The validity of employing A-table methods for material identification in three-energy-bin spectral CT was confirmed through both simulations and experimental studies.

RESULTS

In both noise-free and noisy simulations, the thickness estimation experiments and quantitative imaging results demonstrated high accuracy. In the thickness estimation experiment using the practical spectral CT system, the mean absolute error for the estimated thickness of the decomposed Al basis material was 0.014 ± 0.010 mm, with a mean relative error of 0.66% ± 0.42%. Similarly, for the decomposed polymethyl methacrylate (PMMA) basis material, the mean absolute error in thickness estimation was 0.064 ± 0.058 mm, with a mean relative error of 0.70% ± 0.38%. Additionally, employing the equivalent thickness of the basis material allowed for accurate synthesis of 70 keV virtual monoenergetic images (relative error 1.85% ± 1.26%), electron density (relative error 1.81% ± 0.97%), and effective atomic number (relative error 2.64% ± 1.26%) of the tested materials. In addition, the average synthesis error of the linear attenuation coefficient curves in the energy range from 40 to 150 keV was 1.89% ± 1.07%.

CONCLUSIONS

Both simulation and experimental results demonstrate the accurate generation of 70 keV virtual monoenergetic images, electron density, and effective atomic number images using the A-table method. Quantitative imaging results indicate that the YSO/SiPM-based photon-counting detector is capable of accurately reconstructing virtual monoenergetic images, electron density images, effective atomic number images, and linear attenuation coefficient curves, thereby achieving precise material identification.

An efficient dual-domain deep learning network for sparse-view CT reconstruction

BACKGROUND AND OBJECTIVE

We develop an efficient deep-learning based dual-domain reconstruction method for sparse-view CT reconstruction with small training parameters and comparable running time. We aim to investigate the model's capability and its clinical value by performing objective and subjective quality assessments using clinical CT projection data acquired on commercial scanners.

METHODS

We designed two lightweight networks, namely Sino-Net and Img-Net, to restore the projection and image signal from the DD-Net reconstructed images in the projection and image domains, respectively. The proposed network has small training parameters and comparable running time among dual-domain based reconstruction networks and is easy to train (end-to-end). We prospectively collected clinical thoraco-abdominal CT projection data acquired on a Siemens Biograph 128 Edge CT scanner to train and validate the proposed network. Further, we quantitatively evaluated the CT Hounsfield unit (HU) values on 21 organs and anatomic structures, such as the liver, aorta, and ribcage. We also analyzed the noise properties and compared the signal-to-noise ratio (SNR) and the contrast-to-noise ratio (CNR) of the reconstructed images. Besides, two radiologists conducted the subjective qualitative evaluation including the confidence and conspicuity of anatomic structures, and the overall image quality using a 1-5 likert scoring system.

RESULTS

Objective and subjective evaluation showed that the proposed algorithm achieves competitive results in eliminating noise and artifacts, restoring fine structure details, and recovering edges and contours of anatomic structures using 384 views (1/6 sparse rate). The proposed method exhibited good computational cost performance on clinical projection data.

CONCLUSION

This work presents an efficient dual-domain learning network for sparse-view CT reconstruction on raw projection data from a commercial scanner. The study also provides insights for designing an organ-based image quality assessment pipeline for sparse-view reconstruction tasks, potentially benefiting organ-specific dose reduction by sparse-view imaging.

Source-detector trajectory optimization for CBCT metal artifact reduction based on PICCS reconstruction

Precise instrument placement plays a critical role in all interventional procedures, especially percutaneous procedures such as needle biopsies, to achieve successful tumor targeting and increased diagnostic accuracy. C-arm cone beam computed tomography (CBCT) has the potential to precisely visualize the anatomy in direct vicinity of the needle and evaluate the adequacy of needle placement during the intervention, allowing for instantaneous adjustment in case of misplacement. However, even with the most advanced C-arm CBCT devices, it can be difficult to identify the exact needle position on CBCT images due to the strong metal artifacts around the needle. In this study, we proposed a framework for customized trajectory design in CBCT imaging based on Prior Image Constrained Compressed Sensing (PICCS) reconstruction with the goal of reducing metal artifacts in needle-based procedures. We proposed to optimize out-of-plane rotations in three-dimensional (3D) space and minimize projection views while reducing metal artifacts at specific volume of interests (VOIs). An anthropomorphic thorax phantom with a needle inserted inside and two tumor models as the imaging targets were used to validate the proposed approach. The performance of the proposed approach was also evaluated for CBCT imaging under kinematic constraints by simulating some collision areas on the geometry of the C-arm. We compared the result of optimized 3D trajectories using the PICCS algorithm and 20 projections with the result of a circular trajectory with sparse view using PICCS and Feldkamp, Davis, and Kress (FDK), both using 20 projections, and the circular FDK method with 313 projections. For imaging targets 1 and 2, the highest values of structural similarity index measure (SSIM) and universal quality index (UQI) between the reconstructed image from the optimized trajectories and the initial CBCT image at the VOI was calculated 0.7521, 0.7308 and 0.7308, 0.7248 respectively. These results significantly outperformed the FDK method (with 20 and 313 projections) and the PICCS method (20 projections) both using the circular trajectory. Our results showed that the proposed optimized trajectories not only significantly reduce metal artifacts but also suggest a dose reduction for needle-based CBCT interventions, considering the small number of projections used. Furthermore, our results showed that the optimized trajectories are compatible with spatially constrained situations and enable CBCT imaging under kinematic constraints when the standard circular trajectory is not feasible.

Simultaneous iodine and barium imaging with photon-counting CT

Objective.The objective of this study is to explore the capabilities of photon-counting computed tomography (PCCT) in simultaneously imaging and differentiating materials with close atomic numbers, specifically barium (Z= 56) and iodine (Z= 53), which is challenging for conventional computed tomography (CT).Approach.Experiments were conducted using a bench-top PCCT system equipped with a cadmium zinc telluride detector. Various phantom setups and contrast agent concentrations (1%-5%) were employed, along with a biological sample. Energy thresholds were tuned to the K-edge absorption energies of barium (37.4 keV) and iodine (33.2 keV) to capture multi-energy CT images. K-edge decomposition was performed using K-edge subtraction and principal component analysis (PCA) techniques to differentiate and quantify the contrast agents.Main results.The PCCT system successfully differentiated and accurately quantified barium and iodine in both phantom combinations and a biological sample, achieving high correlations (R2≈1) between true and reconstructed concentrations. PCA outperformed K-edge subtraction, particularly in the presence of calcium, by providing superior differentiation between barium and iodine.Significance.This study demonstrates the potential of PCCT for reliable, detailed imaging in both clinical and research settings, particularly for contrast agents with similar atomic numbers. The results suggest that PCCT could offer significant improvements in imaging quality over conventional CT, especially in applications requiring precise material differentiation.

The Adaption of Recent New Concepts in Neural Radiance Fields and Their Role for High-Fidelity Volume Reconstruction in Medical Images

Volume reconstruction techniques are gaining increasing interest in medical domains due to their potential to learn complex 3D structural information from sparse 2D images. Recently, neural radiance fields (NeRF), which implicitly model continuous radiance fields based on multi-layer perceptrons to enable volume reconstruction of objects at arbitrary resolution, have gained traction in natural image volume reconstruction. However, the direct application of NeRF to medical volume reconstruction presents unique challenges due to differences in imaging principles, internal structure requirements, and boundary delineation. In this paper, we evaluate different NeRF techniques developed for natural images, including sampling strategies, feature encoding, and the use of complimentary features, by applying them to medical images. We evaluate three state-of-the-art NeRF techniques on four datasets of medical images of different complexity. Our goal is to identify the strengths, limitations, and future directions for integrating NeRF into the medical domain.

Investigations into the Geometric Calibration and Systematic Effects of a Micro-CT System

Micro-Computed Tomography (µCT) systems are used for examining the internal structures of various objects, such as material samples, manufactured parts, and natural objects. Resolving fine details or performing accurate geometric measurements in the voxel data critically depends on the precise calibration of the µCT systems geometry. This paper presents a calibration method for µCT systems using projections of a calibration phantom, where the coordinates of the phantom are initially unknown. The approach involves detecting and tracking steel ball bearings and adjusting the unknown system geometry parameters using non-linear least squares optimization. Multiple geometric models are tested to verify their suitability for a self-calibration approach. The implementation is tested using a calibration phantom captured at different magnifications. The results demonstrate the system's capability to determine the geometry model parameters with a remaining error on the detector between 0.27 px and 0.18 px. Systematic errors that remain after calibration, as well as changing parameters due to system instabilities, are investigated. The source code of this work is published to enable further research.

A Novel and Reliable Pixel Response Correction Method (DAC-Shifting) for Spectral Photon-Counting CT Imaging

Spectral photon-counting cone-beam computed tomography (CT) imaging is challenged by individual pixel response behaviours, which lead to noisy projection images and subsequent image artefacts like rings. Existing methods to correct for this either use calibration measurements, like signal-to-thickness calibration (STC), or perform a post-processing ring artefact correction of sinogram data or scan reconstructions without taking the pixel response explicitly into account. Here, we present a novel post-processing method (digital-to-analogue converter (DAC)-shifting) which explicitly measures the current pixel response using flat-field images and subsequently corrects the projection data. The DAC-shifting method was evaluated using a repeat series of the spectral photon-counting imaging (Medipix3) of a phantom with different density inserts and iodine K-edge imaging. The method was also compared against polymethyl methacrylate (PMMA)-based STC. The DAC-shifting method was shown to be effective in correcting individual pixel responses and was robust against detector instability; it led to a 47.4% average reduction in CT-number variation in homogeneous materials, with a range of 40.7-55.6%. On the contrary, the STC correction showed varying results; a 13.7% average reduction in CT-number variation, ranging from a 43.7% increase to a 45.5% reduction. In K-edge imaging, DAC-shifting provides a sharper attenuation peak and more uniform CT values, which are expected to benefit iodine concentration quantifications.

Evaluating spectral performance for quantitative contrast-enhanced breast CT with a GaAs based photon counting detector: a simulation approach

Quantitative contrast-enhanced breast computed tomography (CT) has the potential to improve the diagnosis and management of breast cancer. Traditional CT methods using energy-integrated detectors and dual-exposure images with different incident spectra for material discrimination can increase patient radiation dose and be susceptible to motion artifacts and spectral resolution loss. Photon Counting Detectors (PCDs) offer a promising alternative approach, enabling acquisition of multiple energy levels in a single exposure and potentially better energy resolution. Gallium arsenide (GaAs) is particularly promising for breast PCD-CT due to its high quantum efficiency and reduction of fluorescence x-rays escaping the pixel within the breast imaging energy range. In this study, the spectral performance of a GaAs PCD for quantitative iodine contrast-enhanced breast CT was evaluated. A GaAs detector with a pixel size of 100μm, a thickness of 500μm was simulated. Simulations were performed using cylindrical phantoms of varying diameters (10 cm, 12 cm, and 16 cm) with different concentrations and locations of iodine inserts, using incident spectra of 50, 55, and 60 kVp with 2 mm of added aluminum filtration and and a mean glandular dose of 10 mGy. We accounted for the effects of beam hardening and energy detector response using TIGRE CT open-source software and the publicly available Photon Counting Toolkit (PcTK). Material-specific images of the breast phantom were produced using both projection and image-based material decomposition methods, and iodine component images were used to estimate iodine intake. Accuracy and precision of the proposed methods for estimating iodine concentration in breast CT images were assessed for different material decomposition methods, incident spectra, and breast phantom thicknesses. The results showed that both the beam hardening effect and imperfection in the detector response had a significant impact on performance in terms of Root Mean Squared Error (RMSE), precision, and accuracy of estimating iodine intake in the breast. Furthermore, the study demonstrated the effectiveness of both material decomposition methods in making accurate and precise iodine concentration predictions using a GaAs-based photon counting breast CT system, with better performance when applying the projection-based material decomposition approach. The study highlights the potential of GaAs-based photon counting breast CT systems as viable alternatives to traditional imaging methods in terms of material decomposition and iodine concentration estimation, and proposes phantoms and figures of merit to assess their performance.

Iterative image reconstruction algorithm analysis for optical CT radiochromic gel dosimetry

Background.Modern radiation therapy technologies aim to enhance radiation dose precision to the tumor and utilize hypofractionated treatment regimens. Verifying the dose distributions associated with these advanced radiation therapy treatments remains an active research area due to the complexity of delivery systems and the lack of suitable three-dimensional dosimetry tools. Gel dosimeters are a potential tool for measuring these complex dose distributions. A prototype tabletop solid-tank fan-beam optical CT scanner for readout of gel dosimeters was recently developed. This scanner does not have a straight raypath from source to detector, thus images cannot be reconstructed using filtered backprojection (FBP) and iterative techniques are required.Purpose.To compare a subset of the top performing algorithms in terms of image quality and quantitatively determine the optimal algorithm while accounting for refraction within the optical CT system. The following algorithms were compared: Landweber, superiorized Landweber with the fast gradient projection perturbation routine (S-LAND-FGP), the fast iterative shrinkage/thresholding algorithm with total variation penalty term (FISTA-TV), a monotone version of FISTA-TV (MFISTA-TV), superiorized conjugate gradient with the nonascending perturbation routine (S-CG-NA), superiorized conjugate gradient with the fast gradient projection perturbation routine (S-CG-FGP), superiorized conjugate gradient with with two iterations of CG performed on the current iterate and the nonascending perturbation routine (S-CG-2-NA).Methods.A ray tracing simulator was developed to track the path of light rays as they traverse the different mediums of the optical CT scanner. Two clinical phantoms and several synthetic phantoms were produced and used to evaluate the reconstruction techniques under known conditions. Reconstructed images were analyzed in terms of spatial resolution, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), signal non-uniformity (SNU), mean relative difference (MRD) and reconstruction time. We developed an image quality based method to find the optimal stopping iteration window for each algorithm. Imaging data from the prototype optical CT scanner was reconstructed and analysed to determine the optimal algorithm for this application.Results.The optimal algorithms found through the quantitative scoring metric were FISTA-TV and S-CG-2-NA. MFISTA-TV was found to behave almost identically to FISTA-TV however MFISTA-TV was unable to resolve some of the synthetic phantoms. S-CG-NA showed extreme fluctuations in the SNR and CNR values. S-CG-FGP had large fluctuations in the SNR and CNR values and the algorithm has less noise reduction than FISTA-TV and worse spatial resolution than S-CG-2-NA. S-LAND-FGP had many of the same characteristics as FISTA-TV; high noise reduction and stability from over iterating. However, S-LAND-FGP has worse SNR, CNR and SNU values as well as longer reconstruction time. S-CG-2-NA has superior spatial resolution to all algorithms while still maintaining good noise reduction and is uniquely stable from over iterating.Conclusions.Both optimal algorithms (FISTA-TV and S-CG-2-NA) are stable from over iterating and have excellent edge detection with ESF MTF 50% values of 1.266 mm-1and 0.992 mm-1. FISTA-TV had the greatest noise reduction with SNR, CNR and SNU values of 424, 434 and 0.91 × 10-4, respectively. However, low spatial resolution makes FISTA-TV only viable for large field dosimetry. S-CG-2-NA has better spatial resolution than FISTA-TV with PSF and LSF MTF 50% values of 1.581 mm-1and 0.738 mm-1, but less noise reduction. S-CG-2-NA still maintains good SNR, CNR, and SNU values of 168, 158 and 1.13 × 10-4, respectively. Thus, S-CG-2-NA is a well rounded reconstruction algorithm that would be the preferable choice for small field dosimetry.

XTransCT: ultra-fast volumetric CT reconstruction using two orthogonal x-ray projections for image-guided radiation therapy via a transformer network

Objective.The aim of this study was to reconstruct volumetric computed tomography (CT) images in real-time from ultra-sparse two-dimensional x-ray projections, facilitating easier navigation and positioning during image-guided radiation therapy.Approach.Our approach leverages a voxel-sapce-searching Transformer model to overcome the limitations of conventional CT reconstruction techniques, which require extensive x-ray projections and lead to high radiation doses and equipment constraints.Main results.The proposed XTransCT algorithm demonstrated superior performance in terms of image quality, structural accuracy, and generalizability across different datasets, including a hospital set of 50 patients, the large-scale public LIDC-IDRI dataset, and the LNDb dataset for cross-validation. Notably, the algorithm achieved an approximately 300% improvement in reconstruction speed, with a rate of 44 ms per 3D image reconstruction compared to former 3D convolution-based methods.Significance.The XTransCT architecture has the potential to impact clinical practice by providing high-quality CT images faster and with substantially reduced radiation exposure for patients. The model's generalizability suggests it has the potential applicable in various healthcare settings.

2D antiscatter grid and scatter sampling based CBCT method for online dose calculations during CBCT guided radiation therapy of pelvis

BACKGROUND

Online dose calculations before the delivery of radiation treatments have applications in dose delivery verification, online adaptation of treatment plans, and simulation-free treatment planning. While dose calculations by directly utilizing CBCT images are desired, dosimetric accuracy can be compromised due to relatively lower HU accuracy in CBCT images.

PURPOSE

In this work, we propose a novel CBCT imaging pipeline to enhance the accuracy of CBCT-based dose calculations in the pelvis region. Our approach aims to improve the HU accuracy in CBCT images, thereby improving the overall accuracy of CBCT-based dose calculations prior to radiation treatment delivery.

METHODS

An in-house developed quantitative CBCT pipeline was implemented to address the CBCT raw data contamination problem. The pipeline combines algorithmic data correction strategies and 2D antiscatter grid-based scatter rejection to achieve high CT number accuracy. To evaluate the effect of the quantitative CBCT pipeline on CBCT-based dose calculations, phantoms mimicking pelvis anatomy were scanned using a linac-mounted CBCT system, and a gold standard multidetector CT used for treatment planning (pCT). A total of 20 intensity-modulated treatment plans were generated for five targets, using 6 and 10 MV flattening filter-free beams, and utilizing small and large pelvis phantom images. For each treatment plan, four different dose calculations were performed using pCT images and three CBCT imaging configurations: quantitative CBCT, clinical CBCT protocol, and a high-performance 1D antiscatter grid (1D ASG). Subsequently, dosimetric accuracy was evaluated for both targets and organs at risk as a function of patient size, target location, beam energy, and CBCT imaging configuration.

RESULTS

When compared to the gold-standard pCT, dosimetric errors in quantitative CBCT-based dose calculations were not significant across all phantom sizes, beam energies, and treatment sites. The largest error observed was 0.6% among all dose volume histogram metrics and evaluated dose calculations. In contrast, dosimetric errors reached up to 7% and 97% in clinical CBCT and high-performance ASG CBCT-based treatment plans, respectively. The largest dosimetric errors were observed in bony targets in the large phantom treated with 6 MV beams. The trends of dosimetric errors in organs at risk were similar to those observed in the targets.

CONCLUSIONS

The proposed quantitative CBCT pipeline has the potential to provide comparable dose calculation accuracy to the gold-standard planning CT in photon radiation therapy for the abdomen and pelvis. These robust dose calculations could eliminate the need for density overrides in CBCT images and enable direct utilization of CBCT images for dose delivery monitoring or online treatment plan adaptations before the delivery of radiation treatments.

Fully automated sinogram-based deep learning model for detection and classification of intracranial hemorrhage

PURPOSE

To propose an automated approach for detecting and classifying Intracranial Hemorrhages (ICH) directly from sinograms using a deep learning framework. This method is proposed to overcome the limitations of the conventional diagnosis by eliminating the time-consuming reconstruction step and minimizing the potential noise and artifacts that can occur during the Computed Tomography (CT) reconstruction process.

METHODS

This study proposes a two-stage automated approach for detecting and classifying ICH from sinograms using a deep learning framework. The first stage of the framework is Intensity Transformed Sinogram Sythesizer, which synthesizes sinograms that are equivalent to the intensity-transformed CT images. The second stage comprises of a cascaded Convolutional Neural Network-Recurrent Neural Network (CNN-RNN) model that detects and classifies hemorrhages from the synthesized sinograms. The CNN module extracts high-level features from each input sinogram, while the RNN module provides spatial correlation of the neighborhood regions in the sinograms. The proposed method was evaluated on a publicly available RSNA dataset consisting of a large sample size of 8652 patients.

RESULTS

The results showed that the proposed method had a notable improvement as high as 27% in patient-wise accuracies when compared to state-of-the-art methods like ResNext-101, Inception-v3 and Vision Transformer. Furthermore, the sinogram-based approach was found to be more robust to noise and offset errors in comparison to CT image-based approaches. The proposed model was also subjected to a multi-label classification analysis to determine the hemorrhage type from a given sinogram. The learning patterns of the proposed model were also examined for explainability using the activation maps.

CONCLUSION

The proposed sinogram-based approach can provide an accurate and efficient diagnosis of ICH without the need for the time-consuming reconstruction step and can potentially overcome the limitations of CT image-based approaches. The results show promising outcomes for the use of sinogram-based approaches in detecting hemorrhages, and further research can explore the potential of this approach in clinical settings.

Technical note: Preprocessing of portal images to improve image quality of VMAT-CT

BACKGROUND

The concept of volumetric modulated arc therapy-computed tomography (VMAT-CT) was proposed more than a decade ago. However, its application has been very limited mainly due to the poor image quality. More specifically, the blurred areas in electronic portal imaging device (EPID) images collected during VMAT heavily degrade the image quality of VMAT-CT.

PURPOSE

The goal of this study was to propose systematic methods to preprocess EPID images and improve the image quality of VMAT-CT.

METHODS

Online region-based active contour method was introduced to binarize portal images. Multi-leaf collimator (MLC) motion modeling was developed to remove the MLC motion blur. Outlier filtering was then applied to replace the remaining artifacts with plausible data. To assess the impact of these preprocessing methods on the image quality of VMAT-CT, 44 clinical VMAT plans for several treatment sites (lung, esophagus, and head & neck) were delivered to a Rando phantom, and several real-patient cases were also acquired. VMAT-CT reconstruction was attempted for all the cases, and image quality was evaluated.

RESULTS

All three preprocessing methods could effectively remove the blurred edges of EPID images. The combined preprocessing methods not only saved VMAT-CT from distortions and artifacts, but also increased the percentage of VMAT plans that can be reconstructed.

CONCLUSIONS

The systematic preprocessing of portal images improves the image quality of VMAT-CT significantly, and facilitates the application of VMAT-CT as an effective image guidance tool.

Material decomposition with a prototype photon-counting detector CT system: expanding a stoichiometric dual-energy CT method via energy bin optimization and K-edge imaging

Objective.Computed tomography (CT) has advanced since its inception, with breakthroughs such as dual-energy CT (DECT), which extracts additional information by acquiring two sets of data at different energies. As high-flux photon-counting detectors (PCDs) become available, PCD-CT is also becoming a reality. PCD-CT can acquire multi-energy data sets in a single scan by spectrally binning the incident x-ray beam. With this, K-edge imaging becomes possible, allowing high atomic number (high-Z) contrast materials to be distinguished and quantified. In this study, we demonstrated that DECT methods can be converted to PCD-CT systems by extending the method of Bourqueet al(2014). We optimized the energy bins of the PCD for this purpose and expanded the capabilities by employing K-edge subtraction imaging to separate a high-atomic number contrast material.Approach.The method decomposes materials into their effective atomic number (Zeff) and electron density relative to water (ρe). The model was calibrated and evaluated using tissue-equivalent materials from the RMI Gammex electron density phantom with knownρevalues and elemental compositions. TheoreticalZeffvalues were found for the appropriate energy ranges using the elemental composition of the materials.Zeffvaried slightly with energy but was considered a systematic error. Anex vivobovine tissue sample was decomposed to evaluate the model further and was injected with gold chloride to demonstrate the separation of a K-edge contrast agent.Main results.The mean root mean squared percent errors on the extractedZeffandρefor PCD-CT were 0.76% and 0.72%, respectively and 1.77% and 1.98% for DECT. The tissue types in theex vivobovine tissue sample were also correctly identified after decomposition. Additionally, gold chloride was separated from theex vivotissue sample with K-edge imaging.Significance.PCD-CT offers the ability to employ DECT material decomposition methods, along with providing additional capabilities such as K-edge imaging.

Improving the accuracy of bone mineral density using a multisource CBCT

Multisource cone beam computed tomography CBCT (ms-CBCT) has been shown to overcome some of the inherent limitations of a conventional CBCT. The purpose of this study was to evaluate the accuracy of ms-CBCT for measuring the bone mineral density (BMD) of mandible and maxilla compared to the conventional CBCT. The values measured from a multi-detector CT (MDCT) were used as substitutes for the ground truth. An anthropomorphic adult skull and tissue equivalent head phantom and a homemade calibration phantom containing inserts with varying densities of calcium hydroxyapatite were imaged using the ms-CBCT, the ms-CBCT operating in the conventional single source CBCT mode, and two clinical CBCT scanners at similar imaging doses; and a clinical MDCT. The images of the anthropomorphic head phantom were reconstructed and registered, and the cortical and cancellous bones of the mandible and the maxilla were segmented. The measured CT Hounsfield Unit (HU) and Greyscale Value (GV) at multiple region-of-interests were converted to the BMD using scanner-specific calibration functions. The results from the various CBCT scanners were compared to that from the MDCT. Statistical analysis showed a significant improvement in the agreement between the ms-CBCT and MDCT compared to that between the CBCT and MDCT.

Projection-Angle-Sensor-Assisted X-ray Computed Tomography for Cylindrical Lithium-Ion Batteries

X-ray computed tomography (XCT) has become a powerful technique for studying lithium-ion batteries, allowing non-destructive 3D imaging across multiple spatial scales. Image quality is particularly important for observing the internal structure of lithium-ion batteries. During multiple rotations, the existence of cumulative errors and random errors in the rotary table leads to errors in the projection angle, affecting the imaging quality of XCT. The accuracy of the projection angle is an important factor that directly affects imaging. However, the impact of the projection angle on XCT reconstruction imaging is difficult to quantify. Therefore, the required precision of the projection angle sensor cannot be determined explicitly. In this research, we selected a common 18650 cylindrical lithium-ion battery for experiments. By setting up an XCT scanning platform and installing an angle sensor to calibrate the projection angle, we proceeded with image reconstruction after introducing various angle errors. When comparing the results, we found that projection angle errors lead to the appearance of noise and many stripe artifacts in the image. This is particularly noticeable in the form of many irregular artifacts in the image background. The overall variation and residual projection error in detection indicators can effectively reflect the trend in image quality. This research analyzed the impact of projection angle errors on imaging and improved the quality of XCT imaging by installing angle sensors on a rotary table.

Iterative Reconstruction Algorithms in Magneto-Acousto-Electrical Computed Tomography (MAE-CT) for Image Quality Improvement

Magneto-acousto-electrical computed tomography (MAE-CT) is a recently developed rotational magneto-acousto-electrical tomography (MAET) method, which can map the conductivity parameter of tissues with high spatial resolution. Since the imaging mode of MAE-CT is similar to that of CT, the reconstruction algorithms for CT are possible to be adopted for MAE-CT. Previous studies have demonstrated that the filtered back-projection (FBP) algorithm, which is one of the most common CT reconstruction algorithms, can be used for MAE-CT reconstruction. However, FBP has some inherent shortcomings of being sensitive to noise and non-uniform distribution of views. In this study, we introduced iterative reconstruction (IR) method in MAE-CT reconstruction and compared its performance with that of the FBP. The numerical simulation, the phantom, and in vitro experiments were performed, and several IR algorithms (ART, SART, SIRT) were used for reconstruction. The results show that the images reconstructed by the FBP and IR are similar when the data is noise-free in the simulation. As the noise level increases, the images reconstructed by SART and SIRT are more robust to the noise than FBP. In the phantom experiment, noise and some stripe artifacts caused by the FBP are removed by SART and SIRT algorithms. In conclusion, the IR method used in CT is applicable in MAE-CT, and it performs better than FBP, which indicates that the state-of-the-art achievements in the CT algorithm can also be adopted for the MAE-CT reconstruction in the future.

Non-coplanar lung SABR treatments delivered with a gantry-mounted x-ray tube

Objective.To create two non-coplanar, stereotactic ablative radiotherapy (SABR) lung patient treatment plans compliant with the radiation therapy oncology group (RTOG) 0813 dosimetric criteria using a simple, isocentric, therapy with kilovoltage arcs (SITKA) system designed to provide low cost external radiotherapy treatments for low- and middle-income countries (LMICs).Approach.A treatment machine design has been proposed featuring a 320 kVp x-ray tube mounted on a gantry. A deep learning cone-beam CT (CBCT) to synthetic CT (sCT) method was employed to remove the additional cost of planning CTs. A novel inverse treatment planning approach using GPU backprojection was used to create a highly non-coplanar treatment plan with circular beam shapes generated by an iris collimator. Treatments were planned and simulated using the TOPAS Monte Carlo (MC) code for two lung patients. Dose distributions were compared to 6 MV volumetric modulated arc therapy (VMAT) planned in Eclipse on the same cases for a Truebeam linac as well as obeying the RTOG 0813 protocols for lung SABR treatments with a prescribed dose of 50 Gy.Main results.The low-cost SITKA treatments were compliant with all RTOG 0813 dosimetric criteria. SITKA treatments showed, on average, a 6.7 and 4.9 Gy reduction of the maximum dose in soft tissue organs at risk (OARs) as compared to VMAT, for the two patients respectively. This was accompanied by a small increase in the mean dose of 0.17 and 0.30 Gy in soft tissue OARs.Significance.The proposed SITKA system offers a maximally low-cost, effective alternative to conventional radiotherapy systems for lung cancer patients, particularly in low-income countries. The system's non-coplanar, isocentric approach, coupled with the deep learning CBCT to sCT and GPU backprojection-based inverse treatment planning, offers lower maximum doses in OARs and comparable conformity to VMAT plans at a fraction of the cost of conventional radiotherapy.

An attenuation field network for dedicated cone beam breast CT with short scan and offset detector geometry

The feasibility of full-scan, offset-detector geometry cone-beam CT has been demonstrated for several clinical applications. For full-scan acquisition with offset-detector geometry, data redundancy from complementary views can be exploited during image reconstruction. Envisioning an upright breast CT system, we propose to acquire short-scan data in conjunction with offset-detector geometry. To tackle the resulting incomplete data, we have developed a self-supervised attenuation field network (AFN). AFN leverages the inherent redundancy of cone-beam CT data through coordinate-based representation and known imaging physics. A trained AFN can query attenuation coefficients using their respective coordinates or synthesize projection data including the missing projections. The AFN was evaluated using clinical cone-beam breast CT datasets (n = 50). While conventional analytical and iterative reconstruction methods failed to reconstruct the incomplete data, AFN reconstruction was not statistically different from the reference reconstruction obtained using full-scan, full-detector data in terms of image noise, image contrast, and the full width at half maximum of calcifications. This study indicates the feasibility of a simultaneous short-scan and offset-detector geometry for dedicated breast CT imaging. The proposed AFN technique can potentially be expanded to other cone-beam CT applications.

tomoCAM: fast model-based iterative reconstruction via GPU acceleration and non-uniform fast Fourier transforms

X-ray-based computed tomography is a well established technique for determining the three-dimensional structure of an object from its two-dimensional projections. In the past few decades, there have been significant advancements in the brightness and detector technology of tomography instruments at synchrotron sources. These advancements have led to the emergence of new observations and discoveries, with improved capabilities such as faster frame rates, larger fields of view, higher resolution and higher dimensionality. These advancements have enabled the material science community to expand the scope of tomographic measurements towards increasingly in situ and in operando measurements. In these new experiments, samples can be rapidly evolving, have complex geometries and restrictions on the field of view, limiting the number of projections that can be collected. In such cases, standard filtered back-projection often results in poor quality reconstructions. Iterative reconstruction algorithms, such as model-based iterative reconstructions (MBIR), have demonstrated considerable success in producing high-quality reconstructions under such restrictions, but typically require high-performance computing resources with hundreds of compute nodes to solve the problem in a reasonable time. Here, tomoCAM, is introduced, a new GPU-accelerated implementation of model-based iterative reconstruction that leverages non-uniform fast Fourier transforms to efficiently compute Radon and back-projection operators and asynchronous memory transfers to maximize the throughput to the GPU memory. The resulting code is significantly faster than traditional MBIR codes and delivers the reconstructive improvement offered by MBIR with affordable computing time and resources. tomoCAM has a Python front-end, allowing access from Jupyter-based frameworks, providing straightforward integration into existing workflows at synchrotron facilities.

Feasibility of dual-energy CBCT material decomposition in the human torso with 2D anti-scatter grids and grid-based scatter sampling

BACKGROUND

Dual-energy (DE) imaging techniques in cone-beam computed tomography (CBCT) have potential clinical applications, including material quantification and improved tissue visualization. However, the performance of DE CBCT is limited by the effects of scattered radiation, which restricts its use to small object imaging.

PURPOSE

This study investigates the feasibility of DE CBCT material decomposition by reducing scatter with a 2D anti-scatter grid and a measurement-based scatter correction method. Specifically, the investigation focuses on iodine quantification accuracy and virtual monoenergetic (VME) imaging in phantoms that mimic head, thorax, abdomen, and pelvis anatomies.

METHODS

A 2D anti-scatter grid prototype was utilized with a residual scatter correction method in a linac-mounted CBCT system to investigate the effects of robust scatter suppression in DE CBCT. Scans were acquired at 90 and 140 kVp using phantoms that mimic head, thorax, and abdomen/pelvis anatomies. Iodine vials with varying concentrations were placed in each phantom, and CBCT images were decomposed into iodine and water basis material images. The effect of a 2D anti-scatter grid with and without residual scatter correction on iodine concentration quantification and contrast visualization in VME images was evaluated. To benchmark iodine concentration quantification accuracy, a similar set of experiments and DE processing were also performed with a conventional multidetector CT scanner.

RESULTS

In CBCT images, a 2D grid with or without scatter correction can differentiate iodine and water after DE processing in human torso-sized phantom images. However, iodine quantification errors were up to 10 mg/mL in pelvis phantoms when only the 2D grid was used. Adding scatter correction to 2D-grid CBCT reduced iodine quantification errors below 1.5 mg/mL in pelvis phantoms, comparable to iodine quantification errors in multidetector CT. While a noticeable contrast-to-noise ratio improvement was not observed in VME CBCT images, contrast visualization was substantially better in 40 keV VME images in visual comparisons with 90 and 140 kVp CBCT images across all phantom sizes investigated.

CONCLUSIONS

This study indicates that accurate DE decomposition is potentially feasible in DE CBCT of the human torso if robust scatter suppression is achieved with 2D anti-scatter grids and residual scatter correction. This approach can potentially enable better contrast visualization and tissue and contrast agent quantification in various CBCT applications.

Three-dimensional spine reconstruction from biplane radiographs using convolutional neural networks

PURPOSE

The purpose of this study was to develop and evaluate a deep learning network for three-dimensional reconstruction of the spine from biplanar radiographs.

METHODS

The proposed approach focused on extracting similar features and multiscale features of bone tissue in biplanar radiographs. Bone tissue features were reconstructed for feature representation across dimensions to generate three-dimensional volumes. The number of feature mappings was gradually reduced in the reconstruction to transform the high-dimensional features into the three-dimensional image domain. We produced and made eight public datasets to train and test the proposed network. Two evaluation metrics were proposed and combined with four classical evaluation metrics to measure the performance of the method.

RESULTS

In comparative experiments, the reconstruction results of this method achieved a Hausdorff distance of 1.85 mm, a surface overlap of 0.2 mm, a volume overlap of 0.9664, and an offset distance of only 0.21 mm from the vertebral body centroid. The results of this study indicate that the proposed method is reliable.

Exploration of differentiability in a proton computed tomography simulation framework

Objective.Gradient-based optimization using algorithmic derivatives can be a useful technique to improve engineering designs with respect to a computer-implemented objective function. Likewise, uncertainty quantification through computer simulations can be carried out by means of derivatives of the computer simulation. However, the effectiveness of these techniques depends on how 'well-linearizable' the software is. In this study, we assess how promising derivative information of a typical proton computed tomography (pCT) scan computer simulation is for the aforementioned applications.Approach.This study is mainly based on numerical experiments, in which we repeatedly evaluate three representative computational steps with perturbed input values. We support our observations with a review of the algorithmic steps and arithmetic operations performed by the software, using debugging techniques.Main results.The model-based iterative reconstruction (MBIR) subprocedure (at the end of the software pipeline) and the Monte Carlo (MC) simulation (at the beginning) were piecewise differentiable. However, the observed high density and magnitude of jumps was likely to preclude most meaningful uses of the derivatives. Jumps in the MBIR function arose from the discrete computation of the set of voxels intersected by a proton path, and could be reduced in magnitude by a 'fuzzy voxels' approach. The investigated jumps in the MC function arose from local changes in the control flow that affected the amount of consumed random numbers. The tracking algorithm solves an inherently non-differentiable problem.Significance.Besides the technical challenges of merely applying AD to existing software projects, the MC and MBIR codes must be adapted to compute smoother functions. For the MBIR code, we presented one possible approach for this while for the MC code, this will be subject to further research. For the tracking subprocedure, further research on surrogate models is necessary.

A quantitative CBCT pipeline based on 2D antiscatter grid and grid-based scatter sampling for image-guided radiation therapy

BACKGROUND

Quantitative accuracy is critical for expanding the role of cone beam CT (CBCT) imaging from target localization to quantitative treatment monitoring and plan adaptations in radiation therapy. Despite advances in CBCT image quality improvement methods, quantitative accuracy gap between CBCT and multi-detector CT (MDCT) remains.

PURPOSE

In this work, a physics-driven approach was investigated that combined robust scatter rejection, raw data correction and iterative image reconstruction to further improve CBCT image quality and quantitative accuracy, referred to as quantitative CBCT (qCBCT).

METHODS

QCBCT approach includes tungsten 2D antiscatter grid hardware, residual scatter correction with grid-based scatter sampling, image lag, and beam hardening correction for offset detector geometry linac-mounted CBCT. Images were reconstructed with iterative image reconstruction to reduce image noise. qCBCT was evaluated using a variety of phantoms to investigate the effect of object size and its composition on image quality, and image quality was benchmarked against clinical CBCT and gold standard MDCT images used for treatment planning.

RESULTS

QCBCT provided statistically significant improvement in CT number accuracy and reduced image artifacts when compared to clinical CBCT images. When compared to gold standard MDCT, mean HU errors in qCBCT and clinical CBCT were 17 ± 9 and 38 ± 29 HU, respectively. Magnitude of phantom size dependent HU variations were comparable between MDCT and qCBCT images. With iterative reconstruction, contrast-to-noise ratio improved by 25% when compared to clinical CBCT protocols.

CONCLUSIONS

Combination of novel scatter suppression techniques and other data correction methods in qCBCT provided CT number accuracy comparable to gold standard MDCT used for treatment planning. This approach may potentially improve CBCT's promise in fulfilling the tasks that demand high quantitative accuracy, such as online dose calculations and treatment response assessment, in image guided radiation therapy.

Potential for Dose Reduction in CT-Derived Left Ventricular Ejection Fraction: A Simulation Study

BACKGROUND

Measuring left ventricular ejection fraction (LVEF) is important for detecting heart failure, e.g., in treatment with potentially cardiotoxic chemotherapy. MRI is considered the reference standard for LVEF, but availability may be limited and claustrophobia or metal implants still present challenges. CT has been shown to be accurate and would be advantageous, as LVEF could be measured in conjunction with routine chest-abdomen-pelvis oncology CT. However, the use of CT is not recommended due to the excessive radiation dose. This study aimed to explore the potential for dose reduction using simulation. Using an anthropomorphic heart phantom scanned at 13 dose levels, a noise simulation algorithm was developed to introduce controlled Poisson noise. Filtered backprojection parameters were iteratively tested to minimise differences in myocardium-to-ventricle contrast/noise ratio, as well as structural similarity index (SSIM) differences between real and simulated images at all dose levels. Fifty-one clinical CT coronary angiographies, scanned with full dose through end-systolic and -diastolic phases, were located retrospectively. Using the developed algorithm, noise was introduced corresponding to 25, 10, 5 and 2% of the original dose level. LVEF was measured using clinical software (Syngo.via VB50) with papillary muscles in and excluded from the LV volume. At each dose level, LVEF was compared to the 100% dose level, using Bland-Altman analysis. The effective dose was calculated from DLP using a conversion factor of 0.026 mSv/mGycm.

RESULTS

In the clinical images, mean CTDIvol and DLP were 47.1 mGy and 771.9 mGycm, respectively (effective dose 20.0 mSv). Measurements with papillary muscles excluded did not exhibit statistically significant LVEF bias to full-dose images at 25, 10 and 5% simulated dose. At 2% dose, a significant bias of 4.4% was found. With papillary muscles included, small but significant biases were found at all simulated dose levels.

CONCLUSION

Provided that measurements are performed with papillary muscles excluded from the LV volume, the dose can be reduced by a factor of 20 without significantly affecting LVEF measurements. This corresponds to an effective dose of 1 mSv. CT can potentially be used for LVEF measurement with minimal excessive radiation.

Extension of the open-source TIGRE toolbox for proton imaging

Proton irradiation is a well-established method to treat deep-seated tumors in radio oncology. Usually, an X-ray computed tomography (CT) scan is used for treatment planning. Since proton therapy is based on the precise knowledge of the stopping power describing the energy loss of protons in the patient tissues, the Hounsfield units of the planning CT have to be converted. This conversion introduces range errors in the treatment plan, which could be reduced, if the stopping power values were extracted directly from an image obtained using protons instead of X-rays. Since protons are affected by multiple Coulomb scattering, reconstruction of the 3D stopping power map results in limited image quality if the curved proton path is not considered. This work presents a substantial code extension of the open-source toolbox TIGRE for proton CT (pCT) image reconstruction based on proton radiographs including a curved proton path estimate. The code extension and the reconstruction algorithms are GPU-based, allowing to achieve reconstruction results within minutes. The performance of the pCT code extension was tested with Monte Carlo simulated data using three phantoms (Catphan® high resolution and sensitometry modules and a CIRS patient phantom). In the simulations, ideal and non-ideal conditions for a pCT setup were assumed. The obtained mean absolute percentage error was found to be below 1% and up to 8 lp/cm could be resolved using an idealized setup. These findings demonstrate that the presented code extension to the TIGRE toolbox offers the possibility for other research groups to use a fast and accurate open-source pCT reconstruction.

Target-oriented deep learning-based image registration with individualized test-time adaptation

BACKGROUND

A classic approach in medical image registration is to formulate an optimization problem based on the image pair of interest, and seek a deformation vector field (DVF) to minimize the corresponding objective, often iteratively. It has a clear focus on the targeted pair, but is typically slow. In contrast, more recent deep-learning-based registration offers a much faster alternative and can benefit from data-driven regularization. However, learning is a process to "fit" the training cohort, whose image or motion characteristics or both may differ from the pair of images to be tested, which is the ultimate goal of registration. Therefore, generalization gap poses a high risk with direct inference alone.

PURPOSE

In this study, we propose an individualized adaptation to improve test sample targeting, to achieve a synergy of efficiency and performance in registration.

METHODS

Using a previously developed network with an integrated motion representation prior module as the implementation backbone, we propose to adapt the trained registration network further for image pairs at test time to optimize the individualized performance. The adaptation method was tested against various characteristics shifts caused by cross-protocol, cross-platform, and cross-modality, with test evaluation performed on lung CBCT, cardiac MRI, and lung MRI, respectively.

RESULTS

Landmark-based registration errors and motion-compensated image enhancement results demonstrated significantly improved test registration performance from our method, compared to tuned classic B-spline registration and network solutions without adaptation.

CONCLUSIONS

We have developed a method to synergistically combine the effectiveness of pre-trained deep network and the target-centric perspective of optimization-based registration to improve performance on individual test data.

4D-CT deformable image registration using unsupervised recursive cascaded full-resolution residual networks

A novel recursive cascaded full-resolution residual network (RCFRR-Net) for abdominal four-dimensional computed tomography (4D-CT) image registration was proposed. The entire network was end-to-end and trained in the unsupervised approach, which meant that the deformation vector field, which presented the ground truth, was not needed during training. The network was designed by cascading three full-resolution residual subnetworks with different architectures. The training loss consisted of the image similarity loss and the deformation vector field regularization loss, which were calculated based on the final warped image and the fixed image, allowing all cascades to be trained jointly and perform the progressive registration cooperatively. Extensive network testing was conducted using diverse datasets, including an internal 4D-CT dataset, a public DIRLAB 4D-CT dataset, and a 4D cone-beam CT (4D-CBCT) dataset. Compared with the iteration-based demon method and two deep learning-based methods (VoxelMorph and recursive cascaded network), the RCFRR-Net achieved consistent and significant gains, which demonstrated that the proposed method had superior performance and generalization capability in medical image registration. The proposed RCFRR-Net was a promising tool for various clinical applications.

Volumetric computed tomography with carbon nanotube X-ray source array for improved image quality and accuracy

Cone beam computed tomography (CBCT) is widely used in medical and dental imaging. Compared to a multidetector CT, it provides volumetric images with high isotropic resolution at a reduced radiation dose, cost and footprint without the need for patient translation. The current CBCT has several intrinsic limitations including reduced soft tissue contrast, inaccurate quantification of X-ray attenuation, image distortions and artefacts, which have limited its clinical applications primarily to imaging hard tissues and made quantitative analysis challenging. Here we report a multisource CBCT (ms-CBCT) which overcomes the short-comings of the conventional CBCT by using multiple narrowly collimated and rapidly scanning X-ray beams from a carbon nanotube field emission source array. Phantom imaging studies show that, the ms-CBCT increases the accuracy of the Hounsfield unit values by 60%, eliminates the cone beam artefacts, extends the axial coverage, and improves the soft tissue contrast-to-noise ratio by 30-50%, compared to the CBCT configuration.

A simulation study to evaluate the effect of 2D antiscatter grid primary transmission on flat panel detector based CBCT image quality

Purpose. Two-dimensional antiscatter grids' (2D-ASGs) septal shadows and their impact on primary transmission play a critical role in cone-beam computed tomography (CBCT) image noise and artifact characteristics. Therefore, a numerical simulation platform was developed to evaluate the effect of 2D-ASG's primary transmission on image quality, as a function of grid geometry and CBCT system properties.Methods. To study the effect of 2D-ASG's septal shadows on primary transmission and CBCT image quality, two new methods were introduced; one to simulate projection signal gradients in septal shadows, and the other to simulate septal shadow variations due to gantry flex. Signal gradients in septal shadows were simulated by generating a system point spread function that was directly extracted from projection images of 2D-ASG prototypes in experiments. Variations in septal shadows due to gantry flex were simulated by generating oversampled shadow profiles extracted from experiments. Subsequently, the effect of 2D-ASG's septal shadows on primary transmission and image quality was evaluated.Results.For an apparent septal thickness of 0.15 mm, the primary transmission of 2D-ASG varied between 72%-90% for grid pitches 1-3 mm. In low-contrast phantoms, the effect of 2D-ASG's radiopaque footprint on information loss was subtle. At high spatial frequencies, information loss manifested itself as undersampling artifacts, however, its impact on image quality is subtle when compared to quantum noise. Effects of additive electronic noise and gantry flex induced ring artifacts on image quality varied as a function of grid pitch and septal thickness. Such artifacts were substantially less in lower resolution images.Conclusion. The proposed simulation platform allowed successful evaluation of CBCT image quality variations as a function of 2D-ASG primary transmission properties and CBCT system characteristics. This platform can be potentially used for optimizing 2D-ASG design properties based on the imaging task and properties of the CBCT system.

Convolutional neural network-based model observer for signal known statistically task in breast tomosynthesis images

BACKGROUND

Since human observer studies are resource-intensive, mathematical model observers are frequently used to assess task-based image quality. The most common implementation of these model observers assume that the signal information is exactly known. However, these tasks cannot thoroughly represent situations where the signal information is not exactly known in terms of size and shape.

PURPOSE

Considering the limitations of the tasks for which signal information is exactly known, we proposed a convolutional neural network (CNN)-based model observer for signal known statistically (SKS) and background known statistically (BKS) detection tasks in breast tomosynthesis images.

METHODS

A wide parameter search was conducted from six different acquisition angles (i.e., 10°, 20°, 30°, 40°, 50°, and 60°) within the same dose level (i.e., 2.3 mGy) under two separate acquisition schemes: (1) constant total number of projections, and (2) constant angular separation between projections. Two different types of signals: spherical (i.e., SKE tasks) and spiculated (i.e., SKS tasks) were used. The detection performance of the CNN-based model observer was compared with that of the Hotelling observer (HO) instead of the IO. Pixel-wise gradient-weighted class activation mapping (pGrad-CAM) map was extracted from each reconstructed tomosynthesis image to provide an intuitive understanding of the trained CNN-based model observer.

RESULTS

The CNN-based model observer achieved a higher detection performance compared to that of the HO for all tasks. Moreover, the improvement in its detection performance was greater for SKS tasks compared to that for SKE tasks. These results demonstrated that the addition of nonlinearity improved the detection performance owing to the variation of the background and signal. Interestingly, the pGrad-CAM results effectively localized the class-specific discriminative region, further supporting the quantitative evaluation results of the CNN-based model observer. In addition, we verified that the CNN-based model observer required fewer images to achieve the detection performance of the HO.

CONCLUSIONS

In this work, we proposed a CNN-based model observer for SKS and BKS detection tasks in breast tomosynthesis images. Throughout the study, we demonstrated that the detection performance of the proposed CNN-based model observer was superior to that of the HO.

Image-based scatter correction for cone-beam CT using flip swin transformer U-shape network

BACKGROUND

Cone beam computed tomography (CBCT) plays an increasingly important role in image-guided radiation therapy. However, the image quality of CBCT is severely degraded by excessive scatter contamination, especially in the abdominal region, hindering its further applications in radiation therapy.

PURPOSE

To restore low-quality CBCT images contaminated by scatter signals, a scatter correction algorithm combining the advantages of convolutional neural networks (CNN) and Swin Transformer is proposed.

METHODS

In this paper a scatter correction model for CBCT image, the Flip Swin Transformer U-shape network (FSTUNet) model, is proposed. In this model, the advantages of CNN in texture detail and Swin Transformer in global correlation are used to accurately extract shallow and deep features, respectively. Instead of using the original Swin Transformer tandem structure, we build the Flip Swin Transformer Block to achieve a more powerful inter-window association extraction. The validity and clinical relevance of the method is demonstrated through extensive experiments on a Monte Carlo (MC) simulation dataset and frequency split dataset generated by a validated method, respectively.

RESULT

Experimental results on the MC simulated dataset show that the root mean square error of images corrected by the method is reduced from over 100 HU to about 7 HU. Both the structural similarity index measure (SSIM) and the universal quality index (UQI) are close to 1. Experimental results on the frequency split dataset demonstrate that the method not only corrects shading artifacts but also exhibits a high degree of structural consistency. In addition, comparison experiments show that FSTUNet outperforms UNet, Deep Residual Convolutional Neural Network (DRCNN), DSENet, Pix2pixGAN, and 3DUnet methods in both qualitative and quantitative metrics.

CONCLUSIONS

Accurately capturing the features at different levels is greatly beneficial for reconstructing high-quality scatter-free images. The proposed FSTUNet method is an effective solution to CBCT scatter correction and has the potential to improve the accuracy of CBCT image-guided radiation therapy.

MCR toolkit: A GPU-based toolkit for multi-channel reconstruction of preclinical and clinical x-ray CT data

BACKGROUND

The advancement of x-ray CT into the domains of photon counting spectral imaging and dynamic cardiac and perfusion imaging has created many new challenges and opportunities for clinicians and researchers. To address challenges such as dose constraints and scanning times while capitalizing on opportunities such as multi-contrast imaging and low-dose coronary angiography, these multi-channel imaging applications require a new generation of CT reconstruction tools. These new tools should exploit the relationships between imaging channels during reconstruction to set new image quality standards while serving as a platform for direct translation between the preclinical and clinical domains.

PURPOSE

We outline and demonstrate a new Multi-Channel Reconstruction (MCR) Toolkit for GPU-based analytical and iterative reconstruction of preclinical and clinical multi-energy and dynamic x-ray CT data. To promote open science, open-source distribution of the Toolkit will coincide with the release of this publication (GPL v3; gitlab.oit.duke.edu/dpc18/mcr-toolkit-public).

METHODS

The MCR Toolkit source code is implemented in C/C++ and NVIDIA's CUDA GPU programming interface, with scripting support from MATLAB and Python. The Toolkit implements matched, separable footprint CT reconstruction operators for projection and backprojection in two geometries: planar, cone-beam CT (CBCT) and 3rd generation, cylindrical multi-detector row CT (MDCT). Analytical reconstruction is performed using filtered backprojection (FBP) for circular CBCT, weighted FBP (WFBP) for helical CBCT, and cone-parallel projection rebinning followed by WFBP for MDCT. Arbitrary combinations of energy and temporal channels are iteratively reconstructed under a generalized multi-channel signal model for joint reconstruction. We solve this generalized model algebraically using the split Bregman optimization method and the BiCGSTAB(l) linear solver interchangeably for both CBCT and MDCT data. Rank-sparse kernel regression (RSKR) and patch-based singular value thresholding (pSVT) are used to regularize the energy and time dimensions, respectively. Under a Gaussian noise model, regularization parameters are estimated automatically from the input data, dramatically reducing algorithm complexity for end users. Multi-GPU parallelization of the reconstruction operators is supported to manage reconstruction times.

RESULTS

Denoising with RSKR and pSVT and post-reconstruction material decomposition are illustrated with preclinical and clinical cardiac photon-counting (PC)CT data. A digital MOBY mouse phantom with cardiac motion is used to illustrate single energy (SE), multi-energy (ME), time resolved (TR), and combined multi-energy and time-resolved (METR) helical, CBCT reconstruction. A fixed set of projection data is used across all reconstruction cases to demonstrate the Toolkit's robustness to increasing data dimensionality. Identical reconstruction code is applied to in vivo cardiac PCCT data acquired in a mouse model of atherosclerosis (METR). Clinical cardiac CT reconstruction is illustrated using the XCAT phantom and the DukeSim CT simulator, while dual-source, dual-energy CT reconstruction is illustrated for data acquired with a Siemens Flash scanner. Benchmarking results with NVIDIA RTX 8000 GPU hardware demonstrate 61%-99% efficiency in scaling computation from one to four GPUs for these reconstruction problems.

CONCLUSIONS

The MCR Toolkit provides a robust solution for temporal and spectral x-ray CT reconstruction problems and was built from the ground up to facilitate translation of CT research and development between preclinical and clinical applications.