Abstract:Target audience: Researchers on MR thermometry and phase image reconstruction methods. Purpose: Phase information is a focal point in MR thermometry using proton resonance frequency shift (PRFS) method [1]. In order to shorten the imaging duration and accomplish real-time temperature mapping, acceleration methods can be applied to the reconstruction of phase information. Thus, a method to improve the accuracy of phase reconstruction in dynamic scans is highly needed to get the precise real-time temperature estimation. K-t space FOCal Underdetermined System Solver (k-t FOCUSS) is an effective tool for sparse reconstruction in both k-space and time domain using initial estimation and residue encoding [2]. In order to make better use of the sparsity of phase differential information, Compressed Sensing with nonlinear filters such as median filter [3] is combined in our method for exploiting gradient sparsity with spatiotemporal correlations. This combination of k-t FOCUSS and CS with nonlinear filter in phase mapping is a promising tool for real-time temperature monitoring using PRFS according to our preliminary results. Method: A new Compressed Sensing framework combining nonlinear filter (NFCS) [3] and k-t FOCUSS is proposed in our research (see Fig.1). K-t FOCUSS is implemented in our method to detect and exploit sparse information in k-t domain based on the temporal average prediction and residue encoding. Considering the fact that k-t FOCUSS is designed for dynamic processes with periodic motion and temperature change is nonlinear during the whole experiment, NFCS is additionally used to enhance the fidelity of the reconstruction from the under-sampled data. The implementation of this combined method is generally divided into four steps: (1) Using the latest N frames (N = 15 in this research) for batch k-t FOCUSS to get the preliminary reconstructed frames since the batch k-t FOCUSS does not require fully-sampled reference scans at the end of all frames, which meets the requirement of dynamic temperature mapping. In the algorithm of batch k-t FOCUSS (see Fig.1), ρ Δ represents the difference between prediction 0 ρ and under-sampled data. F is the Fourier Transform along y and t directions. (2) Using NFCS for second reconstruction. In NFCS algorithm, U ∇ represents the gradient sparsity with spatiotemporal correlations. NF is the number of frames and l S l l F u y − represents the data fidelity term. Specifically, the latest three frames are reconstructed using a 3D median filter in image domain. (3) Calculating the phase difference map using the reconstruction result from NFCS; (4) Calculating the final temperature map based on the PRFS method using the value of echo time. Results: 30 frames of fully-sampled 2D Cartesian data were acquired using a 3T Philips Achieva scanner during a heating experiment. The parameters for dynamic images were: TR/TE =70ms/5ms, resolution=1.6×1.6×4 mm, FOV=256×256 and flip angle=40°. The composition of the phantom is: 100% water for the outer ring, 20% fat mixed with 80% water for the central area in order to add boundary information and improve the gradient sparsity. The fully acquired data were artificially down-sampled with a reduction factor of 4 using a variable density acquisition trajectory. Fig.2 displays the temperature change curves for three different phase reconstruction methods (2D FOCUSS, batch k-t FOCUSS and batch k-t FOCUSS with NFCS) at R=4. Fig. 3 shows the 5-time brightened error maps of phase reconstruction for the under-sampled data using batch k-t FOCUSS with NFCS and batch k-t FOCUSS at R=4. Table 1 shows the maximum temperature error after an initialization of 10 frames based on the theory of PRFS [1]. Discussion: As shown in Fig.2, the temperature change curve of our proposed method demonstrates the best agreement with the fully-sampled data among all the three methods. Thus, combining k-t FOCUSS and NFCS enables an accurate prediction as well as a robust correction for phase information. The results shown in Fig.3 also indicate that our proposed method can accurately reconstruct phase information at a reduction factor of 4 and is able to exploit the advantages of both k-t FOCUSS and NFCS in dynamic reconstruction especially when prediction might be inaccurate. Table 1 shows that the temperature error of our proposed method is lower than other methods after the initialization process. Conclusions: Phantom experiments have demonstrated that our new Compressed Sensing framework combining k-t FOCUSS and nonlinear filter is superior to other methods in maintaining accurate phase information at a reduction factor of 4. Further study could be conducted on in-vivo experiments to demonstrate the potential of the proposed method in clinical practice. References: [1] Ishihara Y, MRM 1995; [2] Jung H, MRM 2009; [3] Montefusco L, IEEE Trans. Med. Imag. 2011. Table 1 Maximum Temperature Error after an Initialization Process

A Novel Phase Unwrapping Method Based on a Quality Map and Local Polynomial Surface Fitting with Applications in MR Temperature Imaging

A New Particle Swarm Optimization-Based Method for Phase Unwrapping of MRI Data

Region-growing Phase Unwrapping Method Based on Mask

A Reliable Phase Unwrapping Algorithm Based on the Local Fitting Plane and Quality Map

Sparse-Representation-Based Direct MinimumLp-Norm Algorithm for MRI Phase Unwrapping

New Quality Map Based on Modulation for Phase Unwrapping Algorithm

New 3D phase-unwrapping method based on voxel clustering and local polynomial modeling: application to quantitative susceptibility mapping

A Hybrid Phase-Unwrapping Method For Optical Interferometry Based On New Parameter Map And Local Plane Approximation

Novel method to construct a quality map for phase unwrapping based on modulation and the phase gradient

Quality Weight Based Discrete Cosine Transform Phase Unwrapping Algorithm in Optical Interferogram Processing

New Quality Weight Used for Phase Unwrapping in Color Fringe Reflection Method

Fast Phase-unwrapping Algorithms Used for Structural-light Measurement

3d Phase Unwrapping Using Global Expected Phase As A Reference: Application To Mri Global Shimming

Phase Unwrapping with Graph Cuts Optimization and Dual Decomposition Acceleration for 3D High‐resolution MRI Data

The PHU‐NET: A robust phase unwrapping method for MRI based on deep learning

Accelerated Phase Reconstruction Using Sharable Information from Reference Scan and Its Application for Real-time MR Thermometry

A Fast Phase Unwrapping Method Based on Region-Division

Multilevel Quality-Guided Phase Unwrapping Algorithm for Real-Time Three-Dimensional Shape Reconstruction

Towards a Novel Phase Unwrapping Model for Three Dimensional Measurement

Fast Temperature Estimation from Undersampled K-Space with Fully-Sampled Center for MR Guided Microwave Ablation

Accelerated Real-time MR Thermometry Using a New Compressed Sensing Framework of Nonlinear Filter and Kt FOCUSS