Technical Note: Simulation of 4DCT tumor motion measurement errors

Stereotactic Arrhythmia Radioablation (STAR): Assessment of cardiac and respiratory heart motion in ventricular tachycardia patients - A STOPSTORM.eu consortium review

AIM

To identify the optimal STereotactic Arrhythmia Radioablation (STAR) strategy for individual patients, cardiorespiratory motion of the target volume in combination with different treatment methodologies needs to be evaluated. However, an authoritative overview of the amount of cardiorespiratory motion in ventricular tachycardia (VT) patients is missing.

METHODS

In this STOPSTORM consortium study, we performed a literature review to gain insight into cardiorespiratory motion of target volumes for STAR. Motion data and target volumes were extracted and summarized.

RESULTS

Out of the 232 studies screened, 56 provided data on cardiorespiratory motion, of which 8 provided motion amplitudes in VT patients (n = 94) and 10 described (cardiac/cardiorespiratory) internal target volumes (ITVs) obtained in VT patients (n = 59). Average cardiac motion of target volumes was < 5 mm in all directions, with maximum values of 8.0, 5.2 and 6.5 mm in Superior-Inferior (SI), Left-Right (LR), Anterior-Posterior (AP) direction, respectively. Cardiorespiratory motion of cardiac (sub)structures showed average motion between 5-8 mm in the SI direction, whereas, LR and AP motions were comparable to the cardiac motion of the target volumes. Cardiorespiratory ITVs were on average 120-284% of the gross target volume. Healthy subjects showed average cardiorespiratory motion of 10-17 mm in SI and 2.4-7 mm in the AP direction.

CONCLUSION

This review suggests that despite growing numbers of patients being treated, detailed data on cardiorespiratory motion for STAR is still limited. Moreover, data comparison between studies is difficult due to inconsistency in parameters reported. Cardiorespiratory motion is highly patient-specific even under motion-compensation techniques. Therefore, individual motion management strategies during imaging, planning, and treatment for STAR are highly recommended.

Dosimetric Study of Deep Learning-Guided ITV Prediction in Cone-beam CT for Lung Stereotactic Body Radiotherapy

Purpose

The purpose of this study was to evaluate the accuracy of a lung stereotactic body radiotherapy (SBRT) treatment plan with the target of a newly predicted internal target volume (ITV_predict) and the feasibility of its clinical application. ITV_predict was automatically generated by our in-house deep learning model according to the cone-beam CT (CBCT) image database.

Method

A retrospective study of 45 patients who underwent SBRT was involved, and Mask R-CNN based algorithm model helped to predict the internal target volume (ITV) using the CBCT image database. The geometric accuracy of ITV_predict was verified by the Dice Similarity Coefficient (DSC), 3D Motion Range (R_3D), Relative Volume Index (RVI), and Hausdorff Distance (HD). The PTV_predict was generated by ITV_predict, which was registered and then projected on free-breath CT (FBCT) images. The PTV_FBCT was margined from the GTV on FBCT images gross tumor volume on free-breath CT (GTV_FBCT). Treatment plans with the target of Predict planning target volume on CBCT images (PTV_predict) and planning target volume on free-breath CT (PTV_FBCT) were respectively re-established, and the dosimetric parameters included the ratio of the volume of patients receiving at least the prescribed dose to the volume of PTV (R_100%), the ratio of the volume of patients receiving at least 50% of the prescribed dose to the volume of PTV in the Radiation Therapy Oncology Group (RTOG) 0813 Trial (R_50%), Gradient Index (GI), and the maximum dose 2 cm from the PTV (D_2cm), which were evaluated via Plan_4DCT, plan which based on PTV_predict (Plan_predict), and plan which based on PTV_FBCT (Plan_FBCT).

Result

The geometric results showed that there existed a good correlation between ITV_predict and ITV on the 4-dimensional CT [ITV_4DCT; DSC= 0.83 ±0.18]. However, the average volume of ITV_predict was 10% less than that of ITV_4DCT (p = 0.333). No significant difference in dose coverage was found in V_100% for the ITV with 99.98 ± 0.04% in the ITV_4DCT vs. 97.56 ± 4.71% in the ITV_predict (p = 0.162). Dosimetry parameters of PTV, including R_100%, R_50%, GI and D_2cm showed no statistically significant difference between each plan (p > 0.05).

Conclusion

Dosimetric parameters of Plan_predict are clinically comparable to those of the original Plan_4DCT. This study confirmed that the treatment plan based on ITV_predict produced by our model could automatically meet clinical requirements. Thus, for patients undergoing lung SBRT, the model has great potential for using CBCT images for ITV contouring which can be used in treatment planning.

4DCT and VMAT for lung patients with irregular breathing

PURPOSE

Irregular breathing in lung cancer patients is a common contra-indication to 4D computerized tomography (4DCT), which may then limit radiotherapy treatment options. For irregular breathers, we investigated whether 3DCT or 4DCT (1) better represents tumor motion, (2) better represents average tumor densities, and (3) better allows for volumetric modulated arc threarpy (VMAT) plans delivered with acceptable dosimetric accuracy.

METHODS

Ten clinical breathing traces were identified with irregularities in phase and amplitude, and fed to a programmable moving platform incorporating an anthropomorphic lung tumor phantom. 3DCT and 4DCT data resorted by phase (4DCT-P) and amplitude (4DCT-A) were acquired for each trace. Tumors were delineated by Hounsfield unit (HU) thresholding and apparent motion range assessed. HU profiles were extracted from each image and agreement with calculated expected profiles quantified using area-under-curve (AUC) scoring. Clinically representative VMAT plans were created for each image, delivered to the irregularly moving phantom, and measured with a small-volume ion chamber at the tumor center.

RESULTS

Median difference from expected tumor motion range for 3DCT, 4DCT-P, and 4DCT-A was 2.5 [1.6-3.6] cm, 1.1 [0.1-1.9] cm, and 1.3 [0.4-1.9] cm, respectively (p = 0.005, 4DCT-P vs. 3DCT). Median AUC scores (ideal = 0) for 3DCT, 4DCT-P, and 4DCT-A were 0.25 [0.14-0.49], 0.12 [0.05-0.42], and 0.13 [0.09-0.44], respectively (p = 0.005, 4DCT-P vs. 3DCT). Nine of ten 4DCT-P plans and all 4DCT-A plans measured within 2.5% of expected dose in the treatment planning system (TPS), compared with seven 3DCT plans.

CONCLUSION

For the cases studied tumor motion range and average density was better represented with 4DCT compared with 3DCT, even in the presence of irregular breathing. 4DCT images allowed for delivery of VMAT plans with acceptable dosimetric accuracy. No significant differences were detected between phase and amplitude resorting. In combination with 4D cone beam imaging at treatment, our findings have given us confidence to introduce 4DCT and VMAT for lung radiotherapy patients with irregular breathing.

The impact of respiratory motion and CT pitch on the robustness of radiomics feature extraction in 4DCT lung imaging

PURPOSE/OBJECTIVE(S)

The precise radiomics analysis on thoracic 4DCT data is easily compromised by the respiratory motion and CT scan parameter setting, thus leading to the risk of overfitting and/or misinterpretation of data in AI-enabled therapeutic model building. In this study, we investigated the impact of respiratory amplitudes, frequencies and CT scan pitch settings within the thoracic 4DCT scan on robust radiomics feature selection.

MATERIALS/METHODS

A Three-dimensional QUSAR^TM lung tumor phantom was used to simulate different respiratory amplitudes and frequencies along with different CT scan pitch settings. A total of 43 tumor respiratory patterns extracted from 43 patients with non-small cell lung cancer were used to drive the QUSAR^TM lung tumor phantom to mimic the human tumor motion. The 4DCT images of the QUSAR^TM lung tumor phantom with different respiratory patterns and different CT scan pitch setups were acquired for radiomics feature extraction. A static high-quality CT images of the phantom acquired were also used as a reference for radiomics feature extraction. The range of respiratory amplitudes was mimicked at 3mm at left and right (LR) and anterior and posterior (AP) directions and 3mm - 15 mm at the superior and inferior (SI) direction with an interval of 2 mm. The respiratory frequencies were set at 10, 11, 12, 13, 14, 15 and 20 beats per minute (BPMs), respectively. The CT scan pitches were set at 0.025, 0.048, 0.071, 0.93, 0.108, 0.14, 0.16, 0.18, 0.21, 0.23, and 0.25, respectively, which was based on a procedure described in Med. Phys. 30(1):88-97. The pairwise Concordance Correlation Coefficient (CCC) was used to determine the robustness of radiomics feature extraction via comparing the agreement in feature values between 1766 radiomics features extracted from each image acquired under different combinations of respiratory amplitudes and frequencies and CT scan pitches of 4DCT and those extracted from the static CT images.

RESULTS

(1) When the respiratory amplitudes were at 3, 5, 7, 9, 12 and 15mm in the SI direction, the maximum CCC index could be achieved at the reconstructed 4DCT phase images of 60%, 70%, 30%, 20%, 60%~70% and 10%, respectively. Under these six amplitudes, the maximum intensity projection (MIP) and average intensity projection (AIP) images reconstructed show mean CCC values of 0.778 and 0.609, respectively, in pairwise radiomics feature extraction comparison between 4DCT and static CT. (2) When the respiratory amplitude was set at 12 mm in the SI direction, the maximum CCC index could be consistently achieved at the reconstructed 4DCT phase of 90% for the seven respiratory frequencies of 10, 11, 12, 13, 14, 15 and 20 BPMs, respectively. Under these respiratory states, the MIP and AIP images reconstructed show mean CCC values of 0.702 and 0.562, respectively. (3) When the respiratory amplitude was set at 12 mm and the respiratory frequency was set at 13 BPM, the maximum CCC index could be obtained at the reconstructed 4DCT phase of 90% for all scan pitches used except the 0% phase which was obtained at the pitch setting of 0.048. Under these CT scan pitch settings, the MIP and AIP images reconstructed show mean CCC values of 0.558 and 0.782, respectively. (4) The total number of robust features were 50, 34 and 35 with different respiratory amplitudes and phases and CT scanning pitch used (CCC values ≥ 0.99).

CONCLUSION

In 4DCT, the respiratory amplitude, frequency and CT scan pitch are three limiting factors that greatly affect the robustness of radiomics feature extraction. The reconstructed 4DCT phases with better robustness along with suitable respiratory amplitude, frequency and CT scan pitch determined could be used to guide the breathing training for patients with lung cancer for radiation therapy to improve the robust radiomics feature extraction process.

Quality assurance of four-dimensional computed tomography in a multicentre trial of stereotactic body radiotherapy of centrally located lung tumours

Background and Purpose

Extensive radiation therapy quality assurance (RTQA) programs are needed when advanced radiotherapy treatments are used. As part of the RTQA four dimensional computed tomography (4DCT) imaging performance needs to be assessed. Here we present the RTQA data related to 4DCT procedures used within the context of stereotactic body radiotherapy (SBRT) of centrally located lung tumours. It provides an overview of the 4DCT acquisition methods and achievable accuracy of imaging lung tumour volumes.

Materials and Methods

3DCT and 4DCT images were acquired from a CIRS phantom with spheres of 7.5 and 12.5 mm radius using the institutional scan protocols. Regular asymmetric tumour motion was simulated with varying amplitudes and periods. Target volumes were reconstructed using auto-contouring with scanner specific thresholds. Volume and amplitudes deviations were assessed.

Results

Although acquisition parameters were rather homogeneous over the eleven institutions analysed, volume deviations were observed. Average volume deviations for the 12.5 mm sphere were 15% (−4% to 69%) at end of inspiration, 2% (−2% to 9.0%) at end of expiration and 12% (0% to 36%) at mid-ventilation. For the 7.5 mm sphere deviations were 13% (−99% to 65%), 16% (−34% to 66%) and 1% (−13% to 20%), respectively. The amplitude deviation was generally within 2 mm although underestimations up to 6 mm were observed.

Conclusions

The expiration phase was the most accurate phase to define the tumour volume and should be preferred for GTV delineation of tumours exhibiting large motion causing motion artefacts when using mid-ventilation or tracking techniques. The large variation found among the institutions indicated that further improvements in 4DCT imaging were possible. Recommendations for 4DCT QA have been formulated.

Preoperative Imaging Overestimates the Tumor Size in Pancreatic Neuroendocrine Neoplasms Associated with Multiple Endocrine Neoplasia Type 1

BACKGROUND

Radiological tumor size of non-functioning pancreatic neuroendocrine neoplasms (Nf-pNENs) associated with multiple endocrine neoplasia type 1 (MEN1) is a crucial parameter to indicate surgery. The aim of this study was to compare radiological size (RS) and pathologic size (PS) of MEN1 associated with pNENs.

METHODS

Prospectively collected data of MEN1 patients who underwent pancreatic resections for pNENs were retrospectively analyzed. RS was defined as the largest tumor diameter measured on endoscopic ultrasound (EUS), magnetic resonance imaging (MRI) or computed tomography (CT). PS was defined as the largest tumor diameter on pathological analysis. Student's t test and linear regression analysis were used to compare the median RS and PS. p < 0.05 was considered significant.

RESULTS

Forty-four patients with a median age of 37 (range 10-68) years underwent primary pancreatic resections for pNENs. Overall, the median RS (20 mm, range 3-100 mm) was significantly larger than the PS (13 mm, range 4-110 mm) (p = 0.001). In patients with pNENs < 20 mm (n = 27), the size difference (median RS 15 mm vs PS 12 mm) was also significant (p = 0.003). However, the only modality that significantly overestimated the PS was EUS (median RS 14 mm vs 11 mm; p = 0.0002). RS overestimated the PS in 21 patients (21 of 27 patients, 78%). Five of 11 patients (12%) with a Nf-pNEN and a RS > 20 mm had in reality a PS < 20 mm. MRI was the imaging technique that best correlated with PS in the total cohort (r = 0.8; p < 0.0001), whereas EUS was the best correlating imaging tool in pNENs < 20 mm (r = 0.5; p = 0.0001).

CONCLUSION

Preoperative imaging, especially EUS, frequently overestimates the size of MEN1-pNENs, especially those with a PS < 20 mm. This should be considered when indicating surgery in MEN1 patients with small Nf-pNENs.

Feasibility of proton pencil beam scanning treatment of free-breathing lung cancer patients

BACKGROUND

The interplay effect might degrade the dose of pencil beam scanning proton therapy to a degree that free-breathing treatment might be impossible without further motion mitigation techniques, which complicate and prolong the treatment. We assessed whether treatment of free-breathing patients without motion mitigation is feasible.

MATERIAL AND METHODS

For 40 lung cancer patients, 4DCT datasets and individual breathing patterns were used to simulate 4D dynamic dose distributions of 3D treatment plans over 33 fractions delivered with an IBA universal nozzle. Evaluation was done by assessing under- and overdosage in the target structure using the parameters V90, V95, V98, D98, D2, V107 and V110. The impact of using beam-specific target volumes and the impact of changes in motion and patient anatomy in control 4DCTs were assessed.

RESULTS

Almost half of the patients had tumour motion amplitudes of less than 5 mm. Under- and overdosage was significantly smaller for patients with tumour motion below 5 mm compared to patients with larger motion (2% vs. 13% average absolute reduction of V95, 2% vs. 8% average increase in V107, p < .01). Simulating a 33-fraction treatment, the dose degradation was reduced but persisted for patients with tumour motion above 5 mm (average ΔV95 of <1% vs. 3%, p < .01). Beam-specific target volumes reduced the dose degradation in a fractionated treatment, but were more relevant for large motion. Repeated 4DCT revealed that changes in tumour motion during treatment might result in unexpected large dose degradations.

CONCLUSION

Tumour motion amplitude is an indicator of dose degradation caused by the interplay effect. Fractionation reduces the dose degradation allowing the unmitigated treatment of patients with small tumour motions of less than 5 mm. The beam-specific target approach improves the dose coverage. The tumour motion and position needs to be assessed during treatment for all patients, to quickly react to possible changes, which might require treatment adaptation.

Initial clinical observations of intra- and interfractional motion variation in MR-guided lung SBRT

OBJECTIVE

To evaluate variations in intra- and interfractional tumour motion, and the effect on internal target volume (ITV) contour accuracy, using deformable image registration of real-time two-dimensional-sagittal cine-mode MRI acquired during lung stereotactic body radiation therapy (SBRT) treatments.

METHODS

Five lung tumour patients underwent free-breathing SBRT treatments on the ViewRay system, with dose prescribed to a planning target volume (defined as a 3-6 mm expansion of the 4DCT-ITV). Sagittal slice cine-MR images (3.5 × 3.5 mm² pixels) were acquired through the centre of the tumour at 4 frames per second throughout the treatments (3-4 fractions of 21-32 min). Tumour gross tumour volumes (GTVs) were contoured on the first frame of the MR cine and tracked for the first 20 min of each treatment using offline optical-flow based deformable registration implemented on a GPU cluster. A ground truth ITV (MR-ITV_20 min) was formed by taking the union of tracked GTV contours. Pseudo-ITVs were generated from unions of the GTV contours tracked over 10 s segments of image data (MR-ITV_10 s).

RESULTS

Differences were observed in the magnitude of median tumour displacement between days of treatments. MR-ITV_10 s areas were as small as 46% of the MR-ITV_20 min.

CONCLUSION

An ITV offers a "snapshot" of breathing motion for the brief period of time the tumour is imaged on a specific day. Real-time MRI over prolonged periods of time and over multiple treatment fractions shows that ITV size varies. Further work is required to investigate the dosimetric effect of these results. Advances in knowledge: Five lung tumour patients underwent free-breathing MRI-guided SBRT treatments, and their tumours tracked using deformable registration of cine-mode MRI. The results indicate that variability of both intra- and interfractional breathing amplitude should be taken into account during planning of lung radiotherapy.

Estimation of lung tumor position from multiple anatomical features on 4D-CT using multiple regression analysis

To estimate the lung tumor position from multiple anatomical features on four‐dimensional computed tomography (4D‐CT) data sets using single regression analysis (SRA) and multiple regression analysis (MRA) approach and evaluate an impact of the approach on internal target volume (ITV) for stereotactic body radiotherapy (SBRT) of the lung. Eleven consecutive lung cancer patients (12 cases) underwent 4D‐CT scanning. The three‐dimensional (3D) lung tumor motion exceeded 5 mm. The 3D tumor position and anatomical features, including lung volume, diaphragm, abdominal wall, and chest wall positions, were measured on 4D‐CT images. The tumor position was estimated by SRA using each anatomical feature and MRA using all anatomical features. The difference between the actual and estimated tumor positions was defined as the root‐mean‐square error (RMSE). A standard partial regression coefficient for the MRA was evaluated. The 3D lung tumor position showed a high correlation with the lung volume (R = 0.92 ± 0.10). Additionally, ITVs derived from SRA and MRA approaches were compared with ITV derived from contouring gross tumor volumes on all 10 phases of the 4D‐CT (conventional ITV). The RMSE of the SRA was within 3.7 mm in all directions. Also, the RMSE of the MRA was within 1.6 mm in all directions. The standard partial regression coefficient for the lung volume was the largest and had the most influence on the estimated tumor position. Compared with conventional ITV, average percentage decrease of ITV were 31.9% and 38.3% using SRA and MRA approaches, respectively. The estimation accuracy of lung tumor position was improved by the MRA approach, which provided smaller ITV than conventional ITV.

A method for improved 4D-computed tomography data acquisition

In four-dimensional time-dependent computed tomography (4D-CT) of the lungs, irregularities in breathing movements can cause errors in data acquisition, or even data loss. We present a method based on sending a synthetic, regular breathing signal to the CT instead of the real signal, which ensures 4D-CT data sets without data loss. Subsequent correction of the signal based on the real breathing curve enables an accurate reconstruction of the size and movement of the target volume. This makes it possible to plan radiation treatment based on the obtained data. The method was tested with dynamic thorax phantom measurements using synthetic and real breathing patterns.