Speeding up the EM algorithm for mixture model-based segmentation of magnetic resonance images

Mixture models implemented via the expectation-maximization (EM) algorithm are being increasingly used in a wide range of problems in pattern recognition such as image segmentation. However, the EM algorithm requires considerable computational time in its application to huge data sets such as a three-dimensional magnetic resonance (MR) image of over 10 million voxels. Recently, it was shown that a sparse, incremental version of the EM algorithm could improve its rate of convergence. In this paper, we show how this modified EM algorithm can be speeded up further by adopting a multiresolution kd-tree structure in performing the E-step. The proposed algorithm outperforms some other variants of the EM algorithm for segmenting MR images of the human brain. (C) 2004 Pattern Recognition Society. Published by Elsevier Ltd. All rights reserved.

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.

Big cat scan: magnetic resonance imaging of the tiger

In August 2002, we performed MRI scans on a female juvenile Bengal tiger. We present the clinical course, imaging and autopsy findings, and some comparative anatomy of the tiger brain and skull. Magnetic resonance images of a tiger have not previously been published.

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.

3-D magnetic resonance angiography versus conventional angiography in peripheral arterial disease: Pilot study

Angiography is usually performed as the preoperative road map for those requiring revascularization for lower extremity peripheral arterial disease (PAD). The alternative investigations are ultrasound, 3-D magnetic resonance angiography (3-D MRA) and computed tomography angiography. This pilot study aimed to assess whether 3-D MRA could replace the gold standard angiography in preoperative planning. Eight patients considered for aortoiliac or infrainguinal arterial bypass surgery were recruited. All underwent both imaging modalities within 7 days. A vascular surgeon and a radiologist each reported on the images from both the 3-D MRA and the angiography, with blinding to patient details and each others reports. Comparisons were made between the reports for the angiographic and the 3-D MRA images, and between the reports of the vascular surgeon and the radiologist. Compared to the gold standard angiogram, 3-D MRA had a sensitivity of 77% and specificity of 94% in detecting occlusion, and a sensitivity of 72% and specificity of 90% in differentiating high grade (> 50%) versus low grade (< 50%) stenoses. There was an overall concordance of 78% between the two investigations with a range of 62% in the peroneal artery to 94% in the aorta. 3-D MRA showed flow in 23% of cases where conventional angiography showed no flow. In the present pilot study, 3-D MRA had reasonable concordance with the gold standard angiography, depending on the level of the lesion. At times it showed vessel flow where occlusion was shown on conventional angiogram. 3-D MRA in peripheral vascular disease is challenging the gold standard, but is inconsistent at present.

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.

Magnetic resonance microscopy of the equine hoof wall: a study of resolution and potential

Reasons for performing study: Obtaining magnetic resonance images of the inner hoof wall tissue at the microscopic level would enable early accurate diagnosis of laminitis and therefore more effective therapy. Objectives: To optimise magnetic resonance imaging (MRI) parameters in order to obtain the highest possible resolution of the structures beneath the equine hoof wall. Methods: Magnetic resonance microscopy (MRM) was performed in front feet from 6 cadaver horses using T-2-weighted fast spin echo (FSE-T-2), and T-1-weighted gradient echo (GRE-T-1) sequences. Results: In T-2 weighted FSE images most of the stratum medium showed no signal, however the coronary, terminal and sole papillae were visible. The stratum lamellatum was clearly visible and primary epidermal lamellae could be differentiated from dermal lamellae. Conclusion: Most structures beneath the hoof wall were differentiated. Conventional scanners for diagnostic MRI in horses are low or high field. However this study used ultra-high field scanners currently not available for clinical use. Signal-to-noise ratio (SIN) increases as a function of field strength. An increase of spatial resolution of the image results in a decreased SIN. SIN can also be improved with better coils and the resolution of high field MRI scanners will increase as technology develops and surface array coils become more readily available. Potential relevance: Although MR images with microscopic resolution were obtained ex vivo, this study demonstrates the potential for detection of lamellar pathology as it occurs. Early recognition of the development of laminitis to instigate effective therapy at an earlier stage and may improve the outcome for laminitic horses. Clinical MR is now readily available at 3 T, while 4 T, 7 T and 9 T systems are being used for human whole body applications.

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.