Magnetic resonance-only simulation and dose calculation in external beam radiation therapy: a feasibility study for pelvic cancers

BACKGROUND

The clinical feasibility of using pseudo-computed tomography (pCT) images derived from magnetic resonance (MR) images for external bean radiation therapy (EBRT) planning for prostate cancer patients has been well demonstrated. This paper investigates the feasibility of applying an MR-derived, pCT planning approach to additional types of cancer in the pelvis.

MATERIAL AND METHODS

Fifteen patients (five prostate cancer patients, five rectal cancer patients, and five gynecological cancer patients) receiving EBRT at Turku University Hospital (Turku, Finland) were included in the study. Images from an MRCAT (Magnetic Resonance for Calculating ATtenuation, Philips, Vantaa, Finland) pCT method were generated as a part of a clinical MR-simulation procedure. Dose calculation accuracy was assessed by comparing the pCT-based calculation with a CT-based calculation. In addition, the degree of geometric accuracy was studied.

RESULTS

The median relative difference of PTV mean dose between CT and pCT images was within 0.8% for all tumor types. When assessing the tumor site-specific accuracy, the median [range] relative dose differences to the PTV mean were 0.7 [-0.11;1.05]% for the prostate cases, 0.3 [-0.25;0.57]% for the rectal cases, and 0.09 [-0.69;0.25]% for the gynecological cancer cases. System-induced geometric distortion was measured to be less than 1 mm for all PTV volumes and the effect on the PTV median dose was less than 0.1%.

CONCLUSIONS

According to the comparison, using pCT for clinical EBRT planning and dose calculation in the three investigated types of pelvic cancers is feasible. Further studies are required to demonstrate the applicability to a larger cohort of patients.

Characterization of spatial distortion in a 0.35 T MRI-guided radiotherapy system

Spatial distortion results in image deformation that can degrade accurate targeting and dose calculations in MRI-guided adaptive radiotherapy. The authors present a comprehensive assessment of a 0.35 T MRI-guided radiotherapy system's spatial distortion using two commercially-available phantoms with regularly spaced markers. Images of the spatial integrity phantoms were acquired using five clinical protocols on the MRI-guided radiotherapy machine with the radiotherapy gantry positioned at various angles. Software was developed to identify and localize all phantom markers using a template matching approach. Rotational and translational corrections were implemented to account for imperfect phantom alignment. Measurements were made to assess uncertainties arising from susceptibility artifacts, image noise, and phantom construction accuracy. For a clinical 3D imaging protocol with a 1.5 mm reconstructed slice thickness, 100% of spheres within a 50 mm radius of isocenter had a 3D deviation of 1 mm or less. Of the spheres within 100 mm of isocenter, 99.9% had a 3D deviation less than 1 mm. 94.8% and 100% of the spheres within 175 mm were found to be within 1 mm and 2 mm of the expected positions in 3D respectively. Maximum 3D distortions within 50 mm, 100 mm and 175 mm of isocenter were 0.76 mm, 1.15 mm and 1.88 mm respectively. Distortions present in images acquired using the real-time imaging sequence were less than 1 mm for 98.1% and 95.0% of the cylinders within 50 mm and 100 mm of isocenter. The corresponding maximum distortion in these regions was 1.10 mm and 1.67 mm. These results may be used to inform appropriate planning target volume (PTV) margins for 0.35 T MRI-guided radiotherapy. Observed levels of spatial distortion should be explicitly considered when using PTV margins of 3 mm or less or in the case of targets displaced from isocenter by more than 50 mm.

Stereotactic radiosurgery planning of vestibular schwannomas: Is MRI at 3 Tesla geometrically accurate?

Purpose

MRI is a mandatory requirement to accurately plan Stereotactic Radiosurgery (SRS) for Vestibular Schwannomas. However, MRI may be distorted due not only to inhomogeneity of the static magnetic field and gradients but also due to susceptibility‐induced effects, which are more prominent at higher magnetic fields. We assess geometrical distortions around air spaces and consider MRI protocol requirements for SRS planning at 3 T.

Methods

Hardware‐related distortion and the effect of incorrect shimming were investigated with structured test objects. The magnetic field was mapped over the head on five volunteers to assess susceptibility‐related distortion in the naso‐oro‐pharyngeal cavities (NOPC) and around the internal ear canal (IAC).

Results

Hardware‐related geometric displacements were found to be less than 0.45 mm within the head volume, after distortion correction. Shimming errors can lead to displacements of up to 4 mm, but errors of this magnitude are unlikely to arise in practice. Susceptibility‐related field inhomogeneity was under 3.4 ppm, 2.8 ppm, and 2.7 ppm for the head, NOPC region and IAC region, respectively. For the SRS planning protocol (890 Hz/pixel, approximately 1 mm³ isotropic), susceptibility‐related displacements were less than 0.5 mm (head), and 0.4 mm (IAC and NOPC). Large displacements are possible in MRI examinations undertaken with lower receiver bandwidth values, commonly used in clinical MRI. Higher receiver bandwidth makes the protocol less vulnerable to sub‐optimal shimming. The shimming volume and the CT‐MR co‐registration must be considered jointly.

Conclusion

Geometric displacements can be kept under 1 mm in the vicinity of air spaces within the head at 3 T with appropriate setting of the receiver bandwidth, correct shimming and employing distortion correction.

Investigation of the 4D composite MR image distortion field associated with tumor motion for MR-guided radiotherapy

PURPOSE

Magnetic resonance (MR) images are affected by geometric distortions due to the specifics of the MR scanner and patient anatomy. Quantifying the distortions associated with mobile tumors is particularly challenging due to real anatomical changes in the tumor's volume, shape, and relative location within the MR imaging volume. In this study, the authors investigate the 4D composite distortion field, which combines the effects of the susceptibility-induced and system-related distortion fields, experienced by mobile lung tumors.

METHODS

The susceptibility (χ) effects were numerically simulated for two specific scenarios: (a) a full motion cycle of a lung tumor due to breathing as depicted on ten phases of a 4D CBCT data set and (b) varying the tumor size and location in lung tissue via a synthetically generated sphere with variable diameter (4-80 mm). The χ simulation procedure relied on the segmentation and generation of 3D susceptibility (χ) masks and computation of the magnetic field by means of finite difference methods. A system-related distortion field, determined with a phantom and image processing algorithm, was used as a reference. The 4D composite distortion field was generated as the vector summation of the χ-induced and system-related fields. The analysis was performed for two orientations of the main magnetic field (B0), which correspond to several MRIgRT system configurations. Specifically, B0 was set along the z-axis as in the case of a cylindrical-bore scanner and in the (x,y)-plane as for a biplanar MR. Computations were also performed for a full revolution at 15° increments in the case of a rotating biplanar magnet. Histograms and metrics such as maximum, mean, and range were used to evaluate the characteristics of the 4D distortion field.

RESULTS

The χ-induced field depends on the change in volume and shape of the moving tumor as well as the local surrounding anatomy. In the case of system-related distortions, the tumor experiences increased field perturbations as it moves further away from the MR isocenter. For a mobile lung tumor, the 4D composite field, corresponding to a 1.5 T field and a readout gradient of 5 mT/m, amounts to 3.0 and 2.8 mm for the MRIgRT system designs featuring B0 oriented along the z-axis (cylindrical-bore scanner) and in the (x,y)-plane (biplanar scanner), respectively. For a rotating biplanar scanner, the composite distortion field varied nonlinearly with the rotation angle. Overall, the dominant contribution to the composite field was from the system-related distortion field. The tumor centroid experienced a systematic shift of 2 mm and showed a negligible perturbation for different B0 values. The dependency on the tumor size was also investigated, namely the max values varied from 1.2 to 2.5 mm for spherical volumes with a diameter between 4 and 80 mm.

CONCLUSIONS

The composite distortion field requires adequate quantification for lung radiation therapy applications such as treatment planning, pretreatment patient setup verification, and real-time treatment delivery guidance. For certain scenarios such as small tumor volumes, the spatial distortions may be corrected by applying systematic shifts derived from a single tumor motion phase. In the case of high readout gradients common to fast imaging applications, the χ distortions were found to be less than 1 mm irrespective of scanner configuration.

Anatomic and dosimetric changes in patients with head and neck cancer treated with an integrated MRI-tri-60Co teletherapy device

OBJECTIVE

Prior studies have relied on CT to assess alterations in anatomy among patients undergoing radiation for head and neck cancer. We sought to determine the feasibility of using MRI-based image-guided radiotherapy to quantify these changes and to ascertain their potential dosimetric implications.

METHODS

6 patients with head and neck cancer were treated with intensity-modulated radiotherapy (IMRT) on a novel tri-⁶⁰Co teletherapy system equipped with a 0.35-T MRI (VR, ViewRay Incorporated, Oakwood Village, OH) to 66-70 Gy in 33 fractions (fx). Pre-treatment MRIs on Fx 1, 5, 10, 15, 20, 25, 30 and 33 were imported into a contouring interface, where the primary gross tumour volume (GTV) and parotid glands were delineated. The centre of mass (COM) shifts for these structures were assessed relative to Day 1. Dosimetric data were co-registered with the MRIs, and doses to the GTV and parotid glands were assessed.

RESULTS

Primary GTVs decreased significantly over the course of IMRT (median % volume loss, 38.7%; range, 29.5-72.0%; p < 0.05) at a median rate of 1.2%/fx (range, 0.92-2.2%/fx). Both the ipsilateral and contralateral parotid glands experienced significant volume loss (p < 0.05, for all) and shifted medially during IMRT. Weight loss correlated significantly with parotid gland volume loss and medial COM shift (p < 0.05).

CONCLUSION

Integrated on-board MRI can be used to accurately contour and analyze primary GTVs and parotid glands over the course of IMRT. COM shifts and significant volume reductions were observed, confirming the results of prior CT-based exercises. Advances in knowledge: The superior resolution of on-board MRI may facilitate online adaptive replanning in the future.

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.

Image Guided Radiation Therapy Using Synthetic Computed Tomography Images in Brain Cancer

PURPOSE

The development of synthetic computed tomography (CT) (synCT) derived from magnetic resonance (MR) images supports MR-only treatment planning. We evaluated the accuracy of synCT and synCT-generated digitally reconstructed radiographs (DRRs) relative to CT and determined their performance for image guided radiation therapy (IGRT).

METHODS AND MATERIALS

Magnetic resonance simulation (MR-SIM) and CT simulation (CT-SIM) images were acquired of an anthropomorphic skull phantom and 12 patient brain cancer cases. SynCTs were generated using fluid attenuation inversion recovery, ultrashort echo time, and Dixon data sets through a voxel-based weighted summation of 5 tissue classifications. The DRRs were generated from the phantom synCT, and geometric fidelity was assessed relative to CT-generated DRRs through bounding box and landmark analysis. An offline retrospective analysis was conducted to register cone beam CTs (n=34) to synCTs and CTs using automated rigid registration in the treatment planning system. Planar MV and KV images (n=37) were rigidly registered to synCT and CT DRRs using an in-house script. Planar and volumetric registration reproducibility was assessed and margin differences were characterized by the van Herk formalism.

RESULTS

Bounding box and landmark analysis of phantom synCT DRRs were within 1 mm of CT DRRs. Absolute planar registration shift differences ranged from 0.0 to 0.7 mm for phantom DRRs on all treatment platforms and from 0.0 to 0.4 mm for volumetric registrations. For patient planar registrations, the mean shift differences were 0.4 ± 0.5 mm (range, -0.6 to 1.6 mm), 0.0 ± 0.5 mm (range, -0.9 to 1.2 mm), and 0.1 ± 0.3 mm (range, -0.7 to 0.6 mm) for the superior-inferior (S-I), left-right (L-R), and anterior-posterior (A-P) axes, respectively. The mean shift differences in volumetric registrations were 0.6 ± 0.4 mm (range, -0.2 to 1.6 mm), 0.2 ± 0.4 mm (range, -0.3 to 1.2 mm), and 0.2 ± 0.3 mm (range, -0.2 to 1.2 mm) for the S-I, L-R, and A-P axes, respectively. The CT-SIM and synCT derived margins were <0.3 mm different.

CONCLUSION

DRRs generated by synCT were in close agreement with CT-SIM. Planar and volumetric image registrations to synCT-derived targets were comparable with CT for phantom and patients. This validation is the next step toward MR-only planning for the brain.

Spatial Precision in Magnetic Resonance Imaging-Guided Radiation Therapy: The Role of Geometric Distortion

Because magnetic resonance imaging-guided radiation therapy (MRIgRT) offers exquisite soft tissue contrast and the ability to image tissues in arbitrary planes, the interest in this technology has increased dramatically in recent years. However, intrinsic geometric distortion stemming from both the system hardware and the magnetic properties of the patient affects MR images and compromises the spatial integrity of MRI-based radiation treatment planning, given that for real-time MRIgRT, precision within 2 mm is desired. In this article, we discuss the causes of geometric distortion, describe some well-known distortion correction algorithms, and review geometric distortion measurements from 12 studies, while taking into account relevant imaging parameters. Eleven of the studies reported phantom measurements quantifying system-dependent geometric distortion, while 2 studies reported simulation data quantifying magnetic susceptibility-induced geometric distortion. Of the 11 studies investigating system-dependent geometric distortion, 5 reported maximum measurements less than 2 mm. The simulation studies demonstrated that magnetic susceptibility-induced distortion is typically smaller than system-dependent distortion but still nonnegligible, with maximum distortion ranging from 2.1 to 2.6 mm at a field strength of 1.5 T. As expected, anatomic landmarks containing interfaces between air and soft tissue had the largest distortions. The evidence indicates that geometric distortion reduces the spatial integrity of MRI-based radiation treatment planning and likely diminishes the efficacy of MRIgRT. Better phantom measurement techniques and more effective distortion correction algorithms are needed to achieve the desired spatial precision.

Continuous table acquisition MRI for radiotherapy treatment planning: distortion assessment with a new extended 3D volumetric phantom

PURPOSE

Accurate geometry is required for radiotherapy treatment planning (RTP). When considering the use of magnetic resonance imaging (MRI) for RTP, geometric distortions observed in the acquired images should be considered. While scanner technology and vendor supplied correction algorithms provide some correction, large distortions are still present in images, even when considering considerably smaller scan lengths than those typically acquired with CT in conventional RTP. This study investigates MRI acquisition with a moving table compared with static scans for potential geometric benefits for RTP.

METHODS

A full field of view (FOV) phantom (diameter 500 mm; length 513 mm) was developed for measuring geometric distortions in MR images over volumes pertinent to RTP. The phantom consisted of layers of refined plastic within which vitamin E capsules were inserted. The phantom was scanned on CT to provide the geometric gold standard and on MRI, with differences in capsule location determining the distortion. MRI images were acquired with two techniques. For the first method, standard static table acquisitions were considered. Both 2D and 3D acquisition techniques were investigated. With the second technique, images were acquired with a moving table. The same sequence was acquired with a static table and then with table speeds of 1.1 mm/s and 2 mm/s. All of the MR images acquired were registered to the CT dataset using a deformable B-spline registration with the resulting deformation fields providing the distortion information for each acquisition.

RESULTS

MR images acquired with the moving table enabled imaging of the whole phantom length while images acquired with a static table were only able to image 50%-70% of the phantom length of 513 mm. Maximum distortion values were reduced across a larger volume when imaging with a moving table. Increased table speed resulted in a larger contribution of distortion from gradient nonlinearities in the through-plane direction and an increased blurring of capsule images, resulting in an apparent capsule volume increase by up to 170% in extreme axial FOV regions. Blurring increased with table speed and in the central regions of the phantom, geometric distortion was less for static table acquisitions compared to a table speed of 2 mm/s over the same volume. Overall, the best geometric accuracy was achieved with a table speed of 1.1 mm/s.

CONCLUSIONS

The phantom designed enables full FOV imaging for distortion assessment for the purposes of RTP. MRI acquisition with a moving table extends the imaging volume in the z direction with reduced distortions which could be useful particularly if considering MR-only planning. If utilizing MR images to provide additional soft tissue information to the planning CT, standard acquisition sequences over a smaller volume would avoid introducing additional blurring or distortions from the through-plane table movement.

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.

Magnetic resonance imaging-guided brachytherapy for cervical cancer: initiating a program

Over the past decade, the application of magnetic resonance imaging (MRI) has increased, and there is growing evidence to suggest that improvements in accuracy of target delineation in MRI-guided brachytherapy may improve clinical outcomes in cervical cancer. To implement a high quality image guided brachytherapy program, a multidisciplinary team is required with appropriate expertise as well as an adequate patient load to ensure a sustainable program. It is imperative to know that the most important source of uncertainty in the treatment process is related to target delineation and therefore, the necessity of training and expertise as well as quality assurance should be emphasized. A short review of concepts and techniques that have been developed for implementation and/or improvement of workflow of a MRI-guided brachytherapy program are provided in this document, so that institutions can use and optimize some of them based on their resources to minimize their procedure times.

Extensive pneumatized air cells causing susceptibility artifacts in the petrosus part of the ICA

INTRODUCTION

3D-Time-of-flight magnetic-resonance-angiography (TOF MRA) is an established method in vessel analysis. However, many artifacts that occur may lead to a false diagnosis. This retrospective study evaluates the coherence of MR artifacts to extensive pneumatized air cells surrounding the internal carotid artery (ICA) in the petrosus part of the temporal bone.

MATERIALS AND METHODS

Patients who received 3D-TOF MRA and multidetector helical computed tomography (CT) angiography were registered from April 2012 to April 2013. Of these patients, both ICAs in the petrosus part were analyzed. Vertical maximum intensity projection (MIP) artifacts were graduated as normal, mild to moderate, and severe artifacts. The distinction of the vertical part of the pneumatized air cells was also categorized in three groups, regarding the circumference of the ICA in pneumatization ≤ 90°, between 90° and 180°, and ≥ 180°.

RESULTS

A total of 203 vessels were collected for analysis. The more extensive the pneumatized air cells were present, the more band-like artifacts and pseudostenosis at the vertical portion of the petrosus part of the ICA were registered.

CONCLUSION

Careful examination of the source images and evaluation of the size of the pneumatized air cells with CT scan are essential to avoid false positive diagnosis in the distal petrosus part of the ICA.

Evaluating organ delineation, dose calculation and daily localization in an open-MRI simulation workflow for prostate cancer patients

BACKGROUND

This study describes initial testing and evaluation of a vertical-field open Magnetic Resonance Imaging (MRI) scanner for the purpose of simulation in radiation therapy for prostate cancer. We have evaluated the clinical workflow of using open MRI as a sole modality for simulation and planning. Relevant results related to MRI alignment (vs. CT) reference dataset with Cone-Beam CT (CBCT) for daily localization are presented.

METHODS

Ten patients participated in an IRB approved study utilizing MRI along with CT simulation with the intent of evaluating the MRI-simulation process. Differences in prostate gland volume, seminal vesicles, and penile bulb were assessed with MRI and compared to CT. To evaluate dose calculation accuracy, bulk-density-assignments were mapped onto respective MRI datasets and treated IMRT plans were re-calculated. For image localization purposes, 400 CBCTs were re-evaluated with MRI as the reference dataset and daily shifts compared against CBCT-to-CT registration. Planning margins based on MRI/CBCT shifts were computed using the van Herk formalism.

RESULTS

Significant organ contour differences were noted between MRI and CT. Prostate volumes were on average 39.7% (p = 0.002) larger on CT than MRI. No significant difference was found in seminal vesicle volumes (p = 0.454). Penile bulb volumes were 61.1% higher on CT, without statistical significance (p = 0.074). MRI-based dose calculations with assigned bulk densities produced agreement within 1% with heterogeneity corrected CT calculations. The differences in shift positions for the cohort between CBCT-to-CT registration and CBCT-to-MRI registration are -0.15 ± 0.25 cm (anterior-posterior), 0.05 ± 0.19 cm (superior-inferior), and -0.01 ± 0.14 cm (left-right).

CONCLUSIONS

This study confirms the potential of using an open-field MRI scanner as primary imaging modality for prostate cancer treatment planning simulation, dose calculations and daily image localization.

Magnetic resonance image guided brachytherapy

The application of magnetic resonance image (MRI)-guided brachytherapy has demonstrated significant growth during the past 2 decades. Clinical improvements in cervix cancer outcomes have been linked to the application of repeated MRI for identification of residual tumor volumes during radiotherapy. This has changed clinical practice in the direction of individualized dose administration, and resulted in mounting evidence of improved clinical outcome regarding local control, overall survival as well as morbidity. MRI-guided prostate high-dose-rate and low-dose-rate brachytherapies have improved the accuracy of target and organs-at-risk delineation, and the potential exists for improved dose prescription and reporting for the prostate gland and organs at risk. Furthermore, MRI-guided prostate brachytherapy has significant potential to identify prostate subvolumes and dominant lesions to allow for dose administration reflecting the differential risk of recurrence. MRI-guided brachytherapy involves advanced imaging, target concepts, and dose planning. The key issue for safe dissemination and implementation of high-quality MRI-guided brachytherapy is establishment of qualified multidisciplinary teams and strategies for training and education.

The ViewRay system: magnetic resonance-guided and controlled radiotherapy

A description of the first commercially available magnetic resonance imaging (MRI)-guided radiation therapy (RT) system is provided. The system consists of a split 0.35-T MR scanner straddling 3 (60)Co heads mounted on a ring gantry, each head equipped with independent doubly focused multileaf collimators. The MR and RT systems share a common isocenter, enabling simultaneous and continuous MRI during RT delivery. An on-couch adaptive RT treatment-planning system and integrated MRI-guided RT control system allow for rapid adaptive planning and beam delivery control based on the visualization of soft tissues. Treatment of patients with this system commenced at Washington University in January 2014.

MRI-based treatment planning with pseudo CT generated through atlas registration

PURPOSE

To evaluate the feasibility and accuracy of magnetic resonance imaging (MRI)-based treatment planning using pseudo CTs generated through atlas registration.

METHODS

A pseudo CT, providing electron density information for dose calculation, was generated by deforming atlas CT images previously acquired on other patients. The authors tested 4 schemes of synthesizing a pseudo CT from single or multiple deformed atlas images: use of a single arbitrarily selected atlas, arithmetic mean process using 6 atlases, and pattern recognition with Gaussian process (PRGP) using 6 or 12 atlases. The required deformation for atlas CT images was derived from a nonlinear registration of conjugated atlas MR images to that of the patient of interest. The contrasts of atlas MR images were adjusted by histogram matching to reduce the effect of different sets of acquisition parameters. For comparison, the authors also tested a simple scheme assigning the Hounsfield unit of water to the entire patient volume. All pseudo CT generating schemes were applied to 14 patients with common pediatric brain tumors. The image similarity of real patient-specific CT and pseudo CTs constructed by different schemes was compared. Differences in computation times were also calculated. The real CT in the treatment planning system was replaced with the pseudo CT, and the dose distribution was recalculated to determine the difference.

RESULTS

The atlas approach generally performed better than assigning a bulk CT number to the entire patient volume. Comparing atlas-based schemes, those using multiple atlases outperformed the single atlas scheme. For multiple atlas schemes, the pseudo CTs were similar to the real CTs (correlation coefficient, 0.787-0.819). The calculated dose distribution was in close agreement with the original dose. Nearly the entire patient volume (98.3%-98.7%) satisfied the criteria of chi-evaluation (<2% maximum dose and 2 mm range). The dose to 95% of the volume and the percentage of volume receiving at least 95% of the prescription dose in the planning target volume differed from the original values by less than 2% of the prescription dose (root-mean-square, RMS < 1%). The PRGP scheme did not perform better than the arithmetic mean process with the same number of atlases. Increasing the number of atlases from 6 to 12 often resulted in improvements, but statistical significance was not always found.

CONCLUSIONS

MRI-based treatment planning with pseudo CTs generated through atlas registration is feasible for pediatric brain tumor patients. The doses calculated from pseudo CTs agreed well with those from real CTs, showing dosimetric accuracy within 2% for the PTV when multiple atlases were used. The arithmetic mean process may be a reasonable choice over PRGP for the synthesis scheme considering performance and computational costs.

Feasibility and limitations of bulk density assignment in MRI for head and neck IMRT treatment planning

Head and neck cancers centered at the base of skull are better visualized on MRI than on CT. The purpose of this investigation was to investigate the accuracy of bulk density assignment in head and neck intensity‐modulated radiation therapy (IMRT) treatment plan optimization. Our study investigates dose calculation differences between density‐assigned MRI and CT, and identifies potential limitations related to dental implants and MRI geometrical distortion in the framework of MRI‐only‐based treatment planning. Bulk density assignment was performed and applied onto MRI to generate three MRI image sets with increasing levels of heterogeneity for seven patients: 1) MRIW: all water‐equivalent; 2) MRIW + B: included bone with density of 1.53 g/cm3; and 3) MRIW + B + A: included bone and air. Using identical planning and optimization parameters, MRI‐based IMRT plans were generated and compared to corresponding, forward‐calculated, CT‐based plans on the basis of target coverage, isodose distributions, and dose‐volume histograms (DVHs). Phantom studies were performed to assess the magnitude and spatial dependence of MRI geometrical distortion. MRIW‐based dose calculations overestimated target coverage by 16.1%. Segmentation of bone reduced differences to within 2% of the coverage area on the CT‐based plan. Further segmentation of air improved conformity near air–tissue interfaces. Dental artifacts caused substantial target coverage overestimation even on MRIW + B + A. Geometrical distortion was less than 1 mm in an imaging volume 20 × 20 × 20 cm3 around scanner isocenter, but up to 4 mm at 17 cm lateral to isocenter. Bulk density assignment in the framework of MRI‐only IMRT head and neck treatment planning is a feasible method with certain limitations. Bone and teeth account for the majority of density heterogeneity effects. While soft tissue is well visualized on MRI compared to CT, dental implants may not be visible on MRI and must be identified by other means and assigned appropriate density for accurate dose calculation. Far off‐center geometrical distortion of the body contour near the shoulder region is a potential source of dose calculation inaccuracy.

PACS numbers: 87.61.‐c, 87.55.‐D

A dual model HU conversion from MRI intensity values within and outside of bone segment for MRI-based radiotherapy treatment planning of prostate cancer

PURPOSE

The lack of electron density information in magnetic resonance images (MRI) poses a major challenge for MRI-based radiotherapy treatment planning (RTP). In this study the authors convert MRI intensity values into Hounsfield units (HUs) in the male pelvis and thus enable accurate MRI-based RTP for prostate cancer patients with varying tissue anatomy and body fat contents.

METHODS

T1/T2*-weighted MRI intensity values and standard computed tomography (CT) image HUs in the male pelvis were analyzed using image data of 10 prostate cancer patients. The collected data were utilized to generate a dual model HU conversion technique from MRI intensity values of the single image set separately within and outside of contoured pelvic bones. Within the bone segment local MRI intensity values were converted to HUs by applying a second-order polynomial model. This model was tuned for each patient by two patient-specific adjustments: MR signal normalization to correct shifts in absolute intensity level and application of a cutoff value to accurately represent low density bony tissue HUs. For soft tissues, such as fat and muscle, located outside of the bone contours, a threshold-based segmentation method without requirements for any patient-specific adjustments was introduced to convert MRI intensity values into HUs. The dual model HU conversion technique was implemented by constructing pseudo-CT images for 10 other prostate cancer patients. The feasibility of these images for RTP was evaluated by comparing HUs in the generated pseudo-CT images with those in standard CT images, and by determining deviations in MRI-based dose distributions compared to those in CT images with 7-field intensity modulated radiation therapy (IMRT) with the anisotropic analytical algorithm and 360° volumetric-modulated arc therapy (VMAT) with the Voxel Monte Carlo algorithm.

RESULTS

The average HU differences between the constructed pseudo-CT images and standard CT images of each test patient ranged from -2 to 5 HUs and from 22 to 78 HUs in soft and bony tissues, respectively. The average local absolute value differences were 11 HUs in soft tissues and 99 HUs in bones. The planning target volume doses (volumes 95%, 50%, 5%) in the pseudo-CT images were within 0.8% compared to those in CT images in all of the 20 treatment plans. The average deviation was 0.3%. With all the test patients over 94% (IMRT) and 92% (VMAT) of dose points within body (lower than 10% of maximum dose suppressed) passed the 1 mm and 1% 2D gamma index criterion. The statistical tests (t- and F-tests) showed significantly improved (p ≤ 0.05) HU and dose calculation accuracies with the soft tissue conversion method instead of homogeneous representation of these tissues in MRI-based RTP images.

CONCLUSIONS

This study indicates that it is possible to construct high quality pseudo-CT images by converting the intensity values of a single MRI series into HUs in the male pelvis, and to use these images for accurate MRI-based prostate RTP dose calculations.

Influence of MRI-based bone outline definition errors on external radiotherapy dose calculation accuracy in heterogeneous pseudo-CT images of prostate cancer patients

BACKGROUND

This work evaluates influences of susceptibility-induced bone outline shift and perturbations, and bone segmentation errors on external radiotherapy dose calculation accuracy in magnetic resonance imaging (MRI)-based pseudo-computed tomography (CT) images of the male pelvis.

MATERIAL AND METHODS

T₁/T₂*-weighted fast gradient echo, T₁-weighted spin echo and T₂-weighted fast spin echo images were used in bone detection investigation. Bone edge location and bone diameter in MRI were evaluated by comparing those in the images with actual physical measurements of fresh deer bones positioned in a gelatine phantom. Dose calculation accuracy in pseudo-CT images was investigated for 15 prostate cancer patients. Bone outlines in T₁/T₂*-weighted images were contoured and additional segmentation errors were simulated by expanding and contracting the bone contours with 1 mm spacing. Heterogeneous pseudo-CT images were constructed by adopting a technique transforming the MRI intensity values into Hounsfield units with separate conversion models within and outside of bone segment.

RESULTS

Bone edges and diameter in the phantom were illustrated correctly within a 1 mm-pixel size in MRI. Each 1 mm-sized systematic error in bone segment resulted in roughly 0.4% change to the prostate dose level in the pseudo-CT images. The prostate average (range) dose levels in pseudo-CT images with additional systematic bone segmentation errors of -2 mm, 0 mm and 2 mm were 0.5% (-0.5-1.4%), -0.2% (-1.0-0.7%), and -0.9% (-1.8-0.0%) compared to those in CT images, respectively, in volumetric modulated arc therapy treatment plans calculated by Monte Carlo algorithm.

CONCLUSIONS

Susceptibility-induced bone outline shift and perturbations do not result in substantial uncertainty for MRI-based dose calculation. Dose consistency of 2% can be achieved reliably for the prostate if heterogeneous pseudo-CT images are constructed with ≤± 2 mm systematic error in bone segment.