Bayesian image decomposition applied to relaxographic imaging

T1 relaxographic imaging is a precise and accurate way to characterize tissue. A number of fast MRI acquisition techniques allow both spatial and magnetization recoveries to be well sampled in reasonable imaging times. However, two limitations common to the analysis of relaxographic imaging data are (1) the assumption of single exponential behavior for each image voxel and (2) the treatment of each pixel as an independent entity. The first assumption disregards tissue heterogeneity known to be present and reduces the information content that can be extracted. The latter assumption reduces both the modeling stability and the accuracy of extracted parameters. A new method that overcomes these limitations is presented here. The method, Bayesian Image Decomposition, recovers individual tissue type magnetization recovery curves and their corresponding tissue‐specific relaxographic images (i.e. segmented images) from a series of inversion recovery images. The general form of the decomposition is given together with its specific implementation to longitudinal relaxographic imaging. The method is validated by comparison of the results with those of the standard method and by comparison across data sets. A specific advantage of the new method is the ability to determine fractional contributions of tissue subtypes to each image voxel. © 2005 Wiley Periodicals, Inc. Int J Imaging Syst Technol, 15;2–9, 2005; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/ima.20033

[1]  T. Brown,et al.  A new method for spectral decomposition using a bilinear Bayesian approach. , 1999, Journal of magnetic resonance.

[2]  J. Olesen,et al.  In vivo determination of T1 and T2 in the brain of patients with severe but stable multiple sclerosis , 1988, Magnetic resonance in medicine.

[3]  John Skilling,et al.  Prior Distributions on Measure Space , 1997 .

[4]  P. Allen,et al.  Proton relaxation in acute and subacute ischemic brain edema. , 1990, Advances in neurology.

[5]  Proton-nuclear magnetic resonance relaxation times in brain edema. , 1990, Advances in neurology.

[6]  M. Weiner,et al.  Metabolite 1H relaxation in normal and hyponatremic brain , 1996, Magnetic resonance in medicine.

[7]  J S Fowler,et al.  Intimate combination of low‐ and high‐resolution image data: I. real‐space PET and 1H2O MRI, PETAMRI , 1999, Magnetic resonance in medicine.

[8]  T R Brown,et al.  NMR spectral quantitation by principal-component analysis. II. Determination of frequency and phase shifts. , 1996, Journal of magnetic resonance. Series B.

[9]  M. Osbakken,et al.  Spin‐lattice relaxation (T1) times of cerebral white matter in multiple sclerosis , 1986, Magnetic resonance in medicine.

[10]  Ole Winther,et al.  Mean-Field Approaches to Independent Component Analysis , 2002, Neural Computation.

[11]  Lowell O. Randall,et al.  CHEMICAL TOPOGRAPHY OF THE BRAIN , 1938 .

[12]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[13]  C. Springer,et al.  Relaxographic imaging. , 1994, Journal of magnetic resonance. Series B.

[14]  Truman R. Brown,et al.  Application of Principal-Component Analysis for NMR Spectral Quantitation , 1995 .

[15]  G Johnson,et al.  Precise relaxation time measurements of normal‐appearing white matter in inflammatory central nervous system disease , 1989, Magnetic resonance in medicine.

[16]  Walter R. Gilks,et al.  [Bayesian Computation and Stochastic Systems]: Comment , 1995 .

[17]  Y Horikawa,et al.  Significance of proton relaxation time measurement in brain edema, cerebral infarction and brain tumors. , 1986, Magnetic resonance imaging.

[18]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.