Simulation of spatiotemporal CT data sets using a 4D MRI-based lung motion model

Continuous time-resolved estimated synthetic 4D-CTs for dose reconstruction of lung tumor treatments at a 0.35 T MR-linac

Objective.To experimentally validate a method to create continuous time-resolved estimated synthetic 4D-computed tomography datasets (tresCTs) based on orthogonal cine MRI data for lung cancer treatments at a magnetic resonance imaging (MRI) guided linear accelerator (MR-linac).Approach.A breathing porcine lung phantom was scanned at a CT scanner and 0.35 T MR-linac. Orthogonal cine MRI series (sagittal/coronal orientation) at 7.3 Hz, intersecting tumor-mimicking gelatin nodules, were deformably registered to mid-exhale 3D-CT and 3D-MRI datasets. The time-resolved deformation vector fields were extrapolated to 3D and applied to a reference synthetic 3D-CT image (sCTref), while accounting for breathing phase-dependent lung density variations, to create 82 s long tresCTs at 3.65 Hz. Ten tresCTs were created for ten tracked nodules with different motion patterns in two lungs. For each dataset, a treatment plan was created on the mid-exhale phase of a measured ground truth (GT) respiratory-correlated 4D-CT dataset with the tracked nodule as gross tumor volume (GTV). Each plan was recalculated on the GT 4D-CT, randomly sampled tresCT, and static sCTrefimages. Dose distributions for corresponding breathing phases were compared in gamma (2%/2 mm) and dose-volume histogram (DVH) parameter analyses.Main results.The mean gamma pass rate between all tresCT and GT 4D-CT dose distributions was 98.6%. The mean absolute relative deviations of the tresCT with respect to GT DVH parameters were 1.9%, 1.0%, and 1.4% for the GTVD98%,D50%, andD2%, respectively, 1.0% for the remaining nodulesD50%, and 1.5% for the lungV20Gy. The gamma pass rate for the tresCTs was significantly larger (p< 0.01), and the GTVD50%deviations with respect to the GT were significantly smaller (p< 0.01) than for the sCTref.Significance.The results suggest that tresCTs could be valuable for time-resolved reconstruction and intrafractional accumulation of the dose to the GTV for lung cancer patients treated at MR-linacs in the future.

Evaluating differences in respiratory motion estimates during radiotherapy: a single planning 4DMRI versus daily 4DMRI

Background

In radiotherapy, respiratory-induced tumor motion is typically measured using a single four-dimensional computed tomography acquisition (4DCT). Irregular breathing leads to inaccurate motion estimates, potentially resulting in undertreatment of the tumor and unnecessary dose to healthy tissue. The aim of the research was to determine if a daily pre-treatment 4DMRI-strategy led to a significantly improved motion estimate compared to single planning 4DMRI (with or without outlier rejection).

Methods

4DMRI data sets from 10 healthy volunteers were acquired. The first acquisition simulated a planning MRI, the respiratory motion estimate (constructed from the respiratory signal, i.e. the 1D navigator) was compared to the respiratory signal in the subsequent scans (simulating 5–29 treatment fractions). The same procedure was performed using the first acquisition of each day as an estimate for the subsequent acquisitions that day (2 per day, 4–20 per volunteer), simulating a daily MRI strategy. This was done for three outlier strategies: no outlier rejection (NoOR); excluding 5% of the respiratory signal whilst minimizing the range (Min95) and excluding the datapoints outside the mean end-inhalation and end-exhalation positions (MeanIE).

Results

The planning MRI median motion estimates were 27 mm for NoOR, 18 mm for Min95, and 13 mm for MeanIE. The daily MRI median motion estimates were 29 mm for NoOR, 19 mm for Min95 and 15 mm for MeanIE. The percentage of time outside the motion estimate were for the planning MRI: 2%, 10% and 32% for NoOR, Min95 and MeanIE respectively. These values were reduced with the daily MRI strategy: 0%, 6% and 17%. Applying Min95 accounted for a 30% decrease in motion estimate compared to NoOR.

Conclusion

A daily MRI improved the estimation of respiratory motion as compared to a single 4D (planning) MRI significantly. Combining the Min95 technique with a daily 4DMRI resulted in a decrease of inclusion time of 6% with a 30% decrease of motion. Outlier rejection alone on a planning MRI often led to underestimation of the movement and could potentially lead to an underdosage.

Trial registration: protocol W15_373#16.007

Supplementary Information

The online version contains supplementary material available at 10.1186/s13014-021-01915-1.

Vector-Field dynamic X-ray (VF-DXR) using Optical Flow Method

OBJECTIVES

To explore the feasibility of Vector-Field DXR (VF-DXR) using optical flow method (OFM).

METHODS

Five healthy volunteers and five COPD patients were studied. DXR was performed in the standing position using a prototype X-ray system (Konica Minolta Inc., Tokyo, Japan). During the examination, participants took several tidal breaths and one forced breath. DXR image file was converted to the videos with different frames per second (fps): 15 fps, 7.5 fps, five fps, three fps, and 1.5 fps. Pixel-value gradient was calculated by the serial change of pixel value, which was subsequently converted mathematically to motion vector using OFM. Color-coding map and vector projection into horizontal and vertical components were also tested.

RESULTS

Dynamic motion of lung and thorax was clearly visualized using VF-DXR with an optimal frame rate of 5 fps. Color-coding map and vector projection into horizontal and vertical components were also presented. VF-DXR technique was also applied in COPD patients.

CONCLUSION

The feasibility of VF-DXR was demonstrated with small number of healthy subjects and COPD patients.

ADVANCES IN KNOWLEDGE

A new Vector-Field Dynamic X-ray (VF-DXR) technique is feasible for dynamic visualization of lung, diaphragms, thoracic cage, and cardiac contour.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

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.

T2-Weighted 4D Magnetic Resonance Imaging for Application in Magnetic Resonance-Guided Radiotherapy Treatment Planning

OBJECTIVES

The aim of this study was to develop and verify a method to obtain good temporal resolution T2-weighted 4-dimensional (4D-T2w) magnetic resonance imaging (MRI) by using motion information from T1-weighted 4D (4D-T1w) MRI, to support treatment planning in MR-guided radiotherapy.

MATERIALS AND METHODS

Ten patients with primary non-small cell lung cancer were scanned at 1.5 T axially with a volumetric T2-weighted turbo spin echo sequence gated to exhalation and a volumetric T1-weighted stack-of-stars spoiled gradient echo sequence with golden angle spacing acquired in free breathing. From the latter, 20 respiratory phases were reconstructed using the recently developed 4D joint MoCo-HDTV algorithm based on the self-gating signal obtained from the k-space center. Motion vector fields describing the respiratory cycle were obtained by deformable image registration between the respiratory phases and projected onto the T2-weighted image volume. The resulting 4D-T2w volumes were verified against the 4D-T1w volumes: an edge-detection method was used to measure the diaphragm positions; the locations of anatomical landmarks delineated by a radiation oncologist were compared and normalized mutual information was calculated to evaluate volumetric image similarity.

RESULTS

High-resolution 4D-T2w MRI was obtained. Respiratory motion was preserved on calculated 4D-T2w MRI, with median diaphragm positions being consistent with less than 6.6 mm (2 voxels) for all patients and less than 3.3 mm (1 voxel) for 9 of 10 patients. Geometrical positions were coherent between 4D-T1w and 4D-T2w MRI as Euclidean distances between all corresponding anatomical landmarks agreed to within 7.6 mm (Euclidean distance of 2 voxels) and were below 3.8 mm (Euclidean distance of 1 voxel) for 355 of 470 pairs of anatomical landmarks. Volumetric image similarity was commensurate between 4D-T1w and 4D-T2w MRI, as mean percentage differences in normalized mutual information (calculated over all respiratory phases and patients), between corresponding respiratory phases of 4D-T1w and 4D-T2w MRI and the tie-phase of 4D-T1w and 3-dimensional T2w MRI, were consistent to 0.41% ± 0.37%. Four-dimensional T2w MRI displayed tumor extent, structure, and position more clearly than corresponding 4D-T1w MRI, especially when mobile tumor sites were adjacent to organs at risk.

CONCLUSIONS

A methodology to obtain 4D-T2w MRI that retrospectively applies the motion information from 4D-T1w MRI to 3-dimensional T2w MRI was developed and verified. Four-dimensional T2w MRI can assist clinicians in delineating mobile lesions that are difficult to define on 4D-T1w MRI, because of poor tumor-tissue contrast.

A hybrid approach for fusing 4D-MRI temporal information with 3D-CT for the study of lung and lung tumor motion

PURPOSE

Accurate visualization of lung motion is important in many clinical applications, such as radiotherapy of lung cancer. Advancement in imaging modalities [e.g., computed tomography (CT) and MRI] has allowed dynamic imaging of lung and lung tumor motion. However, each imaging modality has its advantages and disadvantages. The study presented in this paper aims at generating synthetic 4D-CT dataset for lung cancer patients by combining both continuous three-dimensional (3D) motion captured by 4D-MRI and the high spatial resolution captured by CT using the authors' proposed approach.

METHODS

A novel hybrid approach based on deformable image registration (DIR) and finite element method simulation was developed to fuse a static 3D-CT volume (acquired under breath-hold) and the 3D motion information extracted from 4D-MRI dataset, creating a synthetic 4D-CT dataset.

RESULTS

The study focuses on imaging of lung and lung tumor. Comparing the synthetic 4D-CT dataset with the acquired 4D-CT dataset of six lung cancer patients based on 420 landmarks, accurate results (average error <2 mm) were achieved using the authors' proposed approach. Their hybrid approach achieved a 40% error reduction (based on landmarks assessment) over using only DIR techniques.

CONCLUSIONS

The synthetic 4D-CT dataset generated has high spatial resolution, has excellent lung details, and is able to show movement of lung and lung tumor over multiple breathing cycles.

Radiotherapy planning using MRI

The use of magnetic resonance imaging (MRI) in radiotherapy (RT) planning is rapidly expanding. We review the wide range of image contrast mechanisms available to MRI and the way they are exploited for RT planning. However a number of challenges are also considered: the requirements that MR images are acquired in the RT treatment position, that they are geometrically accurate, that effects of patient motion during the scan are minimized, that tissue markers are clearly demonstrated, that an estimate of electron density can be obtained. These issues are discussed in detail, prior to the consideration of a number of specific clinical applications. This is followed by a brief discussion on the development of real-time MRI-guided RT.