Automatic Segmentation of the Eye in 3D Magnetic Resonance Imaging: A novel Statistical Shape Model for treatment planning of Retinoblastoma

Purpose: Proper delineation of ocular anatomy in 3D imaging is a big challenge, particularly when developing treatment plans for ocular diseases. Magnetic Resonance Imaging (MRI) is nowadays utilized in clinical practice for the diagnosis confirmation and treatment planning of retinoblastoma in infants, where it serves as a source of information, complementary to the Fundus or Ultrasound imaging. Here we present a framework to fully automatically segment the eye anatomy in the MRI based on 3D Active Shape Models (ASM), we validate the results and present a proof of concept to automatically segment pathological eyes. Material and Methods: Manual and automatic segmentation were performed on 24 images of healthy children eyes (3.29±2.15 years). Imaging was performed using a 3T MRI scanner. The ASM comprises the lens, the vitreous humor, the sclera and the cornea. The model was fitted by first automatically detecting the position of the eye center, the lens and the optic nerve, then aligning the model and fitting it to the patient. We validated our segmentation method using a leave-one-out cross validation. The segmentation results were evaluated by measuring the overlap using the Dice Similarity Coefficient (DSC) and the mean distance error. Results: We obtained a DSC of 94.90±2.12% for the sclera and the cornea, 94.72±1.89% for the vitreous humor and 85.16±4.91% for the lens. The mean distance error was 0.26±0.09mm. The entire process took 14s on average per eye. Conclusion: We provide a reliable and accurate tool that enables clinicians to automatically segment the sclera, the cornea, the vitreous humor and the lens using MRI. We additionally present a proof of concept for fully automatically segmenting pathological eyes. This tool reduces the time needed for eye shape delineation and thus can help clinicians when planning eye treatment and confirming the extent of the tumor.

Comparison of gadoteric acid and gadobutrol for detection as well as morphologic and dynamic characterization of lesions on breast dynamic contrast-enhanced magnetic resonance imaging

OBJECTIVE In contrast to conventional breast imaging techniques, one major diagnostic benefit of breast magnetic resonance imaging (MRI) is the simultaneous acquisition of morphologic and dynamic enhancement characteristics, which are based on angiogenesis and therefore provide insights into tumor pathophysiology. The aim of this investigation was to intraindividually compare 2 macrocyclic MRI contrast agents, with low risk for nephrogenic systemic fibrosis, in the morphologic and dynamic characterization of histologically verified mass breast lesions, analyzed by blinded human evaluation and a fully automatic computer-assisted diagnosis (CAD) technique. MATERIALS AND METHODS Institutional review board approval and patient informed consent were obtained. In this prospective, single-center study, 45 women with 51 histopathologically verified (41 malignant, 10 benign) mass lesions underwent 2 identical examinations at 1.5 T (mean time interval, 2.1 days) with 0.1-mmol kg doses of gadoteric acid and gadobutrol. All magnetic resonance images were visually evaluated by 2 experienced, blinded breast radiologists in consensus and by an automatic CAD system, whereas the morphologic and dynamic characterization as well as the final human classification of lesions were performed based on the categories of the Breast imaging reporting and data system MRI atlas. Lesions were also classified by defining their probability of malignancy (morpho-dynamic index; 0%-100%) by the CAD system. Imaging results were correlated with histopathology as gold standard. RESULTS The CAD system coded 49 of 51 lesions with gadoteric acid and gadobutrol (detection rate, 96.1%); initial signal increase was significantly higher for gadobutrol than for gadoteric acid for all and the malignant coded lesions (P < 0.05). Gadoteric acid resulted in more postinitial washout curves and fewer continuous increases of all and the malignant lesions compared with gadobutrol (CAD hot spot regions, P < 0.05). Morphologically, the margins of the malignancies were different between the 2 agents, whereas gadobutrol demonstrated more spiculated and fewer smooth margins (P < 0.05). Lesion classifications by the human observers and by the morpho-dynamic index compared with the histopathologic results did not significantly differ between gadoteric acid and gadobutrol. CONCLUSIONS Macrocyclic contrast media can be reliably used for breast dynamic contrast-enhanced MRI. However, gadoteric acid and gadobutrol differed in some dynamic and morphologic characterization of histologically verified breast lesions in an intraindividual, comparison. Besides the standardization of technical parameters and imaging evaluation of breast MRI, the standardization of the applied contrast medium seems to be important to receive best comparable MRI interpretation.

Guidance of treatment decisions in risk-adapted primary radiotherapy for prostate cancer using multiparametric magnetic resonance imaging: a single center experience.

BACKGROUND Magnetic resonance imaging (MRI) of the prostate is considered to be the most precise noninvasive staging modality for localized prostate cancer. Multiparametric MRI (mpMRI) dynamic sequences have recently been shown to further increase the accuracy of staging relative to morphological imaging alone. Correct radiological staging, particularly the detection of extraprostatic disease extension, is of paramount importance for target volume definition and dose prescription in highly-conformal curative radiotherapy (RT); in addition, it may affect the risk-adapted duration of additional antihormonal therapy. The purpose of our study was to analyze the impact of mpMRI-based tumor staging in patients undergoing primary RT for prostate cancer. METHODS A total of 122 patients admitted for primary RT for prostate cancer were retrospectively analyzed regarding initial clinical and computed tomography-based staging in comparison with mpMRI staging. Both tumor stage shifts and overall risk group shifts, including prostate-specific antigen (PSA) level and the Gleason score, were assessed. Potential risk factors for upstaging were tested in a multivariate analysis. Finally, the impact of mpMRI-based staging shift on prostate RT and antihormonal therapy was evaluated. RESULTS Overall, tumor stage shift occurred in 55.7% of patients after mpMRI. Upstaging was most prominent in patients showing high-risk serum PSA levels (73%), but was also substantial in patients presenting with low-risk PSA levels (50%) and low-risk Gleason scores (45.2%). Risk group changes occurred in 28.7% of the patients with consequent treatment adaptations regarding target volume delineation and duration of androgen deprivation therapy. High PSA levels were found to be a significant risk factor for tumor upstaging and newly diagnosed seminal vesicle infiltration assessed using mpMRI. CONCLUSIONS Our findings suggest that mpMRI of the prostate leads to substantial tumor upstaging, and can considerably affect treatment decisions in all patient groups undergoing risk-adapted curative RT for prostate cancer.

Accuracy of cartilage-specific 3-Tesla 3D-DESS magnetic resonance imaging in the diagnosis of chondral lesions: comparison with knee arthroscopy

BACKGROUND Arthroscopy is considered as "the gold standard" for the diagnosis of traumatic intraarticular knee lesions. However, recent developments in magnetic resonance imaging (MRI) now offer good opportunities for the indirect assessment of the integrity and structural changes of the knee articular cartilage. The study was to investigate whether cartilage-specific sequences on a 3-Tesla MRI provide accurate assessment for the detection of cartilage defects. METHODS A 3-Tesla (3-T) MRI combined with three-dimensional double-echo steady-state (3D-DESS) cartilage specific sequences was performed on 210 patients with knee pain prior to knee arthroscopy. Sensitivity, specificity, and positive and negative predictive values of magnetic resonance imaging were calculated and correlated to the arthroscopic findings of cartilaginous lesions. Lesions were classified using the modified Outerbridge classification. RESULTS For the 210 patients (1260 cartilage surfaces: patella, trochlea, medial femoral condyle, medial tibia, lateral femoral condyle, lateral tibia) evaluated, the sensitivities, specificities, positive predictive values, and negative predictive values of 3-T MRI were 83.3, 99.8, 84.4, and 99.8 %, respectively, for the detection of grade IV lesions; 74.1, 99.6, 85.2, and 99.3 %, respectively, for grade III lesions; 67.9, 99.2, 76.6, and 98.2 %, respectively, for grade II lesions; and 8.8, 99.5, 80, and 92 %, respectively, for grade I lesions. CONCLUSIONS For grade III and IV lesions, 3-T MRI combined with 3D-DESS cartilage-specific sequences represents an accurate diagnostic tool. For grade II lesions, the technique demonstrates moderate sensitivity, while for grade I lesions, the sensitivity is limited to provide reliable diagnosis compared to knee arthroscopy.

Comparison between magnetic resonance imaging and computed tomography of the lung in patients with cystic fibrosis with regard to clinical, laboratory, and pulmonary functional parameters

OBJECTIVE To evaluate whether magnetic resonance imaging (MRI) is effective as computed tomography (CT) in determining morphologic and functional pulmonary changes in patients with cystic fibrosis (CF) in association with multiple clinical parameters. MATERIALS AND METHODS Institutional review board approval and patient written informed consent were obtained. In this prospective study, 30 patients with CF (17 men and 13 women; mean (SD) age, 30.2 (9.2) years; range, 19-52 years) were included. Chest CT was acquired by unenhanced low-dose technique for clinical purposes. Lung MRI (1.5 T) comprised T2- and T1-weighted sequences before and after the application of 0.1-mmol·kg gadobutrol, also considering lung perfusion imaging. All CT and MR images were visually evaluated by using 2 different scoring systems: the modified Helbich and the Eichinger scores. Signal intensity of the peribronchial walls and detected mucus on T2-weighted images as well as signal enhancement of the peribronchial walls on contrast-enhanced T1-weighted sequences were additionally assessed on MRI. For the clinical evaluation, the pulmonary exacerbation rate, laboratory, and pulmonary functional parameters were determined. RESULTS The overall modified Helbich CT score had a mean (SD) of 15.3 (4.8) (range, 3-21) and median of 16.0 (interquartile range [IQR], 6.3). The overall modified Helbich MR score showed slightly, not significantly, lower values (Wilcoxon rank sum test and Student t test; P > 0.05): mean (SD) of 14.3 (4.7) (range, 3-20) and median of 15.0 (IQR, 7.3). Without assessment of perfusion, the overall Eichinger score resulted in the following values for CT vs MR examinations: mean (SD), 20.3 (7.2) (range, 4-31); and median, 21.0 (IQR, 9.5) vs mean (SD), 19.5 (7.1) (range, 4-33); and median, 20.0 (IQR, 9.0). All differences between CT and MR examinations were not significant (Wilcoxon rank sum tests and Student t tests; P > 0.05). In general, the correlations of the CT scores (overall and different imaging parameters) to the clinical parameters were slightly higher compared to the MRI scores. However, if all additional MRI parameters were integrated into the scoring systems, the correlations reached the values of the CT scores. The overall image quality was significantly higher for the CT examinations compared to the MRI sequences. CONCLUSIONS One major diagnostic benefit of lung MRI in CF is the possible acquisition of several different morphologic and functional imaging features without the use of any radiation exposure. Lung MRI shows reliable associations with CT and clinical parameters, which suggests its implementation in CF for routine diagnosis, which would be particularly important in follow-up imaging over the long term.

A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings. especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging. (C) 2004 Elsevier Inc. All rights reserved.

Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHz)

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to H-1 resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

Magnetic resonance imaging analysis of the upper cervical spine extensor musculature in an asymptomatic cohort: an index of fat within muscle

AIM: To establish a simple method to quantify muscle/fat constituents in cervical muscles of asymptomatic women using magnetic resonance imaging (MRI), and to determine whether there is an age effect within a defined age range. MATERIALS AND METHODS: MRI of the upper cervical spine was performed for 42 asymptomatic women aged 18-45 years. The muscle and fat signal intensities on axial spin echo T1-weighted images were quantitatively classified by taking a ratio of the pixel intensity profiles of muscle against those of intermuscular fat for the rectus capitis posterior major and minor and inferior obliquus capitis muscles bilaterally. Inter- and intra-examiner agreement was scrutinized. RESULTS: The average relative values of fat within the upper cervical musculature compared with intermuscular fat indicated that there were only slight variations in indices between the three sets of muscles. There was no significant correlation between age and fat indices. There were significant differences for the relative fat within the muscle compared with intermuscular fat and body mass index for the right rectus capitis posterior major and right and left inferior obliquus capitis muscles (p = 0.032). Intraclass correlation coefficients for intraobserver agreement ranged from 0.94 to 0.98. Inter-rater agreement of the measurements ranged from 0.75 to 0.97. CONCLUSION: A quantitative measure of muscle/fat constituents has been developed, and results of this study indicate that relative fatty infiltration is not a feature of age in the upper cervical extensor muscles of women aged 18-45 years. (C) 2005 The Royal College of Radiologists. Published by Elsevier Ltd. All rights reserved.

Noninvasive localization of petroleum-derived spray oil in plants with chemical shift selective magnetic resonance imaging

Nuclear magnetic resonance (NMR) spectroscopy and magnetic resonance imaging (MRI) were used to detect petroleum-derived spray oils (PDSOs) in citrus seedlings and trees. The NMR spectrum of the phantom containing 10% (v/v) of a nC24 agricultural mineral oil (AMO) showed the resonance of the water protons at delta = 5 ppm, while the resonance of the oil protons at delta = 1.3 to 1.7 ppm. The peak resolution and the chemical shift difference of more than 3.3 ppm between water and oil protons effectively differentiated water and the oil. Chemical shift selective imaging (CSSI) was performed to localize the AMO within the stems of Citrus trifoliata L. seedlings after the application of a 4% (v/v) spray. The chemical shift selective images of the oil were acquired by excitation at delta = 1.5 ppm by averaging over 400 transients in each phase-encoding step. Oil was mainly detected in the outer cortex of stems within 10 d of spray application; some oil was also observed in the inner vascular bundle and pith of the stems at this point. CSSI was also applied to investigate the persistence of oil deposits in sprayed mature Washington navel orange (Citrus x aurantium L.) trees in an orchard. The trees were treated with either fourteen 0.25%, fourteen 0.5%, four 1.75%, or single 7% sprays of a nC23 horticultural mineral oil (HMO) 12 to 16 months before examination of plant tissues by CSSI, and were still showing symptoms of chronic phytotoxicity largely manifested as reduced yield. The oil deposits were detected in stems of sprayed flushes and unsprayed flushes produced 4 to 5 months after the last spray was applied, suggesting a potential movement of the oil via phloem and a correlation of the persistence of oil deposit in plants and the phytotoxicity. The results demonstrate that MRI is an effective method to probe the uptake and localization of PDSOs and other xenobiotics in vivo in plants noninvasively and nondestructively.

Method for a detailed measurement of image intensity nonuniformity in magnetic resonance imaging

In magnetic resonance imaging (MRI), the MR signal intensity can vary spatially and this spatial variation is usually referred to as MR intensity nonuniformity. Although the main source of intensity nonuniformity arises from B, inhomogeneity of the coil acting as a receiver and/or transmitter, geometric distortion also alters the MR signal intensity. It is useful on some occasions to have these two different sources be separately measured and analyzed. In this paper, we present a practical method for a detailed measurement of the MR intensity nonuniformity. This method is based on the same three-dimensional geometric phantom that was recently developed for a complete measurement of the geometric distortion in MR systems. In this paper, the contribution to the intensity nonuniformity from the geometric distortion can be estimated and thus, it provides a mechanism for estimation of the intensity nonuniformity that reflects solely the spatial characteristics arising from B-1. Additionally, a comprehensive scheme for characterization of the intensity nonuniformity based on the new measurement method is proposed. To demonstrate the method, the intensity nonuniformity in a 1.5 T Sonata MR system was measured and is used to illustrate the main features of the method. (c) 2005 American Association of Physicists in Medicine.

Salt diffusion and distribution in meat studied by Na-23 nuclear magnetic resonance imaging and relaxometry

This study introduces the use of combined Na-23 magnetic resonance imaging (MRI) and Na-23 NMR relaxometry for the study of meat curing. The diffusion of sodium ions into the meat was measured using Na-23 MRI on a 1 kg meat sample brined in 10% w/w NaCl for 3-100 h. Calculations revealed a diffusion coefficient of 1 x 10(-5) cm(2)/s after 3 h of curing and subsequently decreasing to 8 x 10(-6) cm(2)/s at longer curing times, suggesting that changes occur in the microscopic structure of the meat during curing. The microscopic mobility and distribution of sodium was measured using Na-23 relaxometry. Two sodium populations were observed, and with increasing length of curing time the relaxation times of these changed, reflecting a salt-induced swelling and increase in myofibrillar pore sizes. Accordingly, the present study demonstrated that pore size and thereby salt-induced swelling in meat can be assessed using Na-23 relaxometry.

A new decoupling method for quadrature coils in magnetic resonance imaging

A powerful decoupling method is introduced to obtain decoupled signal voltages from quadrature coils in magnetic resonance imaging (MRI). The new method uses the knowledge of the position of the signal source in MRI, the active slice, to define a new mutual impedance which accurately quantifies the coupling voltages and enables them to be removed almost completely. Results show that by using the new decoupling method, the percentage errors in the decoupled voltages are of the order of 10(-7)% and isolations between two coils are more than 170 dB.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.

Comparison of two-dimensional speckle and tissue velocity based strain and validation with harmonic phase magnetic resonance imaging

Two-dimensional (2-D) strain (epsilon(2-D)) on the basis of speckle tracking is a new technique for strain measurement. This study sought to validate epsilon(2-D) and tissue velocity imaging (TVI)based strain (epsilon(TVI)) with tagged harmonic-phase (HARP) magnetic resonance imaging (MRI). Thirty patients (mean age. 62 +/- 11 years) with known or suspected ischemic heart disease were evaluated. Wall motion (wall motion score index 1.55 +/- 0.46) was assessed by an expert observer. Three apical images were obtained for longitudinal strain (16 segments) and 3 short-axis images for radial and circumferential strain (18 segments). Radial epsilon(TVI) was obtained in the posterior wall. HARP MRI was used to measure principal strain, expressed as maximal length change in each direction. Values for epsilon(2-D), epsilon(TVI), and HARP MRI were comparable for all 3 strain directions and were reduced in dysfunctional segments. The mean difference and correlation between longitudinal epsilon(2-D) and HARP MRI (2.1 +/- 5.5%, r = 0.51, p < 0.001) were similar to those between longitudinal epsilon(TVI), and HARP MRI (1.1 +/- 6.7%, r = 0.40, p < 0.001). The mean difference and correlation were more favorable between radial epsilon(2-D) and HARP MRI (0.4 +/- 10.2%, r = 0.60, p < 0.001) than between radial epsilon(TVI), and HARP MRI (3.4 +/- 10.5%, r = 0.47, p < 0.001). For circumferential strain, the mean difference and correlation between epsilon(2-D) and HARP MRI were 0.7 +/- 5.4% and r = 0.51 (p < 0.001), respectively. In conclusion, the modest correlations of echocardiographic and HARP MRI strain reflect the technical challenges of the 2 techniques. Nonetheless, epsilon(2-D) provides a reliable tool to quantify regional function, with radial measurements being more accurate and feasible than with TVI. Unlike epsilon(TVI), epsilon(2-D) provides circumferential measurements. (c) 2006 Elsevier Inc. All rights reserved.