TIGRE: a MATLAB-GPU toolbox for CBCT image reconstruction

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.

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.

Sinogram Inpainting with Generative Adversarial Networks and Shape Priors

X-ray computed tomography is a widely used, non-destructive imaging technique that computes cross-sectional images of an object from a set of X-ray absorption profiles (the so-called sinogram). The computation of the image from the sinogram is an ill-posed inverse problem, which becomes underdetermined when we are only able to collect insufficiently many X-ray measurements. We are here interested in solving X-ray tomography image reconstruction problems where we are unable to scan the object from all directions, but where we have prior information about the object's shape. We thus propose a method that reduces image artefacts due to limited tomographic measurements by inferring missing measurements using shape priors. Our method uses a Generative Adversarial Network that combines limited acquisition data and shape information. While most existing methods focus on evenly spaced missing scanning angles, we propose an approach that infers a substantial number of consecutive missing acquisitions. We show that our method consistently improves image quality compared to images reconstructed using the previous state-of-the-art sinogram-inpainting techniques. In particular, we demonstrate a 7 dB Peak Signal-to-Noise Ratio improvement compared to other methods.

Characterization of Pharmaceutical Tablets by X-ray Tomography

Solid dosage forms such as tablets are extensively used in drug administration for their simplicity and large-scale manufacturing capabilities. High-resolution X-ray tomography is one of the most valuable non-destructive techniques to investigate the internal structure of the tablets for drug product development as well as for a cost effective production process. In this work, we review the recent developments in high-resolution X-ray microtomography and its application towards different tablet characterizations. The increased availability of powerful laboratory instrumentation, as well as the advent of high brilliance and coherent 3rd generation synchrotron light sources, combined with advanced data processing techniques, are driving the application of X-ray microtomography forward as an indispensable tool in the pharmaceutical industry.

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.

PEPI Lab: a flexible compact multi-modal setup for X-ray phase-contrast and spectral imaging

This paper presents a new flexible compact multi-modal imaging setup referred to as PEPI (Photon-counting Edge-illumination Phase-contrast imaging) Lab, which is based on the edge-illumination (EI) technique and a chromatic detector. The system enables both X-ray phase-contrast (XPCI) and spectral (XSI) imaging of samples on the centimeter scale. This work conceptually follows all the stages in its realization, from the design to the first imaging results. The setup can be operated in four different modes, i.e. photon-counting/conventional, spectral, double-mask EI, and single-mask EI, whereby the switch to any modality is fast, software controlled, and does not require any hardware modification or lengthy re-alignment procedures. The system specifications, ranging from the X-ray tube features to the mask material and aspect ratio, have been quantitatively studied and optimized through a dedicated Geant4 simulation platform, guiding the choice of the instrumentation. The realization of the imaging setup, both in terms of hardware and control software, is detailed and discussed with a focus on practical/experimental aspects. Flexibility and compactness (66 cm source-to-detector distance in EI) are ensured by dedicated motion stages, whereas spectral capabilities are enabled by the Pixirad-1/Pixie-III detector in combination with a tungsten anode X-ray source operating in the range 40-100 kVp. The stability of the system, when operated in EI, has been verified, and drifts leading to mask misalignment of less than 1 [Formula: see text]m have been measured over a period of 54 h. The first imaging results, one for each modality, demonstrate that the system fulfills its design requirements. Specifically, XSI tomographic images of an iodine-based phantom demonstrate the system's quantitativeness and sensibility to concentrations in the order of a few mg/ml. Planar XPCI images of a carpenter bee specimen, both in single and double-mask modes, demonstrate that refraction sensitivity (below 0.6 [Formula: see text]rad in double-mask mode) is comparable with other XPCI systems based on microfocus sources. Phase CT capabilities have also been tested on a dedicated plastic phantom, where the phase channel yielded a 15-fold higher signal-to-noise ratio with respect to attenuation.

A total variation prior unrolling approach for computed tomography reconstruction

BACKGROUND

With the rapid development of deep learning technology, deep neural networks can effectively enhance the performance of computed tomography (CT) reconstructions. One kind of commonly used method to construct CT reconstruction networks is to unroll the conventional iterative reconstruction (IR) methods to convolutional neural networks (CNNs). However, most unrolling methods primarily unroll the fidelity term of IR methods to CNNs, without unrolling the prior terms. The prior terms are always directly replaced by neural networks.

PURPOSE

In conventional IR methods, the prior terms play a vital role in improving the visual quality of reconstructed images. Unrolling the hand-crafted prior terms to CNNs may provide a more specialized unrolling approach to further improve the performance of CT reconstruction. In this work, a primal-dual network (PD-Net) was proposed by unrolling both the data fidelity term and the total variation (TV) prior term, which effectively preserves the image edges and textures in the reconstructed images.

METHODS

By further deriving the Chambolle-Pock (CP) algorithm instance for CT reconstruction, we discovered that the TV prior updates the reconstructed images with its divergences in each iteration of the solution process. Based on this discovery, CNNs were applied to yield the divergences of the feature maps for the reconstructed image generated in each iteration. Additionally, a loss function was applied to the predicted divergences of the reconstructed image to guarantee that the CNNs' results were the divergences of the corresponding feature maps in the iteration. In this manner, the proposed CNNs seem to play the same roles in the PD-Net as the TV prior in the IR methods. Thus, the TV prior in the CP algorithm instance can be directly unrolled to CNNs.

RESULTS

The datasets from the Low-Dose CT Image and Projection Data and the Piglet dataset were employed to assess the effectiveness of our proposed PD-Net. Compared with conventional CT reconstruction methods, our proposed method effectively preserves the structural and textural information in reference to ground truth.

CONCLUSIONS

The experimental results show that our proposed PD-Net framework is feasible for the implementation of CT reconstruction tasks. Owing to the promising results yielded by our proposed neural network, this study is intended to inspire further development of unrolling approaches by enabling the direct unrolling of hand-crafted prior terms to CNNs.

Design and CT imaging of casper, an anthropomorphic breathing thorax phantom

The goal of this work was to build an anthropomorphic thorax phantom capable of breathing motion with materials mimicking human tissues in x-ray imaging applications. The thorax phantom, named Casper, was composed of resin (body), foam (lungs), glow polyactic acid (bones) and natural polyactic acid (tumours placed in the lungs). X-ray attenuation properties of all materials prior to manufacturing were evaluated by means of photon-counting computed tomography (CT) imaging on a table-top system. Breathing motion was achieved by a scotch-yoke mechanism with diaphragm motion frequencies of 10-20 rpm and displacements of 1 to 2 cm. Casper was manufactured by means of 3D printing of moulds and ribs and assembled in a complex process. The final phantom was then scanned using a clinical CT scanner to evaluate material CT numbers and the extent of tumour motion. Casper CT numbers were close to human CT numbers for soft tissue (46 HU), ribs (125 HU), lungs (-840 HU) and tumours (-45 HU). For a 2 cm diaphragm displacement the largest tumour displacement was 0.7 cm. The five tumour volumes were accurately assessed in the static CT images with a mean absolute error of 4.3%. Tumour sizes were either underestimated for smaller tumours or overestimated for larger tumours in dynamic CT images due to motion blurring with a mean absolute difference from true volumes of 10.3%. More Casper information including a motion movie and manufacturing data can be downloaded from http://web.uvic.ca/~bazalova/Casper/.

Machine learning for detection of 3D features using sparse x-ray tomographic reconstruction

In many inertial confinement fusion (ICF) experiments, the neutron yield and other parameters cannot be completely accounted for with one and two dimensional models. This discrepancy suggests that there are three dimensional effects that may be significant. Sources of these effects include defects in the shells and defects in shell interfaces, the fill tube of the capsule, and the joint feature in double shell targets. Due to their ability to penetrate materials, x rays are used to capture the internal structure of objects. Methods such as computational tomography use x-ray radiographs from hundreds of projections, in order to reconstruct a three dimensional model of the object. In experimental environments, such as the National Ignition Facility and Omega-60, the availability of these views is scarce, and in many cases only consists of a single line of sight. Mathematical reconstruction of a 3D object from sparse views is an ill-posed inverse problem. These types of problems are typically solved by utilizing prior information. Neural networks have been used for the task of 3D reconstruction as they are capable of encoding and leveraging this prior information. We utilize half a dozen, different convolutional neural networks to produce different 3D representations of ICF implosions from the experimental data. Deep supervision is utilized to train a neural network to produce high-resolution reconstructions. These representations are used to track 3D features of the capsules, such as the ablator, inner shell, and the joint between shell hemispheres. Machine learning, supplemented by different priors, is a promising method for 3D reconstructions in ICF and x-ray radiography, in general.

Deep 3D reconstruction of synchrotron X-ray computed tomography for intact lungs

Synchrotron X-rays can be used to obtain highly detailed images of parts of the lung. However, micro-motion artifacts induced by such as cardiac motion impede quantitative visualization of the alveoli in the lungs. This paper proposes a method that applies a neural network for synchrotron X-ray Computed Tomography (CT) data to reconstruct the high-quality 3D structure of alveoli in intact mouse lungs at expiration, without needing ground-truth data. Our approach reconstructs the spatial sequence of CT images by using a deep-image prior with interpolated input latent variables, and in this way significantly enhances the images of alveolar structure compared with the prior art. The approach successfully visualizes 3D alveolar units of intact mouse lungs at expiration and enables us to measure the diameter of the alveoli. We believe that our approach helps to accurately visualize other living organs hampered by micro-motion.

3D Reconstruction of Wrist Bones from C-Arm Fluoroscopy Using Planar Markers

In orthopedic surgeries, such as osteotomy and osteosynthesis, an intraoperative 3D reconstruction of the bone would enable surgeons to quickly assess the fracture reduction procedure with preoperative planning. Scanners equipped with such functionality are often more expensive than a conventional C-arm fluoroscopy device. Moreover, a C-arm fluoroscopy device is commonly available in many orthopedic facilities. Based on the widespread use of such equipment, this paper proposes a method to reconstruct the 3D structure of bone with a conventional C-arm fluoroscopy device. We focus on wrist bones as the target of reconstruction in this research as this will facilitate a flexible imaging scheme. Planar markers are attached to the target object and are tracked in the fluoroscopic image for C-arm pose estimation. The initial calibration of the device is conducted using a checkerboard pattern. In general, reconstruction algorithms are sensitive to geometric calibration errors. To assess the practicality of the method for reconstruction, a simulation study demonstrating the effect of checkerboard thickness and spherical marker size on reconstruction quality was conducted.

An empirical method for geometric calibration of a photon counting detector-based cone beam CT system

BACKGROUND

Geometric calibration is essential in developing a reliable computed tomography (CT) system. It involves estimating the geometry under which the angular projections are acquired. Geometric calibration of cone beam CTs employing small area detectors, such as currently available photon counting detectors (PCDs), is challenging when using traditional-based methods due to detectors' limited areas.

OBJECTIVE

This study presented an empirical method for the geometric calibration of small area PCD-based cone beam CT systems.

METHODS

Unlike the traditional methods, we developed an iterative optimization procedure to determine geometric parameters using the reconstructed images of small metal ball bearings (BBs) embedded in a custom-built phantom. An objective function incorporating the sphericities and symmetries of the embedded BBs was defined to assess performance of the reconstruction algorithm with the given initial estimated set of geometric parameters. The optimal parameter values were those which minimized the objective function. The TIGRE toolbox was employed for fast tomographic reconstruction. To evaluate the proposed method, computer simulations were carried out using various numbers of spheres placed in various locations. Furthermore, efficacy of the method was experimentally assessed using a custom-made benchtop PCD-based cone beam CT.

RESULTS

Computer simulations validated the accuracy and reproducibility of the proposed method. The precise estimation of the geometric parameters of the benchtop revealed high-quality imaging in CT reconstruction of a breast phantom. Within the phantom, the cylindrical holes, fibers, and speck groups were imaged in high fidelity. The CNR analysis further revealed the quantitative improvements of the reconstruction performed with the estimated parameters using the proposed method.

CONCLUSION

Apart from the computational cost, we concluded that the method was easy to implement and robust.

Metal artifact correction in photon-counting detector computed tomography: metal trace replacement using high-energy data

BACKGROUND

Metal artifacts have been an outstanding issue in computed tomography (CT) since its first uses in the clinic and continue to interfere. Metal artifact reduction (MAR) methods continue to be proposed and photon-counting detectors (PCDs) have recently been the subject of research toward this purpose. PCDs offer the ability to distinguish the energy of incident x-rays and sort them in a set number of energy bins. High-energy data captured using PCDs have been shown to reduce metal artifacts in reconstructions due to reduced beam hardening.

PURPOSE

High-energy reconstructions using PCD-CT have their drawbacks, such as reduced image contrast and increased noise. Here, we demonstrate a MAR algorithm, trace replacement MAR (TRMAR), in which the data corrupted by metal artifacts in full energy spectrum projections are corrected using the high-energy data captured during the same scan. The resulting reconstructions offer similar MAR to that seen in high-energy reconstructions, but with improved image quality.

METHODS

Experimental data were collected using a bench-top PCD-CT system with a cadmium zinc telluride PCD. Simulations were performed to determine the optimal high-energy threshold and to test TRMAR in simulations using the XCAT phantom and a biological sample. For experiments a 100-mm diameter cylindrical phantom containing vials of water, two screws, various densities of Ca(ClO₄ )₂ , and a spatial resolution phantom was imaged with and without the screws. The screws were segmented in the initial reconstruction and forward projected to identify them in the sinogram space in order to perform TRMAR. The resulting reconstructions were compared to the control and to reconstructions corrected using normalized metal artifact reduction (NMAR). Additionally, a beef short rib was imaged with and without metal to provide a more realistic phantom.

RESULTS

XCAT simulations showed a reduction in the streak artifact from -978 HU in uncorrected images to -10 HU with TRMAR. The magnitude of the metal artifact in uncorrected images of the 100-mm phantom was -442 HU, compared to the desired -81 HU with no metal. TRMAR reduced the magnitude of the artifact to -142 HU, with NMAR reducing the magnitude to -96 HU. Relative image noise was reduced from 176% in the high-energy image to 56% using TRMAR. Density quantification was better with NMAR, with the Ca(ClO₄ )₂ vial affected most by metal artifacts showing 0.8% error compared to 2.1% with TRMAR. Small features were preserved to a greater extent with TRMAR, with the limiting spatial frequency at 20% of the MTF fully maintained at 1.31 lp/mm, while with NMAR it was reduced to 1.22 lp/mm. Images of the beef short rib showed better delineation of the shape of the metal using TRMAR.

CONCLUSIONS

NMAR offers slightly better performance compared to TRMAR in streak reduction and image quality metrics. However, TRMAR is less susceptible to metal segmentation errors and can closely approximate the reduction in the streak metal artifact seen in NMAR at 1/3 the computation time. With the recent introduction of PCD-CT into the clinic, TRMAR offers notable potential for fast, effective MAR.

Multi-material spectral photon-counting micro-CT with minimum residual decomposition and self-supervised deep denoising

Spectral micro-CT imaging with direct-detection energy discriminating photon counting detectors having small pixel size (< 100×100 µm²) is mainly hampered by: i) the limited energy resolution of the imaging device due to charge sharing effects and ii) the unavoidable noise amplification in the images resulting from basis material decomposition. In this work, we present a cone-beam micro-CT setup that includes a CdTe photon counting detector implementing a charge summing hardware solution to correct for the charge-sharing issue and an innovative image processing pipeline based on accurate modeling of the spectral response of the imaging system, an improved basis material decomposition (BMD) algorithm named minimum-residual BMD (MR-BMD), and self-supervised deep convolutional denoising. Experimental tomographic projections having a pixel size of 45×45 µm² of a plastinated mouse sample including I, Ba, and Gd small cuvettes were acquired. Results demonstrate the capability of the combined hardware and software tools to sharply discriminate even between materials having their K-Edge separated by a few keV, such as e.g., I and Ba. By evaluating the quality of the reconstructed decomposed images (water, bone, I, Ba, and Gd), the quantitative performances of the spectral system are here assessed and discussed.

A dual-domain neural network based on sinogram synthesis for sparse-view CT reconstruction

OBJECTIVE

The dual-domain deep learning-based reconstruction techniques have enjoyed many successful applications in the field of medical image reconstruction. Applying the analytical reconstruction based operator to transfer the data from the projection domain to the image domain, the dual-domain techniques may suffer from the insufficient suppression or removal of streak artifacts in areas with the missing view data, when addressing the sparse-view reconstruction problems. In this work, to overcome this problem, an intelligent sinogram synthesis based back-projection network (iSSBP-Net) was proposed for sparse-view computed tomography (CT) reconstruction. In the iSSBP-Net method, a convolutional neural network (CNN) was involved in the dual-domain method to inpaint the missing view data in the sinogram before CT reconstruction.

METHODS

The proposed iSSBP-Net method fused a sinogram synthesis sub-network (SS-Net), a sinogram filter sub-network (SF-Net), a back-projection layer, and a post-CNN into an end-to-end network. Firstly, to inpaint the missing view data, the SS-Net employed a CNN to synthesize the full-view sinogram in the projection domain. Secondly, to improve the visual quality of the sparse-view CT images, the synthesized sinogram was filtered by a CNN. Thirdly, the filtered sinogram was brought into the image domain through the back-projection layer. Finally, to yield images of high visual sensitivity, the post-CNN was applied to restore the desired images from the outputs of the back-projection layer.

RESULTS

The numerical experiments demonstrate that the proposed iSSBP-Net is superior to all competing algorithms under different scanning condintions for sparse-view CT reconstruction. Compared to the competing algorithms, the proposed iSSBP-Net method improved the peak signal-to-noise ratio of the reconstructed images about 1.21 dB, 0.26 dB, 0.01 dB, and 0.37 dB under the scanning conditions of 360, 180, 90, and 60 views, respectively.

CONCLUSION

The promising reconstruction results indicate that involving the SS-Net in the dual-domain method is could be an effective manner to suppress or remove the streak artifacts in sparse-view CT images. Due to the promising results reconstructed by the iSSBP-Net method, this study is intended to inspire the further development of sparse-view CT reconstruction by involving a SS-Net in the dual-domain method.

TIGRE-VarianCBCT for on-board cone-beam computed tomography, an open-source toolkit for imaging, dosimetry and clinical research

We presented TIGRE-VarianCBCT, an open-source toolkit Matlab-GPU for Varian on-board cone-beam CT with particular emphasis to address challenges in raw data preprocessing, artifacts correction, tomographic reconstruction and image post-processing. The aim of this project is to provide not only a tool to bridge the gap between clinical usage of CBCT scan data and research algorithms but also a framework that breaks down the imaging chain into individual processes so that research effort can be focused on a specific part. The entire imaging chain, module-based architecture, data flow and techniques used in the creation of the toolkit are presented. Raw scan data are first decoded to extract X-ray fluoro image series and set up the imaging geometry. Data conditioning operations including scatter correction, normalization, beam-hardening correction, ring removal are performed sequentially. Reconstruction is supported by TIGRE with FDK as well as a variety of iterative algorithms. Pixel-to-HU mapping is calibrated by a Catphan^TM 504 phantom. Imaging dose in CTDIw is calculated in an empirical formula. The performance was validated on real patient scans with good agreement with respect to vendor-designed program. Case studies in scan protocol optimization, low dose imaging and iterative algorithm comparison demonstrated its substantial potential in performing scan data based clinical studies. The toolkit is released under the BSD license, imposing minimal restrictions on its use and distribution. The toolkit is accessible as a module at https://github.com/CERN/TIGRE.

A Fully Automated Method for 3D Individual Tooth Identification and Segmentation in Dental CBCT

Accurate and automatic segmentation of three-dimensional (3D) individual teeth from cone-beam computerized tomography (CBCT) images is a challenging problem because of the difficulty in separating an individual tooth from adjacent teeth and its surrounding alveolar bone. Thus, this paper proposes a fully automated method of identifying and segmenting 3D individual teeth from dental CBCT images. The proposed method addresses the aforementioned difficulty by developing a deep learning-based hierarchical multi-step model. First, it automatically generates upper and lower jaws panoramic images to overcome the computational complexity caused by high-dimensional data and the curse of dimensionality associated with limited training dataset. The obtained 2D panoramic images are then used to identify 2D individual teeth and capture loose- and tight- regions of interest (ROIs) of 3D individual teeth. Finally, accurate 3D individual tooth segmentation is achieved using both loose and tight ROIs. Experimental results showed that the proposed method achieved an F1-score of 93.35 percent for tooth identification and a Dice similarity coefficient of 94.79 percent for individual 3D tooth segmentation. The results demonstrate that the proposed method provides an effective clinical and practical framework for digital dentistry.

Image fusion for neutron tomography of nuclear fuel

Image fusion, the process of combining different images together, can be useful to create a more complete picture. In this work, image fusion is applied to neutron tomography of nuclear fuel with the goal of enhancing the information obtained about the fuel. Different reconstruction methods, such as Feldkamp, Davis and Kress filtered back projection and Simultaneous Reconstruction Technique, were combined to enhance image quality. This methodology was shown to reduce noise and ring artifacts without sacrificing sharp edges, allowing for a more accurate representation of sample geometry. Technique enhancements and future applications for the neutron imaging community are also discussed.

Deep learning method for reducing metal artifacts in dental cone-beam CT using supplementary information from intra-oral scan

Objective. Recently, dental cone-beam computed tomography (CBCT) methods have been improved to significantly reduce radiation dose while maintaining image resolution with minimal equipment cost. In low-dose CBCT environments, metallic inserts such as implants, crowns, and dental fillings cause severe artifacts, which result in a significant loss of morphological structures of teeth in reconstructed images. Such metal artifacts prevent accurate 3D bone-teeth-jaw modeling for diagnosis and treatment planning. However, the performance of existing metal artifact reduction (MAR) methods in handling the loss of the morphological structures of teeth in reconstructed CT images remains relatively limited. In this study, we developed an innovative MAR method to achieve optimal restoration of anatomical details. Approach. The proposed MAR approach is based on a two-stage deep learning-based method. In the first stage, we employ a deep learning network that utilizes intra-oral scan data as side-inputs and performs multi-task learning of auxiliary tooth segmentation. The network is designed to improve the learning ability of capturing teeth-related features effectively while mitigating metal artifacts. In the second stage, a 3D bone-teeth-jaw model is constructed with weighted thresholding, where the weighting region is determined depending on the geometry of the intra-oral scan data. Main results. The results of numerical simulations and clinical experiments are presented to demonstrate the feasibility of the proposed approach. Significance. We propose for the first time a MAR method using radiation-free intra-oral scan data as supplemental information on the tooth morphological structures of teeth, which is designed to perform accurate 3D bone-teeth-jaw modeling in low-dose CBCT environments.

Concurrent kilovoltage CBCT imaging and megavoltage beam delivery: suppression of cross-scatter with 2D antiscatter grids and grid-based scatter sampling

Objective. The concept of using kilovoltage (kV) and megavoltage (MV) beams concurrently has potential applications in cone beam computed tomography (CBCT) guided radiation therapy, such as single breath hold scans, metal artifact reduction, and simultaneous imaging during MV treatment delivery. However, MV cross-scatter generated during MV beam delivery degrades CBCT image quality. To address this, a 2D antiscatter grid and a cross-scatter correction method were investigated in the context of high dose MV treatment delivery.Approach. A 3D printed, tungsten 2D antiscatter grid prototype was utilized in kV CBCT scans to reduce MV cross-scatter fluence during concurrent MV beam delivery. Remaining cross-scatter in projections was corrected by using the 2D grid itself as a cross-scatter intensity sampling device, referred to as grid-based scatter sampling (GSS). To test this approach, kV CBCT acquisitions were performed while delivering 6 and 10 MV beams, mimicking high dose rate treatment delivery scenarios. kV and MV beam deliveries were not synchronized to eliminate MV beam delivery interruption. MV cross-scatter suppression performance of the proposed approach was evaluated in projections and CBCT images of phantoms.Main results. 2D grid reduced the intensity of MV cross-scatter in kV projections by a factor of 3 on the average, when compared to conventional antiscatter grid. Remaining cross scatter as measured by the GSS method was within 7% of measured reference intensity values, and subsequently corrected. CBCT image quality was improved substantially during concurrent kV-MV beam delivery. Median Hounsfield Unit (HU) inaccuracy was up to 182 HU without our methods, and it was reduced to a median 6.5 HU with our 2D grid and scatter correction approach. Our methods provided a factor of 2-6 improvement in contrast-to-noise ratio.Significance. This investigation demonstrates the utility of 2D antiscatter grids and grid-based scatter sampling in suppressing MV cross-scatter. Our approach successfully minimized the effects of MV cross-scatter in concurrent kV CBCT imaging and high dose MV treatment delivery scenarios. Hence, robust MV cross-scatter suppression is potentially feasible without MV beam delivery interruption or compromising kV image acquisition rate.

DIOR: Deep Iterative Optimization-Based Residual-Learning for Limited-Angle CT Reconstruction

Limited-angle CT is a challenging problem in real applications. Incomplete projection data will lead to severe artifacts and distortions in reconstruction images. To tackle this problem, we propose a novel reconstruction framework termed Deep Iterative Optimization-based Residual-learning (DIOR) for limited-angle CT. Instead of directly deploying the regularization term on image space, the DIOR combines iterative optimization and deep learning based on the residual domain, significantly improving the convergence property and generalization ability. Specifically, the asymmetric convolutional modules are adopted to strengthen the feature extraction capacity in smooth regions for deep priors. Besides, in our DIOR method, the information contained in low-frequency and high-frequency components is also evaluated by perceptual loss to improve the performance in tissue preservation. Both simulated and clinical datasets are performed to validate the performance of DIOR. Compared with existing competitive algorithms, quantitative and qualitative results show that the proposed method brings a promising improvement in artifact removal, detail restoration and edge preservation.

A generalized image quality improvement strategy of cone-beam CT using multiple spectral CT labels in Pix2pix GAN

Objective. The quantitative and routine imaging capabilities of cone-beam CT (CBCT) are hindered from clinical applications due to the severe shading artifacts of scatter contamination. The scatter correction methods proposed in the literature only consider the anatomy of the scanned objects while disregarding the impact of incident x-ray energy spectra. The multiple-spectral model is in urgent need for CBCT scatter estimation. Approach. In this work, we incorporate the multiple spectral diagnostic multidetector CT labels into the pixel-to-pixel (Pix2pix) GAN to estimate accurate scatter distributions from CBCT projections acquired at various imaging volume sizes and x-ray energy spectra. The Pix2pix GAN combines the residual network as the generator and the PatchGAN as the discriminator to construct the correspondence between the scatter-contaminated projection and scatter distribution. The network architectures and loss function of Pix2pix GAN are optimized to achieve the best performance on projection-to-scatter transition. Results. The CBCT data of a head phantom and abdominal patients are applied to test the performance of the proposed method. The error of the corrected CBCT image using the proposed method is reduced from over 200 HU to be around 20 HU in both phantom and patient studies. The mean structural similarity index of the CT image is improved from 0.2 to around 0.9 after scatter correction using the proposed method compared with the MC-simulation method, which indicates a high similarity of the anatomy in the images before and after the proposed correction. The proposed method achieves higher accuracy of scatter estimation than using the Pix2pix GAN with the U-net generator. Significance. The proposed scheme is an effective solution to the multiple spectral CBCT scatter correction. The scatter-correction software using the proposed model will be available at: https://github.com/YangkangJiang/Cone-beam-CT-scatter-correction-tool.

Generative Adversarial Network (GAN) for Automatic Reconstruction of the 3D Spine Structure by Using Simulated Bi-Planar X-ray Images

In this study, we modified the previously proposed X2CT-GAN to build a 2Dto3D-GAN of the spine. This study also incorporated the radiologist’s perspective in the adjustment of input signals to prove the feasibility of the automatic production of three-dimensional (3D) structures of the spine from simulated bi-planar two-dimensional (2D) X-ray images. Data from 1012 computed tomography (CT) studies of 984 patients were retrospectively collected. We tested this model under different dataset sizes (333, 666, and 1012) with different bone signal conditions to observe the training performance. A 10-fold cross-validation and five metrics—Dice similarity coefficient (DSC) value, Jaccard similarity coefficient (JSC), overlap volume (OV), and structural similarity index (SSIM)—were applied for model evaluation. The optimal mean values for DSC, JSC, OV, SSIM_anteroposterior (AP), and SSIM_Lateral (Lat) were 0.8192, 0.6984, 0.8624, 0.9261, and 0.9242, respectively. There was a significant improvement in the training performance under empirically enhanced bone signal conditions and with increasing training dataset sizes. These results demonstrate the potential of the clinical implantation of GAN for automatic production of 3D spine images from 2D images. This prototype model can serve as a foundation in future studies applying transfer learning for the development of advanced medical diagnostic techniques.

Cone-beam breast CT using an offset detector: effect of detector offset and image reconstruction algorithm

Objective.A dedicated cone-beam breast computed tomography (BCT) using a high-resolution, low-noise detector operating in offset-detector geometry has been developed. This study investigates the effects of varying detector offsets and image reconstruction algorithms to determine the appropriate combination of detector offset and reconstruction algorithm.Approach.Projection datasets (300 projections in 360°) of 30 breasts containing calcified lesions that were acquired using a prototype cone-beam BCT system comprising a 40 × 30 cm flat-panel detector with 1024 × 768 detector pixels were reconstructed using Feldkamp-Davis-Kress (FDK) algorithm and served as the reference. The projection datasets were retrospectively truncated to emulate cone-beam datasets with sinograms of 768×768 and 640×768 detector pixels, corresponding to 5 cm and 7.5 cm lateral offsets, respectively. These datasets were reconstructed using the FDK algorithm with appropriate weights and an ASD-POCS-based Fast, total variation-Regularized, Iterative, Statistical reconstruction Technique (FRIST), resulting in a total of 4 offset-detector reconstructions (2 detector offsets × 2 reconstruction methods). Signal difference-to-noise ratio (SDNR), variance, and full-width at half-maximum (FWHM) of calcifications in two orthogonal directions were determined from all reconstructions. All quantitative measurements were performed on images in units of linear attenuation coefficient (1/cm).Results.The FWHM of calcifications did not differ (P > 0.262) among reconstruction algorithms and detector formats, implying comparable spatial resolution. For a chosen detector offset, the FRIST algorithm outperformed FDK in terms of variance and SDNR (P < 0.0001). For a given reconstruction method, the 5 cm offset provided better results.Significance.This study indicates the feasibility of using the compressed sensing-based, FRIST algorithm to reconstruct sinograms from offset-detectors. Among the reconstruction methods and detector offsets studied, FRIST reconstructions corresponding to a 30 cm × 30 cm with 5 cm lateral offset, achieved the best performance. A clinical prototype using such an offset geometry has been developed and installed for clinical trials.

Megavoltage cross-scatter rejection and correction using 2D antiscatter grids in kilovoltage CBCT imaging

Simultaneous use of kilovoltage (kV) and megavoltage (MV) beams has numerous potential applications in cone beam computed tomography (CBCT)-guided radiotherapy, such as fast MV+kV CBCT for single breath-hold scan, tumor localization with kV CBCT imaging during MV therapy delivery, and metal artifact suppression. However, the introduction of MV beams results in a large MV-cross scatter fluence incident on the kV Flat Panel Detector (FPD), and thus, deteriorating the low contrast visualization and Hounsfield Unit (HU) accuracy. In this work, we introduced a novel and robust method for reducing the effects of MV cross scatter. First, we implemented a 2D antiscatter grid atop the detector which rejects a large section of MV cross scatter. This hardware-based approach, while effective, allows a fraction of MV cross scatter to be transmitted to the FPD, resulting in artifacts and degraded HU accuracy in CBCT images. We thus introduced a data correction step, which aimed to estimate and correct the remaining MV cross scatter. This approach, referred to as Grid-Based Scatter Sampling, utilized 2D antiscatter grid itself to measure and correct remaining MV cross scatter in projections. We investigated the performance of the proposed approach in experiments by simultaneously acquiring kV CBCT and delivering MV beams with a clinical linac. The results show that the proposed method can substantially reduce HU inaccuracy and increase contrast-to-noise ratio (CNR). Our method does not require synchronization of kV and MV beam pulses, reduction of kV frame acquisition rate, or MV dose rate, and therefore, it is more practical to implement in radiation therapy clinical setting.

Low-dose cone-beam computed tomography reconstruction through a fast three-dimensional compressed sensing method based on the three-dimensional pseudo-polar fourier transform

Background:

Reconstruction of high quality two dimensional images from fan beam computed tomography (CT) with a limited number of projections is already feasible through Fourier based iterative reconstruction method. However, this article is focused on a more complicated reconstruction of three dimensional (3D) images in a sparse view cone beam computed tomography (CBCT) by utilizing Compressive Sensing (CS) based on 3D pseudo polar Fourier transform (PPFT).

Method:

In comparison with the prevalent Cartesian grid, PPFT re gridding is potent to remove rebinning and interpolation errors. Furthermore, using PPFT based radon transform as the measurement matrix, reduced the computational complexity.

Results:

In order to show the computational efficiency of the proposed method, we compare it with an algebraic reconstruction technique and a CS type algorithm. We observed convergence in <20 iterations in our algorithm while others would need at least 50 iterations for reconstructing a qualified phantom image. Furthermore, using a fast composite splitting algorithm solver in each iteration makes it a fast CBCT reconstruction algorithm. The algorithm will minimize a linear combination of three terms corresponding to a least square data fitting, Hessian (HS) Penalty and l1 norm wavelet regularization. We named it PP-based compressed sensing-HS-W. In the reconstruction range of 120 projections around the 360° rotation, the image quality is visually similar to reconstructed images by Feldkamp-Davis-Kress algorithm using 720 projections. This represents a high dose reduction.

Conclusion:

The main achievements of this work are to reduce the radiation dose without degrading the image quality. Its ability in removing the staircase effect, preserving edges and regions with smooth intensity transition, and producing high-resolution, low-noise reconstruction results in low-dose level are also shown.