Systems and methods for statistical reconstruction of magnetic resonance fingerprinting data

Patent metadata

Field	Value
Publication number	US-10241176-B2
Application number	US-201715409378-A
Country	US
Kind code	B2
Filing date	Jan 18, 2017
Priority date	Jan 20, 2016
Publication date	Mar 26, 2019
Grant date	Mar 26, 2019

Abstract

Systems and methods for reconstructing magnetic resonance (MR) tissue parameter maps of a subject from magnetic resonance fingerprinting (MRF) data acquired using a magnetic resonance imaging (MRI) system. The method includes providing MRF data acquired from a subject using an MRI system and performing an iterative, maximum-likelihood reconstruction of the MRF data to create MR tissue parameter maps of the subject.

First claim

The invention claimed is:

1. A magnetic resonance imaging (MRI) system comprising:
a magnet system configured to generate a polarizing magnetic field about at least a portion of a subject arranged in the MRI system;
a plurality of gradient coils configured to apply a gradient field to the polarizing magnetic field;
a radio frequency (RF) system configured to apply an excitation field to the subject and acquire MR image data from a ROI;
a computer system programmed to:
control the plurality of gradient coils and the RF system to acquire magnetic resonance fingerprinting (MRF) data from a subject; and
perform an iterative, maximum-likelihood reconstruction of the MRF data to create magnetic resonance (MR) tissue parameter maps of the subject.

2. The system of claim 1 wherein the iterative maximum-likelihood reconstruction models data, including the Bloch equation based magnetization dynamics and noise characteristics of the MRF data.

3. The system of claim 1 wherein the computer system is programmed to perform the iterative maximum-likelihood reconstruction by estimating tissue parameters directly from the MRF data without reconstructing contrast-weighted images.

4. The system of claim 3 wherein, to perform the maximum-likelihood reconstruction, the computer system is further programmed to analyze the estimated tissue parameters for convergence and, if convergence is not reached, perform an iteration to estimate updated tissue parameters from the MRF data.

5. The system of claim 1 wherein, to perform the maximum-likelihood reconstruction, the computer system is further programmed to represent an image reconstructed from the MRF data as:
I m =Φ m ( T 1 ,T 2 ,f 0 )ρ

wherein T 1 denotes longitudinal relaxation time, T 2 denotes transverse relaxation time, ρ denotes spin density, f 0 and denotes off-resonance frequency, where both I m and ρ are N×1 complex vectors, T 1 , T 2 , and f 0 are all N×1 real vectors, and Φ m (⋅)∈£ N×N is a diagonal matrix with [Φ m ] n,n =ϕ m (T 1 [n],T 2 [n],f 0 [n]).

6. The system of claim 5 wherein, to perform the maximum-likelihood reconstruction, the computer system is programmed to solve a maximum likelihood estimation problem of:

{

T
^

1

,

T
^

2

,

f
^

0

,

ρ
^

}

=

arg

T
1

,

T
2

,

f
0

,
ρ

⁢
min
⁢

∑

c
=
1

N
c

⁢

∑

m
=
1

M

⁢



d

m
,
c

-

F
m

⁢

S
c

⁢

Φ
m

⁡

(

T
1

,

T
2

,

f
0

)

⁢
ρ



2
2

wherein F m ∈£ P×N denotes the undersampled Fourier encoding matrix, S c ∈£ N×N is a diagonal matrix whose diagonal entries contain the coil sensitivities, n m ∈£ P×1 denotes the noise vector, and d m,c =F m S c Φ m (T 1 ,T 2 ,f 0 )ρ+n m,c , assuming that {n m,c } m,c=1 M,N c is Gaussian noise.

7. The system of claim 1 wherein the computer system is further programmed to solve an optimization problem presented by the maximum-likelihood reconstruction using at least one of a variable-splitting technique, an alternating direction method of multipliers (ADMM), and a variable projection (VARPRO) method.

8. The system of claim 7 wherein the variable splitting-technique separates a Bloch equation based physical model from a data consistency constraint.

9. The system of claim 1 wherein performing the maximum-likelihood reconstruction includes finding a maximum-likelihood optimal solution.

10. The system of claim 1 wherein, to create MR tissue parameter maps of the subject, the computer system is further programmed to compute bounding theoretical error-bars for tissue parameters that form the MR tissue parameter maps.

11. The system of claim 10 wherein the computer system is further programmed to calculate a Cramer-Rao bound (CRB) to provide error bars for tissue parameter estimates.

Assignees

• Massachusetts Gen Hospital

Inventors

• Zhao Bo
• Setsompop Kawin
• Wald Lawrence

Classifications

• G01R33/5608 Primary

  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

• G01R33/50

  based on the determination of relaxation times {, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences} · CPC title