Monitoring and noise masking of thermal therapy

1. A method for thermal monitoring of tissue during thermal therapy, comprising:

prior to initiation of thermal therapy,

obtaining, by processing circuitry, a plurality of magnetic resonance (MR) images of a patient, wherein

the plurality of MR images correspond to a same imaging slice captured over a time period, and

each MR image of the plurality of MR images includes respective phase data,

identifying, by the processing circuitry, a volume in the patient corresponding to a) a region of interest at which to apply thermal therapy and b) surrounding tissue,

analyzing, by the processing circuitry, pixel locations of the phase data of the plurality of MR images corresponding to the volume to identify one or more pixel locations demonstrating excessive phase noise,

determining, by the processing circuitry, a noise mask for the region of interest, wherein the noise mask excludes the one or more pixel locations demonstrating excessive phase noise, and

selecting, from the pixels included in the noise mask, a plurality of reference points, wherein

the plurality of reference points identify locations within the surrounding tissue, and

each reference point of the plurality of reference points corresponds to a particular pixel location of the phase data having a particular phase value; and

initiating monitoring of the thermal therapy, wherein monitoring comprises

receiving a current MR image of the patient,

adjusting, by the processing circuitry, the phase data of the current MR image according to the plurality of reference points, wherein adjusting the phase data of the current MR image comprises

determining, based upon at least a subset of the plurality of reference points without considering phase data corresponding to pixel locations within the region of interest, a plurality of phase drift values to correct for phase drift, and

extrapolating the phase drift values into a phase map, wherein at least a portion of the phase map corresponds to the pixel locations within the region of interest,

applying, by the processing circuitry, the phase map to the current phase data to obtain drift-adjusted phase data, and

calculating, by the processing circuitry, a temperature map using the drift-adjusted phase data, wherein the temperature map corresponds to at least the region of interest.

2. The method of claim 1 , further comprising, after determining the noise mask and prior to initiating thermal therapy, displaying, on a graphical user interface, a rendering of the noise mask overlaying the volume in the patient in a first MR image of the plurality of MR images.

3. The method of claim 2 , wherein selecting the plurality of reference points comprises receiving, via interaction with the graphical user interface, a respective user selection of each of the plurality of reference points from the rendering of the noise mask.

4. The method of claim 1 , wherein selecting the plurality of reference points comprises automatically selecting, by the processing circuitry independently of user selection, at least a portion of the plurality of reference points.

5. The method of claim 4 , further comprising receiving user input confirming acceptance of the automatically selected reference points of the plurality of reference points.

6. The method of claim 1 , wherein selecting the plurality of reference points comprises verifying, by the processing circuitry, that each reference point of the plurality of reference points corresponds to a respective tissue region of a same tissue type as the remaining reference points.

7. The method of claim 1 , wherein the plurality of MR images comprises at least three MR images.

8. The method of claim 1 , further comprising, prior to initiation of thermal therapy:

validating, by the processing circuitry, stability of drift-adjusted phase data, wherein

validating the stability of drift-adjusted phase data comprises

obtaining a present MRI image of the volume;

computing, at every location within the volume, a respective temperature difference from a reference tissue temperature, wherein computing the respective temperature difference comprises adjusting the phase data of the present MR image according to the plurality of reference points;

applying, to the respective temperature differences, a temperature variance threshold to identify a plurality of noisy locations; and

illustrating, on a graphical user interface, a rendering of a temperature mask overlaying the volume in the patient, wherein the noisy locations are specified by visual indications.

9. The method of claim 8 , wherein the noisy locations are indicated by transparent portions of a color overlay.

10. A system for thermal monitoring of tissue during thermal therapy, comprising:

a processor in communication with a display and a magnetic resonance (MR) imaging device; and

a memory having instructions stored thereon, wherein the instructions, when executed by the processor, cause the processor to:

prior to initiation of thermal therapy,

obtain, from the MR imaging device, a plurality of magnetic resonance (MR) images of a patient, wherein

the plurality of MR images correspond to a same imaging slice over a time period, and

each MR image of the plurality of MR images includes respective magnetic resonance (MR) phase data,

identify a volume in the patient corresponding to a) a region of interest at which to apply thermal therapy and b) surrounding tissue, wherein the plurality of MR images comprises respective magnetic resonance (MR) image data representing a portion of the volume including a portion of the region of interest, and

determine a plurality of reference points, wherein

the plurality of reference points identify locations within the surrounding tissue, and

each reference point of the plurality of reference points corresponds to a particular pixel location of the phase data having a particular phase value; and

initiate monitoring of thermal therapy, wherein monitoring comprises

receiving, from the MR imaging device, a current MR image of the patient,

calculating a plurality of phase drift values associated with the plurality of reference points, wherein the phase drift values are based upon phase values of each respective reference point of the plurality of reference points over successive magnetic resonance (MR) phase data of the current MR image and a former MR image without considering phase data corresponding to pixel locations within the region of interest,

extrapolating the phase drift values of the plurality of reference points into a phase map, wherein at least a portion of the phase map corresponds to the pixel locations within the region of interest,

applying the phase map to the MR phase data of the current MR image to obtain drift-compensated phase data,

calculating a temperature map based upon the drift-compensated phase data, wherein the temperature map corresponds to at least the region of interest, and

presenting, on the display, a rendering of the MR image including a graphical representation of temperatures in the region of interest according to the temperature map.

11. The system of claim 10 , wherein the graphical representation of temperatures comprises a color overlay.

12. The system of claim 10 , further comprising, prior to determining the plurality of reference points:

analyzing the respective MR phase data of each MR image of the plurality of MR images to identify a plurality of noisy data points; and

marking the plurality of noisy data points as ineligible for selection as reference points.

13. The system of claim 12 , wherein identifying a particular noisy data point of the plurality of noisy data points comprises determining whether a difference of phase values of the particular noisy data point between two MR images of the plurality of MR images exceeds a threshold level.

14. The system of claim 10 , wherein calculating the temperature map comprises multiplying the drift-compensated phase data by a constant and adding a baseline temperature of the volume.

15. The system of claim 14 , wherein the instructions, when executed by the processor, cause the processor to receive, through a user interface, the baseline temperature.

16. The method of claim 8 , further comprising;

receiving, through a user interface responsive to illustrating the rendering of the temperature mask, user input modifying the plurality of reference points, wherein the user input corresponds to at least one of a) a deletion of a first reference point of the plurality of reference points, and b) an addition of a new reference point to the plurality of reference points; and

repeating validation of stability of the drift-adjusted phase data using the modified plurality of reference points.

17. The method of claim 8 , further comprising receiving, through a user interface prior to validating the stability of the drift-adjusted phase data, the reference tissue temperature.

18. The system of claim 10 , wherein the plurality of reference points comprises fewer than twenty-five reference points.

19. The system of claim 10 , wherein:

the plurality of reference points comprises at least eight reference points; and

extrapolating the phase drift values into the phase map comprises applying the phase drift values to the phase data using a polynomial surface fit.

20. The method of claim 6 , wherein the tissue type is brain tissue.