Rabih Hammoud

Characterization of 3D geometric distortion of magnetic resonance imaging scanners commissioned for radiation therapy planning

OBJECTIVE To develop a method for the assessment and characterization of 3D geometric distortion as part of routine quality assurance for MRI scanners commissioned for Radiation Therapy planning. MATERIALS AND METHODS In this study, the in-plane and through-plane geometric distortions on a 1.5T GE MRI-SIM unit are characterized and the 2D and 3D correction algorithms provided by the vendor are evaluated. We used a phantom developed by GE Healthcare that covers a large field of view of 500mm, and consists of layers of foam embedded with a matrix of ellipsoidal markers. An in-house Java-based software module was developed to automatically assess the geometric distortion by calculating the center of each marker using the center of mass method, correcting of gross rotation errors and comparing the corrected positions with a CT gold standard data set. Spatial accuracy of typical pulse sequences used in RT planning was assessed (2D T1/T2 FSE, 3D CUBE, T1 SPGR) using the software. The accuracy of vendor specific geometric distortion correction (GDC) algorithms was quantified by measuring distortions before and after the application of the 2D and 3D correction algorithms. RESULTS Our algorithm was able to accurately calculate geometric distortion with sub-pixel precision. For all typical MR sequences used in Radiotherapy, the vendors GDC was able to substantially reduce the distortions. Our results showed also that the impact of the acquisition produced a maximum variation of 0.2mm over a radial distance of 200mm. It has been shown that while the 2D correction algorithm remarkably reduces the in-plane geometric distortion, 3D geometric distortion further reduced the geometric distortion by correcting both in-plane and through-plane distortions in all acquisitions. CONCLUSION The presented methods represent a valuable tool for routine quality assurance of MR applications that require stringent spatial accuracy assessment such as radiotherapy. The phantom used in this study provides three dimensional arrays of control points. These tools and the detailed results can be also used for developing new geometric distortion correction algorithms or improving the existing ones.

Development and validation of a novel large field of view phantom and a software module for the quality assurance of geometric distortion in magnetic resonance imaging

OBJECTIVE To develop and validate a large field of view phantom and quality assurance software tool for the assessment and characterization of geometric distortion in MRI scanners commissioned for radiation therapy planning. MATERIALS AND METHODS A purpose built phantom was developed consisting of 357 rods (6mm in diameter) of polymethyl-methacrylat separated by 20mm intervals, providing a three dimensional array of control points at known spatial locations covering a large field of view up to a diameter of 420mm. An in-house software module was developed to allow automatic geometric distortion assessment. This software module was validated against a virtual dataset of the phantom that reproduced the exact geometry of the physical phantom, but with known translational and rotational displacements and warping. For validation experiments, clinical MRI sequences were acquired with and without the application of a commercial 3D distortion correction algorithm (Gradwarp™). The software module was used to characterize and assess system-related geometric distortion in the sequences relative to a benchmark CT dataset, and the efficacy of the vendor geometric distortion correction algorithms (GDC) was also assessed. RESULTS Results issued from the validation of the software against virtual images demonstrate the algorithms ability to accurately calculate geometric distortion with sub-pixel precision by the extraction of rods and quantization of displacements. Geometric distortion was assessed for the typical sequences used in radiotherapy applications and over a clinically relevant 420mm field of view (FOV). As expected and towards the edges of the field of view (FOV), distortion increased with increasing FOV. For all assessed sequences, the vendor GDC was able to reduce the mean distortion to below 1mm over a field of view of 5, 10, 15 and 20cm radius respectively. CONCLUSION Results issued from the application of the developed phantoms and algorithms demonstrate a high level of precision. The results indicate that this platform represents an important, robust and objective tool to perform routine quality assurance of MR-guided therapeutic applications, where spatial accuracy is paramount.

Generation of synthetic CT using multi-scale and dual-contrast patches for brain MRI-only external beam radiotherapy

PURPOSE To create a synthetic CT (sCT) from conventional brain MRI using a patch-based method for MRI-only radiotherapy planning and verification. METHODS Conventional T1 and T2-weighted MRI and CT datasets from 13 patients who underwent brain radiotherapy were included in a retrospective study whereas 6 patients were tested prospectively. A new contribution to the Non-local Means Patch-Based Method (NMPBM) framework was done with the use of novel multi-scale and dual-contrast patches. Furthermore, the training dataset was improved by pre-selecting the closest database patients to the target patient for computation time/accuracy balance. sCT and derived DRRs were assessed visually and quantitatively. VMAT planning was performed on CT and sCT for hypothetical PTVs in homogeneous and heterogeneous regions. Dosimetric analysis was done by comparing Dose Volume Histogram (DVH) parameters of PTVs and organs at risk (OARs). Positional accuracy of MRI-only image-guided radiation therapy based on CBCT or kV images was evaluated. RESULTS The retrospective (respectively prospective) evaluation of the proposed Multi-scale and Dual-contrast Patch-Based Method (MDPBM) gave a mean absolute error MAE=99.69±11.07HU (98.95±8.35HU), and a Dice in bones DIbone=83±0.03 (0.82±0.03). Good agreement with conventional planning techniques was obtained; the highest percentage of DVH metric deviations was 0.43% (0.53%) for PTVs and 0.59% (0.75%) for OARs. The accuracy of sCT/CBCT or DRRsCT/kV images registration parameters was <2mm and <2°. Improvements with MDPBM, compared to NMPBM, were significant. CONCLUSION We presented a novel method for sCT generation from T1 and T2-weighted MRI potentially suitable for MRI-only external beam radiotherapy in brain sites.

Voluntary deep inspiration breath-hold reduces the heart dose without compromising the target volume coverage during radiotherapy for left-sided breast cancer

Abstract Background During radiotherapy of left-sided breast cancer, parts of the heart are irradiated, which may lead to late toxicity. We report on the experience of single institution with cardiac-sparing radiotherapy using voluntary deep inspiration breath hold (V-DIBH) and compare its dosimetric outcome with free breathing (FB) technique. Patients and methods Left-sided breast cancer patients, treated at our department with postoperative radiotherapy of breast/chest wall +/– regional lymph nodes between May 2015 and January 2017, were considered for inclusion. FB-computed tomography (CT) was obtained and dose-planning performed. Cases with cardiac V25Gy ≥ 5% or risk factors for heart disease were coached for V-DIBH. Compliant patients were included. They underwent additional CT in V-DIBH for planning, followed by V-DIBH radiotherapy. Dose volume histogram parameters for heart, lung and optimized planning target volume (OPTV) were compared between FB and BH. Treatment setup shifts and systematic and random errors for V-DIBH technique were compared with FB historic control. Results Sixty-three patients were considered for V-DIBH. Nine (14.3%) were non-compliant at coaching, leaving 54 cases for analysis. When compared with FB, V-DIBH resulted in a significant reduction of mean cardiac dose from 6.1 +/– 2.5 to 3.2 +/– 1.4 Gy (p < 0.001), maximum cardiac dose from 51.1 +/– 1.4 to 48.5 +/– 6.8 Gy (p = 0.005) and cardiac V25Gy from 8.5 +/– 4.2 to 3.2 +/– 2.5% (p < 0.001). Heart volumes receiving low (10–20 Gy) and high (30–50 Gy) doses were also significantly reduced. Mean dose to the left anterior coronary artery was 23.0 (+/– 6.7) Gy and 14.8 (+/– 7.6) Gy on FB and V-DIBH, respectively (p < 0.001). Differences between FB- and V-DIBH-derived mean lung dose (11.3 +/– 3.2 vs. 10.6 +/– 2.6 Gy), lung V20Gy (20.5 +/– 7 vs. 19.5 +/– 5.1 Gy) and V95% for the OPTV (95.6 +/– 4.1 vs. 95.2 +/– 6.3%) were non-significant. V-DIBH-derived mean shifts for initial patient setup were ≤ 2.7 mm. Random and systematic errors were ≤ 2.1 mm. These results did not differ significantly from historic FB controls. Conclusions When compared with FB, V-DIBH demonstrated high setup accuracy and enabled significant reduction of cardiac doses without compromising the target volume coverage. Differences in lung doses were non-significant.

MRI reduces variation of contouring for boost clinical target volume in breast cancer patients without surgical clips in the tumour bed

Abstract Background Omitting the placement of clips inside tumour bed during breast cancer surgery poses a challenge for delineation of lumpectomy cavity clinical target volume (CTVLC). We aimed to quantify inter-observer variation and accuracy for CT- and MRI-based segmentation of CTVLC in patients without clips. Patients and methods CT- and MRI-simulator images of 12 breast cancer patients, treated by breast conserving surgery and radiotherapy, were included in this study. Five radiation oncologists recorded the cavity visualization score (CVS) and delineated CTVLC on both modalities. Expert-consensus (EC) contours were delineated by a senior radiation oncologist, respecting opinions of all observers. Inter-observer volumetric variation and generalized conformity index (CIgen) were calculated. Deviations from EC contour were quantified by the accuracy index (AI) and inter-delineation distances (IDD). Results Mean CVS was 3.88 +/− 0.99 and 3.05 +/− 1.07 for MRI and CT, respectively (p = 0.001). Mean volumes of CTVLC were similar: 154 +/− 26 cm3 on CT and 152 +/− 19 cm3 on MRI. Mean CIgen and AI were superior for MRI when compared with CT (CIgen: 0.74 +/− 0.07 vs. 0.67 +/− 0.12, p = 0.007; AI: 0.81 +/− 0.04 vs. 0.76 +/− 0.07; p = 0.004). CIgen and AI increased with increasing CVS. Mean IDD was 3 mm +/− 1.5 mm and 3.6 mm +/− 2.3 mm for MRI and CT, respectively (p = 0.017). Conclusions When compared with CT, MRI improved visualization of post-lumpectomy changes, reduced interobserver variation and improved the accuracy of CTVLC contouring in patients without clips in the tumour bed. Further studies with bigger sample sizes are needed to confirm our findings.

SU‐E‐T‐169: Initial Investigation into the Use of Optically Stimulated Luminescent Dosimeters (OSLDs) for In‐Vivo Dosimetry of TBI Patients

PURPOSE To report on an initial investigation into the use of optically stimulated luminescent dosimeters (OSLDs) for in-vivo dosimetry for total body irradiation (TBI) treatments. Specifically, we report on the determination of angular dependence, sensitivity correction factors and the dose calibration factors. METHODS The OSLD investigated in our work was InLight/OSL nanoDot dosimeters (Landauer Inc.). Nanodots are 5 mm diameter, 0.2 mm thick disk-shaped Carbon-doped Al2O3, and were read using a Landauer InLight microstar reader and associated software.OSLDs were irradiated under two setup conditions: a) typical clinical reference conditions (95cm SSD, 5cm depth in solid water, 10×10 cm field size), and b) TBI conditions (520cm SSD, 5cm depth in solid water, 40×40 cm field size,). The angular dependence was checked for angles ranging ±60 degree from normal incidence. In order to directly compare the sensitivity correction factors, a common dose was delivered to the OSLDs for the two setups. Pre- and post-irradiation readings were acquired. OSLDs were optically annealed under various techniques (1) by keeping over a film view box, (2) Using multiple scan on a flat bed optical scanner and (3) Using natural room light. RESULTS Under reference conditions, the calculated sensitivity correction factors of the OSLDs had a SD of 2.2% and a range of 5%. Under TBI conditions, the SD increased to 3.4% and the range to 6.0%. The variation in sensitivity correction factors between individual OSLDs across the two measurement conditions was up to 10.3%. Angular dependence of less than 1% is observed. The best bleaching method we found is to keep OSLDs for more than 3 hours on a film viewer which will reduce normalized response to less than 1%. CONCLUSIONS In order to obtain the most accurate results when using OSLDs for in-vivo dosimetry for TBI treatments, sensitivity correction factors and dose calibration factors should all be determined under clinical TBI conditions.

Geometric accuracy of the MR imaging techniques in the presence of motion

Abstract Magnetic Resonance Imaging (MRI) is increasingly being used for improving tumor delineation and tumor tracking in the presence of respiratory motion. The purpose of this work is to design and build an MR compatible motion platform and to use it for evaluating the geometric accuracy of MR imaging techniques during respiratory motion. The motion platform presented in this work is composed of a mobile base made up of a flat plate and four wheels. The mobile base is attached from one end and through a rigid rod to a synchrony motion table by Accuray® placed at the end of the MRI table and from the other end to an elastic rod. The geometric accuracy was measured by placing a control point‐based phantom on top of the mobile base. In‐house software module was used to automatically assess the geometric distortion. The blurring artifact was also assessed by measuring the Full Width Half Maximum (FWHM) of each control point. Our results were assessed for 50, 100, and 150 mm radial distances, with a mean geometric distortion during the superior–inferior motion of 0.27, 0.41, and 0.55 mm, respectively. Adding the anterior–posterior motion, the mean geometric distortions increased to 0.4, 0.6, and 0.8 mm. Blurring was observed during motion causing an increase in the FWHM of ≈30%. The platform presented in this work provides a valuable tool for the assessment of the geometric accuracy and blurring artifact for MR during motion. Although the main objective was to test the spatial accuracy of an MR system during motion, the modular aspect of the presented platform enables the use of any commercially available phantom for a full quality control of the MR system during motion.

Off-line magnetic resonance imaging navigation of cervix cancer brachytherapy in patients with risk factors for uterine perforation

Purpose There are no reports on pre-insertion identification of cervix cancer patients at risk for uterine perforation during brachytherapy (BT). Our aim was to assess the incidence of risk factors in our patient cohort, and assess feasibility of a novel technique of magnetic resonance imaging (MRI)-guided navigation for applicator insertion (NAI) in high-risk cases. Material and methods All patients with locally advanced cervical cancer, treated with image guided adaptive BT at our department between October 2013 and June 2017 were considered for analysis. Tumor characteristics on initial MRI (MRIinitial), pre-BT MRI (MRIpre-BT), and BT MRI (MRIBT) were assessed. Frequency of risk factors (age above 60 years, retroverted/retroflected uterus, tumor necrosis, non-visible cervical orifice, distorted cervical canal) was recorded. Patients with two or more factors underwent MRI guided NAI. Time needed for NAI was estimated and procedure feasibility score assigned using a three-tiered scoring system. Results Twenty-seven patients (98 insertions) were included. Mean tumor volume was 70.2 (± 47.9), 17.8 (± 18.9), and 10.3 (± 9.1) cm3 on MRIinitial, MRIpre-BT, and MRIBT1, respectively (p < 0.05). In 16 (59%) cases, ≥ 1 perforation risk factor was found on MRIpre-BT: distorted canal in 12 (44%), necrosis in 9 (33%), retroverted/retroflected uterus in 8 (30%) cases. Nine (33%) patients had ≥ 2 risk factors and underwent MRI guided NAI. Additional time to perform NAI was estimated at 105 minutes, and feasibility score was 1 in all cases. There were no cases of uterine perforation. Conclusions Using pre-insertion MRI, we found ≥ 2 risk factors for uterine perforation in 1/3 of patients. Off-line MRI navigation was feasible and enabled non-complicated insertion in all cases. Further studies with larger sample size are warranted to assess its clinical efficacy.

MO‐F‐CAMPUS‐J‐05: Toward MRI‐Only Radiotherapy: Novel Tissue Segmentation and Pseudo‐CT Generation Techniques Based On T1 MRI Sequences

Purpose: To develop and validate a 4 class tissue segmentation approach (air cavities, background, bone and soft-tissue) on T1 -weighted brain MRI and to create a pseudo-CT for MRI-only radiation therapy verification. Methods: Contrast-enhanced T1-weighted fast-spin-echo sequences (TR = 756ms, TE= 7.152ms), acquired on a 1.5T GE MRI-Simulator, are used.MRIs are firstly pre-processed to correct for non uniformity using the non parametric, non uniformity intensity normalization algorithm. Subsequently, a logarithmic inverse scaling log(1/image) is applied, prior to segmentation, to better differentiate bone and air from soft-tissues. Finally, the following method is enrolled to classify intensities into air cavities, background, bone and soft-tissue:Thresholded region growing with seed points in image corners is applied to get a mask of Air+Bone+Background. The background is, afterward, separated by the scan-line filling algorithm. The air mask is extracted by morphological opening followed by a post-processing based on knowledge about air regions geometry. The remaining rough bone pre-segmentation is refined by applying 3D geodesic active contours; bone segmentation evolves by the sum of internal forces from contour geometry and external force derived from image gradient magnitude.Pseudo-CT is obtained by assigning −1000HU to air and background voxels, performing linear mapping of soft-tissue MR intensities in [-400HU, 200HU] and inverse linear mapping of bone MR intensities in [200HU, 1000HU]. Results: Three brain patients having registered MRI and CT are used for validation. CT intensities classification into 4 classes is performed by thresholding. Dice and misclassification errors are quantified. Correct classifications for soft-tissue, bone, and air are respectively 89.67%, 77.8%, and 64.5%. Dice indices are acceptable for bone (0.74) and soft-tissue (0.91) but low for air regions (0.48). Pseudo-CT produces DRRs with acceptable clinical visual agreement to CT-based DRR. Conclusion: The proposed approach makes it possible to use T1-weighted MRI to generate accurate pseudo-CT from 4-class segmentation.

SU-E-I-59: Virtual Phantom: A First Step to a Comprehensive Automated Quality Control Program for Magnetic Resonance Image Guided Applications

Purpose: To design a virtual phantom for the calculation of quality metrics required for controlling MRI guided applications and to develop software tools to automatically perform the required controls. Methods: The virtual phantom is designed using 3D objects arranged in different dispositions in the space. Six cylinders of 6mm diameter and 500mm length, tilted 30 degrees with respect to the Y axis, are used for slice thickness and location. These cylinders are placed in coronal planes distanced by 60mm. Their width and position allow measuring the slice location and thickness. A 15*15*15 mm3 cube centered at the origin and rotated 3 degrees with respect to the Z axis is used for the calculation of the vertical and horizontal Modulation Transfer Function. Spheres of 10 mm diameter distanced by 20 mm covering a field of view of 500mm are used for estimating the in-plane and through-plane geometric distortion. Sets of cylinders with lengths ranging from 0.5mm to 2.0mm and diameters ranging from 4mm to 10mm are used for the low contrast. The Low Contrast test assesses the number of cylinders that can be detected.In-house Java based software is developed and validated. The software automatically corrects positioning errors before constructing ROIs and acquiring measurements with respect to each control. Results: The virtual phantom was used to produce a large set of 2D DICOM images. Results showed that the software robustly calculated positions and distances with sub-voxel accuracy. For a 500mm FOV and a 0.9mm pixel size, mean errors were in the order of 0.15mm. Conclusion: This virtual phantom is an important step for building a new physical phantom and for validating new medical image processing algorithms. These tools represent a first step toward designing a comprehensive phantom and software for a complete Quality Control program for MRI guided application.