Systems and methods for determining equilibrium CEST effect for multi-slice CEST imaging using quasi-steady state (QUASS) analysis
The devices, systems, and methods can overcome the impact of scan protocols on experimental measurements to reconstruct the CEST equilibrium effect and can therefore improve the quantification of the underlying CEST contrast mechanism. The methods may include determining a spinlock relaxation rate for each RF offset for at least the first signal average so that the normalized spectrum for each RF offset and each slice first signal average corresponds to at least each transient state associated with each saturated slice readout normalized by at least each control state associated with each control slice readout based on the spinlock relaxation rate. The method may further include reconstructing CEST equilibrium effect for each RF offset using the respective spinlock relaxation rate for the RF offset. The method may include determining one or more quantitative information using the CEST equilibrium effect.
1 . A method for determining a chemical exchange saturation transfer (CEST) equilibrium effect, comprising:
acquiring based on one or more parameters, radiofrequency (RF) offset, and the pulse sequence for one or more rounds: (i) two or more control slice readouts of magnetic resonance imaging (MRI) CEST imaging based on parameters and pulse sequence, ii) two or more saturated slice readouts of the MRI CEST imaging based on one or more parameters, RF offset, and the pulse sequence, iii) one or more signal averages for each slice readout and offset, the one or more signal averages including a first signal average for each slice readout and offset, and iv) normalized spectrum for each RF offset, the one or more parameters including relaxation delay;
the two or more control slice readouts including a first control slice readout and one or more additional control slice readouts; and
the two or more saturated slice readouts including a first saturated slice readout and one or more additional saturated slice readouts;
determining a saturation time for each saturated slice readout for each signal average;
determining a control relaxation time for each control slice readout using the respective saturation time and the relaxation delay for each signal average;
determining a saturated relaxation time for each saturated slice readout using the relaxation delay for each signal average;
determining a spinlock relaxation rate for each RF offset for at least the first signal average, using the one or more parameters, each saturated time for each saturated slice readout, each control relaxation time for each control slice readout, and each saturated relaxation time for each saturated slice readout for the first signal average, so that the normalized spectrum for each RF offset and each slice first signal average corresponds to at least each transient state associated with each saturated slice readout normalized by at least each control state associated with each control slice readout based on the spinlock relaxation rate;
reconstructing CEST equilibrium effect for each RF offset using the respective spinlock relaxation rate for the RF offset; and
determining one or more quantitative information using the CEST equilibrium effect.
2 . The method of claim 1 , wherein each transient state includes an equilibrium signal.
3 . The method of claim 1 ,
wherein the one or more signal averages, including one or more additional signal averages for each slice readout and offset; and
wherein the acquiring includes acquiring one or more preparation modules and readout time associated with each signal average;
the method further comprising:
determining a second saturation time for each saturated slice readout for each additional signal average;
determining a second control relaxation time for each control slice readout using the respective saturation time for each additional signal average; and
determining a second saturated relaxation time for each saturated slice readout using the preparation module and the readout time for each additional signal average;
wherein the determining the spinlock relaxation rate for each RF offset further uses the second saturated time for each saturated slice readout of each additional signal average, the second control relaxation time for each control slice readout of each additional signal average, the second saturated relaxation time for each saturated slice readout of each additional signal average, so that the normalized spectrum for each RF offset corresponds to the sum of at least the first transient state associated with the first saturated readout and a second transient state associated with the one or additional saturated readouts normalized by the sum of at least the first control state associated with the first control readout and a second transient state associated with the one or additional control readouts.
4 . The method of claim 1 , wherein the sequence is a multi-slice gradient echo sequence and the one or more parameters includes RF flip angle.
5 . The method of claim 1 , wherein the readouts are for at least one offset.
6 . The method of claim 1 , further comprising:
determining one or more quantitative information based on the CEST equilibrium signal.
7 . The method of claim 6 , further comprising:
reporting the quantification method.
8 . The method of claim 1 , wherein the sequence is a multi-slice spin echo sequence.
9 . A system for determining a chemical exchange saturation transfer (CEST) equilibrium effect, comprising:
at least one computing device comprising a processor and a memory; and
machine-readable instructions stored in the memory that, when executed by the processor, cause the computing device to at least:
acquire based on one or more parameters, radiofrequency (RF) offset, and the pulse sequence for one or more rounds: (i) two or more control slice readouts of magnetic resonance imaging (MRI) CEST imaging based on parameters and pulse sequence, ii) two or more saturated slice readouts of the MRI CEST imaging based on one or more parameters, RF offset, and the pulse sequence, iii) one or more signal averages for each slice readout and offset, the one or more signal averages including a first signal average for each slice readout and offset, and iv) normalized spectrum for each RF offset, the one or more parameters including relaxation delay;
the two or more control slice readouts including a first control slice readout and one or more additional control slice readouts; and
the two or more saturated slice readouts including a first saturated slice readout and one or more additional saturated slice readouts;
determine a saturation time for each saturated slice readout for each signal average;
determine a control relaxation time for each control slice readout using the respective saturation time and the relaxation delay for each signal average;
determine a saturated relaxation time for each saturated slice readout using the relaxation delay for each signal average;
determine a spinlock relaxation rate for each RF offset for at least the first signal average, using the one or more parameters, each saturated time for each saturated slice readout, each control relaxation time for each control slice readout, and each saturated relaxation time for each saturated slice readout for the first signal average, so that the normalized spectrum for each RF offset and each slice first signal average corresponds to at least each transient state associated with each saturated slice readout normalized by at least each control state associated with each control slice readout based on the spinlock relaxation rate;
reconstruct CEST equilibrium effect for each RF offset using the respective spinlock relaxation rate for the RF offset; and
determine one or more quantitative information using the CEST equilibrium effect.
10 . The system of claim 9 , wherein each transient state includes an equilibrium signal.
11 . The system of claim 9 ,
wherein the one or more signal averages, including one or more additional signal averages for each slice readout and offset;
wherein the acquiring includes acquiring one or more preparation modules and readout time associated with each signal average; and
wherein the machine-readable instructions, when executed, further cause the at least one computing device to at least:
determine a second saturation time for each saturated slice readout for each additional signal average;
determine a second control relaxation time for each control slice readout using the respective saturation time for each additional signal average; and
determine a second saturated relaxation time for each saturated slice readout using the preparation module and the readout time for each additional signal average;
wherein the determination of the spinlock relaxation rate for each RF offset further uses the second saturated time for each saturated slice readout of each additional signal average, the second control relaxation time for each control slice readout of each additional signal average, the second saturated relaxation time for each saturated slice readout of each additional signal average, so that the normalized spectrum for each RF offset corresponds to the sum of at least the first transient state associated with the first saturated readout and a second transient state associated with the one or additional saturated readouts normalized by the sum of at least the first control state associated with the first control readout and a second transient state associated with the one or additional control readouts.
12 . The system of claim 9 , wherein the sequence is a multi-slice gradient echo sequence and the one or more parameters includes RF flip angle.
13 . The system of claim 9 , wherein the readouts are for at least one offset.
14 . The system of claim 9 , wherein the machine-readable instructions, when executed, further cause the at least one computing device to at least:
determine one or more quantitative information based on the CEST equilibrium signal.
15 . The system of claim 14 , wherein the machine-readable instructions, when executed, further cause the at least one computing device to at least:
report the quantification method.
16 . The system of claim 9 , wherein the sequence is a multi-slice spin echo sequence.