Advances in 4D radiation therapy for managing respiration: part I - 4D imaging

Radiation therapy for ventricular arrhythmias

Ventricular arrhythmias (VA) can be life-threatening arrhythmias that result in significant morbidity and mortality. Catheter ablation (CA) is an invasive treatment modality that can be effective in the treatment of VA where medications fail. Recurrence occurs commonly following CA due to an inability to deliver lesions of adequate depth to cauterise the electrical circuits that drive VA or reach areas of scar responsible for VA. Stereotactic body radiotherapy is a non-invasive treatment modality that allows volumetric delivery of energy to treat circuits that cannot be reached by CA. It overcomes the weaknesses of CA and has been successfully utilised in small clinical trials to treat refractory VA. This article summarises the current evidence for this novel treatment modality and the steps that will be required to bring it to the forefront of VA treatment.

Musculoskeletal CT Imaging: State-of-the-Art Advancements and Future Directions

CT is one of the most widely used modalities for musculoskeletal imaging. Recent advancements in the field include the introduction of four-dimensional CT, which captures a CT image during motion; cone-beam CT, which uses flat-panel detectors to capture the lower extremities in weight-bearing mode; and dual-energy CT, which operates at two different x-ray potentials to improve the contrast resolution to facilitate the assessment of tissue material compositions such as tophaceous gout deposits and bone marrow edema. Most recently, photon-counting CT (PCCT) has been introduced. PCCT is a technique that uses photon-counting detectors to produce an image with higher spatial and contrast resolution than conventional multidetector CT systems. In addition, postprocessing techniques such as three-dimensional printing and cinematic rendering have used CT data to improve the generation of both physical and digital anatomic models. Last, advancements in the application of artificial intelligence to CT imaging have enabled the automatic evaluation of musculoskeletal pathologies. In this review, the authors discuss the current state of the above CT technologies, their respective advantages and disadvantages, and their projected future directions for various musculoskeletal applications.

Real-time marker-less tumor tracking with TOF PET: in silico feasibility study

Purpose. Although positron emission tomography (PET) can provide a functional image of static tumors for RT guidance, it’s conventionally very challenging for PET to track a moving tumor in real-time with a multiple frame/s sampling rate. In this study, we developed a novel method to enable PET based three-dimension (3D) real-time marker-less tumor tracking (RMTT) and demonstrated its feasibility with a simulation study. Methods. For each line-of-response (LOR) acquired, its positron-electron annihilation position is calculated based on the time difference between the two gamma interactions detected by the TOF PET detectors. The accumulation of these annihilation positions from data acquired within a single sampling frame forms a coarsely measured 3D distribution of positron-emitter radiotracer uptakes of the lung tumor and other organs and tissues (background). With clinically relevant tumor size and sufficient differential radiotracer uptake concentrations between the tumor and background, the high-uptake tumor can be differentiated from the surrounding low-uptake background in the measured distribution of radiotracer uptakes. With a volume-of-interest (VOI) that closely encloses the tumor, the count-weighted centroid of the annihilation positions within the VOI can be calculated as the tumor position. All these data processes can be conducted online. The feasibility of the new method was investigated with a simulated cardiac-torso digital phantom and stationary dual-panel TOF PET detectors to track a 28 mm diameter lung tumor with a 4:1 tumor-to-background 18FDG activity concentration ratio. Results. The initial study shows TOF PET based RMTT can achieve <2.0 mm tumor tracking accuracy with 5 frame s−1 sampling rate under the simulated conditions. In comparison, using reconstructed PET images to track a similar size tumor would require >30 s acquisition time to achieve the same tracking accuracy. Conclusion. With the demonstrated feasibility, the new method may enable TOF PET based RMTT for practical RT applications.

CycN-Net: A Convolutional Neural Network Specialized for 4D CBCT Images Refinement

Four-dimensional cone-beam computed tomography (4D CBCT) has been developed to provide a sequence of phase-resolved reconstructions in image-guided radiation therapy. However, 4D CBCT images are degraded by severe streaking artifacts and noise because the phase-resolved image is an extremely sparse-view CT procedure wherein a few under-sampled projections are used for the reconstruction of each phase. Aiming at improving the overall quality of 4D CBCT images, we proposed two CNN models, named N-Net and CycN-Net, respectively, by fully excavating the inherent property of 4D CBCT. To be specific, the proposed N-Net incorporates the prior image reconstructed from entire projection data based on U-Net to boost the image quality for each phase-resolved image. Based on N-Net, a temporal correlation among the phase-resolved images is also considered by the proposed CycN-Net. Extensive experiments on both XCAT simulation data and real patient 4D CBCT datasets were carried out to verify the feasibility of the proposed CNNs. Both networks can effectively suppress streaking artifacts and noise while restoring the distinct features simultaneously, compared with the existing CNN models and two state-of-the-art iterative algorithms. Moreover, the proposed method is robust in handling complicated tasks of various patient datasets and imaging devices, which implies its excellent generalization ability.

Target motion management in breast cancer radiation therapy

Background

Over the last two decades, breast cancer remains the main cause of cancer deaths in women. To treat this type of cancer, radiation therapy (RT) has proved to be efficient. RT for breast cancer is, however, challenged by intrafractional motion caused by respiration. The problem is more severe for the left-sided breast cancer due to the proximity to the heart as an organ-at-risk. While particle therapy results in superior dose characteristics than conventional RT, due to the physics of particle interactions in the body, particle therapy is more sensitive to target motion.

Conclusions

This review highlights current and emerging strategies for the management of intrafractional target motion in breast cancer treatment with an emphasis on particle therapy, as a modern RT technique. There are major challenges associated with transferring real-time motion monitoring technologies from photon to particles beams. Surface imaging would be the dominant imaging modality for real-time intrafractional motion monitoring for breast cancer. The magnetic resonance imaging (MRI) guidance and ultra high dose rate (FLASH)-RT seem to be state-of-the-art approaches to deal with 4D RT for breast cancer.

Reduction of margin to compensate the respiratory tumor motion by the analysis of dosimetric internal target volume in lung SBRT with nonuniform volume prescription method

PURPOSE

To develop a dosimetric internal target volume (ITV) margin (DIM) for respiratory motion in lung stereotactic body radiotherapy (SBRT) and to evaluate DIM with a nonuniform volume prescription (NVP) and the point prescription (PP).

METHODS

Volumetric modulated arc therapy (VMAT) treatment plans with PP and NVP were created on a heterogeneous programmable respiratory motion phantom, with a tumor (30-mm diameter) inside a cylindrical lung insert. The tumor was defined as the gross tumor volume (GTV), equal to the clinical target volume (CTV). Five-millimeter and 0-mm margins were used for the ITV and setup margins, respectively. The phantom was moved in cranio-caudal direction with a biquadratic sinusoidal waveform with a 4-s cycle and an amplitude of ±5-10 mm. The interplay effect was evaluated by measuring the dose profile with a film in the sagittal plane for different respiratory periods and different initial respiratory phases. DIM was based on the respiratory motion amplitude that satisfied 100% and 95% coverage of the prescribed dose by the minimum dose of the CTV. Moreover, the absolute dose was measured with and without respiratory motion for NVP by a pinpoint chamber.

RESULTS

The dose difference in the tumor region due to the interplay effect was within 1.0%. The gamma passing rate was over 95.1% for different respiratory periods and 98.6% for different initial respiratory phases. DIM with PP was almost equivalent to the margin of the respiratory motion. However, DIM with NVP was 2.0 and 1.8 times larger than the margin of the respiratory motion for the 100% and 95% coverage of the prescribed doses, respectively.

CONCLUSION

The interplay effects experienced between the MLC sequence and tumor motion were negligible for NVP. The DIM analysis revealed that the margin to compensate the respiratory tumor motion could be reduced by more than 44-50% for NVP in SBRT.

Technical Note: Investigating internal-external motion correlation using fast helical CT

PURPOSE

To quantify the use of anterior torso skin surface position measurement as a breathing surrogate.

METHODS

Fourteen patients were scanned 25 times in alternating directions using a free-breathing low-mA fast helical CT protocol. Simultaneously, an abdominal pneumatic bellows was used as a real-time breathing surrogate. The imaged diaphragm dome position was used as a gold standard surrogate, characterized by localizing the most superior points of the diaphragm dome in each lung. These positions were correlated against the bellows signal acquired at the corresponding scan times. The bellows system has been shown to have a slow linear drift, and the bellows-to-CT synchronization process had a small uncertainty, so the drift and time offset were determined by maximizing the correlation coefficient between the craniocaudal diaphragm position and the drift-corrected bellows signal. The corresponding fit was used to model the real-time diaphragm position. To estimate the effectiveness of skin surface positions as surrogates, the anterior torso surface position was measured from the CT scans and correlated against the diaphragm position model. The residual error was defined as the root-mean-square correlation residual with the breathing amplitude normalized to the 5^th to 95^th breathing amplitude percentiles. The fit residual errors were analyzed over the surface for the fourteen studied patients and reported as percentages of the 5^th to 95^th percentile ranges.

RESULTS

A strong correlation was measured between the diaphragm motion and the abdominal bellows signal with an average residual error of 9.21% and standard deviation of 3.77%. In contrast, the correlations between the diaphragm position model and patient surface positions varied throughout the torso and from patient to patient. However, a consistently high correlation was found near the abdomen for each patient, and the average minimum residual error relating the skin surface to the diaphragm was 11.8% with a standard deviation of 4.61%.

CONCLUSIONS

The thoracic patient surface was found to be an accurate surrogate, but the accuracy varied across the surface sufficiently that care would need to be taken to use the surface as an accurate and reliable surrogate. Future studies will use surface imaging to determine surface patch algorithms that utilize the entire chest as well as thoracic and abdominal breathing relationships.

A preliminary investigation of re-evaluating the irradiation dose in hepatocellular carcinoma radiotherapy applying 4D CT and deformable registration

Purpose

To investigate the effect of breathing motion on dose distribution for hepatocellular carcinoma (HCC) patients using four‐dimensional (4D) CT and deformable registration.

Methods

Fifty HCC patients who were going to receive radiotherapy were enrolled in this study. All patients had been treated with transarterial chemoembolization beforehand. Three‐dimensional (3D) and 4D CT scans in free breathing were acquired sequentially. Volumetric modulated arc therapy (VMAT) was planned on the 3D CT images and maximum intensity projection (MIP) images. Thus, the 3D dose (Dose‐_3D) and MIP dose (Dose‐_MIP) were obtained, respectively. Then, the Dose‐_3D and Dose‐_MIP were recalculated on 10 phases of 4D CT images, respectively, in which the end‐inhale and end‐exhale phase doses were defined as Dose‐_3D‐EI, Dose‐_3D‐EE, Dose‐_MIP‐EI, and Dose‐_MIP‐EE. The 4D dose (Dose‐_4D‐3D and Dose‐_4D‐MIP) were obtained by deforming 10 phase doses to the end‐exhale CT to accumulate. The dosimetric difference in Dose‐_3D, Dose‐_EI3D, Dose‐_EE3D, Dose‐_4D‐3D, Dose‐_MIP, Dose‐_EIMIP, Dose‐_EEMIP, and Dose‐_4D‐MIP were compared to evaluate the motion effect on dose delivery to the planning target volume (PTV) and normal liver.

Results

Compared with Dose‐_3D, PTV D99 in Dose‐_EI3D, Dose‐_EE3D and Dose‐_4D‐3D decreased by an average of 6.02%, 1.32%, 2.43%, respectively (P < 0.05); while PTV D95 decreased by an average of 3.34%, 1.51%, 1.93%, respectively (P < 0.05). However, CI and HI of the PTV in Dose‐_3D was superior to the other three distributions (P < 0.05). There was no significant differences for the PTV between Dose‐_EI and Dose‐_EE, and between the two extreme phase doses and Dose‐_4D (P> 0.05). Negligible difference was observed for normal liver in all dose distributions (P> 0.05).

Conclusions

Four‐dimensional dose calculations potentially ensure target volume coverage when breathing motion may affect the dose distribution. Dose escalation can be considered to improve the local control of HCC on the basis of accurately predicting the probability of radiation‐induced liver disease.

Slice-stacking T2-weighted MRI for fast determination of internal target volume for liver tumor

Background

To investigate the feasibility of generating maximum intensity projection (MIP) images to determine internal target volume (ITV) using slice-stacking MRI (SS-MRI) technique.

Methods

Slice-stacking is a technique which applies a multi-slice MRI acquisition to generate a 3D MIP for ITV contouring, without reconstructing 4D-MRI. 4D digital extended cardiac-torso (XCAT) phantom was used to generate MIP images with sequential 2D HASTE sequence, with different tumor diameters (10, 30 and 50 mm) and with simulated regular and irregular (patient) breathing motions. A reference MIP was generated using all acquisition images. Consecutive repetitions were then used to generate MIP to analyze the relationship between Dice's similarity coefficient (DSC) and the number of repetitions, and the relationship between the relative ITV volume difference and the number of repetitions. Images from XCAT phantom and from three hepatic carcinoma patients were collected in this study to demonstrate the feasibility of this technique.

Results

For both regular and irregular breathing motion, the average DSC of ITV is >0.94 and the average relative ITV volume difference is <10% (approximately 0.15 cm³) when using 5 repeated scanning images to reconstruct MIP for tumor diameter of 10 mm. As tumor diameter increases, the DSC of ITV is >0.97 and the relative ITV volume difference is <5% for regular breathing motion, and the DSC of ITV is >0.97 and the relative ITV volume difference is <5.5% for irregular breathing motion when using 5 repeated scanning images to reconstruct MIP. In patient image study, the mean relative ITV volume difference is <3% and the mean DSC is 0.99 when using 5 repeated scanning images to reconstruct MIP.

Conclusions

The number of scans required to generate tumor ITV for slice-stacking method (5-7 repetition) is 3-4 times less than that of 4D-MRI (15-20 repetitions). It is feasible to generate a fast clinically acceptable ITV using slice-stacking method with sequential 2D MR images.

Protecting the Heart: A Practical Approach to Account for the Full Extent of Heart Motion in Radiation Therapy Planning

PURPOSE

Emerging evidence suggests that the heart is more radiosensitive than previously assumed; therefore, accounting for heart motion in radiation therapy planning is becoming more critical. In this study, we determined how much heart delineations based on 3-dimensional (3D) computed tomography (CT), 4-dimensional (4D) average projection (AVG), and maximum intensity projection (MIP) images should be extended to represent the full extent of heart motion during 4D imaging acquisition.

METHODS AND MATERIALS

The 3D and 4D CT scans of 10 lung cancer patients treated with stereotactic ablative radiation therapy were used. Median surfaces were derived from heart delineations of 3 observers on the 3D CT, AVG, MIP, and 25% exhale scans. Per patient, the 25% exhale contour was propagated on every phase of the 4D scan. The union of all 4D phase delineations (U4D) represented the full extent of heart motion during imaging acquisition. Surface distances from U4D to 3D, AVG, and MIP volumes were calculated. Distances in the most extreme surface points (1.5 cm most superoinferior, 10% most right/left/anteroposterior) were used to derive margins accounting only for systematic (delineation) errors.

RESULTS

Heart delineations on the MIP were the closest to the full extent of motion, requiring only ≤2.5-mm margins. Delineations on the AVG and 3D scans required margins up to 3.4 and 7.1 mm, respectively. The largest margins were for the inferior, right, and anterior aspects for the delineations on the 3D, AVG, and MIP scans, respectively.

CONCLUSION

Delineations on 3D, AVG, or MIP scans required extensions for representing the heart's full extent of motion, with the MIP requiring the smallest margins. Research including daily imaging to determine the random components for the margins and dosimetric measurements to determine the relevance of creating a planning organ at risk volume of the heart is required.

Reliability of ITV approach to varying treatment fraction time: a retrospective analysis based on 2D cine MR images

BACKGROUND

Internal Target Volume (ITV) is one of the most common strategies to passively manage tumour motion in Radiotherapy (RT). The reliability of this approach is based on the assumption that the tumour motion estimated during pre-treatment 4D Computed Tomography (CT) acquisition is representative of the motion during the whole RT treatment. With the introduction of Magnetic Resonance-guided RT (MRgRT), it has become possible to monitor tumour motion during the treatment and verify this assumption. Aim of this study was to investigate the reliability of the ITV approach with respect to the treatment fraction time (TFT) in abdominal and thoracic lesions.

METHODS

A total of 12 thoracic and 15 abdominal lesions was analysed. Before treatment, a 10-phase 4DCT was acquired and ITV margins were estimated considering the envelope of the lesion contoured on the different 4DCT phases. All patients underwent MRgRT treatment in free-breathing, monitoring the tumour position on a sagittal plane with 4 frames per second (sec). ITV margins were projected on the tumour trajectory and the percentage of treatment time in which the tumour was inside the ITV (%TT) was measured to varying of TFT. The ITV approach was considered moderately reliable when %TT ≥ 90% and strongly reliable when %TT ≥ 95%. Additional ITV margins required to achieve %TT ≥ 95% were also calculated.

RESULTS

In the analysed cohort of patients, ITV strategy can be considered strongly reliable only for lung lesions with TFT ≤ 7 min (min). The ITV strategy can be considered only moderately reliable for abdominal lesions, and additional margins are required to obtain %TT ≥ 95%. Considering a TFT ≤ 4 min, additional margins of 2 mm in cranio-caudal (CC) and 1 mm in antero-posterior (AP) are suggested for pancreatic lesions, 3 mm in CC and 2 mm in AP for renal and liver ones.

CONCLUSIONS

On the basis of the analysed cases, the ITV approach appears to be reliable in the thorax, while it results more challenging in the abdomen, due to the higher uncertainty in ITV definition and to the observed larger intra and inter-fraction motion variability. The addition of extra margins based on the TFT may represent a valid tool to compensate such limitations.

Non-rigid image registration of 4D-MRI data for improved delineation of moving tumors

Background

To increase the image quality of end-expiratory and end-inspiratory phases of retrospective respiratory self-gated 4D MRI data sets using non-rigid image registration for improved target delineation of moving tumors.

Methods

End-expiratory and end-inspiratory phases of volunteer and patient 4D MRI data sets are used as targets for non-rigid image registration of all other phases using two different registration schemes: In the first, all phases are registered directly (dir-Reg) while next neighbors are successively registered until the target is reached in the second (nn-Reg). Resulting data sets are quantitatively compared using diaphragm and tumor sharpness and the coefficient of variation of regions of interest in the lung, liver, and heart. Qualitative assessment of the patient data regarding noise level, tumor delineation, and overall image quality was performed by blinded reading based on a 4 point Likert scale.

Results

The median coefficient of variation was lower for both registration schemes compared to the target. Median dir-Reg coefficient of variation of all ROIs was 5.6% lower for expiration and 7.0% lower for inspiration compared with nn-Reg. Statistical significant differences between the two schemes were found in all comparisons. Median sharpness in inspiration is lower compared to expiration sharpness in all cases. Registered data sets were rated better compared to the targets in all categories. Over all categories, mean expiration scores were 2.92 ± 0.18 for the target, 3.19 ± 0.22 for nn-Reg and 3.56 ± 0.14 for dir-Reg and mean inspiration scores 2.25 ± 0.12 for the target, 2.72 ± 215 0.04 for nn-Reg and 3.78 ± 0.04 for dir-Reg.

Conclusions

In this work, end-expiratory and inspiratory phases of a 4D MRI data sets are used as targets for non-rigid image registration of all other phases. It is qualitatively and quantitatively shown that image quality of the targets can be significantly enhanced leading to improved target delineation of moving tumors.

High-quality initial image-guided 4D CBCT reconstruction

PURPOSE

Four-dimensional cone-beam computed tomography (4D CBCT) has been developed to provide a sequence of phase-resolved reconstructions in image-guided radiation therapy. However, 4D CBCT images are degraded by severe streaking artifacts because the 4D CBCT reconstruction process is an extreme sparse-view CT procedure wherein only under-sampled projections are used for the reconstruction of each phase. To obtain a set of 4D CBCT images achieving both high spatial and temporal resolution, we propose an algorithm by providing a high-quality initial image at the beginning of the iterative reconstruction process for each phase to guide the final reconstructed result toward its optimal solution.

METHODS

The proposed method consists of three steps to generate the initial image. First, a prior image is obtained by an iterative reconstruction method using the measured projections of the entire set of 4D CBCT images. The prior image clearly shows the appearance of structures in static regions, although it contains blurring artifacts in motion regions. Second, the robust principal component analysis (RPCA) model is adopted to extract the motion components corresponding to each phase-resolved image. Third, a set of initial images are produced by the proposed linear estimation model that combines the prior image and the RPCA-decomposed motion components. The final 4D CBCT images are derived from the simultaneous algebraic reconstruction technique (SART) equipped with the initial images. Qualitative and quantitative evaluations were performed by using two extended cardiac-torso (XCAT) phantoms and two sets of patient data. Several state-of-the-art 4D CBCT algorithms were performed for comparison to validate the performance of the proposed method.

RESULTS

The image quality of phase-resolved images is greatly improved by the proposed method in both phantom and patient studies. The results show an outstanding spatial resolution, in which streaking artifacts are suppressed to a large extent, while detailed structures such as tumors and blood vessels are well restored. Meanwhile, the proposed method depicts a high temporal resolution with a distinct respiratory motion change at different phases. For simulation phantom, quantitative evaluations of the simulation data indicate that an average of 36.72% decrease at EI phase and 42% decrease at EE phase in terms of root-mean-square error (RMSE) are achieved by our method when comparing with PICCS algorithm in Phantom 1 and Phantom 2. In addition, the proposed method has the lowest entropy and the highest normalized mutual information compared with the existing methods in simulation experiments, such as PRI, RPCA-4DCT, SMART, and PICCS. And for real patient cases, the proposed method also achieves the lowest entropy value compared with the competitive method.

CONCLUSIONS

The proposed algorithm can generate an optimal initial image to improve iterative reconstruction performance. The final sequence of phase-resolved volumes guided by the initial image achieves high spatiotemporal resolution by eliminating motion-induced artifacts. This study presents a practical 4D CBCT reconstruction method with leading image quality.

Use of 4D-CT for radiotherapy planning and reality in France: Data from a national survey

PURPOSE

Lung and some digestive tumours move during a respiratory cycle. Four-dimensional scanography (4D-CT) is commonly used in treatment planning to account for respiratory motion. Although many French radiotherapy centres are now equipped, there are no guidelines on this subject to date. We wanted to draw up a description of the use of the 4D-CT for the treatment planning in France.

METHODS AND MATERIAL

We conducted a survey in all French radiotherapy centres between March and April 2017.

RESULTS

One hundred and seventy-two were contacted. The participation rate was 88.37%. The use of the 4D-CT seems to be common and concerned planning for 15.28% of kidney and adrenal cancers, 19.72% of pancreatic cancers, 27.78% of oesophageal cancers and 73.24% of lung cancers in case of normofractionated treatments. The use of the 4D-CT was also widespread in the case of stereotactic body radiation therapy: with 61.11% in the case of pulmonary irradiation and 34.72% in the case of hepatic irradiation. Many centres declared they carried out several 4D-CT for treatment planning (29, 55% in case of stereotactic body radiation therapy for lung tumours and 20% for liver tumours). Private centres tend to repeat 4D-CT more.

CONCLUSION

Although the use of the 4D-CT appears to be developing, it remains very heterogeneous. To date, the repetition of the 4D-CT has been very poorly studied and could be the subject of clinical studies, allowing to define in which indications and for which populations there is a real benefit.

A retrospective 4D-MRI based on 2D diaphragm profiles for lung cancer patients

INTRODUCTION

4D-MRI, compared to 4D-CT, provides better soft-tissue contrast for target delineation. However, motion artefacts are often observed due to residual breathing variations. This study is to present a retrospective 4D-MRI reconstruction method based on 2D diaphragm profiles to improve the quality of 4D-MR images in the presence of significant breathing variations.

METHODS

The proposed 4D-MRI reconstruction method utilized diaphragm profiles (2D cine images on a single sagittal plan at the peak diaphragm) in conjunction with 4D-MR scans (2D-cine images on multiple pre-determined coronal planes along the anterior-posterior direction over a volume of interest). The diaphragm profile images were exploited to sort the 4D-MR scans by matching respiratory amplitude of diaphragm on the 4D-MR scans to the diaphragm profiles. To evaluate reconstructed 4D-MR images (ten 3D-MR images), sagittal images on ten 3D-MR images under free breathing (FB) and respiratory guidance (GB) were compared with diaphragm profile images (reference) from 13 healthy volunteers.

RESULTS

Forty-four 4D-MR scan datasets were successfully reconstructed without distinct respiratory-related motion artefacts even with the presence of breathing variation. The differences in diaphragm profiles between the reference and corresponding reconstructed images in the mean of root mean square were similar between FB (3.5 mm) and GB (3.0 mm), confirming that the 4D-MRI reconstruction method was effective even with significant breathing variation.

CONCLUSIONS

The diaphragm profiles were utilized to reconstruct 4D-MR images with spatial reliability and a fixed scan time under FB and GB. Our method can provide reliable 4D information of thoracic and abdominal regions for MRI-guided radiotherapy.

MRI-guidance for motion management in external beam radiotherapy: current status and future challenges

High precision conformal radiotherapy requires sophisticated imaging techniques to aid in target localisation for planning and treatment, particularly when organ motion due to respiration is involved. X-ray based imaging is a well-established standard for radiotherapy treatments. Over the last few years, the ability of magnetic resonance imaging (MRI) to provide radiation-free images with high-resolution and superb soft tissue contrast has highlighted the potential of this imaging modality for radiotherapy treatment planning and motion management. In addition, these advantageous properties motivated several recent developments towards combined MRI radiation therapy treatment units, enabling in-room MRI-guidance and treatment adaptation. The aim of this review is to provide an overview of the state-of-the-art in MRI-based image guidance for organ motion management in external beam radiotherapy. Methodological aspects of MRI for organ motion management are reviewed and their application in treatment planning, in-room guidance and adaptive radiotherapy described. Finally, a roadmap for an optimal use of MRI-guidance is highlighted and future challenges are discussed.

Evaluation of the combined use of two different respiratory monitoring systems for 4D CT simulation and gated treatment

Purpose

Two different respiratory monitoring systems (Varian's Real‐Time Position Management (RPM) System and Siemens’ ANZAI belt Respiratory Gating System) are compared in the context of respiratory signals and 4D CT images that are accordingly reconstructed. This study aims to evaluate the feasibility of combined use of RPM and ANZAI systems for 4DCT simulation and gated radiotherapy treatment, respectively.

Methods

The RPM infrared reflecting marker and the ANZAI belt pressure sensor were both placed on the patient's abdomen during 4DCT scans. The respiratory signal collected by the two systems was synchronized. Fifteen patients were enrolled for respiratory signal collection and analysis. The discrepancies between the RPM and ANZAI traces can be characterized by phase shift and shape distortion. To reveal the impact of the changes in respiratory signals on 4D images, two sets of 4D images based on the same patient's raw data were reconstructed using the RPM and ANZAI data for phase sorting, respectively. The volume of whole lung and the position of diaphragm apex were measured and compared for each respiratory phase.

Results

The mean phase shift was measured as 0.2 ± 0.1 s averaged over 15 patients. The shape of the breathing trace was found to be in disagreement. For all the patients, the ANZAI trace had a steeper falloff in exhalation than RPM. The inhalation curve, however, was matched for nine patients, steeper in ANZAI for five patients and steeper in RPM for one patient. For 4D image comparison, the difference in whole‐lung volume was about −4% to +4% and the difference in diaphragm position was about −5 mm to +4 mm, compared in each individual phase and averaged over seven patients.

Conclusions

Combined use of one system for 4D CT simulation and the other for gated treatment should be avoided as the resultant gating window would not fully match with each other due to the remarkable discrepancy in breathing traces acquired by the two different surrogate systems.

Impact of deformable image registration on dose accumulation applied electrocardiograph-gated 4DCT in the heart and left ventricular myocardium during esophageal cancer radiotherapy

Background

The deformable image registration (DIR) technique has the potential to realize the dose accumulation during radiotherapy. This study will analyze the feasibility of evaluating dose-volume parameters for the heart and left ventricular myocardium (LVM) by applying DIR.

Methods

The electrocardiograph-gated four-dimensional CT (ECG-gated 4DCT) data of 21 patients were analyzed retrospectively. The heart and LVM were contoured on 20 phases of 4DCT (0%, 5%,…,95%). The heart and LVM in the minimum volume/dice similarity coefficient (DSC) phase (Volume _min/DSC _min) were deformed to the maximum volume/DSC phase (Volume _max/ DSC _max), which used the intensity-based free-form DIR algorithm of MIM software. The dose was deformed according to the deformation vector. The variations in volume, mean dose (D_mean), V₂₀, V₃₀ and V₄₀ for the heart and LVM before and after DIR were compared, and the reference phase was the Volume _max/DSC _max phase.

Results

For the heart, the difference between the pre- and post-registration Volume _min and Volume _max were reduced from 13.87 to 1.72%; the DSC was increased from 0.899 to 0.950 between the pre- and post-registration DSC _min phase relative to the DSC _max phase. The post-registration D_mean, V₂₀, V₃₀ and V₄₀ of the heart were statistically significant compared to those in the Volume _max/DSC _max phase (p < 0.05). For the LVM, the difference between the pre- and post-registration Volume _min and Volume _max were only reduced from 18.77 to 17.38%; the DSC reached only 0.733 in the post-registration DSC _min phase relative to the DSC _max phase. The pre- and post-registration volume, D_mean, V₂₀, V₃₀ and V₄₀ of the LVM were all statistically significant compared to those in the Volume _max/DSC _max phase (p < 0.05).

Conclusions

There was no significant relationship between the variation in dose-volume parameters and the variation in the volume and morphology for the heart; however, the inconsistency of the variation in the volume and morphology for the LVM was a major factor that led to uncertainty in the dose-volume evaluation. In addition, the individualized local deformation registration technology should be applied in dose accumulation for the heart and LVM.

Ultrasound-driven 4D MRI

We present an ultrasound-driven 4D magnetic resonance imaging (US-4DMRI) method for respiratory motion imaging in the thorax and abdomen. The proposed US-4DMRI comes along with a high temporal resolution, and allows for organ motion imaging beyond a single respiratory cycle. With the availability of the US surrogate both inside and outside the MR bore, 4D MR images can be reconstructed for 4D treatment planning and online respiratory motion prediction during radiotherapy. US-4DMRI relies on simultaneously acquired 2D liver US images and abdominal 2D MR multi-slice scans under free respiration. MR volumes are retrospectively composed by grouping the MR slices corresponding to the most similar US images. We present two different US similarity metrics: an intensity-based approach, and a similarity measure relying on predefined fiducials which are being tracked over time. The proposed method is demonstrated on MR liver scans of eight volunteers acquired over a duration of 5.5 min each at a temporal resolution of 2.6 Hz with synchronous US imaging at 14 Hz-17 Hz. Visual inspection of the reconstructed MR volumes revealed satisfactory results in terms of continuity in organ boundaries and blood vessels. In quantitative leave-one-out experiments, both US similarity metrics reach the performance level of state-of-the-art navigator-based approaches.

Quantification of variation in dose-volume parameters for the heart, pericardium and left ventricular myocardium during thoracic tumor radiotherapy

Cardiac activity can induce dose-volume evaluation errors for cardiac structures. The purpose of this study was to quantify the variation in dose-volume parameters for the heart, pericardium and left ventricular myocardium (LVM) throughout the cardiac circle. The heart, pericardium and LVM of 22 patients were contoured on 20 phases of electrocardiography-gated 4D computed tomography (4DCT) images acquired during breath-hold. Radiotherapy plans were designed on 0% phase of the 4DCT images, and the dose distributions of the plans were imported into MIM Maestro and deformed to each phase to generate distributions for all phases. Variations in dose-volume parameters for the heart, pericardium and LVM were compared among different phases. The rates of variation in Dmean for the heart and pericardium were 3.33 ± 1.04% and 2.66 ± 1.15%, respectively. The mean values of the maximum difference in V5, V10, V20, V30 and V40 were all <2% for the heart and pericardium and were not statistically significant (P > 0.05). The rate of variation in Dmean for the LVM reached 87.05 ± 38.34%, and the maximum differences in V5, V10, V20, V30 and V40 were 13.76 ± 4.46%, 13.64 ± 4.33%, 12.84 ± 4.55%, 11.62 ± 4.85% and 3.63 ± 2.56%, respectively; all differences were statistically significant (P < 0.05). Variations in dose-volume parameters were more significant in the LVM than in the heart and pericardium (P < 0.05). The dose-volume parameters for the LVM were significantly influenced by cardiac activity, whereas those for the heart and pericardium were not; therefore, individual dosimetric evaluation and limitation must be performed for the LVM.

Quantification of heart, pericardium, and left ventricular myocardium movements during the cardiac cycle for thoracic tumor radiotherapy

Purpose

The purpose of this study was to quantify variations in the heart, pericardium, and left ventricular myocardium (LVM) caused by cardiac movement using the breath-hold technique.

Patients and methods

In this study, the electrocardiography-gated four-dimensional computed tomography (CT) images of 22 patients were analyzed, which were sorted into 20 phases (0-95%) according to the cardiac cycle. The heart, pericardium, and LVM were contoured on each phase of the CT images. The positions, volume, dice similarity coefficient (DSC) in reference to 0% phase, and morphological parameters (max 3D diameter, roundness, spherical disproportion, sphericity, and surface area) in different phases of the heart, pericardium, and LVM were analyzed, which were presented as mean ± standard deviation.

Results

The mean values of displacements along the X, Y, and Z axes respectively were as follows: 1.2 mm, 0.6 mm, and 0.6 mm for the heart; 0.5 mm, 0.4 mm, and 0.8 mm for the pericardium; and 1.0 mm, 4.1 mm, and 1.9 mm for the LVM. The maximum variations in volume and DSC respectively were 16.49%±3.85% and 10.08%±2.14% for the heart, 12.62%±3.94% and 5.20%±1.54% for the pericardium, and 24.23%±11.35% and 184.33%±128.61% for the LVM. The differences in the morphological parameters between the maximum and minimum DSC phases for the heart and pericardium were not significantly different (p>0.05) but were significantly different for the LVM (p<0.05).

Conclusion

The volumetric and morphological variations of the heart were similar to those of pericardium, and all were significantly smaller than those of the LVM. This inconsistency in the volumetric and morphological variations between the LVM and the heart and pericardium indicates that special protection of the LVM should be considered.

Advanced Radiation DOSimetry phantom (ARDOS): a versatile breathing phantom for 4D radiation therapy and medical imaging

A novel breathing phantom was designed for being used in conventional and ion-beam radiotherapy as well as for medical imaging. Accurate dose delivery and patient safety are aimed to be verified for four-dimensional (4D) treatment techniques compensating for breathing-induced tumor motion. The phantom includes anthropomorphic components representing an average human thorax. It consists of real tissue equivalent materials to fulfill the requirements for dosimetric experiments and imaging purposes. The different parts of the torso (lungs, chest wall, and ribs) and the tumor can move independently. Simple regular movements, as well as more advanced patient-specific breathing cycles are feasible while a reproducible setup can be guaranteed. The phantom provides the flexibility to use different types of dosimetric devices and was designed in a way that it is robust, transportable and easy to handle. Tolerance levels and the reliability of the phantom setup were determined in combination with tests on motion accuracy and reproducibility by using infrared optical tracking technology. Different imaging was performed including positron emission tomography imaging, 4D computed tomography as well as real-time in-room imaging. The initial dosimetric benchmarking studies were performed in a photon beam where dose parameters are predictable and the dosimetric procedures well established.

A review of GPU-based medical image reconstruction

Tomographic image reconstruction is a computationally demanding task, even more so when advanced models are used to describe a more complete and accurate picture of the image formation process. Such advanced modeling and reconstruction algorithms can lead to better images, often with less dose, but at the price of long calculation times that are hardly compatible with clinical workflows. Fortunately, reconstruction tasks can often be executed advantageously on Graphics Processing Units (GPUs), which are exploited as massively parallel computational engines. This review paper focuses on recent developments made in GPU-based medical image reconstruction, from a CT, PET, SPECT, MRI and US perspective. Strategies and approaches to get the most out of GPUs in image reconstruction are presented as well as innovative applications arising from an increased computing capacity. The future of GPU-based image reconstruction is also envisioned, based on current trends in high-performance computing.

Development of a Computerized 4-D MRI Phantom for Liver Motion Study

Purpose:

To develop a 4-dimensional computerized magnetic resonance imaging phantom with image textures extracted from real patient scans for liver motion studies.

Methods:

The proposed phantom was developed based on the current version of 4-dimensional extended cardiac-torso computerized phantom and a clinical magnetic resonance scan. Initially, the extended cardiac-torso phantom was voxelized in abdominal–chest region at the end of exhalation phase. Structures/tissues were classified into 4 categories: (1) Seven key textured organs, including liver, gallbladder, spleen, stomach, heart, kidneys, and pancreas, were mapped from a clinical T1-weighted liver magnetic resonance scan using deformable registration. (2) Large textured soft tissue volumes were simulated via an iterative pattern generation method using the same magnetic resonance scan. (3) Lung and intestine structures were generated by assigning uniform intensity with proper noise modeling. (4) Bony structures were generated by assigning the magnetic resonance values. A spherical hypointensity tumor was inserted into the liver. Other respiratory phases of the 4-dimensional phantom were generated using the backward deformation vector fields exported by the extended cardiac-torso program, except that bony structures were generated separately for each phase. A weighted image filtering process was utilized to improve the overall tissue smoothness at each phase.

Results:

Three 4-dimensional series with different motion amplitudes were generated. The developed motion phantom produced good illustrations of abdominal–chest region with anatomical structures in key organs and texture patterns in large soft tissue volumes. In a standard series, the tumor volume was measured as 13.90 ± 0.11 cm³ in a respiratory cycle and the tumor’s maximum center-of-mass shift was 2.95 cm/1.84 cm on superior–inferior/anterior–posterior directions. The organ motion during the respiratory cycle was well rendered. The developed motion phantom has the flexibility of motion pattern variation, organ geometry variation, and tumor modeling variation.

Conclusions:

A 4-D computerized phantom was developed and could be used to produce image series with synthetic magnetic resonance textures for magnetic resonance imaging research of liver motion.

Motion management strategies and technical issues associated with stereotactic body radiotherapy of thoracic and upper abdominal tumors: A review from NRG oncology

The efficacy of stereotactic body radiotherapy (SBRT) has been well demonstrated. However, it presents unique challenges for accurate planning and delivery especially in the lungs and upper abdomen where respiratory motion can be significantly confounding accurate targeting and avoidance of normal tissues. In this paper, we review the current literature on SBRT for lung and upper abdominal tumors with particular emphasis on addressing respiratory motion and its affects. We provide recommendations on strategies to manage motion for different, patient-specific situations. Some of the recommendations will potentially be adopted to guide clinical trial protocols.

Four-dimensional diffusion-weighted MR imaging (4D-DWI): a feasibility study

PURPOSE

Diffusion-weighted Magnetic Resonance Imaging (DWI) has been shown to be a powerful tool for cancer detection with high tumor-to-tissue contrast. This study aims to investigate the feasibility of developing a four-dimensional DWI technique (4D-DWI) for imaging respiratory motion for radiation therapy applications.

MATERIALS/METHODS

Image acquisition was performed by repeatedly imaging a volume of interest (VOI) using an interleaved multislice single-shot echo-planar imaging (EPI) 2D-DWI sequence in the axial plane. Each 2D-DWI image was acquired with an intermediately low b-value (b = 500 s/mm² ) and with diffusion-encoding gradients in x, y, and z diffusion directions. Respiratory motion was simultaneously recorded using a respiratory bellow, and the synchronized respiratory signal was used to retrospectively sort the 2D images to generate 4D-DWI. Cine MRI using steady-state free precession was also acquired as a motion reference. As a preliminary feasibility study, this technique was implemented on a 4D digital human phantom (XCAT) with a simulated pancreas tumor. The respiratory motion of the phantom was controlled by regular sinusoidal motion profile. 4D-DWI tumor motion trajectories were extracted and compared with the input breathing curve. The mean absolute amplitude differences (D) were calculated in superior-inferior (SI) direction and anterior-posterior (AP) direction. The technique was then evaluated on two healthy volunteers. Finally, the effects of 4D-DWI on apparent diffusion coefficient (ADC) measurements were investigated for hypothetical heterogeneous tumors via simulations.

RESULTS

Tumor trajectories extracted from XCAT 4D-DWI were consistent with the input signal: the average D value was 1.9 mm (SI) and 0.4 mm (AP). The average D value was 2.6 mm (SI) and 1.7 mm (AP) for the two healthy volunteers.

CONCLUSION

A 4D-DWI technique has been developed and evaluated on digital phantom and human subjects. 4D-DWI can lead to more accurate respiratory motion measurement. This has a great potential to improve the visualization and delineation of cancer tumors for radiotherapy.

Accuracy of respiratory motion measurement of 4D-MRI: A comparison between cine and sequential acquisition

PURPOSE

The authors have recently developed a cine-mode T2*/T1-weighted 4D-MRI technique and a sequential-mode T2-weighted 4D-MRI technique for imaging respiratory motion. This study aims at investigating which 4D-MRI image acquisition mode, cine or sequential, provides more accurate measurement of organ motion during respiration.

METHODS

A 4D digital extended cardiac-torso (XCAT) human phantom with a hypothesized tumor was used to simulate the image acquisition and the 4D-MRI reconstruction. The respiratory motion was controlled by the given breathing signal profiles. The tumor was manipulated to move continuously with the surrounding tissue. The motion trajectories were measured from both sequential- and cine-mode 4D-MRI images. The measured trajectories were compared with the average trajectory calculated from the input profiles, which was used as references. The error in 4D-MRI tumor motion trajectory (E) was determined. In addition, the corresponding respiratory motion amplitudes of all the selected 2D images for 4D reconstruction were recorded. Each of the amplitude was compared with the amplitude of its associated bin on the average breathing curve. The mean differences from the average breathing curve across all slice positions (D) were calculated. A total of 500 simulated respiratory profiles with a wide range of irregularity (Ir) were used to investigate the relationship between D and Ir. Furthermore, statistical analysis of E and D using XCAT controlled by 20 cancer patients' breathing profiles was conducted. Wilcoxon Signed Rank test was conducted to compare two modes.

RESULTS

D increased faster for cine-mode (D = 1.17 × Ir + 0.23) than sequential-mode (D = 0.47 × Ir + 0.23) as irregularity increased. For the XCAT study using 20 cancer patients' breathing profiles, the median E values were significantly different: 0.12 and 0.10 cm for cine- and sequential-modes, respectively, with a p-value of 0.02. The median D values were significantly different: 0.47 and 0.24 cm for cine- and sequential-modes, respectively, with a p-value < 0.001.

CONCLUSIONS

Respiratory motion measurement may be more accurate and less susceptible to breathing irregularity in sequential-mode 4D-MRI than that in cine-mode 4D-MRI.

Evaluation of gantry speed on image quality and imaging dose for 4D cone-beam CT acquisition

Background

This study investigates the effect of gantry speed on 4DCBCT image quality and dose for the Varian On-Board Imager®.

Methods

A thoracic 4DCBCT protocol was designed using a 125 kVp spectrum. Image quality parameters were evaluated for 4DCBCT acquisition using Catphan® phantom with real-time position management™ system for gantry speeds varying between 1.0 to 6.0°/s. Superior-inferior motion of the phantom was executed using a sinusoidal waveform with five second period. Scans were retrospectively sorted into 4 phases (CBCT-4 ph) and 10 phases (CBCT-10 ph); average 4DCBCT (CBCT-ave), using all image data from the 4DCBCT acquisitions was also evaluated. The 4DCBCT images were evaluated using the following image quality metrics: spatial resolution, contrast-to-noise ratio (CNR), and uniformity index (UI). Additionally, Hounsfield unit (HU) sensitivity compared to a baseline CBCT and percent differences and RMS errors (RMSE) of excursion were also determined. Imaging dose was evaluated using an IBA CC13 ion chamber placed within CIRS Thorax phantom using the same sinusoidal motion and image acquisition settings as mentioned above.

Results

Spatial resolution decreased linearly from 5.93 to 3.82 lp/cm as gantry speed increased from 1.0 to 6.0°/s. CNR decreased linearly from 4.80 to 1.82 with gantry speed increasing from 1.0 to 6.0°/s, respectively. No noteworthy variations in UI, HU sensitivity, or excursion metrics were observed with changes in gantry speed. Ion chamber dose rates measured ranged from 2.30 (lung) to 5.18 (bone) E-3 cGy/mAs.

Conclusions

A quantitative analysis of the Varian OBI’s 4DCBCT capabilities was explored. Changing gantry speed changes the number of projections used for reconstruction, affecting both image quality and imaging dose if x-ray tube current is held constant. From the results of this study, a gantry speed between 2 and 3°/s was optimal when considering image quality, dose, and reconstruction time. The future of 4DCBCT clinical utility relies on further investigation of image acquisition and reconstruction optimization.

4D respiratory motion-compensated image reconstruction of free-breathing radial MR data with very high undersampling

PURPOSE

To develop four-dimensional (4D) respiratory time-resolved MRI based on free-breathing acquisition of radial MR data with very high undersampling.

METHODS

We propose the 4D joint motion-compensated high-dimensional total variation (4D joint MoCo-HDTV) algorithm, which alternates between motion-compensated image reconstruction and artifact-robust motion estimation at multiple resolution levels. The algorithm is applied to radial MR data of the thorax and upper abdomen of 12 free-breathing subjects with acquisition times between 37 and 41 s and undersampling factors of 16.8. Resulting images are compared with compressed sensing-based 4D motion-adaptive spatio-temporal regularization (MASTeR) and 4D high-dimensional total variation (HDTV) reconstructions.

RESULTS

For all subjects, 4D joint MoCo-HDTV achieves higher similarity in terms of normalized mutual information and cross-correlation than 4D MASTeR and 4D HDTV when compared with reference 4D gated gridding reconstructions with 8.4 ± 1.1 times longer acquisition times. In a qualitative assessment of artifact level and image sharpness by two radiologists, 4D joint MoCo-HDTV reveals higher scores (P < 0.05) than 4D HDTV and 4D MASTeR at the same undersampling factor and the reference 4D gated gridding reconstructions, respectively.

CONCLUSIONS

4D joint MoCo-HDTV enables time-resolved image reconstruction of free-breathing radial MR data with undersampling factors of 16.8 while achieving low-streak artifact levels and high image sharpness. Magn Reson Med 77:1170-1183, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Four dimensional magnetic resonance imaging with retrospective k-space reordering: a feasibility study

PURPOSE

Current four dimensional magnetic resonance imaging (4D-MRI) techniques lack sufficient temporal/spatial resolution and consistent tumor contrast. To overcome these limitations, this study presents the development and initial evaluation of a new strategy for 4D-MRI which is based on retrospective k-space reordering.

METHODS

We simulated a k-space reordered 4D-MRI on a 4D digital extended cardiac-torso (XCAT) human phantom. A 2D echo planar imaging MRI sequence [frame rate (F) = 0.448 Hz; image resolution (R) = 256 × 256; number of k-space segments (NKS) = 4] with sequential image acquisition mode was assumed for the simulation. Image quality of the simulated "4D-MRI" acquired from the XCAT phantom was qualitatively evaluated, and tumor motion trajectories were compared to input signals. In particular, mean absolute amplitude differences (D) and cross correlation coefficients (CC) were calculated. Furthermore, to evaluate the data sufficient condition for the new 4D-MRI technique, a comprehensive simulation study was performed using 30 cancer patients' respiratory profiles to study the relationships between data completeness (Cp) and a number of impacting factors: the number of repeated scans (NR), number of slices (NS), number of respiratory phase bins (NP), NKS, F, R, and initial respiratory phase at image acquisition (P0). As a proof-of-concept, we implemented the proposed k-space reordering 4D-MRI technique on a T2-weighted fast spin echo MR sequence and tested it on a healthy volunteer.

RESULTS

The simulated 4D-MRI acquired from the XCAT phantom matched closely to the original XCAT images. Tumor motion trajectories measured from the simulated 4D-MRI matched well with input signals (D = 0.83 and 0.83 mm, and CC = 0.998 and 0.992 in superior-inferior and anterior-posterior directions, respectively). The relationship between Cp and NR was found best represented by an exponential function (CP=1001-e(-0.18NR) , when NS = 30, NP = 6). At a CP value of 95%, the relative error in tumor volume was 0.66%, indicating that NR at a CP value of 95% (NR,95%) is sufficient. It was found that NR,95% is approximately linearly proportional to NP (r = 0.99), and nearly independent of all other factors. The 4D-MRI images of the healthy volunteer clearly demonstrated respiratory motion in the diaphragm region with minimal motion induced noise or aliasing.

CONCLUSIONS

It is feasible to generate respiratory correlated 4D-MRI by retrospectively reordering k-space based on respiratory phase. This new technology may lead to the next generation 4D-MRI with high spatiotemporal resolution and optimal tumor contrast, holding great promises to improve the motion management in radiotherapy of mobile cancers.

T2-weighted four dimensional magnetic resonance imaging with result-driven phase sorting

PURPOSE

T2-weighted MRI provides excellent tumor-to-tissue contrast for target volume delineation in radiation therapy treatment planning. This study aims at developing a novel T2-weighted retrospective four dimensional magnetic resonance imaging (4D-MRI) phase sorting technique for imaging organ/tumor respiratory motion.

METHODS

A 2D fast T2-weighted half-Fourier acquisition single-shot turbo spin-echo MR sequence was used for image acquisition of 4D-MRI, with a frame rate of 2-3 frames/s. Respiratory motion was measured using an external breathing monitoring device. A phase sorting method was developed to sort the images by their corresponding respiratory phases. Besides, a result-driven strategy was applied to effectively utilize redundant images in the case when multiple images were allocated to a bin. This strategy, selecting the image with minimal amplitude error, will generate the most representative 4D-MRI. Since we are using a different image acquisition mode for 4D imaging (the sequential image acquisition scheme) with the conventionally used cine or helical image acquisition scheme, the 4D dataset sufficient condition was not obviously and directly predictable. An important challenge of the proposed technique was to determine the number of repeated scans (NR) required to obtain sufficient phase information at each slice position. To tackle this challenge, the authors first conducted computer simulations using real-time position management respiratory signals of the 29 cancer patients under an IRB-approved retrospective study to derive the relationships between NR and the following factors: number of slices (NS), number of 4D-MRI respiratory bins (NB), and starting phase at image acquisition (P0). To validate the authors' technique, 4D-MRI acquisition and reconstruction were simulated on a 4D digital extended cardiac-torso (XCAT) human phantom using simulation derived parameters. Twelve healthy volunteers were involved in an IRB-approved study to investigate the feasibility of this technique.

RESULTS

4D data acquisition completeness (Cp) increases as NR increases in an inverse-exponential fashion (Cp = 100 - 99 × exp(-0.18 × NR), when NB = 6, fitted using 29 patients' data). The NR required for 4D-MRI reconstruction (defined as achieving 95% completeness, Cp = 95%, NR = NR,95) is proportional to NB (NR,95 ∼ 2.86 × NB, r = 1.0), but independent of NS and P0. Simulated XCAT 4D-MRI showed a clear pattern of respiratory motion. Tumor motion trajectories measured on 4D-MRI were comparable to the average input signal, with a mean relative amplitude error of 2.7% ± 2.9%. Reconstructed 4D-MRI for healthy volunteers illustrated clear respiratory motion on three orthogonal planes, with minimal image artifacts. The artifacts were presumably caused by breathing irregularity and incompleteness of data acquisition (95% acquired only). The mean relative amplitude error between critical structure trajectory and average breathing curve for 12 healthy volunteers is 2.5 ± 0.3 mm in superior-inferior direction.

CONCLUSIONS

A novel T2-weighted retrospective phase sorting 4D-MRI technique has been developed and successfully applied on digital phantom and healthy volunteers.

Advances in 4D treatment planning for scanned particle beam therapy - report of dedicated workshops

We report on recent progress in the field of mobile tumor treatment with scanned particle beams, as discussed in the latest editions of the 4D treatment planning workshop. The workshop series started in 2009, with about 20 people from 4 research institutes involved, all actively working on particle therapy delivery and development. The first workshop resulted in a summary of recommendations for the treatment of mobile targets, along with a list of requirements to apply these guidelines clinically. The increased interest in the treatment of mobile tumors led to a continuously growing number of attendees: the 2012 edition counted more than 60 participants from 20 institutions and commercial vendors. The focus of research discussions among workshop participants progressively moved from 4D treatment planning to complete 4D treatments, aiming at effective and safe treatment delivery. Current research perspectives on 4D treatments include all critical aspects of time resolved delivery, such as in-room imaging, motion detection, beam application, and quality assurance techniques. This was motivated by the start of first clinical treatments of hepato cellular tumors with a scanned particle beam, relying on gating or abdominal compression for motion mitigation. Up to date research activities emphasize significant efforts in investigating advanced motion mitigation techniques, with a specific interest in the development of dedicated tools for experimental validation. Potential improvements will be made possible in the near future through 4D optimized treatment plans that require upgrades of the currently established therapy control systems for time resolved delivery. But since also these novel optimization techniques rely on the validity of the 4DCT, research focusing on alternative 4D imaging technique, such as MRI based 4DCT generation will continue.

Feasibility and potential benefits of defining the internal gross tumor volume of hepatocellular carcinoma using contrast-enhanced 4D CT images obtained by deformable registration

Objective

To study the feasibility and the potential benefits of defining the internal gross tumor volume (IGTV) of hepatocellular carcinoma (HCC) using contrast-enhanced 4D CT images obtained by combining arterial-phase (AP) contrast-enhanced (CE) 3D CT and non-contrast-enhanced (NCE) 4D CT images using deformable registration (DR).

Methods

Ten HCC patients who had received radiotherapy beforehand were selected for this study. The following CT simulation images were acquired sequentially: NCE 4D CT in free breathing, NCE 3D CT and APCE 3D CT in end-expiration breath holding. All 4D CT images were sorted into ten phases according to breath cycle (CT₀₀ ~ CT₉₀). Gross tumor volumes (GTVs) were contoured on all CT images and the IGTV_-1 was obtained by merging the GTVs in each phase of 4D CT images. The GTV on the APCE 3D CT image was deformably registered to each 4D CT phase image according to liver shape using RayStation^TM 3.99.0.7 version treatment planning system. The IGTV_-DR was obtained by merging the GTVs after DR on the 4D CT images. Volume differences among the GTVs and between the IGTV_-1 and the IGTV_-DR were compared.

Results

The edge of most lesions could be definitively identified using APCE 3D CT images compared to NCE 4D and 3D CT images. The GTV volume on APCE 3D CT images increased by an average of 34.79% (P < 0.05). There was no significant difference among the GTV volumes obtained using NCE 4D and 3D CT images (P > 0.05). The GTV volumes after DR on 4D CT different phase images increased by an average of 36.29% (P < 0.05), as was observed using the APCE 3D CT image (P > 0.05). Lastly, the volume of IGTV_-DR increased by an average of 19.91% compared to that of IGTV_-1 (P < 0.05).

Conclusion

NCE 4D CT imaging alone has the potential risk of missing a partial volume of the HCC. The combination of APCE 3D CT and NCE 4D CT images using the DR technique improved the accuracy of the definition of the IGTV in HCC.

Four-dimensional computed tomography based respiratory-gated radiotherapy with respiratory guidance system: analysis of respiratory signals and dosimetric comparison

PURPOSE

To investigate the effectiveness of respiratory guidance system in 4-dimensional computed tomography (4 DCT) based respiratory-gated radiation therapy (RGRT) by comparing respiratory signals and dosimetric analysis of treatment plans.

METHODS

The respiratory amplitude and period of the free, the audio device-guided, and the complex system-guided breathing were evaluated in eleven patients with lung or liver cancers. The dosimetric parameters were assessed by comparing free breathing CT plan and 4 DCT-based 30-70% maximal intensity projection (MIP) plan.

RESULTS

The use of complex system-guided breathing showed significantly less variation in respiratory amplitude and period compared to the free or audio-guided breathing regarding the root mean square errors (RMSE) of full inspiration (P = 0.031), full expiration (P = 0.007), and period (P = 0.007). The dosimetric parameters including V(5 Gy), V(10 Gy), V(20 Gy), V(30 Gy), V(40 Gy), and V(50 Gy) of normal liver or lung in 4 DCT MIP plan were superior over free breathing CT plan.

CONCLUSIONS

The reproducibility and regularity of respiratory amplitude and period were significantly improved with the complex system-guided breathing compared to the free or the audio-guided breathing. In addition, the treatment plan based on the 4D CT-based MIP images acquired with the complex system guided breathing showed better normal tissue sparing than that on the free breathing CT.

Estimation of lung motion fields in 4D CT data by variational non-linear intensity-based registration: A comparison and evaluation study

Accurate and robust estimation of motion fields in respiration-correlated CT (4D CT) images, usually performed by non-linear registration of the temporal CT frames, is a precondition for the analysis of patient-specific breathing dynamics and subsequent image-supported diagnostics and treatment planning. In this work, we present a comprehensive comparison and evaluation study of non-linear registration variants applied to the task of lung motion estimation in thoracic 4D CT data. In contrast to existing multi-institutional comparison studies (e.g. MIDRAS and EMPIRE10), we focus on the specific but common class of variational intensity-based non-parametric registration and analyze the impact of the different main building blocks of the underlying optimization problem: the distance measure to be minimized, the regularization approach and the transformation space considered during optimization. In total, 90 different combinations of building block instances are compared. Evaluated on proprietary and publicly accessible 4D CT images, landmark-based registration errors (TRE) between 1.14 and 1.20 mm for the most accurate registration variants demonstrate competitive performance of the applied general registration framework compared to other state-of-the-art approaches for lung CT registration. Although some specific trends can be observed, effects of interchanging individual instances of the building blocks on the TRE are in general rather small (no single outstanding registration variant existing); the same level of accuracy is, however, associated with significantly different degrees of motion field smoothness and computational demands. Consequently, the building block combination of choice will depend on application-specific requirements on motion field characteristics.

VMAT techniques for lymph node-positive left sided breast cancer

PURPOSE

To investigate the plan quality of two different volumetric modulated arc therapy (VMAT) techniques for lymph node-positive left-sided breast cancer.

METHODS

Two VMAT plans were generated for 10 lymph node-positive left-sided breast cancer patients: one plan using one single segment of a full rotation, typically an arc segment of 230° (1s-VMAT); and a second plan consisting of 2 small tangential arc segments of about 50° (2s-VMAT). For plan comparison, various dose and dose volume metrics (Dmean, D98%, D2% for target volumes, D2%, Dmean and Vx% for organs at risk (OAR)) were evaluated.

RESULTS

Both techniques fulfilled both clinical target dose and OAR goals. 1s-VMAT achieved a slightly better homogeneity and better target coverage (D2%= 54.2 ± 0.7 Gy, D98%= 30.3 ± 1.8 Gy) compared to 2s-VMAT (D2%= 55.0 ± 1.1 Gy, D98%= 29.9 ± 1.7 Gy). For geometrical reasons, OAR sparing was noticeable but not significant better using 2s-VMAT, particularly heart and contralateral breast. The heart received a mean dose of 4.4 ± 0.8 Gy using 1s-VMAT and 3.3 ± 1.0 Gy using 2s-VMAT; the contralateral breast received 1.5 ± 0.3 Gy and 0.9 ± 0.3 Gy, respectively.

CONCLUSIONS

A VMAT technique based on two small tangential arc segments enables improved OAR sparing; the differences between the two techniques in target coverage and homogeneity are minor. Patient age and -anatomy must be considered for each individual case when deciding which technique to be used.

Actively triggered 4d cone-beam CT acquisition

PURPOSE

4d cone-beam computed tomography (CBCT) scans are usually reconstructed by extracting the motion information from the 2d projections or an external surrogate signal, and binning the individual projections into multiple respiratory phases. In this "after-the-fact" binning approach, however, projections are unevenly distributed over respiratory phases resulting in inefficient utilization of imaging dose. To avoid excess dose in certain respiratory phases, and poor image quality due to a lack of projections in others, the authors have developed a novel 4d CBCT acquisition framework which actively triggers 2d projections based on the forward-predicted position of the tumor.

METHODS

The forward-prediction of the tumor position was independently established using either (i) an electromagnetic (EM) tracking system based on implanted EM-transponders which act as a surrogate for the tumor position, or (ii) an external motion sensor measuring the chest-wall displacement and correlating this external motion to the phase-shifted diaphragm motion derived from the acquired images. In order to avoid EM-induced artifacts in the imaging detector, the authors devised a simple but effective "Faraday" shielding cage. The authors demonstrated the feasibility of their acquisition strategy by scanning an anthropomorphic lung phantom moving on 1d or 2d sinusoidal trajectories.

RESULTS

With both tumor position devices, the authors were able to acquire 4d CBCTs free of motion blurring. For scans based on the EM tracking system, reconstruction artifacts stemming from the presence of the EM-array and the EM-transponders were greatly reduced using newly developed correction algorithms. By tuning the imaging frequency independently for each respiratory phase prior to acquisition, it was possible to harmonize the number of projections over respiratory phases. Depending on the breathing period (3.5 or 5 s) and the gantry rotation time (4 or 5 min), between ∼90 and 145 projections were acquired per respiratory phase resulting in a dose of ∼1.7-2.6 mGy per respiratory phase. Further dose savings and decreases in the scanning time are possible by acquiring only a subset of all respiratory phases, for example, peak-exhale and peak-inhale only scans.

CONCLUSIONS

This study is the first experimental demonstration of a new 4d CBCT acquisition paradigm in which imaging dose is efficiently utilized by actively triggering only those projections that are desired for the reconstruction process.