Dosimetric evaluation of MRI-based treatment planning for prostate cancer

Clinical validation of MR imaging time reduction for substitute/synthetic CT generation for prostate MRI-only treatment planning

Radiotherapy treatment planning based only on magnetic resonance imaging (MRI) has become clinically achievable. Though computed tomography (CT) is the gold standard for radiotherapy imaging, directly providing the electron density values needed for planning calculations, MRI has superior soft tissue visualisation to guide treatment planning decisions and optimisation. MRI-only planning removes the need for the CT scan, but requires generation of a substitute/synthetic/pseudo CT (sCT) for electron density information. Shortening the MRI imaging time would improve patient comfort and reduce the likelihood of motion artefacts. A volunteer study was previously carried out to investigate and optimise faster MRI sequences for a hybrid atlas-voxel conversion to sCT for prostate treatment planning. The aim of this follow-on study was to clinically validate the performance of the new optimised sequence for sCT generation in a treated MRI-only prostate patient cohort. 10 patients undergoing MRI-only treatment were scanned on a Siemens Skyra 3T MRI as part of the MRI-only sub-study of the NINJA clinical trial (ACTRN12618001806257). Two sequences were used, the standard 3D T2-weighted SPACE sequence used for sCT conversion which has been previously validated against CT, and a modified fast SPACE sequence, selected based on the volunteer study. Both were used to generate sCT scans. These were then compared to evaluate the fast sequence conversion for anatomical and dosimetric accuracy against the clinically approved treatment plans. The average Mean Absolute Error (MAE) for the body was 14.98 ± 2.35 HU, and for bone was 40.77 ± 5.51 HU. The external volume contour comparison produced a Dice Similarity Coefficient (DSC) of at least 0.976, and an average of 0.985 ± 0.004, and the bony anatomy contour comparison a DSC of at least 0.907, and an average of 0.950 ± 0.018. The fast SPACE sCT agreed with the gold standard sCT within an isocentre dose of -0.28% ± 0.16% and an average gamma pass rate of 99.66% ± 0.41% for a 1%/1 mm gamma tolerance. In this clinical validation study, the fast sequence, which reduced the required imaging time by approximately a factor of 4, produced an sCT with similar clinical dosimetric results compared to the standard sCT, demonstrating its potential for clinical use for treatment planning.

Dosimetric Accuracy of MR-Guided Online Adaptive Planning for Nasopharyngeal Carcinoma Radiotherapy on 1.5 T MR-Linac

Purpose

The aim of this study is to evaluate the dose accuracy of bulk relative electron density (rED) approach for application in 1.5 T MR-Linac and assess the reliability of this approach in the case of online adaptive MR-guided radiotherapy for nasopharyngeal carcinoma (NPC) patients.

Methods

Ten NPC patients formerly treated on conventional linac were included in this study, with their original planning CT and MRI collected. For each patient, structures such as the targets, organs at risk, bone, and air regions were delineated on the original CT in the Monaco system (v5.40.02). To simulate the online adaptive workflow, firstly all contours were transferred to MRI from the original CT using rigid registration in the Monaco system. Based on the structures, three different types of synthetic CT (sCT) were generated from MRI using the bulk rED assignment approach: the sCT_ICRU uses the rED values recommended by ICRU46, the sCT_tailor uses the patient-specific mean rED values, and the sCT_Homogeneity uses homogeneous water equivalent values. The same treatment plan was calculated on the three sCTs and the original CT. Dose calculation accuracy was investigated in terms of gamma analysis, point dose comparison, and dose volume histogram (DVH) parameters.

Results

Good agreement of dose distribution was observed between sCT_tailor and the original CT, with a gamma passing rate (3%/3 mm) of 97.81% ± 1.06%, higher than that of sCT_ICRU (94.27% ± 1.48%, p = 0.005) and sCT_Homogeneity (96.50% ± 1.02%, p = 0.005). For stricter criteria 1%/1 mm, gamma passing rates for plans on sCT_tailor, sCT_ICRU, and sCT_Homogeneity were 86.79% ± 4.31%, 79.81% ± 3.63%, and 77.56% ± 4.64%, respectively. The mean point dose difference in PTV_nx between sCT_tailor and planning CT was −0.14% ± 1.44%, much lower than that calculated on sCT_ICRU (−8.77% ± 2.33%) and sCT_Homogeneity (1.65% ± 2.57%), all with p < 0.05. The DVH differences for the plan based on sCT_tailor were much smaller than sCT_ICRU and sCT_Homogeneity.

Conclusions

The bulk rED-assigned sCT by adopting the patient-specific rED values can achieve a clinically acceptable level of dose calculation accuracy in the presence of a 1.5 T magnetic field, making it suitable for online adaptive MR-guided radiotherapy for NPC patients.

Safety of stereotactic body radiation therapy for localized prostate cancer without treatment planning MRI

Background

The use of treatment planning prostate MRI for Stereotactic Body Radiation Therapy (SBRT) is largely a standard, yet not all patients can receive MRI for a variety of clinical reasons. Thus, we aim to investigate the safety of patients who received CT alone based SBRT planning for the definitive treatment of localized prostate cancer.

Methods

Our study analyzed 3410 patients with localized prostate cancer who were treated with SBRT at a single academic institution between 2006 and 2020. Acute and late toxicity was evaluated using the Common Terminology Criteria for Adverse Events version 5.0. Expanded Prostate Cancer Index Composite (EPIC) questionnaires evaluated QOL and PSA nadir was evaluated to detect biochemical failures.

Results

A total of 162 patients (4.75%) received CT alone for treatment planning. The CT alone group was older relative to the MRI group (69.9 vs 67.2, p < 0.001) and had higher risk and grade disease (p < 0.001). Additionally, the CT group exhibited a trend in larger CTVs (82.56 cc vs 76.90 cc; p = 0.055), lower total radiation doses (p = 0.048), and more frequent pelvic nodal radiation versus the MRI group (p < 0.001). There were only two reported cases of Grade 3 + toxicity within the CT alone group. Quality of life data within the CT alone group revealed declines in urinary and bowel scores at one month with return to baseline at subsequent follow up. Early biochemical failure data at median time of 2.3 years revealed five failures by Phoenix definition.

Conclusions

While clinical differences existed between the MRI and CT alone group, we observed tolerable toxicity profiles in the CT alone cohort, which was further supported by EPIC questionnaire data. The overall clinical outcomes appear comparable in patients unable to receive MRI for their SBRT treatment plan with early clinical follow up.

Effects of MR imaging time reduction on substitute CT generation for prostate MRI-only treatment planning

The introduction of MRI linear accelerators (MR-linacs) and the increased use of MR imaging in radiotherapy, requires improved approaches to MRI-only radiotherapy. MRI provides excellent soft tissue visualisation but does not provide any electron density information required for radiotherapy dose calculation, instead MRI is registered to CT images to enable dose calculations. MRI-only radiotherapy eliminates registration errors and reduces patient discomfort, workload and cost. Electron density requirements may be addressed in different ways, from manually applying bulk density corrections, to more computationally intensive methods to produce substitute CT datasets (sCT), requiring additional sequences, increasing overall imaging time. Reducing MR imaging time would reduce potential artefacts from intrafraction motion and patient discomfort. The aim of this study was to assess the impact of reducing MR imaging time on a hybrid atlas-voxel sCT conversion for prostate MRI-only treatment planning, considering both anatomical and dosimetric parameters. 10 volunteers were scanned on a Siemens Skyra 3T MRI. Sequences included the 3D T2-weighted (T2-w) SPACE sequence used for sCT conversion as previously validated against CT, along with variations to this sequence in repetition time (TR), turbo factor, and combinations of these to reduce the imaging time. All scans were converted to sCT and were compared to the sCT from the original SPACE sequence, evaluating for anatomical changes and dosimetric differences for a standard prostate VMAT plan. Compared to the previously validated T2-w SPACE sequence, scan times were reduced by up to 80%. The external volume and bony anatomy were compared, with all but one sequence meeting a DICE coefficient of 0.9 or better, with the largest variations occurring at the edges of the external body volume. The generated sCT agreed with the gold standard sCT within an isocentre dose of 1% and a gamma pass rate of 99% for a 1%/1 mm gamma tolerance for all but one sequence. This study demonstrates that the MR imaging sequence time was able to be reduced by approximately 80% with similar dosimetric results.

Development of a hybrid magnetic resonance/computed tomography-compatible phantom for magnetic resonance guided radiotherapy

The purpose of the present study was to develop a hybrid magnetic resonance/computed tomography (MR/CT)-compatible phantom and tissue-equivalent materials for each MR and CT image. Therefore, the essential requirements necessary for the development of a hybrid MR/CT-compatible phantom were determined and the development process is described. A total of 12 different tissue-equivalent materials for each MR and CT image were developed from chemical components. The uniformity of each sample was calculated. The developed phantom was designed to use 14 plugs that contained various tissue-equivalent materials. Measurement using the developed phantom was performed using a 3.0-T scanner with 32 channels and a Somatom Sensation 64. The maximum percentage difference of the signal intensity (SI) value on MR images after adding K2CO3 was 3.31%. Additionally, the uniformity of each tissue was evaluated by calculating the percent image uniformity (%PIU) of the MR image, which was 82.18 ±1.87% with 83% acceptance, and the average circular-shaped regions of interest (ROIs) on CT images for all samples were within ±5 Hounsfield units (HU). Also, dosimetric evaluation was performed. The percentage differences of each tissue-equivalent sample for average dose ranged from -0.76 to 0.21%. A hybrid MR/CT-compatible phantom for MR and CT was investigated as the first trial in this field of radiation oncology and medical physics.

A Multi-center Prospective Study for Implementation of an MRI-Only Prostate Treatment Planning Workflow

Purpose: This project investigates the feasibility of implementation of MRI-only prostate planning in a prospective multi-center study.

Method and Materials: A two-phase implementation model was utilized where centers performed retrospective analysis of MRI-only plans for five patients followed by prospective MRI-only planning for subsequent patients. Feasibility was assessed if at least 23/25 patients recruited to phase 2 received MRI-only treatment workflow. Whole-pelvic MRI scans (T2 weighted, isotropic 1.6 mm voxel 3D sequence) were converted to pseudo-CT using an established atlas-based method. Dose plans were generated using MRI contoured anatomy with pseudo-CT for dose calculation. A conventional CT scan was acquired subsequent to MRI-only plan approval for quality assurance purposes (QA-CT). 3D Gamma evaluation was performed between pseudo-CT calculated plan dose and recalculation on QA-CT. Criteria was 2%, 2 mm criteria with 20% low dose threshold. Gold fiducial marker positions for image guidance were compared between pseudo-CT and QA-CT scan prior to treatment.

Results: All 25 patients recruited to phase 2 were treated using the MRI-only workflow. Isocenter dose differences between pseudo-CT and QA-CT were −0.04 ± 0.93% (mean ± SD). 3D Gamma dose comparison pass-rates were 99.7% ± 0.5% with mean gamma 0.22 ± 0.07. Results were similar for the two centers using two different scanners. All gamma comparisons exceeded the 90% pass-rate tolerance with a minimum gamma pass-rate of 98.0%. In all cases the gold fiducial markers were correctly identified on MRI and the distances of all seeds to centroid were within the tolerance of 1.0 mm of the distances on QA-CT (0.07 ± 0.41 mm), with a root-mean-square difference of 0.42 mm.

Conclusion: The results support the hypothesis that an MRI-only prostate workflow can be implemented safely and accurately with appropriate quality assurance methods.

[Radiotherapy treatment planning of prostate cancer using magnetic resonance imaging]

PURPOSE

Magnetic resonance imaging (MRI) plays an increasing role in radiotherapy dose planning. Indeed, MRI offers superior soft tissue contrast compared to computerized tomography (CT) and therefore could provide a better delineation of target volumes and organs at risk than CT for radiotherapy. Furthermore, an MRI-only radiotherapy workflow would suppress registration errors inherent to the registration of MRI with CT. However, the estimation of the electronic density of tissues using MRI images is still a challenging issue. The purpose of this work was to design and evaluate a pseudo-CT generation method for prostate cancer treatments.

MATERIALS AND METHODS

A pseudo-CT was generated for ten prostate cancer patients using an elastic deformation based method. For each patient, dose delivered to the patient was calculated using both the planning CT and the pseudo-CT. Dose differences between CT and pseudo-CT were investigated.

RESULTS

Mean dose relative difference in the planning target volume is 0.9% on average and ranges from 0.1% to 1.7%. In organs at risks, this value is 1.8%, 0.8%, 0.8% and 1% on average in the rectum, the right and left femoral heads, and the bladder respectively.

CONCLUSION

The dose calculated using the pseudo-CT is very close to the dose calculated using the CT for both organs at risk and PTV. These results confirm that pseudo-CT images generated using the proposed method could be used to calculate radiotherapy treatment doses on MRI images.

Magnetic Resonance Imaging only Workflow for Radiotherapy Simulation and Planning in Prostate Cancer

Magnetic resonance imaging (MRI) is often combined with computed tomography (CT) in prostate radiotherapy to optimise delineation of the target and organs-at-risk (OAR) while maintaining accurate dose calculation. Such a dual-modality workflow requires two separate imaging sessions, and it has some fundamental and logistical drawbacks. Due to the availability of new MRI hardware and software solutions, CT examinations can be omitted for prostate radiotherapy simulations. All information for treatment planning, including electron density maps and bony anatomy, can nowadays be obtained with MRI. Such an MRI-only simulation workflow reduces delineation ambiguities, eases planning logistics, and improves patient comfort; however, careful validation of the complete MRI-only workflow is warranted. The first institutes are now adopting this MRI-only workflow for prostate radiotherapy. In this article, we will review technology and workflow requirements for an MRI-only prostate simulation workflow.

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.

Evaluating the repeatability and set-up sensitivity of a large field of view distortion phantom and software for magnetic resonance-only radiotherapy

BACKGROUND AND PURPOSE

Magnetic Resonance (MR)-only radiotherapy requires geometrically accurate MR images over the full scanner Field of View (FoV). This study aimed to investigate the repeatability of distortion measurements made using a commercial large FoV phantom and analysis software and the sensitivity of these measurements to small set-up errors.

MATERIALS AND METHODS

Geometric distortion was measured using a commercial phantom and software with 2D and 3D acquisition sequences on three different MR scanners. Two sets of repeatability measurements were made: three scans acquired without moving the phantom between scans (single set-up) and five scans acquired with the phantom re-set up in between each scan (repeated set-up). The set-up sensitivity was assessed by scanning the phantom with an intentional 1 mm lateral offset and independently an intentional 1° rotation.

RESULTS

The mean standard deviation of distortion for all phantom markers for the repeated set-up scans was < 0.4 mm for all scanners and sequences. For the 1 mm lateral offset scan 90 % of the markers agreed within two standard deviations of the mean of the repeated set-up scan (median of all scanners and sequences, range 78%-93%). For the 1° rotation scan, 80% of markers agreed within two standard deviations of the mean (range 69%-93%).

CONCLUSIONS

Geometric distortion measurements using a commercial phantom and associated software appear repeatable, although with some sensitivity to set-up errors. This suggests the phantom and software are appropriate for commissioning a MR-only radiotherapy workflow.

MRI-only treatment planning: benefits and challenges

Over the past decade, the application of magnetic resonance imaging (MRI) has increased, and there is growing evidence to suggest that improvements in the accuracy of target delineation in MRI-guided radiation therapy may improve clinical outcomes in a variety of cancer types. However, some considerations should be recognized including patient motion during image acquisition and geometric accuracy of images. Moreover, MR-compatible immobilization devices need to be used when acquiring images in the treatment position while minimizing patient motion during the scan time. Finally, synthetic CT images (i.e. electron density maps) and digitally reconstructed radiograph images should be generated from MRI images for dose calculation and image guidance prior to treatment. A short review of the concepts and techniques that have been developed for implementation of MRI-only workflows in radiation therapy is provided in this document.

Role of Prostate MR Imaging in Radiation Oncology

The use of prostate MR imaging in radiotherapy continues to evolve. This article describes its current application in the selection of treatment regimens, integration in treatment planning or simulation, and assessment of response. An expert consensus statement from the annual MR in RT symposium is presented, as a list of 21 key quality indicators for the practice of MR imaging simulation in prostate cancer. Although imaging requirements generally follow PIRADSv2 guidelines, additional requirements specific to radiotherapy planning are described. MR imaging-only workflows and MR imaging-guided treatment systems are expected to replace conventional computed tomography-based practice, further adding specific requirements for MR imaging in radiotherapy.

Feasibility of MR-only proton dose calculations for prostate cancer radiotherapy using a commercial pseudo-CT generation method

A magnetic resonance (MR)-only radiotherapy workflow can reduce cost, radiation exposure and uncertainties introduced by CT-MRI registration. A crucial prerequisite is generating the so called pseudo-CT (pCT) images for accurate dose calculation and planning. Many pCT generation methods have been proposed in the scope of photon radiotherapy. This work aims at verifying for the first time whether a commercially available photon-oriented pCT generation method can be employed for accurate intensity-modulated proton therapy (IMPT) dose calculation. A retrospective study was conducted on ten prostate cancer patients. For pCT generation from MR images, a commercial solution for creating bulk-assigned pCTs, called MR for Attenuation Correction (MRCAT), was employed. The assigned pseudo-Hounsfield Unit (HU) values were adapted to yield an increased agreement to the reference CT in terms of proton range. Internal air cavities were copied from the CT to minimise inter-scan differences. CT- and MRCAT-based dose calculations for opposing beam IMPT plans were compared by gamma analysis and evaluation of clinically relevant target and organ at risk dose volume histogram (DVH) parameters. The proton range in beam's eye view (BEV) was compared using single field uniform dose (SFUD) plans. On average, a [Formula: see text] mm) gamma pass rate of 98.4% was obtained using a [Formula: see text] dose threshold after adaptation of the pseudo-HU values. Mean differences between CT- and MRCAT-based dose in the DVH parameters were below 1 Gy ([Formula: see text]). The median proton range difference was [Formula: see text] mm, with on average 96% of all BEV dose profiles showing a range agreement better than 3 mm. Results suggest that accurate MR-based proton dose calculation using an automatic commercial bulk-assignment pCT generation method, originally designed for photon radiotherapy, is feasible following adaptation of the assigned pseudo-HU values.

Quantifying the Effect of 3T Magnetic Resonance Imaging Residual System Distortions and Patient-Induced Susceptibility Distortions on Radiation Therapy Treatment Planning for Prostate Cancer

PURPOSE

To investigate the effect of magnetic resonance system- and patient-induced susceptibility distortions from a 3T scanner on dose distributions for prostate cancers.

METHODS AND MATERIALS

Combined displacement fields from the residual system and patient-induced susceptibility distortions were used to distort 17 prostate patient CT images. VMAT dose plans were initially optimized on distorted CT images and the plan parameters transferred to the original patient CT images to calculate a new dose distribution.

RESULTS

Maximum residual mean distortions of 3.19 mm at a radial distance of 25 cm and maximum mean patient-induced susceptibility shifts of 5.8 mm were found using the lowest bandwidth of 122 Hz per pixel. There was a dose difference of <0.5% between distorted and undistorted treatment plans. The 90% confidence intervals of the mean difference between the dCT and CT treatment plans were all within an equivalence interval of (-0.5, 0.5) for all investigated plan quality measures.

CONCLUSIONS

Patient-induced susceptibility distortions at high field strengths in closed bore magnetic resonance scanners are larger than residual system distortions after using vendor-supplied 3-dimensional correction for the delineated regions studied. However, errors in dose due to disturbed patient outline and shifts caused by patient-induced susceptibility effects are below 0.5%.

Image similarity evaluation of the bulk-density-assigned synthetic CT derived from MRI of intracranial regions for radiation treatment

Objective

Various methods for radiation-dose calculation have been investigated over previous decades, focusing on the use of magnetic resonance imaging (MRI) only. The bulk-density-assignment method based on manual segmentation has exhibited promising results compared to dose-calculation with computed tomography (CT). However, this method cannot be easily implemented in clinical practice due to its time-consuming nature. Therefore, we investigated an automatic anatomy segmentation method with the intention of providing the proper methodology to evaluate synthetic CT images for a radiation-dose calculation based on MR images.

Methods

CT images of 20 brain cancer patients were selected, and their MR images including T1-weighted, T2-weighted, and PETRA were retrospectively collected. Eight anatomies of the patients, such as the body, air, eyeball, lens, cavity, ventricle, brainstem, and bone, were segmented for bulk-density-assigned CT image (_BCT) generation. In addition, water-equivalent CT images (_WCT) with only two anatomies—body and air—were generated for a comparison with _BCT. Histogram comparison and gamma analysis were performed by comparison with the original CT images, after the evaluation of automatic segmentation performance with the dice similarity coefficient (DSC), false negative dice (FND) coefficient, and false positive dice (FPD) coefficient.

Results

The highest DSC value was 99.34 for air segmentation, and the lowest DSC value was 73.50 for bone segmentation. For lens segmentation, relatively high FND and FPD values were measured. The cavity and bone were measured as over-segmented anatomies having higher FPD values than FND. The measured histogram comparison results of _BCT were better than those of _WCT in all cases. In gamma analysis, the averaged improvement of _BCT compared to _WCT was measured. All the measured results of _BCT were better than those of _WCT. Therefore, the results of this study show that the introduced methods, such as histogram comparison and gamma analysis, are valid for the evaluation of the synthetic CT generation from MR images.

Conclusions

The image similarity results showed that _BCT has superior results compared to _WCT for all measurements performed in this study. Consequently, more accurate radiation treatment for the intracranial regions can be expected when the proper image similarity evaluation introduced in this study is performed.

Systematic Review of Synthetic Computed Tomography Generation Methodologies for Use in Magnetic Resonance Imaging-Only Radiation Therapy

Magnetic resonance imaging (MRI) offers superior soft-tissue contrast as compared with computed tomography (CT), which is conventionally used for radiation therapy treatment planning (RTP) and patient positioning verification, resulting in improved target definition. The 2 modalities are co-registered for RTP; however, this introduces a systematic error. Implementing an MRI-only radiation therapy workflow would be advantageous because this error would be eliminated, the patient pathway simplified, and patient dose reduced. Unlike CT, in MRI there is no direct relationship between signal intensity and electron density; however, various methodologies for MRI-only RTP have been reported. A systematic review of these methods was undertaken. The PRISMA guidelines were followed. Embase and Medline databases were searched (1996 to March, 2017) for studies that generated synthetic CT scans (sCT)s for MRI-only radiation therapy. Sixty-one articles met the inclusion criteria. This review showed that MRI-only RTP techniques could be grouped into 3 categories: (1) bulk density override; (2) atlas-based; and (3) voxel-based techniques, which all produce an sCT scan from MR images. Bulk density override techniques either used a single homogeneous or multiple tissue override. The former produced large dosimetric errors (>2%) in some cases and the latter frequently required manual bone contouring. Atlas-based techniques used both single and multiple atlases and included methods incorporating pattern recognition techniques. Clinically acceptable sCTs were reported, but atypical anatomy led to erroneous results in some cases. Voxel-based techniques included methods using routine and specialized MRI sequences, namely ultra-short echo time imaging. High-quality sCTs were produced; however, use of multiple sequences led to long scanning times increasing the chances of patient movement. Using nonroutine sequences would currently be problematic in most radiation therapy centers. Atlas-based and voxel-based techniques were found to be the most clinically useful methods, with some studies reporting dosimetric differences of <1% between planning on the sCT and CT and <1-mm deviations when using sCTs for positional verification.

Dosimetric evaluation of synthetic CT relative to bulk density assignment-based magnetic resonance-only approaches for prostate radiotherapy

BACKGROUND

Magnetic resonance imaging (MRI) has been incorporated as an adjunct to CT to take advantage of its excellent soft tissue contrast for contouring. MR-only treatment planning approaches have been developed to avoid errors introduced during the MR-CT registration process. The purpose of this study is to evaluate calculated dose distributions after incorporating a novel synthetic CT (synCT) derived from magnetic resonance simulation images into prostate cancer treatment planning and to compare dose distributions calculated using three previously published MR-only treatment planning methodologies.

METHODS

An IRB-approved retrospective study evaluated 15 prostate cancer patients that underwent IMRT (n = 11) or arc therapy (n = 4) to a total dose of 70.2-79.2 Gy. Original treatment plans were derived from CT simulation images (CT-SIM). T1-weighted, T2-weighted, and balanced turbo field echo images were acquired on a 1.0 T high field open MR simulator with patients immobilized in treatment position. Four MR-derived images were studied: bulk density assignment (10 HU) to water (MRW), bulk density assignments to water and bone with pelvic bone values derived either from literature (491 HU, MRW+B491) or from CT-SIM population average bone values (300 HU, MRW+B300), and synCTs. Plans were recalculated using fixed monitor units, plan dosimetry was evaluated, and local dose differences were characterized using gamma analysis (1 %/1 mm dose difference/distance to agreement).

RESULTS

While synCT provided closest agreement to CT-SIM for D95, D99, and mean dose (<0.7 Gy (1 %)) compared to MRW, MRW+B491, and MRW+B300, pairwise comparisons showed differences were not significant (p < 0.05). Significant improvements were observed for synCT in the bladder, but not for rectum or penile bulb. SynCT gamma analysis pass rates (97.2 %) evaluated at 1 %/1 mm exceeded those from MRW (94.7 %), MRW+B300 (94.0 %), or MRW+B491 (90.4 %). One subject's synCT gamma (1 %/1 mm) results (89.9 %) were lower than MRW (98.7 %) and MRW+B300 (96.7 %) due to increased rectal gas during MR-simulation that did not affect bulk density assignment-based calculations but was reflected in higher rectal doses for synCT.

CONCLUSIONS

SynCT values provided closest dosimetric and gamma analysis agreement to CT-SIM compared to bulk density assignment-based CT surrogates. SynCTs may provide additional clinical value in treatment sites with greater air-to-soft tissue ratio.

Automatic Substitute Computed Tomography Generation and Contouring for Magnetic Resonance Imaging (MRI)-Alone External Beam Radiation Therapy From Standard MRI Sequences

PURPOSE

To validate automatic substitute computed tomography CT (sCT) scans generated from standard T2-weighted (T2w) magnetic resonance (MR) pelvic scans for MR-Sim prostate treatment planning.

PATIENTS AND METHODS

A Siemens Skyra 3T MR imaging (MRI) scanner with laser bridge, flat couch, and pelvic coil mounts was used to scan 39 patients scheduled for external beam radiation therapy for localized prostate cancer. For sCT generation a whole-pelvis MRI scan (1.6 mm 3-dimensional isotropic T2w SPACE [Sampling Perfection with Application optimized Contrasts using different flip angle Evolution] sequence) was acquired. Three additional small field of view scans were acquired: T2w, T2*w, and T1w flip angle 80° for gold fiducials. Patients received a routine planning CT scan. Manual contouring of the prostate, rectum, bladder, and bones was performed independently on the CT and MR scans. Three experienced observers contoured each organ on MRI, allowing interobserver quantification. To generate a training database, each patient CT scan was coregistered to their whole-pelvis T2w using symmetric rigid registration and structure-guided deformable registration. A new multi-atlas local weighted voting method was used to generate automatic contours and sCT results.

RESULTS

The mean error in Hounsfield units between the sCT and corresponding patient CT (within the body contour) was 0.6 ± 14.7 (mean ± 1 SD), with a mean absolute error of 40.5 ± 8.2 Hounsfield units. Automatic contouring results were very close to the expert interobserver level (Dice similarity coefficient): prostate 0.80 ± 0.08, bladder 0.86 ± 0.12, rectum 0.84 ± 0.06, bones 0.91 ± 0.03, and body 1.00 ± 0.003. The change in monitor units between the sCT-based plans relative to the gold standard CT plan for the same dose prescription was found to be 0.3% ± 0.8%. The 3-dimensional γ pass rate was 1.00 ± 0.00 (2 mm/2%).

CONCLUSIONS

The MR-Sim setup and automatic sCT generation methods using standard MR sequences generates realistic contours and electron densities for prostate cancer radiation therapy dose planning and digitally reconstructed radiograph generation.

Feasibility of MRI-based reference images for image-guided radiotherapy of the pelvis with either cone-beam computed tomography or planar localization images

PURPOSE

This study introduces methods to conduct image-guided radiotherapy (IGRT) of the pelvis with either cone-beam computed tomography (CBCT) or planar localization images by relying solely on magnetic resonance imaging (MRI)-based reference images.

MATERIAL AND METHODS

Feasibility of MRI-based reference images for IGRT was evaluated against kV CBCT (50 scans, 5 prostate cancer patients) and kV & MV planar (5 & 5 image pairs and patients) localization images by comparing the achieved patient position corrections to those obtained by standard CT-based reference images. T1/T2*-weighted in-phase MRI, Hounsfield unit conversion-based heterogeneous pseudo-CT, and bulk pseudo-CT images were applied for reference against localization CBCTs, and patient position corrections were obtained by automatic image registration. IGRT with planar localization images was performed manually by 10 observers using reference digitally reconstructed radiographs (DRRs) reconstructed from the pseudo-CTs and standard CTs. Quality of pseudo-DRRs against CT-DRRs was evaluated with image similarity metrics.

RESULTS

The SDs of differences between CBCT-to-MRI and CBCT-to-CT automatic gray-value registrations were ≤1.0 mm & ≤0.8° and ≤2.5 mm & ≤3.6° with 10 cm diameter cubic VOI and prostate-shaped VOI, respectively. The corresponding values for reference heterogeneous pseudo-CT were ≤1.0 mm & ≤0.7° and ≤2.2 mm & ≤3.3°, respectively. Heterogeneous pseudo-CT was the only type of MRI-based reference image working reliably with automatic bone registration (SDs were ≤0.9 mm & ≤0.7°). The differences include possible residual errors from planning CT to MRI registration. The image similarity metrics were significantly (p≤0.01) better in agreement between heterogeneous pseudo-DRRs and CT-DRRs than between bulk pseudo-DRRs and CT-DRRs. The SDs of differences in manual registrations (3D) with planar kV and MV localization images were ≤1.0 mm and ≤1.7 mm, respectively, between heterogeneous pseudo-DRRs and CT-DRRs, and ≤1.4 mm and ≤2.1 mm between bulk pseudo-DRRs and CT-DRRs.

CONCLUSION

This study demonstrated that it is feasible to conduct IGRT of the pelvis with MRI-based reference images.

MRI simulation: end-to-end testing for prostate radiation therapy using geometric pelvic MRI phantoms

To clinically implement MRI simulation or MRI-alone treatment planning requires comprehensive end-to-end testing to ensure an accurate process. The purpose of this study was to design and build a geometric phantom simulating a human male pelvis that is suitable for both CT and MRI scanning and use it to test geometric and dosimetric aspects of MRI simulation including treatment planning and digitally reconstructed radiograph (DRR) generation.A liquid filled pelvic shaped phantom with simulated pelvic organs was scanned in a 3T MRI simulator with dedicated radiotherapy couch-top, laser bridge and pelvic coil mounts. A second phantom with the same external shape but with an internal distortion grid was used to quantify the distortion of the MR image. Both phantoms were also CT scanned as the gold-standard for both geometry and dosimetry. Deformable image registration was used to quantify the MR distortion. Dose comparison was made using a seven-field IMRT plan developed on the CT scan with the fluences copied to the MR image and recalculated using bulk electron densities. Without correction the maximum distortion of the MR compared with the CT scan was 7.5 mm across the pelvis, while this was reduced to 2.6 and 1.7 mm by the vendor's 2D and 3D correction algorithms, respectively. Within the locations of the internal organs of interest, the distortion was <1.5 and <1 mm with 2D and 3D correction algorithms, respectively. The dose at the prostate isocentre calculated on CT and MRI images differed by 0.01% (1.1 cGy). Positioning shifts were within 1 mm when setup was performed using MRI generated DRRs compared to setup using CT DRRs.The MRI pelvic phantom allows end-to-end testing of the MRI simulation workflow with comparison to the gold-standard CT based process. MRI simulation was found to be geometrically accurate with organ dimensions, dose distributions and DRR based setup within acceptable limits compared to CT.

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.

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.

Magnetic resonance imaging in rectal cancer: A surgeon’s perspective

Magnetic resonance imaging (MRI) in rectal cancer was first investigated in 1999 and has become almost mandatory in planning rectal cancer treatment. MRI has a high accuracy in predicting circumferential resection margin involvement and is used to plan neoadjuvant therapy. The accuracy of MRI in assessing mesorectal lymph nodes remains moderate, as there are no reliable criteria to assess nodal involvement. MRI seems to be good in assessing peritoneal involvement in upper rectal cancer; this however has been assessed in only a few studies and needs further research. For low rectal cancers, mesorectum is thin at the level of levator ani especially in relation to prostate; so predicting circumferential resection margin involvement is not easy. However high spatial resolution coronal imaging shows levator muscles, sphincter complex and intersphincteric plane accurately. This is used to stage low rectal tumors and plan plane of surgery (standard surgery, intersphincteric resection, Extralevator abdominoperineal resection). While most centres perform MRI post chemoradiotherapy, its role in accurate staging post neoadjuvant therapy remains debatable. THe role of Diffusion weighted MRI post neoadjuvant therapy is being evaluated in research settings.

Treatment planning of intracranial targets on MRI derived substitute CT data

BACKGROUND

The use of magnetic resonance imaging (MRI) as a complement to computed tomography (CT) in the target definition procedure for radiotherapy is increasing. To eliminate systematic uncertainties due to image registration, a workflow based entirely on MRI may be preferable. In the present pilot study, we investigate dose calculation accuracy for automatically generated substitute CT (s-CT) images of the head based on MRI. We also produce digitally reconstructed radiographs (DRRs) from s-CT data to evaluate the feasibility of patient positioning based on MR images.

METHODS AND MATERIALS

Five patients were included in the study. The dose calculation was performed on CT, s-CT, s-CT data without inhomogeneity correction and bulk density assigned MRI images. Evaluation of the results was performed using point dose and dose volume histogram (DVH) comparisons, and gamma index evaluation.

RESULTS

The results demonstrate that the s-CT images improve the dose calculation accuracy compared to the method of non-inhomogeneity corrected dose calculations (mean improvement 2.0% points) and that it performs almost identically to the method of bulk density assignment. The s-CT based DRRs appear to be adequate for patient positioning of intra-cranial targets, although further investigation is needed on this subject.

CONCLUSION

The s-CT method is very fast and yields data that can be used for treatment planning without sacrificing accuracy.

Voxel-wise uncertainty in CT substitute derived from MRI

PURPOSE

In an earlier work, we demonstrated that substitutes for CT images can be derived from MR images using ultrashort echo time (UTE) sequences, conventional T2 weighted sequences, and Gaussian mixture regression (GMR). In this study, we extend this work by analyzing the uncertainties associated with the GMR model and the information contributions from the individual imaging sequences.

METHODS

An analytical expression for the voxel-wise conditional expected absolute deviation (EAD) in substitute CT (s-CT) images was derived. The expression depends only on MR images and can thus be calculated along with each s-CT image. The uncertainty measure was evaluated by comparing the EAD to the true mean absolute prediction deviation (MAPD) between the s-CT and CT images for 14 patients. Further, the influence of the different MR images included in the GMR model on the generated s-CTs was investigated by removing one or more images and evaluating the MAPD for a spectrum of predicted radiological densities.

RESULTS

The largest EAD was predicted at air-soft tissue and bone-soft tissue interfaces. The EAD agreed with the MAPD in both these regions and in regions with lower EADs, such as the brain. Two of the MR images included in the GMR model were found to be mutually redundant for the purpose of s-CT generation.

CONCLUSIONS

The presented uncertainty estimation method accurately predicts the voxel-wise MAPD in s-CT images. Also, the non-UTE sequence previously used in the model was found to be redundant.

Imaging for target volume delineation in rectal cancer radiotherapy--a systematic review

The global move towards more conformal radiotherapy for rectal cancer requires better imaging modalities that both visualise the disease accurately and are reproducible; to reduce interobserver variation. This review explores the advances in imaging modalities used in target volume delineation, with a view to make recommendations for current clinical practice and to propose future directions for research. A systematic review was conducted using MEDLINE and EMBASE. Articles considered relevant by the authors were included. Planning with orthogonal films is being replaced by computed tomography (CT) simulation. This is now considered the 'gold standard' and allows conformal three-dimensional planning. Magnetic resonance imaging (MRI) has been shown to overcome some of the limitations of CT and can be used either as a diagnostic image to visually aid planning, or as a 'planning' MRI carried out in the treatment position and co-registered with the planning CT. The latter approach has been shown to change the treated volumes compared with CT and in prostate cancer patients has been shown to reduce interobserver variation. There are remaining issues with four-dimensional motion that are yet to be fully appreciated or overcome. 2-[18F] fluoro-2-deoxy-d-glucose positron emission tomography/CT co-registered with planning CT results in smaller volumes than CT alone and also reduces interobserver variation, but requires further validation before routine implementation. Experimental work utilising novel positron emission tomography tracers and diffusion-weighted MRI shows promise and requires further evaluation. Rigorous quality assurance is important with processing of newer imaging modalities. Further work needs to be conducted into both interobserver variation and the formal evaluation of the clinical benefits of newer imaging modalities. Developments in image-guided radiotherapy are also required to ensure that improvements in target definition at the planning stage are reproducible throughout treatment.

Registration accuracy for MR images of the prostate using a subvolume based registration protocol

Background

In recent years, there has been a considerable research effort concerning the integration of magnetic resonance imaging (MRI) into the external radiotherapy workflow motivated by the superior soft tissue contrast as compared to computed tomography. Image registration is a necessary step in many applications, e.g. in patient positioning and therapy response assessment with repeated imaging. In this study, we investigate the dependence between the registration accuracy and the size of the registration volume for a subvolume based rigid registration protocol for MR images of the prostate.

Methods

Ten patients were imaged four times each over the course of radiotherapy treatment using a T2 weighted sequence. The images were registered to each other using a mean square distance metric and a step gradient optimizer for registration volumes of different sizes. The precision of the registrations was evaluated using the center of mass distance between the manually defined prostates in the registered images. The optimal size of the registration volume was determined by minimizing the standard deviation of these distances.

Results

We found that prostate position was most uncertain in the anterior-posterior (AP) direction using traditional full volume registration. The improvement in standard deviation of the mean center of mass distance between the prostate volumes using a registration volume optimized to the prostate was 3.9 mm (p < 0.001) in the AP direction. The optimum registration volume size was 0 mm margin added to the prostate gland as outlined in the first image series.

Conclusions

Repeated MR imaging of the prostate for therapy set-up or therapy assessment will both require high precision tissue registration. With a subvolume based registration the prostate registration uncertainty can be reduced down to the order of 1 mm (1 SD) compared to several millimeters for registration based on the whole pelvis.

MRI-guided prostate radiation therapy planning: Investigation of dosimetric accuracy of MRI-based dose planning

BACKGROUND AND PURPOSE

Dose planning requires a CT scan which provides the electron density distribution for dose calculation. MR provides superior soft tissue contrast compared to CT and the use of MR-alone for prostate planning would provide further benefits such as lower cost to the patient. This study compares the accuracy of MR-alone based dose calculations with bulk electron density assignment to CT-based dose calculations for prostate radiotherapy.

MATERIALS AND METHODS

CT and whole pelvis MR images were contoured for 39 prostate patients. Plans with uniform density and plans with bulk density values assigned to bone and tissue were compared to the patient's gold standard full density CT plan. The optimal bulk density for bone was calculated using effective depth measurements. The plans were evaluated using ICRU point doses, dose volume histograms, and Chi comparisons. Differences in spatial uniformity were investigated for the CT and MR scans.

RESULTS

The calculated dose for CT bulk bone and tissue density plans was 0.1±0.6% (mean±1 SD) higher than the corresponding full density CT plan. MR bulk bone and tissue density plans were 1.3±0.8% lower than the full density CT plan. CT uniform density plans and MR uniform density plans were 1.4±0.9% and 2.6±0.9% lower, respectively. Paired t-tests performed on specific points on the DVH graphs showed that points on DVHs for all bulk electron density plans were equivalent with two exceptions. There was no significant difference between doses calculated on Pinnacle and Eclipse. The dose distributions of six patients produced Chi values outside the acceptable range of values when MR-based plans were compared to the full density plan.

CONCLUSIONS

MR-alone bulk density planning is feasible provided bone is assigned a density, however, manual segmentation of bone on MR images will have to be replaced with automatic methods. The major dose differences for MR bulk density plans are due to differences in patient external contours introduced by the MR couch-top and pelvic coil.

Treatment planning using MRI data: an analysis of the dose calculation accuracy for different treatment regions

BACKGROUND

Because of superior soft tissue contrast, the use of magnetic resonance imaging (MRI) as a complement to computed tomography (CT) in the target definition procedure for radiotherapy is increasing. To keep the workflow simple and cost effective and to reduce patient dose, it is natural to strive for a treatment planning procedure based entirely on MRI. In the present study, we investigate the dose calculation accuracy for different treatment regions when using bulk density assignments on MRI data and compare it to treatment planning that uses CT data.

METHODS

MR and CT data were collected retrospectively for 40 patients with prostate, lung, head and neck, or brain cancers. Comparisons were made between calculations on CT data with and without inhomogeneity corrections and on MRI or CT data with bulk density assignments. The bulk densities were assigned using manual segmentation of tissue, bone, lung, and air cavities.

RESULTS

The deviations between calculations on CT data with inhomogeneity correction and on bulk density assigned MR data were small. The maximum difference in the number of monitor units required to reach the prescribed dose was 1.6%. This result also includes effects of possible geometrical distortions.

CONCLUSIONS

The dose calculation accuracy at the investigated treatment sites is not significantly compromised when using MRI data when adequate bulk density assignments are made. With respect to treatment planning, MRI can replace CT in all steps of the treatment workflow, reducing the radiation exposure to the patient, removing any systematic registration errors that may occur when combining MR and CT, and decreasing time and cost for the extra CT investigation.

Systematisation of spatial uncertainties for comparison between a MR and a CT-based radiotherapy workflow for prostate treatments

Background

In the present work we compared the spatial uncertainties associated with a MR-based workflow for external radiotherapy of prostate cancer to a standard CT-based workflow. The MR-based workflow relies on target definition and patient positioning based on MR imaging. A solution for patient transport between the MR scanner and the treatment units has been developed. For the CT-based workflow, the target is defined on a MR series but then transferred to a CT study through image registration before treatment planning, and a patient positioning using portal imaging and fiducial markers.

Methods

An "open bore" 1.5T MRI scanner, Siemens Espree, has been installed in the radiotherapy department in near proximity to a treatment unit to enable patient transport between the two installations, and hence use the MRI for patient positioning. The spatial uncertainty caused by the transport was added to the uncertainty originating from the target definition process, estimated through a review of the scientific literature. The uncertainty in the CT-based workflow was estimated through a literature review.

Results

The systematic uncertainties, affecting all treatment fractions, are reduced from 3-4 mm (1Sd) with a CT based workflow to 2-3 mm with a MR based workflow. The main contributing factor to this improvement is the exclusion of registration between MR and CT in the planning phase of the treatment.

Conclusion

Treatment planning directly on MR images reduce the spatial uncertainty for prostate treatments.

A simulation of MRI based dose calculations on the basis of radiotherapy planning CT images

BACKGROUND

The advantage of MRI-based radiotherapy planning is the superior soft tissue differentiation. However, for accurate patient dose calculations, a conversion of the MR images into Hounsfield CT maps is necessary. The aim of the present study was to investigate the dose accuracy that can be achieved with segmented MR-images derived from the planning CT images by assigning fixed densities to different classes of tissues.

METHODS

Treatment plans for ten prostate cancer patients were selected. A collapsed cone algorithm was used to calculate patient dose distributions. The dose calculations were based on four different image sets: (1) the original CT-series (DD(DP)), (2) a simulated MR series with all tissue set to a homogenous water equivalent material of density 1.02 g/cm(3) (DD(W)), (3) a simulated MR series with soft tissue set to a water equivalent material with density 1.02 g/cm(3) and the bone set to a density of 1.3 g/cm(3) (DD(W+B1.3)), and (4) a simulated MR series identical to (3) but with a bone density equal to 2.1 g/cm(3) (DD(W+B2.1)). The dose distributions were compared by analysing dose difference histograms as well as through a visual display of spatial dose deviations.

RESULTS

The population based minimum, mean and maximum dose difference between the DD(DP) and DD(W) in the target volume was -2.8, -1.0 and 0.6%, respectively. Corresponding differences between DD(DP) and DD(W+B1.3) were -1.6, 0.2 and 1.5%, respectively, and between DD(DP) and DD(W+B2.1) -4.3, 4.2 and 9.7%, respectively. For the rectum, the differences between CT(DP) and the other image sets were in the range of -19.5 to 8.8%. For the bladder, the differences were in the range of -9.6 to 7.0%.

CONCLUSIONS

A systematic study using segmented MR images was undertaken. To achieve an acceptable accuracy in the CTV dose, the MR images should be segmented into bone and water equivalent tissue. Still, significant dose deviation for the organs at risk may be present. As tissue segmentation in real MR images is introduced, segmentation errors and errors that stem from geometrical non-linearities may further reduce the accuracy.

Accurate heterogeneous dose calculation for lung cancer patients without high-resolution CT densities

The aim of this study was to investigate the relative accuracy of megavoltage photon-beam dose calculations employing either 5 bulk densities or independent voxel densities determined by calibration of the CT Houndsfield number. Full-resolution CT and bulk density treatment plans were generated for 70 lung or esophageal cancer tumors (66 cases) using a commercial treatment planning system with an adaptive convolution dose calculation algorithm (Pinnacle3, Philips Medicals Systems). Bulk densities were applied to segmented regions. Individual and population average densities were compared to the full-resolution plan for each case. Monitor units were kept constant and no normalizations were employed. Dose volume histograms (DVH) and dose difference distributions were examined for all cases. The average densities of the segmented air, lung, fat, soft tissue, and bone for the entire set were found to be 0.14, 0.26, 0.89, 1.02, and 1.12 g/cc, respectively. In all cases, the normal tissue DVH agreed to better than 2% in dose. In 62 of 70 target DVHs, agreement to better than 3% in dose was observed. Six cases demonstrated emphysema, one with bullous formations and one with a hiatus hernia having a large volume of gas. They required the additional assignment of density to the emphysemic lung and inflammatory changes to the lung, the regions of collapsed lung, the bullous formations, and the hernia gas. Bulk tissue density dose calculation provides an accurate method of heterogeneous dose calculation. However, patients with advanced emphysema may require high-resolution CT studies for accurate treatment planning.

MRI-based treatment planning with electron density information mapped from CT images: a preliminary study

Nowadays magnetic resonance imaging (MRI) has been profoundly used in radiotherapy (RT) planning to aid the contouring of targets and critical organs in brain and intracranial cases, which is attributable to its excellent soft tissue contrast and multi-planar imaging capability. However, the lack of electron density information in MRI, together with the image distortion issues, precludes its use as the sole image set for RT planning and dose calculation. The purpose of this preliminary study is to probe the feasibility and evaluate an MRI-based radiation dose calculation process by providing MR images the necessary electron density (ED) information from a patient's readily available diagnostic/staging computed tomography (CT) images using an image registration model. To evaluate the dosimetric accuracy of the proposed approach, three brain and three intracranial cases were selected retrospectively for this study. For each patient, the MR images were registered to the CT images, and the ED information was then mapped onto the MR images by in-house developed software generating a modified set of MR images. Another set of MR images with voxel values assigned with the density of water was also generated. The original intensity modulated radiation treatment (IMRT) plan was then applied to the two sets of MR images and the doses were calculated. The dose distributions from the MRI-based calculations were compared to that of the original CT-based calculation. In all cases, the MRI-based calculations with mapped ED yielded dose values very close (within 2%) to that of the CT-based calculations. The MRI-based calculations with voxel values assigned with water density indicated a dosimetric error of 3-5%, depending on the treatment site. The present approach offers a means of utilizing MR images for accurate dose calculation and affords a potential to eliminate the redundant simulation CT by planning a patient's treatment with only simulation MRI and any available diagnostic/staging CT data.

Concurrent multimodality image segmentation by active contours for radiotherapy treatment planning

Multimodality imaging information is regularly used now in radiotherapy treatment planning for cancer patients. The authors are investigating methods to take advantage of all the imaging information available for joint target registration and segmentation, including multimodality images or multiple image sets from the same modality. In particular, the authors have developed variational methods based on multivalued level set deformable models for simultaneous 2D or 3D segmentation of multimodality images consisting of combinations of coregistered PET, CT, or MR data sets. The combined information is integrated to define the overall biophysical structure volume. The authors demonstrate the methods on three patient data sets, including a nonsmall cell lung cancer case with PET/CT, a cervix cancer case with PET/CT, and a prostate patient case with CT and MRI. CT, PET, and MR phantom data were also used for quantitative validation of the proposed multimodality segmentation approach. The corresponding Dice similarity coefficient (DSC) was 0.90 +/- 0.02 (p < 0.0001) with an estimated target volume error of 1.28 +/- 1.23% volume. Preliminary results indicate that concurrent multimodality segmentation methods can provide a feasible and accurate framework for combining imaging data from different modalities and are potentially useful tools for the delineation of biophysical structure volumes in radiotherapy treatment planning.

[Dosimetric evaluation of an automatic segmentation tool of pelvic structures from MRI images for prostate cancer radiotherapy]

PURPOSE

An automatic segmentation tool of pelvic structures from MRI images for prostate cancer radiotherapy was developed and dosimetric evaluation of differences of delineation (automatic versus human) is presented here.

MATERIALS AND METHODS

CTV, rectum and bladder were defined automatically and by a physician in 20 patients. Treatment plans based on "automatic" volumes were transferred on "manual" volumes and reciprocally. Dosimetric characteristics of PTV (V(95), minimal, maximal and mean doses), rectum (V(50), V(70), maximal and mean doses) and bladder (V(70), maximal and mean doses) were compared.

RESULTS

Automatic delineation of CTV did not significantly influence dosimetric characteristics of "manual" PTV. Rectal V(50) and V(70) were not significantly different; mean rectal dose is slightly superior (43.2 versus 44.4Gy, p=0.02, Student test). Bladder V(70) was significantly superior too (19.3 versus 21.6, p=0.004). Organ-at-risk (OAR) automatic delineation had little influence on their dosimetric characteristics; rectal V(70) was slightly underestimated (20 versus 18.5Gy, p=0.001).

CONCLUSION

CTV and OAR automatic delineation had little influence on dosimetric characteristics. Software developments are ongoing to enable routine use and interobserver evaluation is needed.