# Patent

## Title
Simultaneous acquisitions of spin- and stimulated-echo planar imaging

## Abstract

The disclosure provides echo planar imaging (EPI) based single-shot imaging techniques for acquiring spin-EPI (SEPI) and stimulated-EPI after a single RF excitation. In certain embodiments, the SEPI and STEPI acquired in a singleshot are used to compute a T1 map in realtime, which can be used for realtime monitoring of concentrations of paramagnetic-ion based contrast agent in dynamic contrast enhanced MRI. In certain embodiments, B1 field inhomogeneity correction is provided for the STEPI. In certain embodiments, a gradient-EPI (GEPI) is also acquired after the single RF excitation providing singleshot SEPI, GEPI and STEPI acquisition. In certain embodiments, the phase difference between the SEPI and GEPI is used to compute temperature-dependent chemical shift in a subject in realtime. The temperature-dependent chemical shift can be used to monitor temperature changes in the subject in realtime, for example, during a heat treatment to control the dosage of heat energy.

## Background
FIELD
This application relates to magnetic resonance imaging.
BACKGROUND
Spin-echo sequences, gradient-echo pulse sequences, and their hybrid are commonly used for data acquisition in magnetic resonance imaging (MRI). The stimulated-echo (STE) pulse sequence is often used to investigate the transport phenomena, such as molecular translational diffusion and displacement-encoded imaging (DENSE). STE is also used in localized single-voxel spectroscopy for studying short T2 species. In STE imaging, half of the prepared magnetization is restored into the longitudinal space by a second 90° RF (tipup) pulse, and other half of the prepared magnetization, which remains on the transverse plane after the tipup RF pulse, is discarded. This results in wasting of a quadrature component of the magnetization and quadruples the acquisition time to achieve a signal-to-noise ratio (SNR) comparable to spin-echo MRI. Considering the low sensitivity of NMR measurements, stimulated echo NMR/MRI is an inefficient measurement technique.
SUMMARY
The disclosure provides echo planar imaging (EPI) based singleshot imaging techniques for acquiring spin-EPI (SEPI) and stimulated-EPI after a single RF excitation. In certain embodiments, the techniques uses half of the prepared magnetization that is discarded in conventional stimulated echo imaging to acquire the SEPI. Because the SEPI and STEPI are acquired substantially simultaneously, the SEPI and STEPI can be used to perform quantitative analysis (e.g., T1 estimation) in realtime that is free from motion-related artifacts.
In certain embodiments, the SEPI and STEPI acquired in a singleshot are used to compute a T1 map within the imaging field of view (FOV) in realtime. Because the SEPI and STEPI are acquired in a singleshot, the T1 map is free from motion-related artifacts. In certain embodiments, the T1 map is used for realtime monitoring of concentrations of paramagnetic-ion based contrast agent in dynamic contrast enhanced MRI.
In certain embodiments, the STEPI is modified to correct for inhomogeneity in the B1 field. This is done by acquiring a STEPI with minimum mixing time (TM) and generating a B1 correction map using the STEPI acquired with minimum TM and the SEPI. The STEPI is then modified based on the B1 correction map to correct for B1 field inhomogeneity. This correction enables accurate quantitative analysis using the SEPI and STEP that is free from error due to B1 field inhomogeneity.
In certain embodiments, a gradient-EPI (GEPI) is also acquired after the single RF excitation providing singleshot SEPI, GEPI and STEPI acquisition. In certain embodiments, the phase difference between the SEPI and GEPI is used to compute temperature-dependent chemical shift in a subject in realtime. The temperature-dependent chemical shift can be used to monitor temperature changes in the subject in realtime. For example, realtime monitoring of temperature changes can be used in a heat treatment to control the dosage of heat energy during the heat treatment.
Additional features and advantages of the invention will be set forth in the description below, and in part will be apparent from the description, or may be learned by practice of the invention. The advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
It is to be understood that both the foregoing summary of the invention and the following detailed description are exemplary and explanatory and are intended to provide further explanation of the invention as claimed.



BRIEF DESCRIPTION OF THE DRAWINGS
The accompanying drawings, which are included to provide further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate aspects of the invention and together with the description serve to explain the principles of the invention.
FIG. 1 shows a pulse sequence diagram for acquiring a singleshot SEPI and STEPI (2D ss-SESTEPI) according to an embodiment of the invention.
FIG. 2 shows a plot of estimated T1 without B1 correction for different flipangles and mixing times.
FIG. 3 shows a pulse sequence diagram for acquiring SEPI with minimum mixing time according to an embodiment of the invention.
FIG. 4 shows a pulse sequence diagram for acquiring a singleshot SEPI, GEPI and STEPI (2D ss-SGSTEPI) according to an embodiment of the invention.
FIG. 5 illustrates experimental results for a B1 correction process according to an embodiment of the invention.
FIG. 6 shows plots of the effects of flip angle variation on T1 measurements according to an embodiment of the invention.
FIG. 7 shows plots of the effects of mixing time (TM) variation on T1 measurements according to an embodiment of the invention.
FIG. 8 shows a coil assembly, a coil interface and gradient-echo images according to an embodiment of the invention.
FIG. 9 show illustrates experimental results for T1 mapping of tumor-bearing mice according to an embodiment of the invention.
FIG. 10 shows a forward transfer map and a extracellular fractional map according to an embodiment of the invention.
FIG. 11 shows singleshot T1 and phase-difference mapping on a human brain according to an embodiment of the invention.
FIG. 12 shows plots of signal intensities of 2D ss-SGSTEPI, 2D turbo FLASH and T1 vs. concentration of [Gd-BOPTA].
FIG. 13 shows a diagram of an imaging system according to an embodiment of the invention.



DETAILED DESCRIPTION
FIG. 1 shows a pulse sequence diagram for acquiring a singleshot SEPI and STEPI (2D ss-SESTEPI) according to an embodiment of the invention. The diagram in FIG. 1 shows both RF pulses and data acquisition windows of the MR signal. Gcr and Gsp indicate the gradient pulses for the crusher and the spoiler, respectively. Both spin-echo and stimulated-echo are formed at the same TE, where the effect from the field inhomogeniety is removed and the T2 decay time is the same.
Referring to FIG. 1, a first 90°x RF pulse is applied to excite an initial transverse magnetization. Three reference echoes are acquired immediately after the first 90°x RF pulse for EPI phase-correction. The transverse magnetization excited by the first RF pulse is dephased by the dephasing crusher Gcd, and then split into two quadrature components MI and MQ with equal magnitudes by a perfect tipup 90°−x, RF pulse. The in-phase (I) magnetization MI component is rephased by the rephasing crusher Gcr1, and a spin-echo is formed at position “SEPI.” After SEPI acquisition, the MI magnetization is spoiled by Gsp1. The MI magnetization is given as follows:




  
    
      
        
          
            M
            +
            I
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            1
            2
          
          ⁢
          
            
              
                M
                +
                0
              
              (
              
                r
                ->
              
              )
            
            ·
            
              e
              
                
                  -
                  TE
                
                /
                
                  T
                  2
                
              
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




where {right arrow over (r)} is a position vector, and M+0({right arrow over (r)}) is the initial transverse magnetization at position {right arrow over (r)}. The “+” notation indicates transverse magnetization.

The Quadrature (Q) magnetization component is stored in the longitudinal space during the mixing time TM, during which time the fresh longitudinal magnetization is also recovered. The longitudinal magnetization just before a third RF pulse is the sum of T1 decayed prepared magnetization and freshly recovered magnetization during the mixing time TM. The longitudinal magnetization is given as follows:




  
    
      
        
          
            M
            z
          
          (
          
            
              r
              ->
            
            ,
            TM
          
          )
        
        =
        
          
            
              
                M
                0
              
              (
              
                r
                ->
              
              )
            
            ⁢
            
              (
              
                1
                -
                
                  e
                  
                    
                      -
                      TM
                    
                    /
                    
                      T
                      1
                    
                  
                
              
              )
            
          
          +
          
            
              1
              2
            
            ⁢
            
              
                
                  M
                  +
                  0
                
                (
                
                  r
                  ->
                
                )
              
              ·
              
                e
                
                  
                    -
                    TM
                  
                  /
                  
                    T
                    1
                  
                
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




where Mz({right arrow over (r)}, TM) is the longitudinal magnetization, the first term represents the freshly recovered magnetization and the second term represents the prepared magnetization. The “z” notation indicates longitudinal magnetization.

The third 90° RF pulse tips this longitudinal magnetization down to the transverse plane. Gcr2 rephases only the prepared magnetization previously dephased by Gcd, while spoiling the freshly recovered magnetization as well as all other spin-echoes that are mirrored by the third 90° RF pulse. The stimulated echo is formed at the position “STEPI,” where the center of phase-encoding (ky=0) is also positioned. Therefore, the STEPI signal is only from the purely prepared magnetization, without being contaminated by the fresh magnetization. The prepared magnetization at the position “STEPI” is given by:




  
    
      
        
          
            M
            +
            Q
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            1
            2
          
          ⁢
          
            
              
                M
                +
                0
              
              (
              
                r
                ->
              
              )
            
            ·
            
              e
              
                
                  -
                  TE
                
                /
                
                  T
                  2
                
              
            
            ·
            
              e
              
                
                  -
                  TM
                
                /
                
                  T
                  1
                
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




where M+Q({right arrow over (r)}) is the prepared magnetization, which corresponds to the Q magnetization component after the initial transverse magnetization is split into the I and Q magnetization components.

Thus, the pulse sequence in FIG. 1 acquires a SEPI at position “SEPI” and a STEPI at position “STEPI” in a singleshot following the same RF excitation. The two EPI images may have the same resolution and spatial dimensions. Further, the two EPI images are acquired substantially simultaneously with different T1 decay time. For perfect 90° pulses, half the initial magnetization is stored in the longitudinal space and the other half remains in the transverse plane after the tipup 90° RF pulse. The magnetization remaining in the transverse plane is used for the SEPI and the magnetization stored in the longitudinal space is used for the STEPI. Crusher gradient Gcr2 selectively rephases the prepared magnetization, while spoiling the fresh longitudinal magnetization that recovers during TM.
The pulse sequence after the initial transverse magnetization is spilt into the quadrature components MI and MQ with equal magnitudes by a perfect tipup 90°−x, RF pulse (second RF pulse), may be summarized as follows:
(1) MI component is rephased by rephasing crusher Gcr1, and refocused at the spin-echo position “SEPI.”
(2) MI is spoiled by Gsp1 after SEPI acquisition.
(3) During the mixing time TM, the MQ component undergoes T1 decay and is mixed with the freshly recovered longitudinal magnetization.
(4) The third 90° RF pulse tips the longitudinal magnetization down to the transverse plane, which is the sum of T1 decayed stored magnetization and the freshly recovered magnetization.
(5) Gcr2 rephases only prepared magnetization previously dephased by Gcd, while spoiling the freshly recovered magnetization.
(6) The stimulated echo is formed at the position “STEPI,” where the center of phase-encoding ky is also positioned.
Assuming perfect 90° RF pulses, the magnetizations for the SEPI and the STEPI are given by:




  
    
      
        
          
            M
            SE
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            1
            2
          
          ⁢
          
            
              
                M
                +
                0
              
              (
              
                r
                ->
              
              )
            
            ·
            
              e
              
                
                  -
                  TE
                
                /
                
                  T
                  2
                
              
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
            M
            STE
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            1
            2
          
          ⁢
          
            
              
                M
                +
                0
              
              (
              
                r
                ->
              
              )
            
            ·
            
              e
              
                
                  -
                  TE
                
                /
                
                  T
                  2
                
              
            
            ·
            
              e
              
                
                  -
                  TM
                
                /
                
                  T
                  1
                
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




For the ideal case of perfect 90° RF pulses, the magnetizations for the SEPI and the STEPI only differ by the T1 decay (e−TM/t1).

However, the fligangles of the RF pulses deviate from the perfect 90° RF pulses due to B1 field inhomogeneity. The effect of to B1 field inhomogeneity has to be taken into account to accurately compute T1. The magnetizations for the SEPI and STEPI may be modified as follows to include the effect from spatial variation of the B1 field:




  
    
      
        
          
            M
            SE
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            
              
                f
                I
              
              (
              
                r
                ->
              
              )
            
            ·
            
              1
              2
            
          
          ⁢
          
            
              
                M
                +
                0
              
              (
              
                r
                ->
              
              )
            
            ·
            
              e
              
                
                  -
                  TE
                
                /
                
                  T
                  2
                
              
            
          
        
      
    
    
      
        (
        6
        )
      
    
  
  
    
      
        
          
            
              
                
                  M
                  STE
                
                (
                
                  r
                  ->
                
                )
              
              =
              
              ⁢
              
                
                  
                    
                      f
                      Q
                    
                    (
                    
                      r
                      ->
                    
                    )
                  
                  ·
                  
                    1
                    2
                  
                
                ⁢
                
                  
                    
                      M
                      +
                      0
                    
                    (
                    
                      r
                      ->
                    
                    )
                  
                  ·
                  
                    e
                    
                      
                        -
                        TE
                      
                      /
                      
                        T
                        2
                      
                    
                  
                  ·
                  
                    e
                    
                      
                        -
                        TM
                      
                      /
                      
                        T
                        1
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  g
                  (
                  
                    r
                    ->
                  
                  )
                
                ·
                
                  
                    M
                    STE
                  
                  (
                  
                    r
                    ->
                  
                  )
                
                ·
                
                  e
                  
                    
                      -
                      TM
                    
                    /
                    
                      T
                      1
                    
                  
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




where fI({right arrow over (F)}) and fQ({right arrow over (r)}) are position-dependent fractions for spin-echo and stimulated-echo, respectively, and g({right arrow over (r)}) is defined as g({right arrow over (r)})=fQ(f)/fI({right arrow over (r)}). Because the effect in the flipangle deviation of the first RF is identical for both SEPI and STEPI, the RF functions fI({right arrow over (r)}) and fQ({right arrow over (r)}) are caused by B1 field distribution of the second and third RF pulses. Note that fI({right arrow over (r)}) and fQ({right arrow over (r)}) are smoothly varying and position dependent functions. It is necessary to obtain a map of these B1 functions to use the magnitude images of SEPI and STEPI for quantitative analysis, such as T1 or diffusion measurement.

Assuming the actual flipangle α for the second and third RF pulses are a at position {right arrow over (r)}, Eqs. 5 and 6 for the magnetization components measured at SE and STE are given as




  
    
      
        
          
            M
            SE
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            
              M
              +
              0
            
            (
            
              r
              ->
            
            )
          
          ⁢
          
            sin
            2
          
          ⁢
          
            
              α
              2
            
            ·
            
              e
              
                
                  -
                  TE
                
                /
                
                  T
                  2
                
              
            
          
        
      
    
    
      
        (
        8
        )
      
    
  
  
    
      
        
          
            
              
                
                  M
                  STE
                
                (
                
                  r
                  ->
                
                )
              
              =
              
              ⁢
              
                
                  1
                  2
                
                ⁢
                
                  
                    M
                    +
                    0
                  
                  (
                  
                    r
                    ->
                  
                  )
                
                ⁢
                
                  sin
                  2
                
                ⁢
                
                  α
                  ·
                  
                    e
                    
                      
                        -
                        TE
                      
                      /
                      
                        T
                        2
                      
                    
                  
                  ·
                  
                    e
                    
                      
                        -
                        TM
                      
                      /
                      
                        T
                        1
                      
                    
                  
                
              
            
          
        
        
          
            
              =
              
              ⁢
              
                
                  g
                  (
                  
                    r
                    ->
                  
                  )
                
                ·
                
                  
                    M
                    STE
                  
                  (
                  
                    r
                    ->
                  
                  )
                
                ·
                
                  e
                  
                    
                      -
                      TM
                    
                    /
                    
                      T
                      1
                    
                  
                
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




Here, g({right arrow over (r)}) indicates the error in magnitude from the case with the perfect 90° RF due to B1 inhomogeneity. For perfect 90° RF pulses, Eqs. 8 and 9 become the same as Eqs. 4 and 5 respectively. From these two equations, the theoretical B1 correction function g({right arrow over (r)}) is obtained as follows





  
    
      
        
          g
          (
          
            r
            ->
          
          )
        
        =
        
          
            
              
                sin
                2
              
              ⁢
              β
            
            
              2
              ⁢
              
                
              
              ⁢
              
                sin
                2
              
              ⁢
              
                α
                2
              
            
          
          =
          
            2
            ⁢
            
              
            
            ⁢
            
              cos
              2
            
            ⁢
            
              α
              2
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




The measured spin-relaxation time T1meas of the measurement without B1 correction can be given as follows:





  
    
      
        
          T
          1
          meas
        
        =
        
          
            
              {
              
                
                  
                    ln
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        S
                        SE
                      
                      (
                      
                        r
                        ->
                      
                      )
                    
                  
                  -
                  
                    ln
                    ⁢
                    
                      
                    
                    ⁢
                    
                      
                        S
                        STE
                      
                      (
                      
                        r
                        ->
                      
                      )
                    
                  
                
                TM
              
              }
            
            
              -
              1
            
          
          =
          
            
              {
              
                
                  1
                  
                    T
                    1
                    true
                  
                
                -
                
                  
                    1
                    TM
                  
                  ⁢
                  
                    (
                    
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      
                        cos
                        2
                      
                      ⁢
                      
                        α
                        2
                      
                    
                    )
                  
                
              
              }
            
            
              -
              1
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




where T1meas is the measured T1 and T1true is the true T1, ln is a natural log function and SSE and SSTE are the signal intensities of the SEPI and STEPI, respectively. Eq. 11 indicates that the error in measured T1 will differ for different mixing times as well as different actual flipangle. Eq. 11 is only for the transmit function without including the difference in signal reception.

Eq. 11 was used to simulate the expected error on the T1 estimation with respect to the mixing time and actual flipangle. One second was used as the true T1 in the simulation. FIG. 2 shows a plot of estimated T1 without B1 correction for different flipangles and mixing times TM=0.50, 0.75, 1.0, 1.25 and 1.5 T1 from top to bottom. T1 tends to be prolonged with smaller flipangle than 90° and short mixing time TM, while it is decreased with larger flipangle and mixing time TM. As indicated in Eq. 11, the estimated T1 agrees with the true T1 for perfect 90° flipangle, as indicated by vertical arrow.
It is not necessary to accurately measure the actual position-dependent flipangle α({right arrow over (r)}) in order to correct for B1 field inhomogeniety. Any difference between SEPI and STEPI with when the mixing time TM approaches zero (i.e., no or negligible T1 decay in STEPI) is directly related to the B1 correction function. STEPI with TM approaching zero is denoted STEPI0. Therefore, the correction function g({right arrow over (r)}) can be experimentally measured as a part of the data acquisition by setting the mixing time TM→0. FIG. 3 shows a pulse sequence in which the data acquisition window for SEPI is removed to acquire the STEPIo image with minimal mixing time (TM ̃0). In practice, the minimal TM is about 5 ms and the T1 decay for typical T1 at 3 T is less than 1%. The B1 correction map g({right arrow over (r)}) is then obtained by dividing STEPIo by SEPIo as follows:




  
    
      
        
          g
          (
          
            r
            ->
          
          )
        
        =
        
          
            
              
                f
                Q
              
              (
              
                r
                ->
              
              )
            
            
              
                f
                I
              
              (
              
                r
                ->
              
              )
            
          
          =
          
            
              STEPI
              0
            
            
              SEPI
              0
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




where STEPIo is acquired using the pulse sequence in FIG. 3 with minimum TM and SEPIo may be the same as the SEPI acquired using 2D ss-SESPTEPI. In the pulse sequence shown in FIG. 3, STEPI is acquired after three 90° RF pulses with the minimum mixing time TM to remove the T1 decay. A fourth 90° RF pulse is applied after STEPIo to achieve exactly the same recovery time as that of ss-SESTEPI.

The B1 correction function g({right arrow over (r)}) may be computed for each position {right arrow over (r)} or pixel within the imaging FOV by dividing STEPIo by SEPIo pixel-by-pixel. The inverse correction map g−1 ({right arrow over (r)}) may then be multiplied by the STEPI signal pixel-by-pixel to remove the difference caused by B1 inhomogeneity. This eliminates error in the resultant quantification that is caused by the inhomogeneous B1 field.
In practice, STEPIo is measured only once by pre-imaging with increased signal averaging. The same number of averages may be used for STEPIo and SEPIo for the B1 correction function g({right arrow over (r)}). SEPIo is the average SEPI from pre-dynamic images and is not measured separately. Two to four magnitude averaging may be used for STEPIo and SEPIo. It is worth noting that no change in T1 should be induced during B1 correction mapping acquisition, i.e., 2*Navgs, where Navgs indicates the number of average for B1 correction mapping.
The STEPI corrected for B1 inhomogeneity may then be used to accurately compute a T1 map for the positions {right arrow over (r)} or pixels within the imaging FOV. The T1 map may be determine by computing the logarithm of the signal intensities of the SEPI and the B1 corrected STEPI and performing linear fitting on the logarithms assuming the SEPI and the B1 corrected STEPI differ by the T1 decay (e−TM/T1). The B1 correction may be applied to the SEPI instead. Since the SEPI and STEPI are acquired in a singleshot, the T1 map can be computed in substantially realtime and is free from motion-related artifacts.
The realtime T1 map can be used to monitor concentrations of contrast agents in dynamic contrast enhanced MRI. For example, the concentration of paramagnetic-ion based contrast agents are approximately linearly related to T1 allowing the concentrations of these contrast agents in the patient to be accurately monitored in realtime using the T1 map.
The 2D ss-SESTEPI may also be used to measure the phase-change between SEPI and STEPI. As mentioned previously, both echoes are formed at the same TE and therefore, they are supposed to have the same phase map. The pulse sequence can easily be modified to encode the displacement of an imaging voxel during the mixing time by adjusting the zeroth moment of the crusher gradient. The spatial distribution of this phase may be used to extract the displacement and the average velocity during TM, and change of local magnetic field during the mixing time. Unlike multishot imaging for which imaging phase is mixed with pseudo-random phase error caused by the physiological motion, SEPI and STEPI share the identical phase-error and their difference must be insensitive to the phase-error.
FIG. 4 shows a pulse sequence diagram for acquiring a GEPI in addition to the SEPI and the STEPI in a singleshot according to an embodiment of the invention. The diagram in FIG. 4 is similar to the pulse sequence in FIG. 1 with the addition of GEPI acquisition. In this embodiment, the GEPI is acquired following SEPI acquisition using the transverse magnetization MI. The transverse magnetization is further dephased and forms a gradient echo at the position “GEPI.” The phase difference between GEPI and SEPI can be used to estimate the local field variation, as discussed further below. The SEPI and STEPI acquisitions are the same as before.
After completion of SEPI measurement, all imaging gradients are rewound. The transverse magnetization is further evolved with T2* decay by e−ΔTE/T*2 and the GEPI is measured, as indicated by “GEPI” in FIG. 4. The GEPI provides the information about phase-accumulation during ΔTE caused by local field variation, such as by change of chemical shift or susceptibility field. Therefore, the phase difference between SEPI and GEPI can be used to estimate the variation in Bo field induced by, for instance, temperature-dependent chemical shift.
In this embodiment, a pair of magnitude and phase measurements are acquired for each of the SEPI, GEPI, and STEPI. Thus, three pairs of magnitude and phase images are acquired in a singleshot enabling simultaneous measurement of T1 and phase difference. The SEPI and STEPI have the same TE but different T1 decay, while the SEPI and GEPI have different T2* decay and phase with the same T1 weighting. The signal intensity of the GEPI becomes:




  
    
      
        
          
            S
            GRE
          
          (
          
            r
            ->
          
          )
        
        =
        
          
            
              S
              SE
            
            (
            
              r
              ->
            
            )
          
          ·
          
            e
            
              -
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  TE
                
                
                  T
                  2
                  *
                
              
            
          
        
      
    
    
      
        (
        13
        )
      
    
  




where SGRE({right arrow over (r)}) and SSE({right arrow over (r)}) are the signal intensities of the GEPI and SEPI, respectively.

The phase difference between SEPI and STEPI can be used to compute temperature-dependent chemical shift (e.g., of water proton) as follows:




  
    
      
        
          Δ
          ⁢
          
            
          
          ⁢
          
            σ
            ppm
          
        
        =
        
          
            
              
                θ
                GRE
              
              -
              
                θ
                SE
              
            
            
              2
              ⁢
              
                
              
              ⁢
              
                π
                ·
                Δ
              
              ⁢
              
                
              
              ⁢
              TE
            
          
          ·
          
            1
            
              γ
              ⁢
              
                
              
              ⁢
              
                B
                o
              
            
          
        
      
    
    
      
        (
        14
        )
      
    
  




where Δσppm is the temperature-dependent chemical shift in parts per million, θGRE and θSE are the phases of the GEPI and SEPI, respectively, B0 is the main field of the MRI (e.g., 3.0 T) and γ is the gyromagentic ratio. The term γB0 represents an MRI frequency and the term





  
    
      θ
      GRE
    
    -
    
      θ
      SE
    
  
  
    2
    ⁢
    
      
    
    ⁢
    
      π
      ·
      Δ
    
    ⁢
    
      
    
    ⁢
    TE
  




represents a frequency distribution based on the phase difference between the GEPI and SEPI. Since B0 and γ are known, Eq. 14 can be used to compute the temperature-dependent chemical shift in realtime based on the phase difference between SEPI and GEPI. Eq. 14 can also be used to determine the local variation in the magnetic field due to the chemical shift by multiplying the chemical shift by the field B0. Eq. (14) can be used to compute a chemical shift distribution for the positions {right arrow over (r)} or pixels within the imaging FOV by determining the phase difference between SEPI and GEPI pixel-by-pixel and using Eq. (14) to compute the chemical shift for each position {right arrow over (r)} or pixel.

The computed chemical shift can then be used to determine changes in temperature within a patient. Because the SEPI and GEPI are acquired in a singleshot, the chemical shift and therefore the temperature change can be monitored in realtime. Realtime temperature monitoring may be used in a heat treatment to control the dosage of heat energy during the heat treatment. For example, realtime temperature monitoring can be used in high intensity ultrasound (HIFU) treatments to monitor the temperature of normal tissue surrounding a tumor being ablated to make sure that the temperature of the normal tissue does not exceed a threshold, thereby avoiding heat damage.
Thus, the novel 2D ss-SGSTEPI imaging technique according to embodiments of the invention provides for simultaneous acquisitions of spin-/gradient-/stimulated-EP images (SEPI, GEPI, STEPI) of localized volumes. Unlike the conventional stimulated-echo imaging techniques which discard half the magnetization, all prepared magnetization is measured to produce three EP images in 2D ss-SGSTEPI.
Because the SEPI, GEPI and STEPI are acquired in a singleshot after a single excitation according to embodiments of the invention, they are free from motion-related error. The technique can be used to measure T1 and phase in a singleshot. One immediate application of the singleshot T1 mapping is to determine the concentration of paramagnetic-ion (such as gadolinium) based contrast agent for studying pharmacokinetics in drug-development research and for evaluation of perfusion kinetics in cancer imaging. Using the imaging technique, T1 and phase-difference maps can be achieved almost in realtime and concentration of the paramagnetic-ion based contrast agent can be accurately obtained to evaluate tumor physiology and tumor response to, for example, heat ablation.
MRI Experiments
MRI studies were performed on a Siemens Trio 3 Tesla MRI system (Siemens Medical Solutions, Erlangen, Germany) with Avanto gradients (45 mT/m strength and 150 T/m/s slew rate). The pulse sequence was developed based on a singleshot spin-echo EPI sequence, using the IDEA pulse sequence development environment. Realtime calculation of T1 mapping was implemented into a custom-made ICE program (Image Construction Environment). The image construction program produced two magnitude images, two phase images, and a T1 map for each slice.
A phantom was constructed with four 20 cc plastic bottles filled with 0.05, 0.1, 0.2, and 0.4 mM MnCl2/water solutions, respectively, and immersed in jello in a rectangular acrylic container with dimensions 8×8×13 cm3. The T1 of each concentration was 1.3, 0.9, 0.55, 0.35 sec for 0.05, 0.1, 0.2, 0.4 mM, respectively. T1 of the phantom was measured using saturation-recovery spin-echo.
The MRI experiment was performed to investigate the effect of flipangle and mixing time TM for B1 correction mapping on the accuracy of T1 measurement, using 128×32 imaging matrix with 62.5% phase-encoding asymmetry, 2×2×2 mm3 spatial resolution, 1.39 kHz/pixel receiver bandwidth, and TR=5.0 s, using a transmit/receive head RF coil. The echo times were 15 ms for both SEPI and STEPI.
2D ss-SESTEPI was applied to in-vivo dynamic T1 mapping MRI to measure in realtime the T1 map of two tumor-bearing female nude mice (15), using a home-made quadrature transmit/receive saddle coil with 1.5 inch inner-diameter and 3 inch length and imaging parameters: in-plane resolution of 1×1 mm2 with 2 mm slice thickness, 128×32 imaging matrix, 5.0 s TR, 120 total repetitions (time frames) for 10 min, TEs=15 msec for both SEPI and STEPI. The first 4 measurements were used to acquire STEPIo images for B1 correction mapping. T1 maps were constructed from the time frame 5. In this application for realtime concentration mapping in DMT-MRI, contrast agent was administered after 8 repetitions, i.e., acquisitions of four STEPIo and four pre-contrast imaging. Multi-Hance (Gadobenate dimeglumine: Gd-BOPTA, Bracco Diagnostics, Inc.) was injected via the catheterized tail vein to a dose of 0.1 mM Gd/kg. Animals were cared for under the guidelines of a protocol approved by the University of Utah Institutional Animal Care and Use Committee.
Results of the MRI experiment are discussed below.
B1 Correction Process
A MRI experiment was performed to demonstrate the accuracy of T1 maps with respect to flipangle, mixing time TM, and number of averages for B1 correction map. A phantom was constructed with 4 plastic bottles filled with MnCl2/water solutions with concentrations of 0.05, 0.1, 0.2, and 0.4 mM, which were surrounded by jello. T1 mapping experiment was conducted using inversion-recovery spin-echo (IR-SE) imaging and 2D ss-SGSTEPI. T1 of each concentration of MnCl2/water solution were measured as 1.5, 0.89, 0.55, 0.35 sec for 0.05, 0.1, 0.2, 0.5 mM, respectively. FIG. 5 shows the T1 calculation process for 4-averaged magnitude images (SEPI and STEPI0) in (a) and (b), respectively, for the B1 correction map. FIG. 5 also shows the B1 correction map g({right arrow over (r)}) in (c), and single average singleshot images (SEPI) and (STEPI) in (d) and (e), respectively. FIG. 5 also shows the T1 maps before and after B1 correction in (f) and (g), respectively.
The T1 map (g) after B1 correction indicates remarkable improvement in uniformity compared to that before correction in (f). T1 values measured using 2D ss-SGSTEPI agree very well with those measured using IR-SE, of which T1 for 0.1 and 0.2 mM are compared in FIG. 6(a). T1 values along the horizontal dotted lines on T1 maps before and after corrections are plotted in (h) and (i), respectively, of FIG. 5. The T1 map (g) and plot (i) clearly indicate substantial improvement by B1 inhomogeneity correction.
Referring to FIG. 5, the B1 correction process is summarized as follows: four averaged (a) SEPI and (b) STEPIo images were used to calculate (c) the B1 correction map, which is used to correct the singleshot (d) SEPI and (e) STEPI. The T1 maps before and after B1 correction are shown in (f) and (g), respectively. T1 values of jello along the horizontal dotted lines in the water proton are presented in plots (h) and (i) for before and after B1 correction, respectively. The T1 map (g) and plot (i) clearly show substantial improvement as a result of B1 inhomogeneity correction.
Flip Angle Variation
2D ss-SESTEPI was applied to investigate the effects of flip angle variation on the resultant T1 value. A T1 map was acquired using 2D ss-SESTEPI with nominal flip angles 60, 75, 90, 105, and 120° for the second and third RF pulses in FIG. 1, using a 0.4 sec mixing time, and four-averaged B1 correction map.
Mean T1 values on regions-of-interest (ROIs) of 0.1 and 0.2 mM solutions are plotted in (a) and (b) of FIG. 6. Closed and open symbols represent T1 measurements before and after B1 correction, respectively and circles and square symbols indicate 0.1 and 0.2 mM MnCl2 concentrations, respectively. Dotted and broken lines are guidelines to the T1 plots without B1 correction. Solid horizontal lines indicate T1 values measured using SR-SE imaging. These experiments demonstrated significant improvement in accuracy by B1 correction for all values of T1 (0.2 ̃1.3 sec.). Plots in (b) are T1 values in (a) normalized to T1=1.0 s to compare the empirical data to the simulated plot in FIG. 3. A TM of 0.4 s corresponds to 0.45 T1 for T1=0.89 s and 0.74 T1 for T1=0.55 s for 0.1 and 0.2 mM MnCl2 concentrations, respectively.
Plots in (b) of FIG. 6 show good agreement with the theoretical simulation in FIG. 2 and robustness of the B1 inhomogeneity correction. In the B1 correction, no modeling or user-defined input is used. Therefore, the resultant measurement is free from any user-dependent bias. SNR is the key element to determine the accuracy of the measurement. Considering that the measurement was made in a singleshot scan with a few seconds acquisition time, the accuracy of the T1 measurement may be considered remarkable.
Referring to FIG. 6, the results for flip angle variation are summarized as follows: T1 measurement before and after B1 correction are indicted by (●,▪) and (◯, □), respectively. The data in (a) are actually measured T1 values, and those in (b) were normalized T1 to 1.0 s to compare with the theoretical plots in FIG. 2. Dotted and broken lines indicate guidelines. Horizontal solid-lines represent T1 using spin-echo MRI.
TM Variation
To accurately measure T1 decay using two signal intensities, the signal difference between SEPI and STEPI must be significantly larger than the noise amplitude, otherwise the resultant measurement will be dominated by noise, rather than T1 decayed signal change. The T1 map was acquired using 2D ss-SESTEPI with a nominal flip angle of 90° and a four-averaged B1 correction map for variable mixing times TM=0.1, 0.2, 0.4, and 0.8 sec. The mean and standard deviation were measured on ROIs selected on each concentration and plotted with respect to the mixing time in FIG. 7. Error bars on the plots indicate the standard deviation of each measurement.
These plots indicate that the uncertainty of the measurements greatly depends on the selection of TM. T1 measurement using the increased TM produces increased accuracy and reduced errors; however, it limits the number of slices in an interleaved multi-slice scan. For T1 below 1.0 sec, the difference in uncertainty for measurements with 0.4 sec and 0.8 sec TM is not significant.
In-Vivo Measurement of Singleshot T1 of a Mouse
To demonstrate the performance of 2D ss-SGSTEPI in in-vivo application, the technique was used to acquire, (1) dynamic T1 mapping in DCE-MRI on two tumor-bearing mice, and (2) rapid T1 mapping of human brain. To increase the SNR, a dedicated transmit/receive mouse RF coil was developed and used for mice MRI experiment.
Because typical clinical RF coil is over-coupled to accommodate the wide range of imaging subject, SNR of the resultant images are relatively low. A transmit/receive quadrature volume RF coil with interface (T/R switch, preamplifier, quadrature hybrid mixer/splitter, and coil selector switch) was developed for imaging of mice at 3T wholebody clinical MRI system. The dimension of this coil was 1.5′′ inner diameter and 3′′ length to image mice and small rats. FIG. 8(a) shows the coil assembly and 8(b) shows the coil-interface.
Using this RF probe, SNR was improved 3 times for the same imaging protocol, compared to that using Siemens wrist RF coil. This corresponds to 9 times efficiency in imaging time. Images (c) and (d) in FIG. 8 are 3D gradient-echo images of pre- and post-contrast, respectively, after Gd-chelate injection with (0.4 mm)3 isotropic resolution acquired in 1.5 min imaging time.
In-Vivo Dynamic T1 mapping (DTM) of Tumor-Bearing Mice
2D ss-SGSTEPI was applied to measure dynamic T1 maps from two in-vivo tumor-bearing mice, using the coil in FIG. 8 and imaging parameters of in-plane resolution of 1×1 mm2 with 2 mm silce thickness, 128×32 imaging matrix, 10 slices, TR 5.0 s, TEs=15, 40, and 15 msec for SEPI, GEPI, and STEPI. The first four measurements and next four pre-dynamic images were used for B1 correction mapping. MultiHance (Gd-BOPTA, Bracco Diagnostics, Inc.) was injected via the catheterized tail vein with the dose of 0.1 mM Gd/kg.
Images in FIG. 9 are (a) T1 and (b) ΔR1 maps of selected time points (every 30 seconds) from 120 total time frames acquired in 10 min. Note that 10 slices were obtained every 5 seconds in this DTM experiment. Images in Figure (c) are T1 weighted spin-echo images of the same location as T1 and ΔR1 maps in (a, b) for pre-contrast and 10 min after the injection. The brightness in (a) and (b) represents the larger values in T1 and ΔR1 maps.
ΔR1 was plotted for selected ROIs with 9 pixels and all time points of 10 minutes, as numeric labels in FIG. 9(e). Regions of 3×3 pixels (18 μL=9×0.1×0.1×0.2 cm3 volume) were selected at 1/kidney, 2/tumor periphery, 3/tumor core, and 4/background noise, and mean ΔR1 was plotted in FIG. 9(d) over the entire duration (10 minutes). Unlike other DCE-MRI method that uses the signal intensity as the concentration, the values ΔR1({right arrow over (r)}) are directly proportional to the PBCA concentration. Therefore, these time-course ΔR1 maps can be used for further evaluation of pharmacokinetics in pixel-by-pixel or organ-based.
ΔR1 maps were used as the concentration maps to calculate pharmacokinetics parametric maps. Two compartment model and standard arterial input function were used to extract the forward transfer rate (ktrans) and the extracellular fractional volume maps, as shown in FIG. 10. FIG. 10(a) shows the ktrans map and FIG. 10(b) shows the extracellular fractional volume map. These maps reflect the capillary permeability and the volume of the extracellular fluid.
Singleshot Measurements of T1 of Human Brain
To demonstrate clinical applicability, 2D ss-SGSTEPI was applied to a human brain with TR 5 sec, 10 slices with 2×2×2 mm3 spatial resolution, 2-average B1 correction mapping, and TM 0.4 sec. After B1 correction mapping acquisition, T1 and phase-difference maps were obtained in every 5 sec. The phase-difference image was also constructed on-line by simply subtracting the phase image of SEPI from that of GEPI. FIG. 11 shows singleshot T1 and phase-difference mappings on a human brain with TR 5 s, 2×2×2 mm3 resolution. These maps were constructed on-line at the MRI system.
Applications for Estimating Concentration of Paramagnetic-Ion Based Contrast Agents
Pharmacokinetic studies for drug development, and for evaluation of tumor perfusion require accurate estimation of the concentration of pharmaceuticals. Paramagnetic-ion based chelates, such as Gd-DTPA, Gd-DOTA, and Mn-EDTA, are MRI contrast agents (PBCA) commonly used to model pharmaceuticals, because they can be monitored approximately quantitatively by MR. Temporal changes in the MRI signal intensity of the water proton is then interpreted to extract kinetic parameters. Upon the introduction of a paramagnetic contrast agent, such as Gd-DTPA, the relaxation times (T1, T2, and T2*) of the water proton change depending on the local concentration of the contrast agent. Therefore, dynamic changes of the local contrast agent concentration may be obtained by measuring either T1 or T2.
Dynamic contrast enhanced MRI (DCE-MRI) has shown promise in detection and characterization of tumors, by providing information about important perfusion parameters such as vessel permeability, microvessel density, and washout, which may allow more accurate characterization of benign versus malignant neoplasms.
Physiologic variables that characterize tissue microcirculation (capillary permeability and volume of the extracellular fluid) are altered in pathological states. Dynamic contrast enhanced MRI (DCE-MRI) is frequently performed for the serial acquisition of MR images of interest to evaluate these parameters. By fitting the concentration-time curve to an appropriate pharmacokinetic model function, these physiologic parameters can be extracted. The accuracy of these parameters depends on the accuracy of the concentration mapping. Therefore, there are continuous demands to develop new and improved methods for precise measurement of kinetic parameters. For instance, conventional MRI of breast provides less specific information about the tumor than DCE-MRI and diffusion-weighted imaging (DWI), which can be used to make quantitative measurements of tissue properties that are relevant to the tumor progression and response to therapy.
Drug delivery research has revealed a clear need for an imaging method to improve quantitative drug measurements. Because of the rapid change in signal intensity and the time constraint of the available imaging techniques after injection of the Gd-Based contrast agent, a fast gradient echo (FGRE or turboFLASH) sequence has been used to dynamically monitor fast signal changes or T1 weighted 2D spin-echo to dynamically monitor slow signal changes. Then, the change in signal intensity has been used to provide a semi-quantitative estimation of the drug quantity. However, the signal intensity of FGRE imaging is weighted by both T1 and T2 (or T2*). This T2* shortening induces an error in the estimation of concentration. To improve the accuracy of drug concentration measurement, T2 or T2* effects must be removed from the GRE signal in DCE-MRI.
Increasing the concentration of the PBCA decreases both T1 and T2 relaxation times, and leads to an increase in the T1WI signal intensity but a decrease in the T2WI signal intensity. T2* is strongly influenced by external factors, such as the location of the imaging FOV relative to the magnet center, and is much shorter than T2. Particularly, T2* shortening happens with high local concentrations of PBCA, as in the tumor periphery. Therefore, the signal intensity does not directly convert to the concentration of the PBCA, although it is commonly taken to do so in DCE-MRI assuming a linear relationship at low concentrations, i.e., changes in the relaxation rate ΔR1(R1=1/T1) is linearly related to the local concentration of Gd(III)-based substance over a range of low concentrations in the solution phantom, and can be used to estimate Gd-DTPA concentration. Pharmacokinetic information can be obtained from the time-variation of ΔR1(r).
To demonstrate the non-linearity of SI vs. concentration, an MRI experiment was performed to measure T1 and signal intensity of turbo-FLASH on a phantom with different concentrations of Gd-BOPTA (MultiHance, Bracco Diagnostics Inc.) using TR/TE=7.2/2.7 ms and flipangle 30°. The results are plotted in FIG. 12. FIG. 12(a) indicates that signal intensities of both turboFLASH (temporal resolution of 6 sec) and spin-EPI are non-linear even at low concentration. The dosage 0.1 ̃0.2 mM [PBCA]/kg is the routine dosage for DCE-MRI and local concentrations may easily exceed 0.2 ̃0.4 mM/kg in organs with high-blood supply or at the periphery of a tumor. The signal behavior of the gradient-echo MRI varies widely with the selection of TR, TE, and flipangle, and for T1. Although a set of these imaging parameters may be linear for a small range of PBCA concentrations and specific T1, it will not be linear for other T1 values. In contrast, the T1 relaxation rate (1/T1) of the water proton was measured to be very linear with the concentration, as shown in FIG. 12(b). Thus, direct conversion of SI to the concentration is subject to large error, which will be propagated into the subsequent kinetic calculation. In contrast, the realtime T1 mapping provided by embodiments of the invention described above is linearly related to the PBCA concentration. Thus, the realtime T1 mapping provides accurate realtime estimation of the PBCA concentration resulting in improvement of the accuracy of drug-delivery or tumor evaluation.
A spin-echo sequence with either inversion-recovery (IR) or saturation-recovery (SR) is widely used to measure T1. Application of these methods is limited to non-dynamic imaging because of their long imaging times of a few tens of minutes. Several alternative methods for rapid T1 mapping have been presented including 2D singleshot-EPI (2D ss-EPI) and its variants (Look-Locker EPI, IR-EPI) which achieve better temporal resolution. However, the spatial resolution using these EPI-based methods is greatly limited by susceptibility artifacts. A multishot spin-echo EPI (ms-SEPI) sequence may be used to trade off imaging speed with the severity of the susceptibility artifact. Gradient-echo imaging with flip angle variation has been used with reasonable temporal and spatial resolution. However, the T1 estimate obtained by this method is very sensitive to B1 inhomogeneity and subject to large error unless B1 inhomogeneity is corrected. All these techniques complete the data acquisition in multiple imagings for multiple recovery times and are subject to mis-registration or phase-error due to the subject's physiologic or voluntary motion. The T1 mapping according to embodiments of the invention overcome these problems by calculating T1 in a singleshot using SEPI and STEPI acquired in a singleshot acquisition using 2D ss-SESTEPI.
Reducing Geometric Distortion in EPI
Advantages of singleshot imaging using 2D singleshot EPI (2D ss-EPI) include rapid imaging, lack of motion-artifact, and high signal-to-noise ratio (SNR) per imaging time. These are crucial for quantitative MRI applications, such as diffusion MRI and phase measurement. However, 2D ss-EPI is very sensitive to the local variation of non-linear background gradient created by magnetic susceptibility differences at/near the tissue/air and tissue/bone interfaces. The resultant images suffer greatly from large geometric distortion that often renders the images useless. This distortion also limits the spatial resolution of EPI. Thus, its application is limited mostly to diffusion and functional MRI of the intracranial human brain, far from sinus and temporal bone, for which there is no alternative technique.
The degree of geometric distortion in 2D ss-EPI is inversely proportional to the increment of phase-encoding Δkytravel. Δkytravel can be increased by reducing the imaging FOV, however, reduced FOV imaging within a large volume induces aliasing in the phase-encoding direction. Parallel imaging is commonly used to remove this aliasing artifact. It requires a multi-channel phased-array RF coil and even then, residual aliasing may greatly reduce the accuracy of the quantitative MRI. Inner-volume-imaging or reduced FOV eliminate the aliasing artifact by applying two RF pulses with slice selection gradients applied perpendicular to each other. Spatial saturation is often used to eliminate the signal contribution from spins outside the reduced FOV. The success of reduced FOV imaging using spatial saturation relies greatly on the spatial homogeneity of the RF field, which is never perfect. ZOOM-EPI was the first interleaved multi inner slice imaging technique. However, it is very time-inefficient because the technique measures images of the adjacent slices in different acquisitions. 2D spectro-spatial excitation was used to limit the RF excitation both in slice and phase-encoding directions. Generally, the spatial profile of the coverage is very poor using the 2D excitation RF pulse.
A more efficient interleaved multi-inner-volume (IMIV) imaging technique may be used, which is disclosed in Jeong E K, Kim S E, Kholmovski E G, Parker D L “High-resolution DTI of a localized volume using 3D single-shot diffusion-weighted STimulated echo-planar imaging (3D ss-DWSTEPI)” Magn Reson Med 2006; 56(6):1173-1181, and in Jeong E K, Kim S E, Guo J, Kholmovski E G, Parker D L “High-resolution DTI with 2D interleaved multislice reduced FOV single-shot diffusion-weighted EPI (2D ss-rFOV-DWEPI)” Magn Reson Med 2005; 54(6):1575-1579, both of which are incorporated herein by reference. Interleaved multi-inner volume/slice imaging (IMIV) can be implemented in 2D ss-SESTEPI by applying a pair of adiabatic 180° pulses sandwiched by crusher gradient pulses with a slice-selection gradient along the phase-encoding direction to spoil transverse magnetization outside the phase-FOV. Thus, the 2D ss-SGSTEPI according to embodiments of the invention can use interleaved multi inner volume/slice (IMIV) imaging capability to reduce the geometric distortion in the EPI images and image a small region within a large surrounding without aliasing
FIG. 13 shows an exemplary imaging system 1300 that can be used to perform the 2D ss-SESTEPI and SGSTEPI according to embodiments of the invention. The imaging system 1300 includes an MRI machine 1310, a processor 1320 and a display console 1330. The MRI machine may comprise, for example, a Siemans Trio 3 Telsa MRI system, which was used in the MRI experiments described above. The processor 1320 may instruct the MRI machine 1310 to acquire SEPI/GEPI/STEPI in a singleshot, for example, by instructing the MRI machine 1310 to perform the pulses sequences in FIGS. 1, 3 and/or 4. The processor 1320 may then receive the resulting SEPI/GEPI/STEPI data from the MRI machine 1310 and process the data, for example using the techniques described above. For example, the processor 1320 may compute a T1 map based on the SEPI and STEPI data according to embodiments of the invention. The processor 1320 may also compute PBCA concentrations in a subject based on the computed T1 map. In another example, the processor 1320 may use the GEPI and STEPI to compute change in temperature based on the phase difference between GEPI and STEPI according to embodiments of the invention. The processor 1320 may output the processed data to a user (physician) as an image on the display console 1330 or a report, such as a T1 map or report for generation or review by a radiologist or other health care professional. The SEPI, STEPI and GEPI images may or may not be displayed to the user. For example, the processor 1320 may use the signal intensities of the SEPI and STEPI to compute the T1 map and display the T1 map on the display console 1330 without displaying the SEPI or STEPI.
Those skilled in the art will appreciate that the systems and methods described in connection with the embodiments disclosed herein may be implemented or performed with general purpose processors, a digital signal processors (DSPs), application specific integrated circuit (ASICs), field programmable gate array (FPGAs), logic programmable devices, discrete hardware components, or any combination thereof. Further, the algorithms performed by the beamformers and controllers may be embodied in software executed by a processor, directly in hardware, or a combination thereof. The software may be stored in memory or any computer-readable storage medium known in the art.
The memories described in the disclosure may comprise machine readable media. Readable media may include storage integrated into a processor, such as might be the case with an ASIC, and/or storage external to a processor. By way of illustration, and not limitation, readable media may include one or more of volatile memory, nonvolatile memory, a Random Access Memory (RAM), a flash memory, a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable PROM (EPROM), a register, a hard disk, a removable disk, a CD-ROM, a DVD, or any other suitable storage device. In addition, readable media may include a transmission line or a carrier wave that encodes a data signal. A readable medium may be a computer-readable medium encoded or stored with a computer program or instructions. The computer program or instructions may be executable by a transmitter or receiver device or by a processing system of a transmitter or receiver device.
Those of skill in the art would appreciate that the various illustrative blocks, modules, elements, components, methods, and algorithms described herein may be implemented as electronic hardware, computer software, or combinations of both. Furthermore, these may be partitioned differently than what is described. To illustrate this interchangeability of hardware and software, various illustrative blocks, modules, elements, components, methods, and algorithms have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application.
It is understood that the specific order or hierarchy of steps or blocks in the processes disclosed is an illustration of exemplary approaches. Based upon design preferences, it is understood that the specific order or hierarchy of steps or blocks in the processes may be rearranged. The accompanying method claims present elements of the various steps in a sample order, and are not meant to be limited to the specific order or hierarchy presented.
The previous description is provided to enable any person skilled in the art to practice the various aspects described herein. Various modifications to these aspects will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other aspects. Thus, the claims are not intended to be limited to the aspects shown herein, but is to be accorded the full scope consistent with the language claims, wherein reference to an element in the singular is not intended to mean “one and only one” unless specifically so stated, but rather “one or more.” All structural and functional equivalents to the elements of the various aspects described throughout this invention that are known or later come to be known to those of ordinary skill in the art are expressly incorporated herein by reference and are intended to be encompassed by the claims. Moreover, nothing disclosed herein is intended to be dedicated to the public regardless of whether such invention is explicitly recited in the claims. No claim element is to be construed under the provisions of 35 U.S.C. §112, sixth paragraph, unless the element is expressly recited using the phrase “means for” or, in the case of a method claim, the element is recited using the phrase “step for.”

## Claims

1. A method of acquiring both magnetic resonance spin echo planar image (SEPI) data and stimulated echo planar image (STEPI) data, in response to a single radio frequency (RF) excitation pulse sequence, the method comprising the activities of:
exciting, with a first radio frequency (RF) pulse, nuclei to produce an initial transverse magnetization in a region of a patient;
splitting, with a second RF pulse, the initial transverse magnetization into a first magnetization component and a second magnetization component, said first magnetization component and said second magnetization component comprising quadrature components;
acquiring, with an RF detector, spin echo planar image (SEPI) data of an image slice of said region using the first magnetization component; and
acquiring, with said RF detector, stimulated echo planar image (STEPI) data of said image slice of said region using the second magnetization component.

2. The method of claim 1, wherein the SEPI data and STEPI data are of approximately a same image slice volume of the patient.

3. The method of claim 1, further comprising outputting an output data set to a display, based on the acquired SEPI data of said image slice and the acquired STEPI data of said image slice, the output data set comprising at least one of a magnetic resonance image and a report.

4. The method of claim 1, further comprising:
after the transverse second magnetization has returned to a longitudinal plane, tipping, with an applied third RF pulse, the second magnetization into the transverse plane, the transverse plane being transverse to the longitudinal plane; and
acquiring the STEPI data with the second magnetization in the transverse plane, in response to the applied third RF pulse.

5. The method of claim 1, wherein the first RF pulse is a 90°x pulse and the second RF pulse is a 90°−x pulse.

6. The method of claim 1, further comprising:
dephasing the initial transverse magnetization using a dephasing crusher gradient pulse; and
rephasing the first magnetization using a rephasing crusher gradient pulse.

7. The method of claim 1, further comprising generating a T1 map of the slice of the patient using the acquired STEPI data of said image slice and the acquired SEPI data of said image slice.

8. The method of claim 7, wherein generating the T1 map comprises:
computing logarithms of signal intensities of the acquired STEPI data of said image slice and the acquired SEPI data of said image slice; and
performing a linear fitting on the logarithms of the signal intensities of the acquired STEPI data of said image slice and acquired SEPI data of said image slice.

9. The method of claim 7, wherein generating the T1 map comprises:
generating a correction map that corrects for B1 field inhomogeneity; and
modifying the acquired STEPI data of said image slice using the correction map, wherein the T1 map is generated using the modified acquired STEPI data of said image slice and the acquired SEPI data of said image slice.

10. The method of claim 9, wherein generating the correction map comprises:
after the transverse second magnetization has returned to a longitudinal plane, tipping, with an applied third RF pulse, the second magnetization into the transverse plane, the transverse plane being transverse to the longitudinal plane; and
acquiring the STEPI data with the second magnetization in the transverse plane, in response to the applied third RF pulse;
repeating the first, second and third RF pulses; and
acquiring a second STEPI dataset of said image slice with a minimum mixing time between the repeated second PF pulse and the repeated third RF pulse, wherein the correction map is generated using the second STEPI dataset of said image slice and the acquired SEPI data of said image slice.

11. The method of claim 10, wherein the minimum mixing time is substantially 10 milliseconds or less.

12. The method of claim 10, wherein generating the correction map comprises dividing the second STEPI dataset of said image slice by the initially acquired SEPI data of said image slice.

13. The method of claim 7, further comprising estimating concentrations levels of a contrast agent in the patient using the generated T1 map.

14. The method of claim 13, wherein the contrast agent is a paramagnetic-ion based contrast agent.

15. The method of claim 1, further comprising acquiring gradient EPI (GEPI) data of said image slice using the first magnetization and determining a phase difference between the acquired SEPI and the GEPI data of said image slice.

16. The method of claim 15, further comprising determining a temperature-dependent chemical shift based on phase difference between the acquired SEPI data of said image slice and the acquired GEPI data of said image slice.

17. The method of claim 15, further comprising determining a temperature change in the patient during a heat treatment procedure based temperature-dependent chemical shift from the acquired magnetic resonance imaging SEPT, SEPI or GEDI data.

18. The method of claim 16, wherein determining the temperature-dependent chemical shift comprises:
dividing the phase difference between the acquired SEPI data and the acquired GEPI data by a time between acquisitions of the acquired SEPI data and the acquired GEPI data in order to obtain a frequency distribution; and
dividing the frequency distribution by a magnetic resonance frequency of an imaged substance.

19. An apparatus configured for processing magnetic resonance spin echo planar image (SEPI) data and stimulated echo planar image (STEPI) data, in response to a single radio frequency (RF) excitation pulse sequence, the apparatus comprising:
a processing module that receives (i) SEPI data, acquired with an RF detector from of an image slice of a region of a patient, in response to a first magnetization component, and (ii) STEPI data, acquired with said RF detector from of said image slice of said region of a patient, in response to a second magnetization component;
wherein the first and second magnetizations result from exciting, with a first RF pulse, nuclei to produce an initial transverse magnetization in said region of said patient, and then splitting, with a second RF pulse, the initial transverse magnetization into said first magnetization component and said second magnetization component, said first magnetization component and said second magnetization component being quadrature components;
wherein the processing module processes the acquired SEPI data and the acquired STEPI data into an output data set provided to a display, the output data set configured to produce at least one of a magnetic resonance image of said patient and a report.

20. The apparatus of claim 19, wherein the processing module sends instructions to a magnetic resonance imaging apparatus in order to perform a single RF pulse sequence for acquiring the SEPI and STEPI data.