Patient-specific respiratory models using dynamic 3D MRI: Preliminary volunteer results

AAPM Task Group Report 303 endorsed by the ABS: MRI Implementation in HDR Brachytherapy-Considerations from Simulation to Treatment

The Task Group (TG) on Magnetic Resonance Imaging (MRI) Implementation in High Dose Rate (HDR) Brachytherapy - Considerations from Simulation to Treatment, TG 303, was constituted by the American Association of Physicists in Medicine's (AAPM's) Science Council under the direction of the Therapy Physics Committee, the Brachytherapy Subcommittee, and the Working Group on Brachytherapy Clinical Applications. The TG was charged with developing recommendations for commissioning, clinical implementation, and on-going quality assurance (QA). Additionally, the TG was charged with describing HDR brachytherapy (BT) workflows and evaluating practical consideration that arise when implementing MR imaging. For brevity, the report is focused on the treatment of gynecologic and prostate cancer. The TG report provides an introduction and rationale for MRI implementation in BT, a review of previous publications on topics including available applicators, clinical trials, previously published BT related TG reports, and new image guided recommendations beyond CT based practices. The report describes MRI protocols and methodologies, including recommendations for the clinical implementation and logical considerations for MR imaging for HDR BT. Given the evolution from prescriptive to risk-based QA,¹ an example of a risk-based analysis using MRI-based, prostate HDR BT is presented. In summary, the TG report is intended to provide clear and comprehensive guidelines and recommendations for commissioning, clinical implementation, and QA for MRI-based HDR BT that may be utilized by the medical physics community to streamline this process. This report is endorsed by the American Brachytherapy Society (ABS). This article is protected by copyright. All rights reserved.

A super-resolution framework for the reconstruction of T2-weighted (T2w) time-resolved (TR) 4DMRI using T1w TR-4DMRI as the guidance

PURPOSE

The purpose of this study was to develop T2-weighted (T2w) time-resolved (TR) four-dimensional magnetic resonance imaging (4DMRI) reconstruction technique with higher soft-tissue contrast for multiple breathing cycle motion assessment by building a super-resolution (SR) framework using the T1w TR-4DMRI reconstruction as guidance.

METHODS

The multi-breath T1w TR-4DMRI was reconstructed by deforming a high-resolution (HR: 2 × 2 × 2 mm³ ) volumetric breath-hold (BH, 20s) three-dimensional magnetic resonance imaging (3DMRI) image to a series of low-resolution (LR: 5 × 5 × 5 mm³ ) 3D cine images at a 2Hz frame rate in free-breathing (FB, 40 s) using an enhanced Demons algorithm, namely [T1^BH →FB] reconstruction. Within the same imaging session, respiratory-correlated (RC) T2w 4DMRI (2 × 2 × 2 mm³ ) was acquired based on an internal navigator to gain HR T2w (T2^HR ) in three states (full exhalation and mid and full inhalation) in ~5 min. Minor binning artifacts in the RC-4DMRI were automatically identified based on voxel intensity correlation (VIC) between consecutive slices as outliers (VIC < VIC^mean -σ) and corrected by deforming the artifact slices to interpolated slices from the adjacent slices iteratively until no outliers were identified. A T2^HR image with minimal deformation (<1 cm at the diaphragm) from the T1^BH image was selected for multi-modal B-Spline deformable image registration (DIR) to establish the T2^HR -T1^BH voxel correspondence. Two approaches to reconstruct T2w TR-4DMRI were investigated: (A) T2^HR →[T1^BH →FB]: to deform T2w HR to T1w BH only as T1w TR-4DMRI was reconstructed, and combine the two displacement vector fields (DVFs) to reconstruct T2w TR-4DMRI, and (B) [T2^HR ←T1^BH ]→FB: to deform T1w BH to T2w HR first and apply the deformed T1w BH to reconstruct T2w TR-4DMRI. The reconstruction times were similar, 8-12 min per volume. To validate the two methods, T2w- and T1w-mapped 4D XCAT digital phantoms were utilized with three synthetic spherical tumors (ϕ = 2.0, 3.0, and 4.0 cm) in the lower or mid lobes as the ground truth to evaluate the tumor location (the center of mass, COM), size (volume ratio, %V), and shape (Dice index). Six lung cancer patients were scanned under an IRB-approved protocol and the T2w TR-4DMRI images reconstructed from the two methods were compared based on the preservation of the three tumor characteristics. The local tumor-contained image quality was also characterized using the VIC and structure similarity (SSIM) indexes.

RESULTS

In the 4D digital phantom, excellent tumor alignment after T2^HR -T1^HR DIR is achieved: ∆COM = 0.8 ± 0.5 mm, %V = 1.06 ± 0.02, and Dice = 0.91 ± 0.03, in both deformation directions using the DIR-target image as the reference. In patients, binning artifacts are corrected with improved image quality: average VIC increases from 0.92 ± 0.03 to 0.95 ± 0.01. Both T2w TR-4DMRI reconstruction methods produce similar tumor alignment errors ∆COM = 2.9 ± 0.6 mm. However, method B ([T2^HR ←T1^BH ]→FB) produces superior results in preserving more T2w tumor features with a higher %V = 0.99 ± 0.03, Dice = 0.81 ± 0.06, VIC = 0.85 ± 0.06, and SSIM = 0.65 ± 0.10 in the T2w TR-4DMRI images.

CONCLUSIONS

This study has demonstrated the feasibility of T2w TR-4DMRI reconstruction with high soft-tissue contrast and adequately-preserved tumor position, size, and shape in multiple breathing cycles. The T2w-centric DIR (method B) produces a superior solution for the SR-based framework of T2w TR-4DMRI reconstruction with highly preserved tumor characteristics and local image features, which are useful for tumor delineation and motion management in radiation therapy.

Novel Super-Resolution Approach to Time-Resolved Volumetric 4-Dimensional Magnetic Resonance Imaging With High Spatiotemporal Resolution for Multi-Breathing Cycle Motion Assessment

PURPOSE

To develop and evaluate a super-resolution approach to reconstruct time-resolved 4-dimensional magnetic resonance imaging (TR-4DMRI) with a high spatiotemporal resolution for multi-breathing cycle motion assessment.

METHODS AND MATERIALS

A super-resolution approach was developed to combine fast 3-dimensional (3D) cine MRI with low resolution during free breathing (FB) and high-resolution 3D static MRI during breath hold (BH) using deformable image registration. A T1-weighted, turbo field echo sequence, coronal 3D cine acquisition, partial Fourier approximation, and SENSitivity Encoding parallel acceleration were used. The same MRI pulse sequence, field of view, and acceleration techniques were applied in both FB and BH acquisitions; the intensity-based Demons deformable image registration method was used. Under an institutional review board-approved protocol, 7 volunteers were studied with 3D cine FB scan (voxel size: 5 × 5 × 5 mm³) at 2 Hz for 40 seconds and a 3D static BH scan (2 × 2 × 2 mm³). To examine the image fidelity of 3D cine and super-resolution TR-4DMRI, a mobile gel phantom with multi-internal targets was scanned at 3 speeds and compared with the 3D static image. Image similarity among 3D cine, 4DMRI, and 3D static was evaluated visually using difference image and quantitatively using voxel intensity correlation and Dice index (phantom only). Multi-breathing-cycle waveforms were extracted and compared in both phantom and volunteer images using the 3D cine as the references.

RESULTS

Mild imaging artifacts were found in the 3D cine and TR-4DMRI of the mobile gel phantom with a Dice index of >0.95. Among 7 volunteers, the super-resolution TR-4DMRI yielded high voxel-intensity correlation (0.92 ± 0.05) and low voxel-intensity difference (<0.05). The detected motion differences between TR-4DMRI and 3D cine were -0.2 ± 0.5 mm (phantom) and -0.2 ± 1.9 mm (diaphragms).

CONCLUSION

Super-resolution TR-4DMRI has been reconstructed with adequate temporal (2 Hz) and spatial (2 × 2 × 2 mm³) resolutions. Further TR-4DMRI characterization and improvement are necessary before clinical applications. Multi-breathing cycles can be examined, providing patient-specific breathing irregularities and motion statistics for future 4D radiation therapy.

Review of potential improvements using MRI in the radiotherapy workflow

The goal of modern radiotherapy is to deliver a lethal amount of dose to tissue volumes that contain a significant amount of tumour cells while sparing surrounding unaffected or healthy tissue. Online image guided radiotherapy with stereotactic ultrasound, fiducial-based planar X-ray imaging or helical/conebeam CT has dramatically improved the precision of radiotherapy, with moving targets still posing some methodical problems regarding positioning. Therefore, requirements for precise target delineation and identification of functional body structures to be spared by high doses become more evident. The identification of areas of relatively radioresistant cells or areas of high tumor cell density is currently under development. This review outlines the state of the art of MRI integration into treatment planning and its importance in follow up and the quantification of biological effects. Finally the current state of the art of online imaging for patient positioning will be outlined and indications will be given what the potential of integrated radiotherapy/online MRI systems is.

Joint surface reconstruction and 4D deformation estimation from sparse data and prior knowledge for marker-less Respiratory motion tracking

PURPOSE

The intraprocedural tracking of respiratory motion has the potential to substantially improve image-guided diagnosis and interventions. The authors have developed a sparse-to-dense registration approach that is capable of recovering the patient's external 3D body surface and estimating a 4D (3D + time) surface motion field from sparse sampling data and patient-specific prior shape knowledge.

METHODS

The system utilizes an emerging marker-less and laser-based active triangulation (AT) sensor that delivers sparse but highly accurate 3D measurements in real-time. These sparse position measurements are registered with a dense reference surface extracted from planning data. Thereby a dense displacement field is recovered, which describes the spatio-temporal 4D deformation of the complete patient body surface, depending on the type and state of respiration. It yields both a reconstruction of the instantaneous patient shape and a high-dimensional respiratory surrogate for respiratory motion tracking. The method is validated on a 4D CT respiration phantom and evaluated on both real data from an AT prototype and synthetic data sampled from dense surface scans acquired with a structured-light scanner.

RESULTS

In the experiments, the authors estimated surface motion fields with the proposed algorithm on 256 datasets from 16 subjects and in different respiration states, achieving a mean surface reconstruction accuracy of ± 0.23 mm with respect to ground truth data-down from a mean initial surface mismatch of 5.66 mm. The 95th percentile of the local residual mesh-to-mesh distance after registration did not exceed 1.17 mm for any subject. On average, the total runtime of our proof of concept CPU implementation is 2.3 s per frame, outperforming related work substantially.

CONCLUSIONS

In external beam radiation therapy, the approach holds potential for patient monitoring during treatment using the reconstructed surface, and for motion-compensated dose delivery using the estimated 4D surface motion field in combination with external-internal correlation models.

Adequacy of inhale/exhale breathhold CT based ITV margins and image-guided registration for free-breathing pancreas and liver SBRT

Purpose

To evaluate use of breath-hold CTs and implanted fiducials for definition of the internal target volume (ITV) margin for upper abdominal stereotactic body radiation therapy (SBRT). To study the statistics of inter- and intra-fractional motion information.

Methods and materials

11 patients treated with SBRT for locally advanced pancreatic cancer (LAPC) or liver cancer were included in the study. Patients underwent fiducial implantation, free-breathing CT and breath-hold CTs at end inhalation/exhalation. All patients were planned and treated with SBRT using volumetric modulated arc therapy (VMAT). Two margin strategies were studied: Strategy I uses PTV = ITV + 3 mm; Strategy II uses PTV = GTV + 1.5 cm. Both CBCT and kV orthogonal images were taken and analyzed for setup before patient treatments. Tumor motion statistics based on skeletal registration and on fiducial registration were analyzed by fitting to Gaussian functions.

Results

All 11 patients met SBRT planning dose constraints using strategy I. Average ITV margins for the 11 patients were 2 mm RL, 6 mm AP, and 6 mm SI. Skeletal registration resulted in high probability (RL = 69%, AP = 4.6%, SI = 39%) that part of the tumor will be outside the ITV. With the 3 mm ITV expansion (Strategy 1), the probability reduced to RL 32%, AP 0.3%, SI 20% for skeletal registration; and RL 1.2%, AP 0%, SI 7% for fiducial registration. All 7 pancreatic patients and 2 liver patients failed to meet SBRT dose constraints using strategy II. The liver dose was increased by 36% for the other 2 liver patients that met the SBRT dose constraints with strategy II.

Conclusions

Image guidance matching to skeletal anatomy is inadequate for SBRT positioning in the upper abdomen and usage of fiducials is highly recommended. Even with fiducial implantation and definition of an ITV, a minimal 3 mm planning margin around the ITV is needed to accommodate intra-fractional uncertainties.