Joshua Kim

Technical Note: Characterization and correction of gradient nonlinearity induced distortion on a 1.0 T open bore MRâ SIM

PURPOSE Distortions in magnetic resonance imaging (MRI) compromise spatial fidelity, potentially impacting delineation and dose calculation. We characterized 2D and 3D large field of view (FOV), sequence-independent distortion at various positions in a 1.0 T high-field open MR simulator (MR-SIM) to implement correction maps for MRI treatment planning. METHODS A 36 × 43 × 2 cm(3) phantom with 255 known landmarks (∼1 mm(3)) was scanned using 1.0 T high-field open MR-SIM at isocenter in the transverse, sagittal, and coronal axes, and a 465 × 350 × 168 mm(3) 3D phantom was scanned by stepping in the superior-inferior direction in three overlapping positions to achieve a total 465 × 350 × 400 mm(3) sampled FOV yielding >13 800 landmarks (3D Gradient-Echo, TE/TR/α = 5.54 ms/30 ms/28°, voxel size = 1 × 1 × 2 mm(3)). A binary template (reference) was generated from a phantom schematic. An automated program converted MR images to binary via masking, thresholding, and testing for connectivity to identify landmarks. Distortion maps were generated by centroid mapping. Images were corrected via warping with inverse distortion maps, and temporal stability was assessed. RESULTS Over the sampled FOV, non-negligible residual gradient distortions existed as close as 9.5 cm from isocenter, with a maximum distortion of 7.4 mm as close as 23 cm from isocenter. Over six months, average gradient distortions were -0.07 ± 1.10 mm and 0.10 ± 1.10 mm in the x and y directions for the transverse plane, 0.03 ± 0.64 and -0.09 ± 0.70 mm in the sagittal plane, and 0.4 ± 1.16 and 0.04 ± 0.40 mm in the coronal plane. After implementing 3D correction maps, distortions were reduced to <1 pixel width (1 mm) for all voxels up to 25 cm from magnet isocenter. CONCLUSIONS Inherent distortion due to gradient nonlinearity was found to be non-negligible even with vendor corrections applied, and further corrections are required to obtain 1 mm accuracy for large FOVs. Statistical analysis of temporal stability shows that sequence independent distortion maps are consistent within six months of characterization.

Implementation of a Novel Algorithm For Generating Synthetic CT Images From Magnetic Resonance Imaging Data Sets for Prostate Cancer Radiation Therapy

PURPOSE To describe and evaluate a method for generating synthetic computed tomography (synCT) images from magnetic resonance simulation (MR-SIM) data for accurate digitally reconstructed radiograph (DRR) generation and dose calculations in prostate cancer radiation therapy. METHODS AND MATERIALS A retrospective evaluation was performed in 9 prostate cancer patients who had undergone MR-SIM in addition to CT simulation (CT-SIM). MR-SIM data were used to generate synCT images by using a novel, voxel-based weighted summation approach. A subset of patients was used for weight optimization, and the number of patients to use during optimization was determined. Hounsfield unit (HU) differences between CT-SIM and synCT images were analyzed via mean absolute error (MAE). Original, CT-based treatment plans were mapped onto synCTs. DRRs were generated, and agreement between CT and synCT-generated DRRs was evaluated via Dice similarity coefficient (DSC). Dose was recalculated, and dose-volume metrics and gamma analysis were used to evaluate resulting treatment plans. RESULTS Full field-of-view synCT MAE across all patients was 74.3 ± 10.9 HU with differences from CTs of 2.0 ± 8.1 HU and 11.9 ± 46.7 HU for soft tissue structures (prostate, bladder, and rectum) and femoral bones, respectively. Calculated DSCs for anterior-posterior and lateral DRRs were 0.90 ± 0.04 and 0.92 ± 0.05, respectively. Differences in D99%, mean dose, and maximum dose to the clinical target volume from CT-SIM dose calculations were 0.75% ± 0.35%, 0.63% ± 0.34%, and 0.54% ± 0.33%, respectively, for synCT-generated plans. Gamma analysis (2%/2 mm dose difference/distance to agreement) revealed pass rates of 99.9% ± 0.1% (range, 99.7%-100%). CONCLUSION Generated synCTs enabled accurate DRR generation and dose computation for prostate MR-only simulation. Dose recalculated on synCTs agreed well with original planning distributions. Further validation using a larger patient cohort is warranted.

Initial clinical experience with a radiation oncology dedicated open 1.0T MR-simulation

The purpose of this study was to describe our experience with 1.0T MR‐SIM including characterization, quality assurance (QA) program, and features necessary for treatment planning. Staffing, safety, and patient screening procedures were developed. Utilization of an external laser positioning system (ELPS) and MR‐compatible couchtop were illustrated. Spatial and volumetric analyses were conducted between CT‐SIM and MR‐SIM using a stereotactic QA phantom with known landmarks and volumes. Magnetic field inhomogeneity was determined using phase difference analysis. System‐related, in‐plane distortion was evaluated and temporal changes were assessed. 3D distortion was characterized for regions of interest (ROIs) 5–20 cm away from isocenter. American College of Radiology (ACR) recommended tests and impact of ELPS on image quality were analyzed. Combined ultrashort echotime Dixon (UTE/Dixon) sequence was evaluated. Amplitude‐triggered 4D MRI was implemented using a motion phantom (2–10 phases, ~2 cm excursion, 3–5 s periods) and a liver cancer patient. Duty cycle, acquisition time, and excursion were evaluated between maximum intensity projection (MIP) datasets. Less than 2% difference from expected was obtained between CT‐SIM and MR‐SIM volumes, with a mean distance of <0.2 mm between landmarks. Magnetic field inhomogeneity was <2 ppm. 2D distortion was <2 mm over 28.6–33.6 mm of isocenter. Within 5 cm radius of isocenter, mean 3D geometric distortion was 0.59±0.32 mm (maximum=1.65 mm) and increased 10–15 cm from isocenter (mean=1.57±1.06 mm, maximum=6.26 mm). ELPS interference was within the operating frequency of the scanner and was characterized by line patterns and a reduction in signal‐to‐noise ratio (4.6–12.6% for TE=50−150 ms). Image quality checks were within ACR recommendations. UTE/Dixon sequences yielded detectability between bone and air. For 4D MRI, faster breathing periods had higher duty cycles than slow (50.4% (3 s) and 39.4% (5 s), p<0.001) and ~ fourfold acquisition time increase was measured for ten‐phase versus two‐phase. Superior–inferior object extent was underestimated 8% (6 mm) for two‐phase as compared to ten‐phase MIPs, although <2% difference was obtained for ≥4 phases. 4D MRI for a patient demonstrated acceptable image quality in ~7 min. MR‐SIM was integrated into our workflow and QA procedures were developed. Clinical applicability was demonstrated for 4D MRI and UTE imaging to support MR‐SIM for single modality treatment planning. PACS numbers: 87.56.Fc, 87.61.‐c, 87.57.cp

Technical Note: Characterization and correction of gradient nonlinearity induced distortion on a 1.0 T open bore MR‐SIM

PURPOSE Distortions in magnetic resonance imaging (MRI) compromise spatial fidelity, potentially impacting delineation and dose calculation. We characterized 2D and 3D large field of view (FOV), sequence-independent distortion at various positions in a 1.0 T high-field open MR simulator (MR-SIM) to implement correction maps for MRI treatment planning. METHODS A 36 × 43 × 2 cm(3) phantom with 255 known landmarks (∼1 mm(3)) was scanned using 1.0 T high-field open MR-SIM at isocenter in the transverse, sagittal, and coronal axes, and a 465 × 350 × 168 mm(3) 3D phantom was scanned by stepping in the superior-inferior direction in three overlapping positions to achieve a total 465 × 350 × 400 mm(3) sampled FOV yielding >13 800 landmarks (3D Gradient-Echo, TE/TR/α = 5.54 ms/30 ms/28°, voxel size = 1 × 1 × 2 mm(3)). A binary template (reference) was generated from a phantom schematic. An automated program converted MR images to binary via masking, thresholding, and testing for connectivity to identify landmarks. Distortion maps were generated by centroid mapping. Images were corrected via warping with inverse distortion maps, and temporal stability was assessed. RESULTS Over the sampled FOV, non-negligible residual gradient distortions existed as close as 9.5 cm from isocenter, with a maximum distortion of 7.4 mm as close as 23 cm from isocenter. Over six months, average gradient distortions were -0.07 ± 1.10 mm and 0.10 ± 1.10 mm in the x and y directions for the transverse plane, 0.03 ± 0.64 and -0.09 ± 0.70 mm in the sagittal plane, and 0.4 ± 1.16 and 0.04 ± 0.40 mm in the coronal plane. After implementing 3D correction maps, distortions were reduced to <1 pixel width (1 mm) for all voxels up to 25 cm from magnet isocenter. CONCLUSIONS Inherent distortion due to gradient nonlinearity was found to be non-negligible even with vendor corrections applied, and further corrections are required to obtain 1 mm accuracy for large FOVs. Statistical analysis of temporal stability shows that sequence independent distortion maps are consistent within six months of characterization.

Efficacy of suicide gene therapy in hypoxic rat 9L glioma cells.

Viral vector mediated suicide gene therapy (SGT) involving thymidine kinase (TK) or cytosine deaminase (CD) have considerable promise in the treatment of malignant brain tumors. An unresolved issue is to what extent tumor hypoxia influences the outcome of SGT since brain tumors characterized by regions of hypoxia have potentially reduced cellular metabolism and SGT’s cytotoxicity is manifest through cellular metabolism. We studied in vitro and in vivo, the effect of hypoxia on the cytotoxicity of SGT in rat 9L glioma cells. Neither acute nor chronic hypoxia affected the cell killing of SGT by TK or CD. In vivo confirmation that SGT efficacy was not adversely affected by tumor hypoxia using the hypoxic cell marker pimonidazole was shown by the absence of a change in tumor hypoxia by SGT. These studies support the use of SGT utilizing either TK or CD gene strategies even when tumors are characterized by a hypoxic microenvironment.

Evaluation of a magnetic resonance guided linear accelerator for stereotactic radiosurgery treatment

INTRODUCTION The purpose of this study was to investigate the systematic localization accuracy, treatment planning capability, and delivery accuracy of an integrated magnetic resonance imaging guided Linear Accelerator (MR-Linac) platform for stereotactic radiosurgery. MATERIALS AND METHODS The phantom for the end-to-end test comprises three different compartments: a rectangular MR/CT target phantom, a Winston-Lutz cube, and a rectangular MR/CT isocenter phantom. Hidden target tests were performed at gantry angles of 0, 90, 180, and 270 degrees to quantify the systematic accuracy. Five patient plans with a total of eleven lesions were used to evaluate the dosimetric accuracy. Single-isocenter IMRT treatment plans using 10-15 coplanar beams were generated to treat the multiple metastases. RESULTS The end-to-end localization accuracy of the system was 1.0 ± 0.1 mm. The conformity index, homogeneity index and gradient index of the plans were 1.26 ± 0.22, 1.22 ± 0.10, and 5.38 ± 1.44, respectively. The average absolute point dose difference between measured and calculated dose was 1.64 ± 1.90%, and the mean percentage of points passing the 3%/1 mm gamma criteria was 96.87%. CONCLUSIONS Our experience demonstrates that excellent plan quality and delivery accuracy was achievable on the MR-Linac for treating multiple brain metastases with a single isocenter.

Technical Note: Evaluation of plastic scintillator detector for small field stereotactic patient-specific quality assurance

Purpose: To evaluate the performance of a commercial plastic scintillator detector (PSD) for small‐field stereotactic patient‐specific quality assurance (QA) measurements using flattening‐filter‐free beam. Methods: A total of 10 spherical targets [volume range: (0.03 cc–2 cc)] were planned with two techniques: (a) dynamic conformal arc (DCA‐10 plans) and (b) volumetric modulated arc therapy (VMAT‐10 plans). All plans were generated using Varian Eclipse treatment planning system, and AcurosXB v.13 algorithm in 1.0 mm grid size. Additionally, 14 previously treated cranial and spine SRS plans were evaluated [6 DCA, 8 VMAT, volume range: (0.04 cc–119.02 cc)]. Plan modulation was quantified via two metrics: MU per prescription dose (MU/Rx) and Average Leaf Pair Opening (ALPO). QA was performed on the Varian Edge linear accelerator equipped with HDMLC. Three detectors were used: (a) PinPoint ion chamber (PTW; active volume 0.015 cc), (b) Exradin W1 PSD (Standard Imaging; active volume 0.002 cc), and (c) Gafchromic EBT3 film (Ashland). PinPoint chamber and PSD were positioned perpendicular to beam axis in a Lucy phantom (Standard Imaging); films were placed horizontally capturing the coronal plane. Results: PSD, film, and PinPoint chamber measured average differences of 1.00 ± 1.54%, 1.30 ± 1.69%, and −0.66 ± 2.36%, respectively, compared to AcurosXB dose calculation. As the target volume decreased, PinPoint chamber measured lower doses (maximum −5.07% at 0.07 cc target), while PSD and film measured higher doses (2.87% and 2.54% at 0.03 cc target) than AcurosXB. Film agreed with the benchmark detector PSD by an average difference of 0.31 ± 1.20%, but suffered from larger uncertainty; PinPoint chamber underestimated dose by more than 4% for targets smaller than 0.2 cc. Taking PSD as the measurement standard, DCA plans achieved good QA results across all volumes studied, with an average of −0.07 ± 0.89%; for VMAT plans, PSD measured consistently higher dose (1.95 ± 1.36%) than AcurosXB. Correlation study revealed that plan modulation quantified by both MU/Rx and ALPO correlated significantly with QA results. Conclusion: Among all three detectors, PSD demonstrated superior performances in plans with small fields and heavy modulation. High consistency and low uncertainty made PSD a suitable detector for clinical routine SRS QA. PinPoint chamber should be avoided for targets smaller than 0.2 cc; film dosimetry can be utilized with careful evaluation of its uncertainty bracket. Compared to PSD measurements, AcurosXB calculation demonstrated high accuracy for nonmodulated small fields. The positive correlation between plan modulation and QA discrepancy calls for our attention for clinical SRS plans with high modulation.

Using Synthetic CT for Partial Brain Radiation Therapy: Impact on Image Guidance

PURPOSE Recent advancements in synthetic computed tomography (synCT) from magnetic resonance (MR) imaging data have made MRI-only treatment planning feasible in the brain, although synCT performance for image guided radiation therapy (IGRT) is not well understood. This work compares geometric equivalence of digitally reconstructed radiographs (DRRs) from CTs and synCTs for brain cancer patients and quantifies performance for partial brain IGRT. METHODS AND MATERIALS Ten brain cancer patients (12 lesions, 7 postsurgical) underwent MR-SIM and CT-SIM. SynCTs were generated by combining ultra-short echo time, T1, T2, and fluid attenuation inversion recovery datasets using voxel-based weighted summation. SynCT and CT DRRs were compared using patient-specific thresholding and assessed via overlap index, Dice similarity coefficient, and Jaccard index. Planar IGRT images for 22 fractions were evaluated to quantify differences between CT-generated DRRs and synCT-generated DRRs in 6 quadrants. Previously validated software was implemented to perform 2-dimensional (2D)-2D rigid registrations using normalized mutual information. Absolute (planar image/DRR registration) and relative (differences between synCT and CT DRR registrations) shifts were calculated for each axis and 3-dimensional vector difference. A total of 1490 rigid registrations were assessed. RESULTS DRR agreements in anteroposterior and lateral views for overlap index, Dice similarity coefficient, and Jaccard index were >0.95. Normalized mutual information results were equivalent in 75% of quadrants. Rotational registration results were negligible (<0.07°). Statistically significant differences between CT and synCT registrations were observed in 9/18 matched quadrants/axes (P < .05). The population average absolute shifts were 0.77 ± 0.58 and 0.76 ± 0.59 mm for CT and synCT, respectively, for all axes/quadrants. Three-dimensional vectors were <2 mm in 77.7 ± 10.8% and 76.5 ± 7.2% of CT and synCT registrations, respectively. SynCT DRRs were sensitive in postsurgical cases (vector displacements >2 mm in affected quadrants). CONCLUSIONS DRR synCT geometry was robust. Although statistically significant differences were observed between CT and synCT registrations, results were not clinically significant. Future work will address synCT generation in postsurgical settings.

FMEA of MR-only Treatment Planning in the Pelvis

Purpose To evaluate the implementation of a magnetic resonance (MR)-only workflow (ie, implementing MR simulation as the primary planning modality) using failure mode and effects analysis (FMEA) in comparison with a conventional multimodality (MR simulation in conjunction with computed tomography simulation) workflow for pelvis external beam planning. Methods and Materials To perform the FMEA, a multidisciplinary 9-member team was assembled and developed process maps, identified potential failure modes (FMs), and assigned numerical values to the severity (S), frequency of occurrence (O), and detectability (D) of those FMs. Risk priority numbers (RPNs) were calculated via the product of S, O, and D as a metric for evaluating relative patient risk. An alternative 3-digit composite number (SOD) was computed to emphasize high-severity FMs. Fault tree analysis identified the causality chain leading to the highest-severity FM. Results Seven processes were identified, 3 of which were shared between workflows. Image fusion and target delineation subprocesses using the conventional workflow added 9 and 10 FMs, respectively, with 6 RPNs >100. By contrast, synthetic computed tomography generation introduced 3 major subprocesses and propagated 46 unique FMs, 15 with RPNs >100. For the conventional workflow, the largest RPN scores were introduced by image fusion (RPN range, 120-192). For the MR-only workflow, the highest RPN scores were from inaccuracies in target delineation resulting from misinterpretation of MR images (RPN = 240) and insufficient management of patient- and system-level distortions (RPN = 210 and 168, respectively). Underestimation (RPN = 140) or overestimation (RPN = 192) of bone volume produced higher RPN scores. The highest SODs for both workflows were related to changes in target location because of internal anatomy changes (conventional = 961, MR-only = 822). Conclusions FMEA identified areas for mitigating risk in MR-only pelvis RTP, and SODs identified high-severity process modes. Efforts to develop a quality management program to mitigate high FMs are underway.

WE‐G‐141‐07: Feasibility of Using a Grid to Detect and Correct the Geometric Variations in Flat‐Panel Based Cone Beam CT

PURPOSE To study the feasibility of using a grid, shown previously to mitigate the scatter problem, to detect and correct geometric variations in flat-panel-based cone beam CT (CBCT) imaging. METHODS The study was performed on the CBCT system of a Varian Trilogy linac. The grid was located before the imaging object and attached to the bowtie-filter with a source-to-grid distance (SGD) of 40-50 cm. We first studied the correlations in geometric variation between the grid, bowtie-filter, source and imager. The grid position was represented by BBs in a plastic-plate replacing the grid. The bowtie-filter position was determined by a BB placed on its surface. The source and imager geometries were derived from the projections of 16 BBs in an IsoCal phantom using an in-house-developed computer program. Multiple scans were performed with different SGD in a 2-month period. Correlations were analyzed and used to predict the geometric variations of the bowtie-filter, source and imager from the grid variation obtained using a grid in a CatPhan scan. RESULTS The grid and bowtie-filter showed exactly same variations (up to 5 mm) in Y-direction during a 360° rotation, with a sine-like pattern superimposed by 1-2 mm random vibrations. The variation patterns were also sine-like curves in the X-direction. However, they were smooth and repeatable, with differences in phase and amplitude between the two, and the amplitude varying with SGD. The source and imager showed similar and repeatable 0.3-0.5 mm sine-like variation patterns. Based on these results, a model was developed to predict and correct the geometric variations in a CBCT scan with the grid. CONCLUSION The preliminary results suggest that it is able to predict geometric variations of the CBCT system. Future study is underway to verify whether correction of these variations will improve the spatial resolution and mitigate the bowtie-filter variation-induced crescent artifact. The study is supported by NIH Grant Number 5 R01 CA166948-02.