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.

Impact of scar thickness on the assessment of viability using dobutamine echocardiography and thallium single-photon emission computed tomography - A comparison with contrast-enhanced magnetic resonance imaging

OBJECTIVES We sought to determine whether the transmural extent of scar (TES) explains discordances between dobutamine echocardiography (DbE) and thallium single-photon emission computed tomography (Tl-SPECT) in the detection of viable myocardium (VM). BACKGROUND Discrepancies between DbE and Tl-SPECT are often attributed to differences between contractile reserve and membrane integrity, but may also reflect a disproportionate influence of nontransmural scar on thickening at DbE. METHODS Sixty patients (age 62 +/- 12 years; 10 women and 50 men) with postinfarction left ventricular dysfunction underwent standard rest-late redistribution Tl-SPECT and DbE. Viable myocardium was identified when dysfunctional segments showed Tl activity >60% on the late-redistribution image or by low-dose augmentation at DbE. Contrast-enhanced magnetic resonance imaging (ceMRI) was used to divide TES into five groups: 0%, 75% of the wall thickness replaced by scar. RESULTS As TES increased, both the mean Tl uptake and change in wall motion score decreased significantly (both p < 0.001). However, the presence of subendocardial scar was insufficient to prevent thickening; >50% of segments still showed contractile function with TES of 25% to 75%, although residual function was uncommon with TES >75%. The relationship of both tests to increasing TES was similar, but Tl-SPECT identified VM more frequently than DbE in all groups. Among segments without scar or with small amounts of scar (50% were viable by SPECT. CONCLUSIONS Both contractile reserve and perfusion are sensitive to the extent of scar. However, contractile reserve may be impaired in the face of no or minor scar, and thickening may still occur with extensive scar. (C) 2004 by the American College of Cardiology Foundation.

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.

A proposed scheme for comprehensive characterization of the measured geometric distortion in magnetic resonance imaging using a three-dimensional phantom

Recently, a 3-dimensional phantom that can provide a comprehensive, accurate and complete measurement of the geometric distortion in MRI has been developed. In this paper, a scheme for characterizing the measured geometric distortion using the 3-D phantom is described. In the proposed scheme, a number of quantitative measures are developed and used to characterize the geometric distortion. These measures encompass the overall and spatial aspects of the geometric distortion. Two specific types of volume of interest, rectangular parallelepipeds (including cubes) and spheres are considered in the proposed scheme. As an illustration, characterization of the geometric distortion in a Siemens 1.5T Sonata MRI system using the proposed scheme is presented. As shown, the proposed scheme provides a comprehensive assessment of the geometric distortion. The scheme can be potentially used as a standard procedure for the assessment of geometric distortion in MRI. (C) 2004 American Association of Physicists in Medicine.

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.

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.

The use of simultaneous H-1 & P-31 magic angle spinning nuclear magnetic resonance measurements to characterize energy metabolism during the conversion of muscle to meat

In order to investigate the potential of magic angle spinning nuclear magnetic resonance (MAS NMR) in the elucidation of post-mortem metabolism in muscle biopsies, simultaneous H-1 and (31)p MAS NMR measurements were made continuously on postmortem (20 min to 24 h) muscle longissimus samples from rabbits. The animals had either been or not been given adrenaline (0.5 mg kg(-1) 4 h pre-slaughter) to deplete stores of muscle glycogen. The intracellular pH was calculated from H-1 spectra, and the post-mortem rate of formation of lactate was followed and quantified. Comparison of measurements made on muscle samples from rabbits treated with adrenaline with measurements made on muscle samples from untreated' rabbits revealed significant effects of adrenaline treatment on both pH (pH24 h = 6.42 vs. pH24 It = 5.60) and formation of lactate (16 mmol g(-1) vs. 65 mmol g(-1)). The P-31 NMR spectra were used to follow the rate of degradation of ATP and phosphocreatine. The present study clearly shows that MAS NMR has potential for the study of post-mortem energy metabolism.

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.

Screening with magnetic resonance imaging and mammography of a UK population at high familial risk of breast cancer: a prospective multicentre cohort study (MARIBS)

Background Women genetically predisposed to breast cancer often develop the disease at a young age when dense breast tissue reduces the sensitivity of X-ray mammography. Our aim was, therefore, to compare contrast enhanced magnetic resonance imaging (CE MRI) with mammography for screening. Methods We did a prospective multicentre cohort study in 649 women aged 35-49 years with a strong family history of breast cancer or a high probability of a BRCA1, BRCA2, or TP53 mutation. We recruited participants from 22 centres in the UK, and offered the women annual screening with CE MRI and mammography for 2-7 years. Findings We diagnosed 35 cancers in the 649 women screened with both mammography and CE MRI (1881 screens): 19 by CE MRI only, six by mammography only, and eight by both, with two interval cases. Sensitivity was significantly higher for CE MRI (77%, 95% CI 60-90) than for mammography (40%, 24-58; p=0.01), and was 94% (81-99) when both methods were used. Specificity was 93% (92-95) for mammography, 81% (80-83) for CE MRI (p

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.

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.

Synthesis of gadolinium-labeled shell-crosslinked nanoparticles for magnetic resonance imaging applications

Shell-crosslinked knedel-like nanoparticles (SCKs; knedel is a Polish term for dumplings) were derivatized with gadolinium Shell chelates and studied as robust magnetic-resonance-imaging-active structures with hydrodynamic diameters of 40 +/- 3 nm. SCKs possessing an amphiphilic core-shell morphology were produced from the aqueous assembly of diblock copolymers of poly(acrylic acid) (PAA) and poly(methyl acrylate) (PMA), PAA(52)-b-PMA(128), and subsequent covalent crosslinking by amidation upon reaction with 2,2'-(ethylenedioxy)bis(ethylamine) throughout the shell layer. The properties of these materials, including non-toxicity towards mammalian cells, non-immunogenicity within mice, and capability for polyvalent targeting, make them ideal candidates for utilization within biological systems. The synthesis of SCKs derivatized with Gd-III and designed for potential use as a unique nanometer-scale contrast agent for MRI applications is described herein. Utilization of an amino-functionalized diethylenetriaminepentaacetic acid-Gd analogue allowed for direct covalent conjugation throughout the hydrophilic shell layer of the SCKs and served to increase the rotational correlation lifetime of the Gd. In addition, the highly hydrated nature of the shell layer in which the Gd was located allowed for rapid water exchange; thus, the resulting material demonstrated large ionic relaxivities (39 s(-1) mM(-1)) in an applied magnetic field of 0.47 T at 40 degrees C and, as a result of the large loading capacity of the material, also demonstrated high molecular relaxivities (20 000 s(-1) mM(-1)).