METHODS AND SYSTEMS FOR CHEMICAL EXCHANGE SATURATION TRANSFER SIGNAL MATCHING

• US 20190324104A1
• Filed: 04/23/2018
• Published: 10/24/2019
• Est. Priority Date: 04/23/2018
• Status: Active Grant

• Current Assignee: Cedars-Sinai Medical Center (Cedars Sinai Health System)
• Original Assignee: Cedars-Sinai Medical Center (Cedars Sinai Health System)
• Inventors: Li, Debiao, Zhou, Zhengwei
• Granted Patent: US 10,859,656 B2
• Time in Patent Office: Days
• CPC Class Codes: G01N 24/08 by using nuclear magnetic resonance G01N24/12 takes precedence
• CPC Class Codes: G01R 33/4838 using spatially selective suppression or saturation of MR signals
• CPC Class Codes: G01R 33/5605 by transferring coherence or polarization from a spin species to another, e.g. creating magnetization transfer contrast [MTC], polarization transfer using nuclear Overhauser enhancement [NOE]

Abstract

A method of determining an exchange rate between a solute material and a solvent material in chemical exchange with each other in a sample comprises performing a plurality of CEST scans on the sample to determine a reference data signal and a label data signal; simulating a plurality of CEST scans on the sample to obtain a reference dictionary signal and a plurality of label dictionary signals, each label dictionary signal corresponding to a candidate exchange rate; correcting the reference dictionary signal and each label dictionary signal for a magnetization transfer effect; determining a final data signal based on the reference data signal and the label data signal; determining a plurality of final dictionary signals based on the reference dictionary signal and the plurality of label dictionary signals; and determining the exchange rate by matching the final data signal to one of the plurality of final dictionary signals.

Claims

1. A method of determining an exchange rate between a solute material and a solvent material in chemical exchange with each other in a sample, the method comprising:

receiving data indicative of a plurality of reference CEST scans performed on the sample and a plurality of label CEST scans performed on the sample;

determining a reference data signal and a label data signal based on the received data;

simulating a plurality of reference CEST scans on the sample to obtain a reference dictionary signal;

simulating a plurality of groups of label CEST scans on the sample to obtain a plurality of label dictionary signals, each of the plurality of label dictionary signals corresponding to a respective one of a plurality of candidate exchange rates between the solvent material and the solute material;

correcting the reference dictionary signal and each of the plurality of label dictionary signals for a magnetization transfer effect between the solvent material and the solute material;

determining a final data signal based on the reference data signal and the label data signal;

determining a plurality of final dictionary signals based on the reference dictionary signal and the plurality of label dictionary signals; and

determining the exchange rate between the solute material and the solvent material in the sample by matching the final data signal to one of the plurality of final dictionary signals.

2. The method of claim 1, further comprising creating an acquisition index containing a plurality of entries n, each of the n entries of the acquisition index specifying a respective saturation power and a respective saturation time.

3. The method of claim 2, wherein an nth one of the plurality of performed reference CEST scans and an nth one of the plurality of simulated reference CEST scans each correspond to the nth one of the plurality of entries in the acquisition index such that the nth one of the plurality of performed reference CEST scans is performed at the saturation power and the saturation time specified by the nth one of the plurality of entries in the acquisition index and the nth one of the plurality of simulated reference CEST scans is simulated at the saturation power and the saturation time specified by the nth one of the plurality of entries in the acquisition index.

4. The method of claim 3, wherein the reference data signal and the reference dictionary signal each contain n data points, the nth data point of the reference dictionary signal being a simulation of the nth data point of the reference data signal.

5. The method of claim 2, wherein an nth one of the plurality of performed label CEST scans and an nth one of the plurality of simulated label CEST scans of each of the plurality of groups of simulated label CEST scans each correspond to the nth one of the plurality of entries in the acquisition index such that the nth one of the plurality of performed label CEST scans is performed at the saturation power and the saturation time specified by the nth one of the plurality of entries in the acquisition index and the nth one of the plurality of simulated label CEST scans of each of the plurality of groups of simulated label CEST scans is simulated at the saturation power and the saturation time specified by the nth one of the plurality of entries in the acquisition index.

6. The method of claim 5, wherein the label data signal and each of the plurality of label dictionary signals each contain n data points, the nth data point of each of the plurality of label dictionary signals being a simulation of the nth data point of the label data signal at the respective one of the plurality of candidate exchange rates between the solute material and the solvent material.

7. The method of claim 2, wherein the saturation power of each of the n entries in the acquisition index is between about 0.0 μ

T and about 2.0 μ

T.

8. The method of claim 2, wherein the saturation time of each of the n entries in the acquisition index is between about 0.0 seconds and about 1.5 seconds.

9. The method of claim 2, wherein the saturation power and the saturation time of each of the n entries are selected on a pseudorandom basis.

10. The method of claim 1, wherein each of the plurality of performed reference CEST scans is performed at a saturation frequency offset of about −

2.0 ppm, and wherein each of the plurality of performed label CEST scans is performed at a saturation frequency offset of about 2.0 ppm.

11. The method of claim 1, wherein each of the plurality of simulated reference CEST scans is simulated at a saturation frequency offset of about −

2.0 ppm, and wherein each of the plurality of simulated label CEST scans is performed at a saturation frequency offset of about 2.0 ppm.

12. The method of claim 1, wherein each of the plurality of performed reference CEST scans is performed at one of a plurality of saturation frequency offsets to obtain a plurality of partial reference data signals, wherein the plurality of partial reference data signals are interpolated to obtain the reference data signal.

13. The method of claim 1, further comprising:

normalizing the reference data signal and the label data signal to the reference data signal; and

normalizing the reference dictionary signal and each of the plurality of label dictionary signals to the reference dictionary signal.

14. The method of claim 1, wherein correcting the reference dictionary signal and each of the plurality of label dictionary signals for the magnetization transfer effect comprises:

comparing the reference data signal and the reference dictionary signal to obtain a magnetization transfer signal;

subtracting the magnetization transfer signal from the reference dictionary signal; and

subtracting the magnetization transfer signal from each of the plurality of label dictionary signals.

15. The method of claim 1, wherein the final data signal is determined by subtracting the label data signal from reference data signal.

16. The method of claim 1, wherein each of final dictionary signals are determined by subtracting each of label dictionary signals from the reference dictionary signal.

17. The method of claim 1, wherein the determining the exchange rate between the solute material and the solvent material comprises:

determining a value of a vector dot product between the final data signal and each of the plurality of final dictionary signals; and

selecting the one of the plurality of final dictionary signals having the largest value of the vector dot product.

18. A method of creating an exchange rate-weighted cross-sectional image of a sample containing a solute material and a solvent material that are in chemical exchange with each other, the method comprising:

defining a cross-sectional area of the sample, the cross-sectional area containing a plurality of pixels;

determining an exchange rate between the solute material and the solvent material at each of the plurality of pixels, the determining the exchange rate at a respective pixel further comprising;

receiving data indicative of a plurality of CEST scans performed on the sample;

determining a data signal based on the received data;

simulating a plurality of groups of CEST scans to obtain a plurality of dictionary signals, each of the plurality of dictionary signals corresponding to a respective one of a plurality of candidate exchange rates between the solute material and the solvent material; and

matching the data signal to one of the plurality of dictionary signals to determine an exchange rate between the solute material and the solvent material at the respective pixel; and

mapping the exchange rate of the solute material to the cross-sectional area by correlating the exchange rate at each of the plurality of pixels to a color of each of the plurality of pixels.

19. The method of claim 18, wherein each of the plurality of performed CEST scans is performed at a pseudorandom saturation power and a pseudorandom saturation time, and wherein each of the plurality of simulated CEST scans is simulated at a pseudorandom saturation power and a pseudorandom saturation time.

20. A system for determining an exchange rate between a solute material and a solvent material in chemical exchange with each other in a sample, the system comprising:

at least one processor; and

at least one memory containing processor-executable instructions, which when executed by the at least one processor, cause the at least one processor to;

receive data indicative of a plurality of reference CEST scans performed on the sample and a plurality of label CEST scans performed on the sample;

determine a reference data signal and a label data signal based on the received data;

simulate a plurality of reference CEST scans on the sample to obtain a reference dictionary signal;

simulate a plurality of groups of label CEST scans on the sample to obtain a plurality of label dictionary signals, each of the plurality of label dictionary signals corresponding to a respective one of a plurality of candidate exchange rates between the solute material and the solvent material;

correct the reference dictionary signal and each of the plurality of label dictionary signals for a magnetization transfer effect between the solute material and the solvent material;

determine a final data signal based on the reference data signal and the label data signal;

determine a plurality of final dictionary signals based on the reference dictionary signal and the plurality of label dictionary signals; and

determine the exchange rate between the solute material and the solvent material in the sample by matching the final data signal to one of the plurality of final dictionary signals.