Multi-contrast image reconstruction with joint bayesian compressed sensing

Patent metadata

Field	Value
Publication number	US-9336611-B2
Application number	US-201113822102-A
Country	US
Kind code	B2
Filing date	Sep 13, 2011
Priority date	Sep 14, 2010
Publication date	May 10, 2016
Grant date	May 10, 2016

Abstract

A method for reconstructing multiple images of a subject depicting multiple different contrast characteristics from medical image data acquired with a medical imaging system is provided. Multiple image data sets are acquired with one or more medical imaging systems and the image data sets used to estimate hyperparameters drawn from a prior distribution, such as a prior distribution of image gradient coefficients. These hyperparameters and the acquired image data sets are utilized to produce a posterior distribution, such as a posterior distribution of image gradients. From this posterior distribution, multiple images with the different contrast characteristics are reconstructed. The medical imaging system may be a magnetic resonance imaging system, an x-ray computed tomography imaging system, an ultrasound system, and so on.

First claim

The invention claimed is:

1. A method for reconstructing images of a subject with a magnetic resonance imaging (MRI) system, the steps of the method comprising:
a) acquiring with the MRI system, multiple image data sets with different contrast characteristics;
b) estimating from the multiple image data sets, hyperparameters related to a prior distribution;
c) producing a posterior distribution of image gradient coefficients using the estimated hyperparameters; and
d) producing multiple images using the posterior distribution and the multiple image data sets, each of the multiple images having an image contrast associated with a contrast characteristic of the respective one of the multiple image data sets.

2. The method as recited in claim 1 in which step d) includes deriving image gradients from the posterior distribution.

3. The method as recited in claim 2 in which the image gradients are derived in step d) as the mean of the posterior distribution.

4. The method as recited in claim 2 in which step d) includes producing transformed image gradients by performing a Fourier transform on the image gradients.

5. The method as recited in claim 4 in which step d) includes producing multiple refined image data sets using the transformed image gradients.

6. The method as recited in claim 5 in which the multiple images are reconstructed by performing a Fourier transform on the multiple refined image data sets.

7. The method as recited in claim 1 in which step c) includes using Bayes' theorem to produce the posterior distribution using the hyperparameters and the acquired multiple image data sets.

8. The method as recited in claim 1 in which the hyperparameters are estimated in step b) by maximizing a log likelihood.

9. The method as recited in claim 1 in which the prior distribution related to the hyperparameters is a Gamma distribution.

10. The method as recited in claim 1 in which the hyperparameters are estimated in step b) using an expectation-maximization algorithm.

11. The method as recited in claim 10 in which step a) includes acquiring at least one image data set that samples k-space in satisfaction of a Nyquist criterion, and reconstructing a prior image from the at least one image data set.

12. The method as recited in claim 11 in which the prior image is used to initialize the expectation-maximization algorithm.

13. A method for reconstructing images of a subject from medical image data acquired with a medical imaging system, the steps of the method comprising:
a) acquiring at least two medical image data sets with at least one medical imaging system, each of the at least two medical image data sets being acquired with a selected contrast characteristic;
b) estimating from the at least two medical image data sets, hyperparameters related to a prior distribution of image gradient coefficients;
c) producing a posterior distribution of image gradients using the estimated hyperparameters and the at least two medical image data sets; and
d) reconstructing images of the subject using the produced posterior distribution of image gradients.

14. The method as recited in claim 13 in which the at least one medical imaging system includes at least one of a magnetic resonance imaging system, a computed tomography system, an x-ray imaging system, a positron emission tomography system, an ultrasound system, an optical imaging system, an electroencephalography system, and a magnetoencephalography system.

15. The method as recited in claim 14 in which the at least one medical imaging system is a magnetic resonance imaging system and the selected contrast characteristic is at least one of a T 1 -weighting, a T 2 -weighting, a T* 2 -weighting, a T 1ρ -weighting, aT 2ρ -weighting, a proton density weighting, and a diffusion weighting.

16. The method as recited in claim 14 in which the at least one medical imaging system is at least one of a computed tomography system and an x-ray imaging system, and the selected contrast characteristic is at least one of an x-ray energy level and an x-ray tube current.

17. The method as recited in claim 13 in which the hyperparameters are estimated in step b) by maximizing a log likelihood.

18. The method as recited in claim 13 in which the prior distribution of image gradient coefficients is a Gamma distribution.

19. The method as recited in claim 13 in which step d) includes deriving image gradients from the posterior distribution of image gradients.

20. The method as recited in claim 19 in which step d) includes producing transformed image gradients by performing a Fourier transform on the derived image gradients.

21. The method as recited in claim 20 in which step d) includes producing at least two refined image data sets using the transformed image gradients.

22. The method as recited in claim 21 in which the images of the subject are reconstructed by performing a Fourier transform on the at least two refined image data sets.

23. The method as recited in claim 13 in which the hyperparameters are estimated in step b) using an expectation-maximization algorithm.

24. A method for reconstructing images of a subject from image data acquired with a magnetic resonance imaging (MRI) system, the steps of the method comprising:
a) acquiring at least two image data sets with the MRI system, wherein the at least two medical image data sets are selected from the group consisting of a structural data set, a spectroscopic data set, a functional MRI (fMRI) data set, a contrast-enhanced angiography data set, a non-contrast enhanced angiography data set, a velocity-encoded data set, a phase data set, and a magnitude data set;
b) estimating from the at least two image data sets, hyperparameters related to a prior distribution of image gradient coefficients;
c) producing a posterior distribution of image gradients using the estimated hyperparameters and the at least two image data sets; and
d) reconstructing images of the subject using the produced posterior distribution of image gradients.

Assignees

• Bilgic Berkin
• Adalsteinsson Elfar
• Massachusetts Inst Technology

Inventors

• Bilgic Berkin
• Adalsteinsson Elfar

Classifications

• G01R33/5602

  by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse · CPC title

• G01R33/5635

  Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA] · CPC title

• G01R33/4806

  Functional imaging of brain activation · CPC title

• G16H50/20

  for computer-aided diagnosis, e.g. based on medical expert systems · CPC title

• G01R33/5608

  Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels (image data processing or generation, in general G06T) · CPC title