Double offsets and powers magnetization transfer ratio for specific magnetization transfer magnetic resonance imaging
Imaging methods for assessing the macromolecular content, such as myelin, are of great interest for understanding brain tissue microstructure, and have shown potentials in diagnosing and prognosing demyelinating diseases. for example. Magnetization transfer (MT) is a MRI contrast mechanism that enables detection of macromolecules. Previously, the MT effect has been analyzed by a semi-quantitative method termed magnetization transfer ratio (MTR) or by a quantitative magnetization transfer (qMT) method. However, because MTR does not have enough sensitivity and specificity to myelin, and qMT takes a very long scan time, their translation into clinical scenarios has been limited. This disclosure describes a MT data analysis metric using double saturation pulse offsets and powers (dopMTR). Simulations and experiments using the systems and methods described in this disclosure show that the dopMTR yields much better sensitivity and specificity to MT effect than the conventional MTR, and requires much less scan time than the qMT.
1. A method for measuring magnetization transfer (MT) of magnetic resonance (MR) signals, the method comprising:
providing a first MT sequence and a second MT sequence, each comprising a radiofrequency (RF) saturation block, wherein:
the first MT sequence comprises a predetermined first RF saturation power and a predetermined first RF saturation frequency offset from water,
the second MT sequence comprises a second RF saturation power and a second RF saturation frequency offset from water,
the second RF saturation power is greater than the first RF saturation power by a factor c, and
the second RF saturation frequency offset is greater than the first saturation frequency offset by the factor c;
acquiring a first MR signal and a second MR signal using the first MT sequence and the second MT sequence; and
measuring the MT via a difference between the first MR signal and the second MR signal.
2. The method of claim 1 , wherein the RF saturation block comprises either
(1) a single continuous wave (CW) saturation pulse followed by imaging readout;
(2) a train of shorter saturation pulses followed by imaging readout;
(3) a train of shorter saturation pulses with interleaved imaging readout; or
(4) a train of shorter saturation pulses for MT saturation followed by another train of shorter saturation pulses with interleaved imaging readout for both MT saturation and data acquisition.
3. The method of claim 1 , wherein the factor c is a rational number.
4. The method of claim 1 , further comprising:
subtracting or inverse subtracting the first MR signal and the second MR signal; and
determining an MT effect using the subtracting or inverse subtracting of the first MR signal and the second MR signal.
5. The method of claim 4 , wherein determining the MT effect comprises using normalization by water longitudinal relaxation time.
6. The method of claim 1 , wherein the first and second MT sequence have a duration between 50 ms and 20 s, and the RF saturation block of the first and second MT sequence is an off-resonance CW-RF saturation block followed by an imaging readout.
7. The method of claim 1 , wherein the first and second MT sequence have a duration between 5 ms and 2 s, and comprise a train of off-resonance saturation pulses followed by an imaging readout.
8. The method of claim 7 , wherein the duration between two adjacent saturation pulses is between 0 s and 1 s.
9. The method of claim 7 , wherein the number of saturation pulses is between 1 and 1000.
10. The method of claim 7 , wherein the shape of the saturation pulses is rectangular, Gaussian, SINC, or any variation of rectangular, Gaussian, or SINC.
11. The method of claim 1 , wherein the first and second MT sequence
have a duration between 5 ms to 2 s, and
comprise a train of off-resonance saturation pulses with interleaved imaging readout between two adjacent saturation pulses.
12. The method of claim 11 , wherein the duration between two adjacent saturation pulses is between 0 s and 1 s.
13. The method of claim 11 , wherein the number of saturation pulses is between 1 and 1000.
14. The method of claim 11 , wherein the shape of the saturation pulses is rectangular, Gaussian, SINC, or any variation of rectangular, Gaussian, and SINC.
15. The method of claim 1 , wherein the first and second MT sequence
have a duration between 5 ms and 2 s, and
comprise a first train of off-resonance saturation pulses followed by a second train of off-resonance saturation pulses with a duration between 5 ms and 2 s with interleaved imaging readout between two adjacent saturation pulses.
16. The method of claim 15 , wherein the duration between the two adjacent saturation pulses is between 0 s and 1 s.
17. The method of claim 15 , wherein the number of saturation pulses is between 1 and 1000.
18. The method of claim 15 , wherein the shape of the saturation pulses is rectangular, Gaussian, SINC, or any variation of rectangular, Gaussian, and SINC.
19. A magnetic resonance imaging scanner that measures the MT of MR signals according to the method of any of the previous claims.