Faster and lower dose imaging: evaluating adaptive, constant gantry velocity and angular separation in fast low-dose 4D cone beam CT imaging

BACKGROUND

The adoption of four-dimensional cone beam computed tomography (4DCBCT) for image-guided lung cancer radiotherapy is increasing, especially for hypofractionated treatments. However, the drawbacks of 4DCBCT include long scan times (∼240 s), inconsistent image quality, higher imaging dose than necessary, and streaking artifacts. With the emergence of linear accelerators that can acquire 4DCBCT scans in a short period of time (9.2 s) there is a need to examine the impact that these very fast gantry rotations have on 4DCBCT image quality.

PURPOSE

This study investigates the impact of gantry velocity and angular separation between x-ray projections on image quality and its implication for fast low-dose 4DCBCT with emerging systems, such as the Varian Halcyon that provide fast gantry rotation and imaging. Large and uneven angular separation between x-ray projections is known to reduce 4DCBCT image quality through increased streaking artifacts. However, it is not known when angular separation starts degrading image quality. The study assesses the impact of constant and adaptive gantry velocity and determines the level when angular gaps impair image quality using state-of-the-art reconstruction methods.

METHODS

This study considers fast low-dose 4DCBCT acquisitions (60-80 s, 200-projection scans). To assess the impact of adaptive gantry rotations, the angular position of x-ray projections from adaptive 4DCBCT acquisitions from a 30-patient clinical trial were analyzed (referred to as patient angular gaps). To assess the impact of angular gaps, variable and static angular gaps (20°, 30°, 40°) were introduced into evenly separated 200 projections (ideal angular separation). To simulate fast gantry rotations, which are on emerging linacs, constant gantry velocity acquisitions (9.2 s, 60 s, 120 s, 240 s) were simulated by sampling x-ray projections at constant intervals using the patient breathing traces from the ADAPT clinical trial (ACTRN12618001440213). The 4D Extended Cardiac-Torso (XCAT) digital phantom was used to simulate projections to remove patient-specific image quality variables. Image reconstruction was performed using Feldkamp-Davis-Kress (FDK), McKinnon-Bates (MKB), and Motion-Compensated-MKB (MCMKB) algorithms. Image quality was assessed using Structural Similarity-Index-Measure (SSIM), Contrast-to-Noise-Ratio (CNR), Signal-to-Noise-Ratio (SNR), Tissue-Interface-Width-Diaphragm (TIW-D), and Tissue-Interface-Width-Tumor (TIW-T).

RESULTS

Patient angular gaps and variable angular gap reconstructions produced similar results to ideal angular separation reconstructions, while static angular gap reconstructions produced lower image quality metrics. For MCMKB-reconstructions, average patient angular gaps produced SSIM-0.98, CNR-13.6, SNR-34.8, TIW-D-1.5 mm, and TIW-T-2.0 mm, static angular gap 40° produced SSIM-0.92, CNR-6.8, SNR-6.7, TIW-D-5.7 mm, and TIW-T-5.9 mm and ideal produced SSIM-1.00, CNR-13.6, SNR-34.8, TIW-D-1.5 mm, and TIW-T-2.0 mm. All constant gantry velocity reconstructions produced lower image quality metrics than ideal angular separation reconstructions regardless of the acquisition time. Motion compensated reconstruction (MCMKB) produced the highest contrast images with low streaking artifacts.

CONCLUSION

Very fast 4DCBCT scans can be acquired provided that the entire scan range is adaptively sampled, and motion-compensated reconstruction is performed. Importantly, the angular separation between x-ray projections within each individual respiratory bin had minimal effect on the image quality of fast low-dose 4DCBCT imaging. The results will assist the development of future 4DCBCT acquisition protocols that can now be achieved in very short time frames with emerging linear accelerators.

Inverse radon transform with deep learning: an application in cardiac motion correction

Objective. This paper addresses performing inverse radon transform (IRT) with artificial neural network (ANN) or deep learning, simultaneously with cardiac motion correction (MC). The suggested application domain is cardiac image reconstruction in emission or transmission tomography where IRT is relevant. Our main contribution is in proposing an ANN architecture that is particularly suitable for this purpose.Approach. We validate our approach with two types of datasets. First, we use an abstract object that looks like a heart to simulate motion-blurred radon transform. With the known ground truth in hand, we then train our proposed ANN architecture and validate its effectiveness in MC. Second, we used human cardiac gated datasets for training and validation of our approach. The gating mechanism bins data over time using the electro-cardiogram (ECG) signals for cardiac motion correction.Main results. We have shown that trained ANNs can perform motion-corrected image reconstruction directly from a motion-corrupted sinogram. We have compared our model against two other known ANN-based approaches.Significance. Our method paves the way for eliminating any need for hardware gating in medical imaging.

Future trends for patient-specific dosimetry methodology in molecular radiotherapy

Molecular radiotherapy is rapidly expanding, and new radiotherapeutics are emerging. The majority of treatments is still performed using empirical fixed activities and not tailored for individual patients. Molecular radiotherapy dosimetry is often seen as a promising candidate that would allow personalisation of treatments as outcome should ultimately depend on the absorbed doses delivered and not the activities administered. The field of molecular radiotherapy dosimetry has made considerable progress towards the feasibility of routine clinical dosimetry with reasonably accurate absorbed-dose estimates for a range of molecular radiotherapy dosimetry applications. A range of challenges remain with respect to the accurate quantification, assessment of time-integrated activity and absorbed dose estimation. In this review, we summarise a range of technological and methodological advancements, mainly focussed on beta-emitting molecular radiotherapeutics, that aim to improve molecular radiotherapy dosimetry to achieve accurate, reproducible, and streamlined dosimetry. We describe how these new technologies can potentially improve the often time-consuming considered process of dosimetry and provide suggestions as to what further developments might be required.

Method for comparison of data driven gating algorithms in emission tomography

Objective.Multiple algorithms have been proposed for data driven gating (DDG) in single photon emission computed tomography (SPECT) and have successfully been applied to myocardial perfusion imaging (MPI). Application of DDG to acquisition types other than SPECT MPI has not been demonstrated so far, as limitations and pitfalls of current methods are unknown.Approach.We create a comprehensive set of phantoms simulating the influence of different motion artifacts, view angles, moving objects, contrast, and count levels in SPECT. We perform Monte Carlo simulation of the phantoms, allowing the characterization of DDG algorithms using quantitative metrics derived from the data and evaluate the Center of Light (COL) and Laplacian Eigenmaps methods as sample DDG algorithms.Main results.View angle, object size, count rate density, and contrast influence the accuracy of both DDG methods. Moreover, the ability to extract the respiratory motion in the phantom was shown to correlate with the contrast of the moving feature to the background, the signal to noise ratio, and the noise in the data.Significance.We showed that reporting the average correlation to an external physical reference signal per acquisition is not sufficient to characterize DDG methods. Assessing DDG methods on a view-by-view basis using the simulations and metrics from this work could enable the identification of pitfalls of current methods, and extend their application to acquisitions beyond SPECT MPI.

SPECT and CT misregistration reduction in [99mTc]Tc-MAA SPECT/CT for precision liver radioembolization treatment planning

PURPOSE

Respiration and body movement induce misregistration between static [99mTc]Tc-MAA SPECT and CT, causing lung shunting fraction (LSF) and tumor-to-normal liver ratio (TNR) errors for 90Y radioembolization planning. We aim to alleviate the misregistration between [99mTc]Tc-MAA SPECT and CT using two registration schemes on simulation and clinical data.

METHODS

In the simulation study, 70 XCAT phantoms were modeled. The SIMIND Monte Carlo program and OS-EM algorithm were used for projection generation and reconstruction, respectively. Low-dose CT (LDCT) at end-inspiration was simulated for attenuation correction (AC), lungs and liver segmentation, while contrast-enhanced CT (CECT) was simulated for tumor and perfused liver segmentation. In the clinical study, 16 patient data including [99mTc]Tc-MAA SPECT/LDCT and CECT with observed SPECT and CT mismatch were analyzed. Two liver-based registration schemes were studied: SPECT registered to LDCT/CECT and vice versa. Mean count density (MCD) of different volumes-of-interest (VOIs), normalized mutual information (NMI), LSF, TNR, and maximum injected activity (MIA) based on the partition model before and after registration were compared. Wilcoxon signed-rank test was performed.

RESULTS

In the simulation study, compared to before registration, registrations significantly reduced estimation errors of MCD of all VOIs, LSF (Scheme 1: - 100.28%, Scheme 2: - 101.59%), and TNR (Scheme 1: - 7.00%, Scheme 2: - 5.67%), as well as MIA (Scheme 1: - 3.22%, Scheme 2: - 2.40%). In the clinical study, Scheme 1 reduced 33.68% LSF and increased 14.75% TNR, while Scheme 2 reduced 38.88% LSF and increased 6.28% TNR compared to before registration. One patient may change from 90Y radioembolization untreatable to treatable and other patients may change the MIA up to 25% after registration. NMI between SPECT and CT was significantly increased after registrations in both studies.

CONCLUSION

Registration between static [99mTc]Tc-MAA SPECT and corresponding CTs is feasible to reduce their spatial mismatch and improve dosimetric estimation. The improvement of LSF is larger than TNR. Our method can potentially improve patient selection and personalized treatment planning for liver radioembolization.

Dosimetric impact of 3D motion-compensated SPECT reconstruction for SIRT planning

BACKGROUND

In selective internal radiation therapy, ^99mTc SPECT images are used to optimize patient treatment planning, but they are affected by respiratory motion. In this study, we evaluated on patient data the dosimetric impact of motion-compensated SPECT reconstruction on several volumes of interest (VOI), on the tumor-to-normal liver (TN) ratio and on the activity to be injected.

METHODS

Twenty-nine patients with liver cancer or hepatic metastases treated by radioembolization were included in this study. The biodistribution of ⁹⁰Y is assumed to be the same as that of ^99mTc when predictive dosimetry is implemented. A total of 31 ^99mTc SPECT images were acquired and reconstructed with two methods: conventional OSEM (3D) and motion-compensated OSEM (3Dcomp). Seven VOI (liver, lungs, tumors, perfused liver, hepatic reserve, healthy perfused liver and healthy liver) were delineated on the CT or obtained by thresholding SPECT images followed by Boolean operations. Absorbed doses were calculated for each reconstruction using Monte Carlo simulations. Percentages of dose difference (PDD) between 3Dcomp and 3D reconstructions were estimated as well as the relative differences for TN ratio and activities to be injected. The amplitude of movement was determined with local rigid registration of the liver between the 3Dcomp reconstructions of the extreme phases of breathing.

RESULTS

The mean amplitude of the liver was 9.5 ± 2.7 mm. Medians of PDD were closed to zero for all VOI except for lungs (6.4%) which means that the motion compensation overestimates the absorbed dose to the lungs compared to the 3D reconstruction. The smallest lesions had higher PDD than the largest ones. Between 3D and 3Dcomp reconstructions, means of differences in lung dose and TN ratio were not statistically significant, but in some cases these differences exceed 1 Gy (4/31) and 8% (2/31). The absolute differences in activity were on average 3.1% ± 5.1% and can reach 22.8%.

CONCLUSION

The correction of respiratory motion mainly impacts the lung and tumor doses but only for some patients. The largest dose differences are observed for the smallest lesions.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

Technical Note: Respiratory impacts on static and respiratory gated 99m Tc-MAA SPECT/CT for liver radioembolization- A simulation study

PURPOSE

We aimed to evaluate respiratory impacts on static and respiratory gated (RG) ^99m Tc-MAA SPECT in terms of respiratory motion (RM) blur, attenuation correction (AC) and volume-of-interest (VOI) segmentation on lung shunt faction (LSF) and tumor-to-normal liver ratio (TNR) estimation for liver radioembolization therapy planning.

METHODS

The XCAT phantom was used to simulate a population of 300 phantoms, modelling various anatomical variations, tumor characteristics, respiratory motion amplitudes, LSFs and TNRs. One hundred and twenty noisy projections of average activity maps near end-expiration (End-EX) and whole respiratory cycle were simulated analytically, modeling attenuation and geometric collimator-detector-response (GCDR). The OS-EM algorithm was employed for reconstruction, modeling AC and GCDR. RM effect was evaluated for static SPECT, while AC and VOI mismatch effects were investigated independently and together for static and RG SPECT utilizing one gate, i.e., End-EX. LSF and TNR errors were measured based on the ground truth. Lesions with different characteristics were also investigated for static and RG SPECT.

RESULTS

RM overestimates LSF and underestimates TNR. The VOI mismatch caused the largest errors in both RG and static SPECT for LSF and TNR estimation, reaching 160% and -52% correspondingly with extremely mismatched VOIs for RG SPECT, even larger than those for static SPECT. With matched AC and VOIs, RG SPECT has better performance than static SPECT. Larger TNR errors are associated with tumors of smaller sizes and higher TNR for static SPECT.

CONCLUSIONS

The VOI segmentation mismatch has a stronger impact, followed by RM and AC in static ^99m Tc-MAA SPECT/CT. This effect is more pronounced for RG SPECT. When VOI masks are derived from a matched CT, RG SPECT is generally superior to static SPECT for LSF and TNR estimation. The performance of RG SPECT could be worse than static SPECT when a mismatched CT is used for segmentation. This article is protected by copyright. All rights reserved.