Jason Dowling

An Atlas-Based Electron Density Mapping Method for Magnetic Resonance Imaging (MRI)-Alone Treatment Planning and Adaptive MRI-Based Prostate Radiation Therapy

PURPOSE Prostate radiation therapy dose planning directly on magnetic resonance imaging (MRI) scans would reduce costs and uncertainties due to multimodality image registration. Adaptive planning using a combined MRI-linear accelerator approach will also require dose calculations to be performed using MRI data. The aim of this work was to develop an atlas-based method to map realistic electron densities to MRI scans for dose calculations and digitally reconstructed radiograph (DRR) generation. METHODS AND MATERIALS Whole-pelvis MRI and CT scan data were collected from 39 prostate patients. Scans from 2 patients showed significantly different anatomy from that of the remaining patient population, and these patients were excluded. A whole-pelvis MRI atlas was generated based on the manually delineated MRI scans. In addition, a conjugate electron-density atlas was generated from the coregistered computed tomography (CT)-MRI scans. Pseudo-CT scans for each patient were automatically generated by global and nonrigid registration of the MRI atlas to the patient MRI scan, followed by application of the same transformations to the electron-density atlas. Comparisons were made between organ segmentations by using the Dice similarity coefficient (DSC) and point dose calculations for 26 patients on planning CT and pseudo-CT scans. RESULTS The agreement between pseudo-CT and planning CT was quantified by differences in the point dose at isocenter and distance to agreement in corresponding voxels. Dose differences were found to be less than 2%. Chi-squared values indicated that the planning CT and pseudo-CT dose distributions were equivalent. No significant differences (p > 0.9) were found between CT and pseudo-CT Hounsfield units for organs of interest. Mean ± standard deviation DSC scores for the atlas-based segmentation of the pelvic bones were 0.79 ± 0.12, 0.70 ± 0.14 for the prostate, 0.64 ± 0.16 for the bladder, and 0.63 ± 0.16 for the rectum. CONCLUSIONS The electron-density atlas method provides the ability to automatically define organs and map realistic electron densities to MRI scans for radiotherapy dose planning and DRR generation. This method provides the necessary tools for MRI-alone treatment planning and adaptive MRI-based prostate radiation therapy.

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 patients 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.

Evaluation of prostate segmentation algorithms for MRI: the PROMISE12 challenge

Prostate MRI image segmentation has been an area of intense research due to the increased use of MRI as a modality for the clinical workup of prostate cancer. Segmentation is useful for various tasks, e.g. to accurately localize prostate boundaries for radiotherapy or to initialize multi-modal registration algorithms. In the past, it has been difficult for research groups to evaluate prostate segmentation algorithms on multi-center, multi-vendor and multi-protocol data. Especially because we are dealing with MR images, image appearance, resolution and the presence of artifacts are affected by differences in scanners and/or protocols, which in turn can have a large influence on algorithm accuracy. The Prostate MR Image Segmentation (PROMISE12) challenge was setup to allow a fair and meaningful comparison of segmentation methods on the basis of performance and robustness. In this work we will discuss the initial results of the online PROMISE12 challenge, and the results obtained in the live challenge workshop hosted by the MICCAI2012 conference. In the challenge, 100 prostate MR cases from 4 different centers were included, with differences in scanner manufacturer, field strength and protocol. A total of 11 teams from academic research groups and industry participated. Algorithms showed a wide variety in methods and implementation, including active appearance models, atlas registration and level sets. Evaluation was performed using boundary and volume based metrics which were combined into a single score relating the metrics to human expert performance. The winners of the challenge where the algorithms by teams Imorphics and ScrAutoProstate, with scores of 85.72 and 84.29 overall. Both algorithms where significantly better than all other algorithms in the challenge (p<0.05) and had an efficient implementation with a run time of 8min and 3s per case respectively. Overall, active appearance model based approaches seemed to outperform other approaches like multi-atlas registration, both on accuracy and computation time. Although average algorithm performance was good to excellent and the Imorphics algorithm outperformed the second observer on average, we showed that algorithm combination might lead to further improvement, indicating that optimal performance for prostate segmentation is not yet obtained. All results are available online at http://promise12.grand-challenge.org/.

Artificial human vision.

Can vision be restored to the blind? As early as 1929 it was discovered that stimulating the visual cortex of an individual led to the perception of spots of light, known as phosphenes [1]. The aim of artificial human vision systems is to attempt to utilize the perception of phosphenes to provide a useful substitute for normal vision. Currently, four locations for electrical stimulation are being investigated; behind the retina (subretinal), in front of the retina (epiretinal), the optic nerve and the visual cortex (using intra- and surface electrodes). This review discusses artificial human vision technology and requirements, and reviews the current development projects.

Current and future prospects for optoelectronic retinal prostheses

There has been accelerating progress in the development of retinal prosthesis systems designed to partially restore vision to people with retinitis pigmentosa or macular degeneration. Current retinal prostheses can be divided into two types: those that receive power and information from passive or active photodiodes (optoelectronic systems) and those based on multielectrode arrays powered by cables or transcutaneous telemetry systems. Currently, four research groups have ongoing chronic implantation clinical studies, and two of these groups plan to have commercial retinal prosthesis systems available within the next 2 years. This paper reviews the development and current status of the most significant international retinal prosthesis research groups and discusses future prospects and issues associated with their retinal prosthesis systems.

Patient Specific Prostate Segmentation in 3-D Magnetic Resonance Images

Accurate localization of the prostate and its surrounding tissue is essential in the treatment of prostate cancer. This paper presents a novel approach to fully automatically segment the prostate, including its seminal vesicles, within a few minutes of a magnetic resonance (MR) scan acquired without an endorectal coil. Such MR images are important in external beam radiation therapy, where using an endorectal coil is highly undesirable. The segmentation is obtained using a deformable model that is trained on-the-fly so that it is specific to the patients scan. This case specific deformable model consists of a patient specific initialized triangulated surface and image feature model that are trained during its initialization. The image feature model is used to deform the initialized surface by template matching image features (via normalized cross-correlation) to the features of the scan. The resulting deformations are regularized over the surface via well established simple surface smoothing algorithms, which is then made anatomically valid via an optimized shape model. Mean and median Dices similarity coefficients (DSCs) of 0.85 and 0.87 were achieved when segmenting 3T MR clinical scans of 50 patients. The median DSC result was equal to the inter-rater DSC and had a mean absolute surface error of 1.85 mm. The approach is showed to perform well near the apex and seminal vesicles of the prostate.

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.

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.

Robust inverse-consistent affine CT-MR registration in MRI-assisted and MRI-alone prostate radiation therapy

BACKGROUND CT-MR registration is a critical component of many radiation oncology protocols. In prostate external beam radiation therapy, it allows the propagation of MR-derived contours to reference CT images at the planning stage, and it enables dose mapping during dosimetry studies. The use of carefully registered CT-MR atlases allows the estimation of patient specific electron density maps from MRI scans, enabling MRI-alone radiation therapy planning and treatment adaptation. In all cases, the precision and accuracy achieved by registration influences the quality of the entire process. PROBLEM Most current registration algorithms do not robustly generalize and lack inverse-consistency, increasing the risk of human error and acting as a source of bias in studies where information is propagated in a particular direction, e.g. CT to MR or vice versa. In MRI-based treatment planning where both CT and MR scans serve as spatial references, inverse-consistency is critical, if under-acknowledged. PURPOSE A robust, inverse-consistent, rigid/affine registration algorithm that is well suited to CT-MR alignment in prostate radiation therapy is presented. METHOD The presented method is based on a robust block-matching optimization process that utilises a half-way space definition to maintain inverse-consistency. Inverse-consistency substantially reduces the influence of the order of input images, simplifying analysis, and increasing robustness. An open source implementation is available online at http://aehrc.github.io/Mirorr/. RESULTS Experimental results on a challenging 35 CT-MR pelvis dataset demonstrate that the proposed method is more accurate than other popular registration packages and is at least as accurate as the state of the art, while being more robust and having an order of magnitude higher inverse-consistency than competing approaches. CONCLUSION The presented results demonstrate that the proposed registration algorithm is readily applicable to prostate radiation therapy planning.

Fast automatic multi-atlas segmentation of the prostate from 3D MR images

A fast fully automatic method of segmenting the prostate from 3D MR scans is presented, incorporating dynamic multi-atlas label fusion. The diffeomorphic demons method is used for non-rigid registration and a comparison of alternate metrics for atlas selection is presented. A comparison of results from an average shape atlas and the multi-atlas approach is provided. Using the same clinical dataset and manual contours from 50 clinical scans as Klein et al. (2008) a median Dice similarity coefficient of 0.86 was achieved with an average surface error of 2.00mm using the multi-atlas segmentation method.