B-value calculation and correction using a linear segment gradient waveform model

• US 6,670,812 B1
• Filed: 09/13/2002
• Issued: 12/30/2003
• Est. Priority Date: 09/13/2002
• Status: Expired due to Term

• Current Assignee: GE Medical Systems Global Technology Company LLC (GE Healthcare Technologies, Inc.)
• Original Assignee: GE Medical Systems Global Technology Company LLC (GE Healthcare Technologies, Inc.)
• Inventors: Michelich, Charles R., Mock, Bryan J.
• Primary Examiner(s): Lefkowitz, Edward
• Assistant Examiner(s): Shrivastav, Brij B.
• Application Number:
  US10/243,162
• Time in Patent Office: 473 Days
• Field of Search: 324/309, 324/307, 324/311, 600/410, 128/653.2
• US Class Current: 324/309
• CPC Class Codes: G01R 33/543 Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription G01R33/546 takes precedence
• CPC Class Codes: G01R 33/5608 Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels image data processing or generation, in general G06T
• CPC Class Codes: G01R 33/56341 Diffusion imaging

Abstract

A technique for calculating actual b-values in a diffusion weighted magnetic resonance imaging system and of adjusting the observed b-value until it is within a tolerance range of the prescribed b-value. A diffusion weighted image with the desired diffusion sensitivity may then be produced. Calculation of the actual b-value is accomplished by integrating the linear segments of the gradient waveform. Adjustment of observed b-value is accomplished by initially scaling the diffusion lobe amplitude and comparing a recalculated b-value to the prescribed b-value. Additional adjustment is provided by means of a numeric search iteratively performed until the recalculated b-value is within tolerance of the prescribed b-value.

Claims

1. A method for calculating a b-value comprising:

generating a gradient waveform comprising a plurality of timing data points and one or more waveform data points using a magnetic resonance imaging technique;

calculating two or more corner points associated with the gradient waveform;

determining an integration start point such that the start point corresponds to an excitation pulse generated by the magnetic resonance imaging technique;

determining an integration end point such that the end point corresponds to an echo time resulting from the excitation pulse;

generating an area sum by integrating one or more linear waveform segments between the start point and the end point wherein each linear waveform segment is defined by a line connecting adjacent corner points; and

calculating the b-value using the area sum.

2. The method of claim 1, further comprising the act of determining one or more inversion points corresponding to one or more refocusing pulses.

3. The method of claim 1, wherein calculating the two or more corner points comprises taking a derivative of the gradient waveform.

4. The method of claim 3, wherein each corner point corresponds to a change in slope.

5. The method of claim 1, wherein determining the integration start point comprises designating the corner point which corresponds to the excitation pulse as the start point.

6. The method of claim 1, wherein determining the integration start point comprises interpolating between two adjacent corner points to create an interpolated point which corresponds to the excitation pulse and designating the interpolated point as the start point.

7. The method of claim 1, wherein determining the integration end point comprises designating the corner point which corresponds to the echo time as the end point.

8. The method of claim 1, wherein determining the integration end point comprises interpolating between two adjacent corner points to create an interpolated point which corresponds to the echo time and designating the interpolated point as the end point.

9. The method of claim 1, further comprising determining one or more sign change points such that each sign change point corresponds to a refocusing pulse generated by the magnetic resonance imaging technique.

10. The method of claim 9, wherein determining each sign change point comprises designating the corner point which corresponds to the refocusing pulse as the sign change point.

11. The method of claim 9, wherein determining each sign change point comprises interpolating between two adjacent corner points to create an interpolated point which corresponds to the refocusing pulse and designating the interpolated point as the sign change point.

12. The method of claim 1, wherein generating the area sum is accomplished by numerical integration.

13. The method of claim 1, further comprising discarding a plurality of non-corner points.

14. The method of claim 1, wherein generating the area sum comprises determining a k-space trajectory at each corner point using numerical integration on each linear waveform segment.

15. The method of claim 14, further comprising determining one or more sign change points such that each sign change point corresponds to a refocusing pulse generated by the magnetic resonance imaging technique.

16. The method of claim 15, wherein determining each sign change point comprises designating the corner point which corresponds to the refocusing pulse as the sign change point.

17. The method of claim 15, wherein determining each sign change point comprises interpolating between two adjacent corner points to create an interpolated point which corresponds to the refocusing pulse and designating the interpolated point as the sign change point.

18. The method of claim 15, wherein generating the area sum comprises changing signs at each sign change point during numerical integration.

19. The method of claim 18, wherein calculating the b-value comprises calculating a cumulative area of the square of the k-space trajectory.

20. The method of claim 14, wherein calculating the b-value comprises calculating a cumulative area of the square of the k-space trajectory.

21. A method for correcting one or more b-values comprising:

scaling an initial amplitude of a gradient diffusion lobe on a gradient axis to produce a scaled amplitude;

calculating a b-value for the axis using the scaled amplitude;

comparing the b-value with a prescribed b-value; and

adjusting the scaled amplitude if the b-value is outside a given tolerance of the prescribed b-value such that the adjustment is performed iteratively until the b-value is not outside the given tolerance of the prescribed b-value.

22. The method of claim 21, wherein the initial amplitude is scaled by a square root of a ratio of the prescribed b-value to an initial b-value.

23. The method of claim 22, wherein the initial b-value is calculated by integrating one or more linear waveform segments defined by two or more corner points of a gradient waveform.

24. The method of claim 21, wherein adjusting the scaled amplitude comprises a numeric search.

25. The method of claim 24, further comprising determining a minimum diffusion amplitude point which is initially less than the scaled amplitude and a maximum diffusion amplitude point which is initially greater than the scaled amplitude, wherein an intermediate point initially equals the scaled amplitude.

26. The method of claim 25, wherein the minimum diffusion amplitude point and the maximum diffusion amplitude point are equidistant from the intermediate point.

27. The method of claim 25, further comprising setting the minimum diffusion amplitude equal to the intermediate point if the b-value is less than the prescribed b-value minus the given tolerance but setting the maximum diffusion amplitude equal to the intermediate point if the b-value is greater that the prescribed b-value plus the given tolerance.

28. The method of claim 27, further comprising updating the intermediate point according to a mathematical function after setting the minimum diffusion amplitude or the maximum diffusion amplitude.

29. The method of claim 28, wherein the mathematical function updates the intermediate point to equal the (minimum amplitude point)+(0.5)(maximum amplitude point−

minimum amplitude point).

30. The method of claim 28, further comprising recalculating the b-value based upon the intermediate point after updating.

31. A magnetic resonance imaging system which calculates a b-matrix consisting of one or more b-values comprising:

a magnetic resonance imaging scanner configured to collect one or more gradient waveforms comprising a plurality of timing data points and one or more waveform data points;

an analysis circuit configured to generate an output, to receive the gradient waveforms, and to calculate one or more b-values from the gradient waveforms such that the analysis circuit integrates one or more linear waveform segments comprising a gradient waveform to generate an area sum such that the actual b-values are calculated based upon the area sums;

a display circuit configured to receive the output and to preferably display the output via a display device; and

a control circuit configured to control the operation of the magnetic resonance imaging scanner, the analysis circuit, and the display circuit.

32. The system of claim 31, wherein the analysis circuit determines two or more corner points, an integration start point, an integration end point, one or more sign change points, and one or more linear waveform segments.

33. The system of claim 32, wherein the analysis circuit determines the corner points by taking the derivative of the gradient waveform.

34. The system of claim 33, wherein the analysis circuit integrates using numerical integration.

35. The system of claim 31, further comprising a memory circuit configured to receive the gradient waveform such that the analysis circuit receives the gradient waveform from the memory circuit.

36. The system of claim 31, further comprising a memory circuit configured to receive the output such that the display circuit receives the output from the memory circuit.

37. The system of claim 31, wherein the b-matrix is a 3×

3 matrix.

38. The system of claim 37, wherein the b-values of the b-matrix are the one or more actual b-values calculated by the analysis circuit.

39. The system of claim 31, wherein the b-matrix is a 3×

1 matrix.

40. The system of claim 39, wherein the b-value of the b-matrix is the actual b-value calculated by the analysis circuit.

41. A device for correcting one or more b-values comprising:

a magnetic resonance imaging scanner configured to collect a gradient waveform comprising a plurality of timing data points and one or more waveform data points;

an analysis circuit configured to form an output and to generate a scaled amplitude from an initial amplitude of a gradient diffusion lobe and to iteratively adjust the scaled amplitude if a calculated b-value is outside a given tolerance to a prescribed b-value;

a display circuit configured to receive the output and to preferably display the output via a display device; and

a control circuit configured to control the operation of the magnetic resonance imaging scanner, the analysis circuit, and the display circuit.

42. The device of claim 41, further comprising a memory circuit configured to receive the gradient waveform such that the analysis circuit receives the gradient waveform from the memory circuit.

43. The device of claim 41, further comprising a memory circuit configured to receive the output such that the display circuit receives the output from the memory circuit.

44. The device of claim 41, wherein the analysis circuit scales the initial amplitude by a square root of a ratio of the prescribed b-value to an initial b-value.

45. The device of claim 41, wherein the analysis circuit iteratively adjusts the scaled amplitude by performing one or more numeric searches.

46. A method for determining a b-matrix of one or more b-values comprising:

calculating each b-value, the act of calculating comprising;

generating a gradient waveform comprising a plurality of timing data points and one or more waveform data points using a magnetic resonance imaging technique;

calculating two or more corner points associated with the gradient waveform;

determining an integration start point such that the start point corresponds to an excitation pulse generated by the magnetic resonance imaging technique;

determining an integration end point such that the end point corresponds to an echo time resulting from the excitation pulse;

generating an area sum by integrating one or more linear waveform segments between the start point and the end point wherein each linear waveform segment is defined by a line connecting adjacent corner points; and

calculating the b-value using the area sum.

47. The method of claim 46, wherein the b-matrix is a 3×

3 matrix.

48. The method of claim 46, wherein the b-matrix is a 3×

1 matrix.

49. The method of claim 46, wherein the act of calculating a b-value further comprises the act of determining one or more inversion points corresponding to one or more refocusing pulses.

50. The method of claim 46, wherein calculating the two or more corner points comprises taking a derivative of the gradient waveform.

51. The method of claim 46, wherein the act of calculating a b-value further comprises determining one or more sign change points such that each sign change point corresponds to a refocusing pulse generated by the magnetic resonance imaging technique.

52. The method of claim 51, wherein generating the area sum comprises changing signs at each sign change point during numerical integration.

53. The method of claim 46, wherein generating the area sum is accomplished by numerical integration.

54. The method of claim 46, wherein generating the area sum comprises determining a k-space trajectory at each corner point using numerical integration on each linear waveform segment.

55. The method of claim 54, wherein calculating the b-value comprises calculating a cumulative area of the square of the k-space trajectory.