Image registration and principal component analysis based multi-baseline phase correction method for proton resonance frequency thermometry
A method for phase correction in proton resonance frequency (PRF) thermometry application includes acquiring a series of magnetic resonance (MR) images comprising a first MR image and plurality of subsequent MR images depicting an anatomical area of interest. The MR images are acquired while tissue in the anatomical area of interest is undergoing a temperature change. Each subsequent MR image is registered to the first MR image to yield a plurality of registered images. A plurality of basis images are computed from the registered images using Principal Component Analysis (PCA). The basis images are used to remove motion-related phase changes from a second series of MR images, thereby yielding a motion corrected second series of MR images. One or more temperature difference maps are generated that depict a relative temperature change for the tissue in the anatomical area of interest based on the motion corrected second series.
1. A method for phase correction in proton resonance frequency (PRF) thermometry applications, the method comprising:
acquiring a series of magnetic resonance (MR) images comprising a first MR image and plurality of subsequent MR images depicting an anatomical area of interest, wherein the MR images are acquired while tissue in the anatomical area of interest is undergoing a temperature change;
registering each subsequent MR image to the first MR image to yield a plurality of registered images;
computing a plurality of basis images from the registered images using Principal Component Analysis (PCA);
using the basis images to remove motion-related phase changes from a second series of MR images, thereby yielding a motion corrected second series of MR images; and
generating one or more temperature difference maps depicting a relative temperature change for the tissue in the anatomical area of interest based on the motion corrected second series.
2. The method of claim 1 , wherein the series of MR images are acquired using a 2D or 3D single echo gradient-recalled echo (GRE) sequence.
3. The method of claim 1 , wherein the series of MR images are acquired using a 3D segmented echo-planar imaging (EPI) sequence.
4. The method of claim 1 , wherein the series of MR images are single shot EPI images with or without navigation.
5. The method of claim 1 , wherein the basis images are computed using a process comprising:
calculating a complex mean of the registered images;
computing a signal-to-noise ratio (SNR)-based spatial mask on the registered images using the complex mean;
performing PCA over the spatial mask to yield a plurality of eigen images; and
retaining a subset of the eigen images as the basis images.
6. The method of claim 5 , wherein the eigen images are computed by performing PCA on complex phase difference between the registered images and the complex mean.
7. The method of claim 1 , wherein the basis images are used to remove the motion-related phase changes from the second MR images according to a process comprising:
computing a phase difference between each registered image and a corresponding previous image in the second series of MR images;
generating the motion corrected second series of MR images by removing a projection of each complex phase difference onto a subspace spanned by the plurality of basis images.
8. A system for phase correction in proton resonance frequency (PRF) thermometry applications, the system comprising:
a plurality of RF coils configured to acquire a series of magnetic resonance (MR) images comprising a first MR image and plurality of subsequent MR images depicting an anatomical area of interest, wherein the MR images are acquired while tissue in the anatomical area of interest is undergoing a temperature change;
one or more displays;
one or more computers connected to the RF coils, wherein the computers are configured to:
register each subsequent MR image to the first MR image to yield a plurality of registered images;
compute a plurality of basis images from the registered images using Principal Component Analysis (PCA);
use the basis images to remove motion-related phase changes from a second series of MR images, thereby yielding a motion corrected second series of MR images; and
generate one or more temperature difference maps for presentation on the display, wherein the temperature difference maps depict a relative temperature change for the tissue in the anatomical area of interest based on the motion corrected second series.
9. The system of claim 8 , wherein the series of MR images are acquired using a 2D or 3D single echo gradient-recalled echo (GRE) sequence.
10. The system of claim 8 , wherein the series of MR images are acquired using a 3D segmented echo-planar imaging (EPI) sequence.
11. The system of claim 8 , wherein the series of MR images are single shot EPI images with or without navigation.
12. The system of claim 8 , wherein the basis images are computed using a process comprising:
calculating a complex mean of the registered images;
computing a signal-to-noise ratio (SNR)-based spatial mask on the registered images using the complex mean;
performing PCA over the spatial mask to yield a plurality of eigen images; and
retaining a subset of the eigen images as the basis images.
13. The system of claim 12 , wherein the eigen images are computed by performing PCA on complex phase difference between the registered images and the complex mean.
14. The system of claim 8 , wherein the basis images are used to remove the motion-related phase changes from the second MR images according to a process comprising:
computing a phase difference between each registered image and a corresponding previous image in the second series of MR images;
generating the motion corrected second series of MR images by removing a projection of each complex phase difference onto a subspace spanned by the plurality of basis images.
15. An article of manufacture for phase correction in proton resonance frequency (PRF) thermometry applications, the article of manufacture comprising a non-transitory, tangible computer-readable medium holding computer-executable instructions for performing a method comprising:
receiving a series of magnetic resonance (MR) images comprising a first MR image and plurality of subsequent MR images depicting an anatomical area of interest, wherein the MR images are acquired while tissue in the anatomical area of interest is undergoing a temperature change;
registering each subsequent MR image to the first MR image to yield a plurality of registered images;
computing a plurality of basis images from the registered images using Principal Component Analysis (PCA);
using the basis images to remove motion-related phase changes from a second series of MR images, thereby yielding a motion corrected second series of MR images; and
generating one or more temperature difference maps depicting a relative temperature change for the tissue in the anatomical area of interest based on the motion corrected second series.
16. The article of manufacture of claim 15 , wherein the basis images are computed using a process comprising:
calculating a complex mean of the registered images;
computing a signal-to-noise ratio (SNR)-based spatial mask on the registered images using the complex mean;
performing PCA over the spatial mask to yield a plurality of eigen images; and
retaining a subset of the eigen images as the basis images.
17. The article of manufacture of claim 16 , wherein the eigen images are computed by performing PCA on complex phase difference between the registered images and the complex mean.
18. The article of manufacture of claim 15 , wherein the basis images are used to remove the motion-related phase changes from the second MR images according to a process comprising:
computing a phase difference between each registered image and a corresponding previous image in the second series of MR images;
generating the motion corrected second series of MR images by removing a projection of each complex phase difference onto a subspace spanned by the plurality of basis images.