Method of magnetic resonance reconstruction imaging from projections using partial data collected in k-space

• US 5,243,284 A
• Filed: 07/24/1991
• Issued: 09/07/1993
• Est. Priority Date: 07/24/1991
• Status: Expired due to Term

• Current Assignee: Board of Trustees of the Leland Stanford Junior University (Stanford University)
• Original Assignee: Board of Trustees of the Leland Stanford Junior University (Stanford University)
• Inventors: Noll, Douglas C.
• Primary Examiner(s): Tokar, Michael J.
• Assistant Examiner(s): MAH, RAYMOND
• Application Number:
  US07/736,268
• Time in Patent Office: 776 Days
• Field of Search: 382/47, 382/68, 324/300, 324/307, 324/309, 324/310, 324/312
• US Class Current: 324/309
• CPC Class Codes: G01R 33/4824 using a non-Cartesian trajectory
• CPC Class Codes: G01R 33/561 by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences

Abstract

The acquisition of radial lines in 2-D Fourier space (k-space) allows the realization of the extremely short echo times useful in the imaging of short-T₂ species and flow using Magnetic Resonance (MR). The disadvantage of this imaging method, which is also known as projection reconstruction, is that in order to prevent aliasing at a given resolution, π times as many excitations are required as the conventional 2-D Fourier transform or spin-warp acquisition. Disclosed is an acquisition and reconstruction method that reduces the required number of radial lines and excitations by using the properties that projections are real or have slowly varying phase and the at radial acquisition methods are oversampled at the center of k-space. This method preserves resolution while reducing imaging time at the expense of signal-to-noise ratio (SNR). This has been verified with both phantom and human subjects.

Claims

1. In magnetic resonance imaging of one object based on projection reconstruction, a method of obtaining data for k-space trajectories comprising the steps of

a) establishing a set of radial k-space trajectories all originating at the center k-space for use in projection reconstruction of an image of said object, each trajectory having a low frequency part and a high frequency part,

b) sequentially exciting said object and detecting data samples along less than all of said set of radial k-space trajectories, and

c) obtaining data samples for low frequency parts of unsampled trajectories based on interpolation of detected data samples of low frequency parts of adjacent sampled trajectories.

2. The method as defined by claim 1 wherein all unsampled trajectories are opposite to and collinear extensions of sampled trajectories, and further including d) obtaining data samples for high frequency parts of unsampled trajectories based on the sampled high frequency parts of linearly extending sampled trajectories.

3. The method as defined by claim 2 wherein step d) includes establishing a phase reference for said data samples for high frequency parts, said phase reference being for use in determining high frequency parts of unsampled trajectories based on a reconstructed low frequency projection.

4. The method as defined by claim 3 and further including doubling said sampled high frequency parts.

5. The method as defined by claim 3 wherein a sampled trajectory has a low frequency part, a'"'"', and a high frequency parts, a;

a collinear unsampled trajectory has a low frequency a part, b'"'"'; and

a high frequency part, b; and

a projection of the line (a,a'"'"',b,b'"'"') is defined as

space="preserve" listing-type="equation">p.sub.θ

(r)e.sup.iφ

.sbsp.θ

.sup.(r) =p.sub.l (r)e.sup.iφ

.sbsp.l.sup.(r) +p.sub.h (r)e.sup.iφ

.sbsp.h.sup.(r)

where p_l (r)eⁱφ

.sbsp.l.sup.(r) and p_h (r)eⁱφ

.sbsp.h.sup.(r) are defined as the inverse Fourier Transform of (a'"'"',b'"'"') and (a,b), respectively, and

p.sub.θ

(r)=projection magnitude

φ

.sub.θ

(r)=projection phase

p_l (r)=projection magnitude low frequency parts

φ

_l (r)=projection low frequency phase

p_h (r)=projection magnitude high frequency parts

φ

_h (r)=projection high frequency phase.

6. The method as defined by claim 5 wherein ##EQU3## where φ

_l (r) is a phase reference based on a low frequency image, and assuming that φ

_h (r)=φ

.sub.θ

(r)=φ

_l (r), and

F₁^-1_D =sampled and interpolated low frequency parts,

F₁^-1_D =sampled high frequency parts.

7. The method as defined by claim 6 wherein the projection for a k-space trajectory is based on the in-phase portion, p.sub.θ

(r)eⁱφ

.sbsp.θ

.sup.(r).

8. A method of magnetic resonance imaging using projection reconstruction comprising the steps of

a) obtaining data samples for a first plurality of radial k-space trajectories all originating at the center of k-space, each trajectory having a low frequency part and a high frequency part,

b) obtaining data for a second plurality of radial k-space trajectories having low frequency parts and high frequency parts each of said second plurality of trajectories being opposite to and collinear extensions of one of said first plurality of trajectories, including interpolating low frequency parts from adjacent sampled trajectories,

c) Fourier transforming collinear trajectories of said first plurality of trajectories and said second plurality of trajectories to obtain projections, and

d) reconstructing an image based on said projections.

9. The method as defined by claim 8 wherein a sampled trajectory has a low frequency parts, a'"'"', and a high frequency part, a;

a collinear unsampled trajectory has a low frequency apart, b'"'"', and a high frequency part, b; and

a projection of the line (a,a'"'"',b'"'"',b) is defined as

space="preserve" listing-type="equation">p.sub.θ

(r)e.sup.iφ

.sbsp.θ

.sup.(r) =p.sub.l (r)e.sup.iφ

.sbsp.l.sup.(r) +p.sub.h (r)e.sup.iφ

.sbsp.h.sup.(r)

where p_l (r)eⁱφ

.sbsp.l.sup.(r) and p_h (r)eⁱφ

.sbsp.h.sup.(r) are defined as the inverse Fourier Transform of (a'"'"',b'"'"') and (a,b), respectively, and

p.sub.θ

(r)=projection magnitude

φ

.sub.θ

(r)=projection phase

p_l (r)=projection magnitude low frequency parts

φ

_l (r)=projection low frequency phase

p_h (r)=projection magnitude high frequency parts

φ

_h (r)=projection high frequency phase.

10. The method as defined by claim 9 and further including doubling said sampled high frequency parts, a, in obtaining high frequency unsampled parts, b.

11. The method as defined by claim 10 wherein ##EQU4## wherein φ

_l (r) is a phase reference based on a low frequency image, and assuming that φ

_h (r)=φ

.sub.θ

(r)=φ

_l (r), and

F₁^-1_D =sampled and interpolated low frequency parts,

R₁^-1_D =sampled high frequency parts.

12. The method as defined by claim 11 wherein the projection for a k-space trajectory is based on the in-phase portion, p.sub.θ

(r)eⁱφ

.sbsp.θ

.sup.(r).