US7310548B2 - Method of magnetic resonance perfusion imaging - Google Patents
Method of magnetic resonance perfusion imaging Download PDFInfo
- Publication number
- US7310548B2 US7310548B2 US10/551,068 US55106805A US7310548B2 US 7310548 B2 US7310548 B2 US 7310548B2 US 55106805 A US55106805 A US 55106805A US 7310548 B2 US7310548 B2 US 7310548B2
- Authority
- US
- United States
- Prior art keywords
- magnetic resonance
- perfusion
- resonance data
- sensitivity
- tensor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related, expires
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56341—Diffusion imaging
Definitions
- the present invention relates to the field of magnetic resonance imaging, and more particularly to perfusion tensor imaging.
- Diffusion weighted magnetic resonance imaging is a known prior art technique.
- diffusion weighted magnetic resonance imaging the diffusion tensor is obtained from, the magnetic resonance measurement signal for a defined body region of interest.
- the diffusion tensor is then visualised by known imaging methods (U.S. Pat. No. 5,539,310; La Bilhan et al, Diffusion Tensor Imaging: Concepts and Applications, Journal of Medical Resonance Imaging 13:534-46 (2001); Lawrence R. Frank, Anisotropy in High Angular Resolution Diffusion-Weighted MRI, Magnetic Resonance in Medicine 45:935-939 (2001),).
- IVIM intravoxel incoherent motion
- Classical perfusion is a measure of the blood delivered to and used by a specified mass of tissue. It is often measured using spin labelling techniques in MRI.
- IVIM measures quasi-random blood movement within a single imaging voxel and results in a bi-exponential signal attenuation in a standard pulsed gradient spin echo (PGSE) experiment (cf. R. J. Moore et al, in vivo intravoxel incoherent motion measurements in the human placenta using echo-planar imaging at 0.5 T, Magnetic Resonance in Medicine, 43;295-302 (2000)).
- PGSE pulsed gradient spin echo
- the present invention aims to provide an improved imaging method which enables imaging and visualisation of the perfusion anisotropy as well as a corresponding computer program product and perfusion imaging apparatus.
- the present invention provides for a method of perfusion imaging which can be used for medical purposes such as abdominal imaging.
- the present invention is based on the discovery that an isotropic perfusion can be described by a tensor which is structurally equivalent to a diffusion tensor used for diffusion tensor imaging (DTI).
- DTI diffusion tensor imaging
- the slope values form the basis for calculating the perfusion tensor, especially its eigenvectors. This is based on the discovery, that for lower b sensitivity values the signal decay is mostly due to perfusion and not diffusion effects. This is why diffusion can be neglected at low b sensitivity values in order to extract the information out of the decay signal which forms the basis for determining the shape of the perfusion tensor. At least, the sorting of the eigenvalues and the direction of the eigenvectors of the perfusion tensor does not depend much on the diffusion-related signal decay. This directional information derived from the perfusion tensor can be visualized. For this purpose software algorithms which are used in the prior art for diffusion visualization can be reused due to the same mathematics governing the visualization of diffusion and perfusion tensors.
- two additional magnetic resonance data acquisitions are performed for higher b sensitivity values.
- one additional magnetic resonance data acquisition is performed for an intermediate b sensitivity value of around 200 and a further magnetic resonance data acquisition is performed at a high b sensitivity value of around 800.
- the decay signal is mostly governed by its diffusion signal component, see e.g. [Petra Mürtz, et al. Abdomen: Diffusion-weighted MR Imaging with Pulse-triggered Single-Shot Sequences, 258-264, Radiology July 2002].
- the diffusion is isotropic, as demonstrated in the cited reference.
- the diffusion coefficient can be obtained from just two higher b-values measurements.
- For brain imaging it is preferred to perform at least six data acquisitions with gradient encodings for random motion in different directions as diffusion is anisotropic in this case.
- only one additional magnetic resonance data acquisition is performed for a high b sensitivity value, leaving out the intermediary measurement.
- one of the magnetic resonance data acquisitions which have been performed for the low b sensitivity values having the highest signal value is selected as a replacement for the intermediate measurement.
- FIG. 1 is illustrative of a logarithmic diagram showing magnetic resonance imaging signals which form the basis for perfusion tensor determination and imaging
- FIG. 2 is illustrative of a flow chart for determining the perfusion tensor based on the signals shown in FIG. 1 ,
- FIG. 3 is an enlarged view of the decay curves of FIG. 1 and is illustrative of a further preferred embodiment of the invention.
- FIG. 4 is a block diagram of an imaging system.
- IVIM intra-voxel incoherent motion imaging
- a 1 +A 2 1.
- a 1 has been designated to represent f, the fraction of flowing material (i.e. blood) in the voxel.
- f denotes the blood fraction
- a 2 1 ⁇ f, which in other words means that the non-flowing part of the signal contributes to the diffusion signal.
- the fraction f being representative of the blood content, is an isotropic quantity for each voxel. This knowledge is used to constrain the data analysis in the rest of this section.
- Equation 2 The first term of equation 2 will be referred to as the perfusion signal component and the second term in equation 2 will be referred to as the diffusion signal component.
- the perfusion signal component is decaying much more rapidly than the diffusion signal component and thus primarily determines the slope of the signal for low b values.
- the perfusion is described by a perfusion tensor, which has the same structure as the diffusion tensor, i.e. the perfusion tensor is of second rank and symmetric.
- the perfusion tensor can be determined on the basis of the bi-exponential signal of equation 2, using the signal decays for low b values whereby the diffusion signal component is neglected. This is illustrated by means of FIG. 1 .
- the magnetic resonance data acquisition for point A is performed without gradient encoding and for the points B 1 , B 2 , . . . B 6 at six different gradient directions, respectively.
- the slopes m 1 , m 2 , . . . m 6 of the signal decay curves for low b-values are determined as illustrated in FIG. 1 by linear approximation.
- the perfusion tensor can be calculated by means of the same mathematics as used for diffusion tensor calculation (c.f. for example M. T. Vlaadringerbroek and J. A. den Boer. Magnetic Resonance Imaging. Springer-Verlag Berlin Heidelberg New York, 1999, in particular section 7.7).
- FIG. 2 is illustrative of the corresponding flow chart.
- magnetic resonance data acquisition is performed for B i for a low b value with gradient encoding into a first direction.
- the i is incremented and step 202 is performed again for the next gradient encoding direction.
- Step 202 is repeated at least 6 times for magnetic resonance data acquisitions into at least six different gradient encoding directions.
- step 206 the slope values m i are obtained from the measurements points A and B i , B 2 , . . . B 6 by linear approximation, i.e. from the slopes of the straight lines connecting A to B 1 , A to B 2 . . . A to B 6 .
- These slope values provide full mapping of the perfusion tensor and are input into the perfusion tensor calculation routine of step 208 .
- FIG. 3 shows the signal decay curves 100 and 102 of FIG. 1 for higher b-values. As apparent from FIG. 3 the signal decay curves become approximately a straight line for higher b values.
- the slope of the straight line between points X and Y is the diffusion coefficient D.
- the perfusion signal component is isolated from the measured data values; as f is also known the perfusion coefficient P can be obtained this way.
- this enables a pre-processing of the signal decay curves 100 , 102 , 104 , . . . of FIG. 1 by subtracting the diffusion signal component before the determination of the slopes m i .
- the perfusion coefficient P is calculated as explained above for every diffusion-sensitisation gradient direction.
- the resulting set of at least six perfusion coefficients is used to calculate the perfusion tensor map.
- the perfusion tensor's eigenvalues, eigenvectors, and rotationally invariant quantities like the trace or the fractional anisotropy are derived and visualized.
- voxels with substantially co-linear main eigenvectors can be connected (starting from a user defined seed position), and the connecting line can be displayed (tractography).
- the highest value of the B 1 , B 2 , . . . B 6 measurements for low b values is B (cf. FIG. 1 ).
- a straight line between Y and B 1 can be used instead of a straight line between X and Y.
- each acquisition prolongs the breath hold duration by typically 1 to 3 seconds.
- Completely avoiding higher b-value acquisitions enables the use of short echo times, because such high b-value acquisitions require much gradient area, thus having a relatively long echo time (typically 60-90 ms, whereas T 2 relaxation in the abdomen is only 100 ms).
- the long echo time relative to T 2 relaxation significantly lowers the SNR of the measurements.
- Using only small b values ( ⁇ 50) significantly reduces the required gradient area, and thus enables lower echo times (typically 20-30 ms). This significantly improves the SNR.
- FIG. 4 shows a block diagram of perfusion imaging system 400 .
- Perfusion imaging system 400 has magnetic resonance data acquisition device 402 which provides magnetic resonance data to computer system 404 .
- the magnetic resonance data acquisition is performed by means of a sequence of single-shot echo-planner MR imaging sequences.
- control program 406 which controls MR data acquisition device 402 accordingly.
- the acquired MR data is stored in storage 408 .
- Program 410 analyses the acquired MR data which are stored in storage 408 in accordance with the principles as explained above with reference to FIGS. 1 to 3 . This way a perfusion tensor is obtained which is stored in storage 412 .
- the perfusion tensors stored in storage 412 are processed by imaging program 414 which generates a perfusion tensor image which is stored in frame buffer 416 for display on display unit 418 connected to computer system 404 .
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Vascular Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The present invention relates to a method of perfusion imaging comprising: performing a first magnetic resonance data acquisition (A) at a first sensitivity (b) value, performing a set of at least six second magnetic resonance data acquisitions (B1, B2, . . . B6) with gradiant encodings in different directions at second sensitivity (b) values, determining a perfusion tensor based on the magnetic resonance data acquisitions, performing a perfusion tensor visualitation step.
Description
The present invention relates to the field of magnetic resonance imaging, and more particularly to perfusion tensor imaging.
Diffusion weighted magnetic resonance imaging is a known prior art technique. In diffusion weighted magnetic resonance imaging the diffusion tensor is obtained from, the magnetic resonance measurement signal for a defined body region of interest. The diffusion tensor is then visualised by known imaging methods (U.S. Pat. No. 5,539,310; La Bilhan et al, Diffusion Tensor Imaging: Concepts and Applications, Journal of Medical Resonance Imaging 13:534-46 (2001); Lawrence R. Frank, Anisotropy in High Angular Resolution Diffusion-Weighted MRI, Magnetic Resonance in Medicine 45:935-939 (2001),). An overview of known diffusion-weighted magnetic resonance techniques is provided by Gray L. MacFall J. Overview of diffusion Imaging. Magnetic Resonance Imaging Clin N Am 1998; 6:125-138.
In addition in vivo intravoxel incoherent motion (IVIM) magnetic resonance imaging is a known technique (La Bilhan D et al, MR imaging of intravoxel inchoherent motions: application to diffusion and perfusion in neurological disorders. Radiology 1986:161:401-407.) In particular, the interpretation of IVIM measurements with respect to classical perfusion has been discussed in the prior art (La Bilhan D et al., Separation of diffusion and perfusion in intravoxel incoherent motion MR Imaging. Radiology 1998:168:497-505; Henkelman R M. Does IVIM measure classical perfusion? Magn Reson Med 1990:16:470-475; Le Bilhan D, Turner R. The capillary network: a link between IVIM and classical perfusion. Magn Reson Med Med 1992:27:171-178). This imaging technique uses gradient encodings for random motion.
Classical perfusion is a measure of the blood delivered to and used by a specified mass of tissue. It is often measured using spin labelling techniques in MRI. In contrast, IVIM measures quasi-random blood movement within a single imaging voxel and results in a bi-exponential signal attenuation in a standard pulsed gradient spin echo (PGSE) experiment (cf. R. J. Moore et al, in vivo intravoxel incoherent motion measurements in the human placenta using echo-planar imaging at 0.5 T, Magnetic Resonance in Medicine, 43;295-302 (2000)).
It is a common disadvantage of prior measurement techniques that only a scalar value for the perfusion is obtained even though perfusion usually is anisotropic. A scalar perfusion value therefore only gives a limited amount of information on the actual nature of the perfusion.
Therefore the present invention aims to provide an improved imaging method which enables imaging and visualisation of the perfusion anisotropy as well as a corresponding computer program product and perfusion imaging apparatus.
The present invention provides for a method of perfusion imaging which can be used for medical purposes such as abdominal imaging.
The present invention is based on the discovery that an isotropic perfusion can be described by a tensor which is structurally equivalent to a diffusion tensor used for diffusion tensor imaging (DTI). In order to provide the data for determination of a perfusion tensor a magnetic resonance data acquisition is performed at a low b sensitivity value, preferably b=0, as well as at least six magnetic resonance data acquisitions with gradient encodings in different directions at b sensitivity values of for example below 50 [s/mm2], preferably between 5 and 15.
In accordance with a preferred embodiment of the invention six slopes are determined from the 7 measurement points. The slope values form the basis for calculating the perfusion tensor, especially its eigenvectors. This is based on the discovery, that for lower b sensitivity values the signal decay is mostly due to perfusion and not diffusion effects. This is why diffusion can be neglected at low b sensitivity values in order to extract the information out of the decay signal which forms the basis for determining the shape of the perfusion tensor. At least, the sorting of the eigenvalues and the direction of the eigenvectors of the perfusion tensor does not depend much on the diffusion-related signal decay. This directional information derived from the perfusion tensor can be visualized. For this purpose software algorithms which are used in the prior art for diffusion visualization can be reused due to the same mathematics governing the visualization of diffusion and perfusion tensors.
In accordance with a further preferred embodiment of the invention, two additional magnetic resonance data acquisitions are performed for higher b sensitivity values. For example one additional magnetic resonance data acquisition is performed for an intermediate b sensitivity value of around 200 and a further magnetic resonance data acquisition is performed at a high b sensitivity value of around 800.
For such high b sensitivity values the decay signal is mostly governed by its diffusion signal component, see e.g. [Petra Mürtz, et al. Abdomen: Diffusion-weighted MR Imaging with Pulse-triggered Single-Shot Sequences, 258-264, Radiology July 2002]. Especially in the abdomen, the diffusion is isotropic, as demonstrated in the cited reference. Thus, the diffusion coefficient can be obtained from just two higher b-values measurements. For brain imaging it is preferred to perform at least six data acquisitions with gradient encodings for random motion in different directions as diffusion is anisotropic in this case.
Mathematical analysis also provides the blood fraction value. This fraction is known to be isotropic for each measured signal, thus constraining the analysis. Based on the diffusion coefficient and the blood fraction value the diffusion component of the decay signal can be estimated. By subtracting the diffusion signal component from the measurement signal the perfusion signal component is obtained. The perfusion coefficient is obtained from the perfusion signal component by determining the signal slope of the perfusion signal component, thereto the b=0 and low b-value measurements are analyzed.
In accordance with a further preferred embodiment of the invention only one additional magnetic resonance data acquisition is performed for a high b sensitivity value, leaving out the intermediary measurement. In this case one of the magnetic resonance data acquisitions which have been performed for the low b sensitivity values having the highest signal value is selected as a replacement for the intermediate measurement.
It is to be noted that for the intermediate and high b data acquisitions no gradient encodings in different directions are necessary as diffusion is isotropic for such high b values.
It is a particular advantage of the present invention that due to the limited number of data acquisitions, all of the required data acquisitions can be performed for a couple of slices through the body (typically 5-10) within a single breath hold, for example in less than 16 seconds.
In the following preferred embodiments of the invention will be described in greater detail by making reference to the drawings in which:
The general idea of intra-voxel incoherent motion imaging (IVIM) is that the observed MR signal decay has a bi-exponential behaviour as function of the diffusion weighting factor b, i.e.
S/S 0 =A1*exp(−b*P)+A2*exp(−b*D), (1)
where P is the perfusion constant (in mm2/sec, typically 0.05-0.08), and D the diffusion constant (in mm2/sec, typically 0.002).
S/S 0 =A1*exp(−b*P)+A2*exp(−b*D), (1)
where P is the perfusion constant (in mm2/sec, typically 0.05-0.08), and D the diffusion constant (in mm2/sec, typically 0.002).
It can readily be appreciated that (for b=0) A1+A2=1. Thus, recognizing that the signal related to perfusion can be attributed to (randomly) flowing material in the voxel, A1 has been designated to represent f, the fraction of flowing material (i.e. blood) in the voxel. Hence f denotes the blood fraction and A2=1−f, which in other words means that the non-flowing part of the signal contributes to the diffusion signal. The fraction f, being representative of the blood content, is an isotropic quantity for each voxel. This knowledge is used to constrain the data analysis in the rest of this section.
This leads to
S/S 0 =f*exp(−bP)+(1−f)*exp(−bD). (2)
The first term of equation 2 will be referred to as the perfusion signal component and the second term in equation 2 will be referred to as the diffusion signal component. The perfusion signal component is decaying much more rapidly than the diffusion signal component and thus primarily determines the slope of the signal for low b values.
S/S 0 =f*exp(−bP)+(1−f)*exp(−bD). (2)
The first term of equation 2 will be referred to as the perfusion signal component and the second term in equation 2 will be referred to as the diffusion signal component. The perfusion signal component is decaying much more rapidly than the diffusion signal component and thus primarily determines the slope of the signal for low b values.
The perfusion is described by a perfusion tensor, which has the same structure as the diffusion tensor, i.e. the perfusion tensor is of second rank and symmetric.
The perfusion tensor can be determined on the basis of the bi-exponential signal of equation 2, using the signal decays for low b values whereby the diffusion signal component is neglected. This is illustrated by means of FIG. 1 .
As it can be seen from the signal decay curves of FIG. 1 the curves are approximately linear for low b values.
For the purpose of determining the perfusion tensor the magnetic resonance data acquisition for point A is performed without gradient encoding and for the points B1, B2, . . . B6 at six different gradient directions, respectively. The slopes m1, m2, . . . m6 of the signal decay curves for low b-values are determined as illustrated in FIG. 1 by linear approximation.
By means of the six slope values the perfusion tensor can be calculated by means of the same mathematics as used for diffusion tensor calculation (c.f. for example M. T. Vlaadringerbroek and J. A. den Boer. Magnetic Resonance Imaging. Springer-Verlag Berlin Heidelberg New York, 1999, in particular section 7.7).
In step 206 the slope values mi are obtained from the measurements points A and Bi, B2, . . . B6 by linear approximation, i.e. from the slopes of the straight lines connecting A to B1, A to B2 . . . A to B6. These slope values provide full mapping of the perfusion tensor and are input into the perfusion tensor calculation routine of step 208.
For increased precision the diffusion signal component of equation (2) is also taken into consideration for the perfusion tensor determination. One way of accomplishing this is illustrated by making reference to FIG. 3 .
For determination of the diffusion coefficient D two magnetic resonance data acquisitions are performed for b values within the linear portion of the signal decay curves. For example this can be done for b=200 (point X) and for b=800 (point Y).
The slope of the straight line between points X and Y is the diffusion coefficient D. The extrapolation of the straight line between X and Y to b=0 provides 1−f as illustrated in FIG. 3 . This way the diffusion signal component of equation (2) is obtained. It is to be noted that the extrapolation by means of the straight line is performed for convenience of explanation; in a practical application it is preferred to perform the analysis of the bi-exponential decay using well known mathematical routines such as available from the LAPACK mathematical library.
The perfusion signal component is analysed by subtracting the estimated diffusion signal component, i.e. Sest=(1−f)*exp(−bD), from the measured value S. For the b=0 value of the curve this provides a value S0*f, denoted as S0′. For any signal related to non-zero b-values, a value S′ is derived. A new curve can now be drawn to show ln(S0′/S′). The slope of this curve provides P.
In other words, by subtracting the estimated diffusion signal component, the perfusion signal component is isolated from the measured data values; as f is also known the perfusion coefficient P can be obtained this way.
For example, this enables a pre-processing of the signal decay curves 100, 102, 104, . . . of FIG. 1 by subtracting the diffusion signal component before the determination of the slopes mi.
Alternatively the perfusion coefficient P is calculated as explained above for every diffusion-sensitisation gradient direction. The resulting set of at least six perfusion coefficients is used to calculate the perfusion tensor map. Specifically, the perfusion tensor's eigenvalues, eigenvectors, and rotationally invariant quantities like the trace or the fractional anisotropy are derived and visualized. Also, voxels with substantially co-linear main eigenvectors can be connected (starting from a user defined seed position), and the connecting line can be displayed (tractography).
In order to avoid the measurement for point X at the intermediate b value it is also possible to use the highest value of the B1, B2, . . . B6 measurements for low b values as an approximation. In the example considered here the highest value is B (cf. FIG. 1 ). In order to obtain an approximation for D and 1−f, a straight line between Y and B1 can be used instead of a straight line between X and Y.
One of the advantages of keeping the number of high b-value acquisitions at a minimum is that each acquisition prolongs the breath hold duration by typically 1 to 3 seconds. Completely avoiding higher b-value acquisitions enables the use of short echo times, because such high b-value acquisitions require much gradient area, thus having a relatively long echo time (typically 60-90 ms, whereas T2 relaxation in the abdomen is only 100 ms). The long echo time relative to T2 relaxation significantly lowers the SNR of the measurements. Using only small b values (<50) significantly reduces the required gradient area, and thus enables lower echo times (typically 20-30 ms). This significantly improves the SNR.
These sequences are determined by control program 406 which controls MR data acquisition device 402 accordingly. The acquired MR data is stored in storage 408. Program 410 analyses the acquired MR data which are stored in storage 408 in accordance with the principles as explained above with reference to FIGS. 1 to 3 . This way a perfusion tensor is obtained which is stored in storage 412. The perfusion tensors stored in storage 412 are processed by imaging program 414 which generates a perfusion tensor image which is stored in frame buffer 416 for display on display unit 418 connected to computer system 404.
- 100 signal decay curve
- 102 signal decay curve
- 104 signal decay curve
- 400 performance imaging system
- 402 MR data acquisition device
- 404 computer system
- 406 control program
- 408 storage
- 410 program
- 412 storage
- 414 imaging program
- 416 frame buffer
- 418 display unit
Claims (20)
1. A method of perfusion imaging comprising:
performing a first magnetic resonance data acquisition with gradient encodings for random motion at a first low b sensitivity value;
performing a set of at least six second magnetic resonance data acquisitions with gradient encodings for random motion in different directions at second sensitivity values larger than the first low b sensitivity value; and
determining a perfusion tensor based on the magnetic resonance data acquisitions.
2. The method of perfusion imaging of claim 1 , the second sensitivity values being below 50 s/mm2.
3. The method of perfusion imaging of claim 1 , wherein the first low b sensitivity value is substantially zero and the second sensitivity values being between five and thirty.
4. The method of perfusion imaging of claim 1 , the magnetic resonance data acquisitions being performed by a series of single-shot echo-planar magnetic resonance data acquisitions.
5. The method of perfusion imaging of claim 1 , further comprising:
performing a perfusion tensor visualisation.
6. The method of perfusion imaging of claim 5 , further comprising:
deriving directional information from the perfusion tensor, the performing of a perfusion tensor visualization including visualizing the derived directional information.
7. The method of perfusion imaging of claim 1 , further comprising:
determining first slope values between each one of the set of magnetic resonance data acquisitions and the first magnetic resonance data acquisition; and
determining the perfusion tensor based on the first slope values.
8. The method of perfusion imaging of claim 1 , further comprising:
performing a third magnetic resonance data acquisition at a third sensitivity value;
performing a fourth magnetic resonance data acquisition at a fourth sensitivity value, the third sensitivity value being higher than the second sensitivity values, and the fourth sensitivity value being higher than the third sensitivity value;
determining a diffusion coefficient and a fraction value based on the third and the fourth magnetic resonance data acquisitions to provide a diffusion signal component;
eliminating the diffusion signal component from the magnetic resonance data acquisitions to provide a perfusion signal component; and
determining a perfusion tensor from the perfusion signal components.
9. The method of perfusion imaging of claim 8 , wherein a set of at least six third magnetic resonance data acquisitions with gradient encodings for random motion in different directions at third sensitivity values is performed, and a set of at least six fourth magnetic resonance data acquisitions with gradient encodings for random motion in different directions at fourth sensitivity values is performed, and the diffusion tensor is determined based on the third and the fourth magnetic resonance data acquisitions to provide a diffusion signal component.
10. The method of perfusion imaging of claim 8 , the third sensitivity value being between 100 and 400, and the second sensitivity value being between 600 and 1200.
11. A method of perfusion imaging comprising:
performing a first magnetic resonance data acquisition with gradient encodings for random motion at a first low b sensitivity value;
performing a set of at least six second magnetic resonance data acquisitions with gradient encodings for random motion in different directions at second low b sensitivity values;
selecting one of the second magnetic resonance data acquisitions having the strongest measured signal decay;
performing a third magnetic resonance data acquisition at a third sensitivity value, the third sensitivity value being higher than the second sensitivity values;
determining a diffusion coefficient and a fraction value based on the selected second and third magnetic resonance data acquisitions to provide a diffusion signal component;
eliminating the diffusion signal component from the magnetic resonance data acquisitions to provide a perfusion signal component; and
determining a perfusion tensor from the perfusion signal components.
12. A computer program product comprising a digital storage medium storing a perfusion imaging program executable to perform a method including determining a perfusion tensor based on a first magnetic resonance data acquisition and a set of at least six second magnetic resonance data acquisitions, the first magnetic resonance data acquisition being performed at a first low b sensitivity value and the second magnetic data resonance data acquisitions being performed at a second low b sensitivity value with gradient encodings in different directions, the first low b sensitivity value being below the second low b sensitivity values, and further executable to perform perfusion tensor imaging.
13. The computer program product of claim 12 , wherein the perfusion imaging program is further executable to determine first slope values for each one of the second magnetic resonance data acquisitions based on the first magnetic resonance data acquisition and to determine the perfusion tensor based on the first slope values.
14. The computer program product of claim 12 , wherein the perfusion imaging program is further executable to determine a diffusion coefficient and a fraction value based on third and fourth magnetic resonance data acquisitions to provide a diffusion signal component, to eliminate the diffusion signal component from the magnetic resonance data acquisitions to provide a perfusion signal component, and to determine a perfusion tensor from the perfusion signal component.
15. The computer program product of claim 14 , wherein the perfusion imaging program is further executable to process a set of at least six third magnetic resonance data acquisitions with gradient encodings for random motion in different directions at third sensitivity values, a set of at least six fourth magnetic resonance data acquisitions with gradient encodings for random motion in different directions at fourth sensitivity values, and to determine the diffusion tensor based on the third and the fourth magnetic resonance data acquisitions to provide a diffusion signal component.
16. The computer program product of claim 12 , wherein the perfusion imaging program is further executable to select one of the second magnetic resonance data acquisitions having the highest data value, determining of a diffusion coefficient and a fraction value based on the selected second magnetic resonance data acquisition and a third magnetic resonance data acquisition being performed at a third sensitivity value, the third sensitivity value being above the second sensitivity values, providing a diffusion signal component based on the diffusion coefficient and the blood fraction value, eliminating of the diffusion signal component from the magnetic resonance data acquisitions to provide a perfusion signal component, and to determine of a perfusion tensor from the perfusion signal components.
17. A perfusion imaging apparatus comprising:
a magnetic resonance data acquisition device; and
a computer system programmed to cause the magnetic resonance data acquisition device to perform a first magnetic resonance data acquisition at a first low b sensitivity value and to perform a set of at least six second magnetic resonance data acquisitions with gradient encodings in different directions at second low b sensitivity values below 50 s/mm2, the first sensitivity value being smaller than the second sensitivity values, and further programmed to determine a perfusion tensor based on the magnetic resonance data acquisitions.
18. The perfusion imaging apparatus as set forth in claim 17 , further comprising:
a display unit cooperating with the computer system to perform perfusion tensor imaging based on the determined perfusion tensor.
19. The perfusion imaging apparatus as set forth in claim 17 , wherein the first low b sensitivity value is substantially zero.
20. The perfusion imaging apparatus as set forth in claim 17 , wherein the computer system is further programmed to cause the magnetic resonance data acquisition device to perform at least one additional magnetic resonance data acquisition at a higher sensitivity value than the second low b sensitivity values, and is further programmed to determine a diffusion signal component of the first and second magnetic resonance data acquisitions based on the at least one additional magnetic resonance data acquisition.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP03100848 | 2003-03-31 | ||
EP03100848.5 | 2003-03-31 | ||
PCT/IB2004/050322 WO2004088345A1 (en) | 2003-03-31 | 2004-03-23 | A method of magnetic resonance perfusion imaging |
Publications (2)
Publication Number | Publication Date |
---|---|
US20060241375A1 US20060241375A1 (en) | 2006-10-26 |
US7310548B2 true US7310548B2 (en) | 2007-12-18 |
Family
ID=33104158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US10/551,068 Expired - Fee Related US7310548B2 (en) | 2003-03-31 | 2004-03-23 | Method of magnetic resonance perfusion imaging |
Country Status (4)
Country | Link |
---|---|
US (1) | US7310548B2 (en) |
EP (1) | EP1611452A1 (en) |
JP (1) | JP2006521863A (en) |
WO (1) | WO2004088345A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090096448A1 (en) * | 2007-10-15 | 2009-04-16 | Siemens Corporate Research, Inc. | b-Value Optimization for Diffusion Weighted Magnetic Resonance Imaging |
US9513358B2 (en) | 2013-03-12 | 2016-12-06 | Vaposun Inc. | Method and apparatus for magnetic resonance imaging |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8381140B2 (en) * | 2011-02-11 | 2013-02-19 | Tokyo Electron Limited | Wide process range library for metrology |
WO2013156222A2 (en) * | 2012-04-16 | 2013-10-24 | Centre Hospitalier Universitaire Vaudois (Chuv) | Method for generating perfusion images |
SE537064C2 (en) | 2012-05-04 | 2014-12-23 | Cr Dev Ab | Analysis for quantification of microscopic anisotropic diffusion |
JP2014195532A (en) * | 2013-03-29 | 2014-10-16 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Estimation device, magnetic resonance device, program, and estimation method |
EP4383276A3 (en) * | 2014-10-20 | 2024-07-17 | New York University | Method to select radiation dosage for tumor treatment based on cellular imaging |
JP6697261B2 (en) | 2014-12-26 | 2020-05-20 | キヤノンメディカルシステムズ株式会社 | Magnetic resonance imaging apparatus, diffusion weighted image generation method, and image processing apparatus |
WO2018088954A1 (en) | 2016-11-09 | 2018-05-17 | Cr Development Ab | A method of performing diffusion weighted magnetic resonance measurements on a sample |
WO2021096946A1 (en) * | 2019-11-15 | 2021-05-20 | The University Of Chicago | Cancer detection based on four quadrant mapping and matrix analysis of image data |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5539310A (en) * | 1993-08-06 | 1996-07-23 | The United States Of America As Represented By The Department Of Health And Human Services | Method and system for measuring the diffusion tensor and for diffusion tensor imaging |
US20010039377A1 (en) * | 2000-03-31 | 2001-11-08 | Maier Stephan E. | Method of fast and reliable tissue differentiation using diffusion-weighted magnetic resonance imaging |
US6452390B1 (en) * | 1999-11-16 | 2002-09-17 | Wollin Ventures, Inc. | Magnetic resonance analyzing flow meter and flow measuring method |
US20040106864A1 (en) * | 2001-03-07 | 2004-06-03 | Rose Stephen Edward | Method of predicting stroke evolution utilising mri |
US7020578B2 (en) * | 2000-02-02 | 2006-03-28 | The General Hospital Corporation | Method for evaluating novel, stroke treatments using a tissue risk map |
US20060098876A1 (en) * | 2002-03-11 | 2006-05-11 | Buscema Paolo M | Method for encoding image pixels a method for processing images and a method for processing images aimed at qualitative recognition of the object reproduced by one or more image pixels |
US7078897B2 (en) * | 2002-01-16 | 2006-07-18 | Washington University | Magnetic resonance method and system for quantification of anisotropic diffusion |
-
2004
- 2004-03-23 JP JP2006506761A patent/JP2006521863A/en not_active Withdrawn
- 2004-03-23 EP EP04722617A patent/EP1611452A1/en not_active Withdrawn
- 2004-03-23 US US10/551,068 patent/US7310548B2/en not_active Expired - Fee Related
- 2004-03-23 WO PCT/IB2004/050322 patent/WO2004088345A1/en active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5539310A (en) * | 1993-08-06 | 1996-07-23 | The United States Of America As Represented By The Department Of Health And Human Services | Method and system for measuring the diffusion tensor and for diffusion tensor imaging |
US6452390B1 (en) * | 1999-11-16 | 2002-09-17 | Wollin Ventures, Inc. | Magnetic resonance analyzing flow meter and flow measuring method |
US7020578B2 (en) * | 2000-02-02 | 2006-03-28 | The General Hospital Corporation | Method for evaluating novel, stroke treatments using a tissue risk map |
US20010039377A1 (en) * | 2000-03-31 | 2001-11-08 | Maier Stephan E. | Method of fast and reliable tissue differentiation using diffusion-weighted magnetic resonance imaging |
US20040106864A1 (en) * | 2001-03-07 | 2004-06-03 | Rose Stephen Edward | Method of predicting stroke evolution utilising mri |
US7078897B2 (en) * | 2002-01-16 | 2006-07-18 | Washington University | Magnetic resonance method and system for quantification of anisotropic diffusion |
US20060098876A1 (en) * | 2002-03-11 | 2006-05-11 | Buscema Paolo M | Method for encoding image pixels a method for processing images and a method for processing images aimed at qualitative recognition of the object reproduced by one or more image pixels |
Non-Patent Citations (14)
Title |
---|
Basser, P.J., et al.; Estimation of the Effective Self-Diffusion Tensor from the NMR Spin Echo; 1994; J. of MR; B:103:247-254. |
Callot, V., et al.; An IVIM based MRI method to study microcirculation in the heart; 2003; Proc. of SPIE; 5031:214-220. |
Frank, L.R.; Anisotropy in High Angular Resolution Diffusion-Weighted MRI; 2001; MR in Med.; 45:935-939. |
Gray, L., et al.; Overview of Diffusion Imaging; 1998; Magnetic Resonance Imaging Clin N Am; 6:(1)125-137. |
Henkelman, R.M.; Does IVIM Measure Classical Perfusion?; 1990; MRM; 16:470-475. |
Le Bihan, D., et al.; Diffusion Tensor Imaging: Concepts and Applications; 2001; J. of MRI; 13:534-546. |
Le Bihan, D., et al.; Intravoxel Incoherent Motion Imaging Using Spin Echoes; 1991; MR in Med.; 19:221-227. |
Le Bihan, D., et al.; MR Imaging of Intravoxel Incoherent Motions: Application to Diffusion and Perfusion in Neurologic Disorders; 1986; Radiology; 161:401-407. |
Le Bihan, D., et al.; Separation of Diffusion and Perfusion in Intravoxel Incoherent Motion MR Imaging; 1988; Radiology; 168:497-505. |
Le Bihan, D., et al.; The Capillary Network: A link between IVIM and classical perfusion; 1992; MRM; 27:171-178. |
Le Bihan, D.; Diffusion & Perfusion in MRI;1996; in Encyclopedia of Nuclear Magnetic Resonance; Grant-Harris; vol. 3; Con-F; pp. 1645-1656. |
Moore, R.J., et al.; In Vivo Intravoxel Incoherent Motion Measurements in the Human Placenta Using Echo-Planar Imaging at 0.5T; 2000; MRM; 43:295-302. |
Murtz, P., et al.; Abdomen: Diffusion-weighted MR Imaging with Pulse-triggered Single-shot Sequences; 2002: Radiology; 224:258-264. |
Vlaadringerbroek, M.T., et al.; 1999; Magnetic Resonance Imaging; Springer-Vertag, Berlin; Section 7.7; pp. 354-360. |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090096448A1 (en) * | 2007-10-15 | 2009-04-16 | Siemens Corporate Research, Inc. | b-Value Optimization for Diffusion Weighted Magnetic Resonance Imaging |
US7560926B2 (en) * | 2007-10-15 | 2009-07-14 | Siemens Aktiengesellschaft | b-value optimization for diffusion weighted magnetic resonance imaging |
US9513358B2 (en) | 2013-03-12 | 2016-12-06 | Vaposun Inc. | Method and apparatus for magnetic resonance imaging |
Also Published As
Publication number | Publication date |
---|---|
JP2006521863A (en) | 2006-09-28 |
EP1611452A1 (en) | 2006-01-04 |
WO2004088345A1 (en) | 2004-10-14 |
US20060241375A1 (en) | 2006-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8406849B2 (en) | Systems and methods for magnetic resonance imaging | |
US6516210B1 (en) | Signal analysis for navigated magnetic resonance imaging | |
US6922580B2 (en) | Blood flow gated MRI | |
EP2798364B1 (en) | Mr imaging with suppression of flow artefacts | |
US10379189B2 (en) | Simultaneous magnetic resonance angiography and perfusion with nuclear magnetic resonance fingerprinting | |
US20080114234A1 (en) | Method, apparatus and user interface for determining an arterial input function used for calculating hemodynamic parameters | |
JP5235283B2 (en) | Method and apparatus for determining the coefficient of diffusion tensor by magnetic resonance | |
RU2683605C1 (en) | Parallel mri with sensitivity mapping using rf coil | |
EP3155442B1 (en) | Mri method to quantify iron amount in tissues using diffusion magnetic resonance imaging | |
US10502801B2 (en) | Method and magnetic resonance apparatus for generating a weighting matrix for reducing artifacts with parallel imaging | |
US20160116560A1 (en) | Mr imaging with enhanced susceptibility contrast | |
US7310548B2 (en) | Method of magnetic resonance perfusion imaging | |
US11313931B2 (en) | System and method for quantifying T1, T2 and resonance frequency using rosette trajectory acquisition and read segmented reconstruction | |
US8995738B2 (en) | System and method for magnetic resonance imaging parametric mapping using confidence maps | |
US20200033432A1 (en) | Method and apparatus for identifying at least one material comprised in a voxel | |
CN109716155A (en) | With Dickson type water/fat separation MR imaging | |
WO2014111266A2 (en) | Synthetic parameter-weigthed images based on relaxometry (spare) | |
WO2004013648A1 (en) | Measurement of the macromolecular proton density involved in magnetization transfer by means of magnetic resonance | |
US11519987B2 (en) | Magnetic resonance fingerprinting thermometry | |
DE112018006309T5 (en) | FAT SUPPRESSION DIFFUSION MR IMAGING | |
Nakashima et al. | Detecting mild lower-limb skeletal muscle fatigue with stimulated-echo q-space imaging | |
US11768264B2 (en) | System and method for magnetic resonance fingerprinting with relaxation and diffusion data acquisition | |
CN111521962A (en) | System and method for separating water and fat signals in spin lock magnetic resonance imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: KONINKLIJKE PHILIPS ELECTRONICS N.V., NETHERLANDS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:VAN DEN BRINK, JOHAN SAMUEL;REEL/FRAME:017851/0479 Effective date: 20041028 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20111218 |