Diffusion imaging of tissues

• US 6,614,226 B2
• Filed: 04/02/2001
• Issued: 09/02/2003
• Est. Priority Date: 03/31/2000
• Status: Expired due to Term

• Current Assignee: The General Hospital Corporation (Mass General Brigham, Inc.)
• Original Assignee: The General Hospital Corporation (Mass General Brigham, Inc.)
• Inventors: Wedeen, Van
• Primary Examiner(s): Lefkowitz, Edward
• Assistant Examiner(s): Shrivastav, Brij B.
• Application Number:
  US09/825,469
• Publication Number:
  US 20020042569A1
• Time in Patent Office: 883 Days
• Field of Search: 324/307, 324/309, 324/303, 324/306, 324/312, 324/314, 324/318, 324/320, 324/322, 600/410
• US Class Current: 324/309
• CPC Class Codes: A61B 5/055 involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
• CPC Class Codes: A61B 5/4064 Evaluating the brain A61B5/031, A61B5/369, A61B5/14553 take precedence
• CPC Class Codes: A61B 5/4519 Muscles A61B5/389, A61B5/224 take precedence
• CPC Class Codes: A61B 5/7257 using Fourier transforms
• CPC Class Codes: G01R 33/56341 Diffusion imaging

Abstract

The invention relates to the use of diffusion spectrum magnetic resonance imaging (MRI) to map complex fiber architectures in tissues. The new methods can be used to resolve intravoxel heterogeneity of diffusion in vivo with MRI of diffusion density spectra.

Claims

1. A method of constructing an image representative of structure within a tissue, the method comprising

(a) inducing a population of spins in the tissue to produce a set of nuclear magnetic resonance (NMR) signals, wherein the signals comprise a set of complex Fourier-encodings of a distribution of three-dimensional displacements of the spins in the population;

(b) converting each of the NMR signals in the set of complex Fourier-encodings into a positive number to form a set of positive numbers to exclude phase information;

(c) reconstructing from the set of positive numbers a function that approximates the distribution of three-dimensional displacements of the spins in the population in multiple directions; and

(d) constructing an image that represents the function, whereby the image represents structure within the tissue.

2. The method of claim 1, wherein each spin in the population is within a three-dimensional voxel.

3. The method of claim 1, wherein the population of spins is induced to produce the set of NMR signals by the application of a set of magnetic gradient pulses.

4. The method of claim 3, wherein the magnetic gradient pulses are applied in a pulse train whose time-intensity integral is zero.

5. The method of claim 4, wherein the pulses in the pulse train are bipolar gradient pulses.

6. The method of claim 5, wherein the gradient pulses are transected by one or more 180°

radio frequency (RF) pulses and wherein a gradient sign is reversed following each 180°

pulse.

7. The method of claim 1, wherein the NMR signals are converted into positive numbers by determining the modulus (z->

|z|).

8. The method of claim 1, wherein the NMR signals are converted into positive numbers by determining the squared modulus (z->

zz*).

9. The method of claim 1, wherein the function is reconstructed by determining the discrete Fourier transform.

10. The method of claim 1, wherein the function is reconstructed by interpolation and regridding followed by determining the discrete Fourier transform.

11. The method of claim 1, wherein the method is performed for multiple contiguous locations, and further constructing a curve that represents fiber tracts in the tissue that conform to the orientation of directions of maximum displacement.

12. The method of claim 1, wherein the image is a three-dimensional image.

13. The method of claim 12, wherein the three-dimensional image represents the three-dimensional distribution of spin displacement for a voxel.

14. The method of claim 12, wherein the image is a three-dimensional graphic of the amplitude of spin displacement in multiple directions.

15. The method of claim 12, wherein the polar plot is colored to represent the amplitude and orientation of spin displacement.

16. The method of claim 15, wherein the color is coded to assign red, green, and blue to the amplitude of spin displacement in each of three orthogonal coordinates.

17. The method of claim 16, wherein the amplitude of spin displacement is the relative probability of spins displacing a constant distance in any direction.

18. The method of claim 12, wherein the image is a density plot of the spin density in a position-angle space.

19. The method of claim 18, wherein the image is a slice through a position-angle space.

20. The method of claim 18, wherein the image is a projection through a position-angle space.

21. The method of claim 1, wherein the tissue is a heterogeneous tissue.

22. The method of claim 21, wherein the tissue is brain tissue.

23. The method of claim 22, wherein the tissue is neural white matter.

24. The method of claim 23, wherein the neural white matter has multiple fiber orientations.

25. The method of claim 21, wherein the tissue comprises normal and pathologic tissue.

26. The method of claim 25, wherein the pathologic tissue comprises cerebral edema, cerebral hematoma, cerebral neoplasm, cerebral metastasis, or ischemic tissue.

27. The method of claim 25, wherein the pathologic tissue comprises a neurodegenerative disease.

28. The method of claim 27, wherein the neurodegenerative disease is Huntington'"'"'s chorea, multiple sclerosis, or stroke.

29. The method of claim 1, wherein the tissue is muscle.

30. The method of claim 29, wherein the muscle is heart or tongue.

31. The method of claim 1, further comprising diagnosing a disorder in a tissue using the image.

32. The method of claim 31, wherein the tissue is brain tissue.

33. The method of claim 32, further comprising constructing a model of fiber tracts in the brain based on the image.

34. The method of claim 1, further comprising mapping a surgical site in the tissue using the image.

35. The method of claim 34, wherein the surgical site is in the brain.

36. The method of claim 34, wherein the surgical site is in the heart.

37. The method of claim 1, further comprising combining the imaging method with one or more additional magnetic resonance imaging (MRI) contrast parameters.

38. The method of claim 37, wherein the contrast parameter is T1, T2, magnetization transfer contrast, or blood oxygen level dependent contrast (BOLD).

40. The method of claim 1, where the structure is one or more fiber orientations within the tissue.

41. The method of claim 14, wherein the graphic is a polar plot.

42. The method of claim 1, wherein the set of complex Fourier-encodings comprises a 3D grid.

43. The method of claim 1, wherein the set of complex Fourier-encodings comprises a 3D spherical shell.

39. A computer-implemented program for constructing an image representative of structure within a tissue, the program comprising a plurality of program instructions stored on a electronic apparatus-readable medium for implementing the steps of:

(a) inducing a population of spins in the tissue to produce nuclear magnetic resonance (NMR) signals, wherein the signals comprise a set of complex Fourier-encodings of a distribution of three-dimensional displacements of the spins in the population;

(b) converting each of the NMR signals in the set of complex Fourier-encodings into a positive number to form a set of positive numbers to exclude phase information;

(c) reconstructing from the set of positive numbers a function that approximates the distribution of three-dimensional displacements of the spins in the population in multiple directions; and

(d) constructing an image that represents the function, whereby the image represents structure within the tissue.