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

The Phase Study of PSF Model in MR.

Improved <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:math>-<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:math> PCA Algorithm Using Artificial Sparsity in Dynamic MRI

Sparsity-constrained SENSE reconstruction: An efficient implementation using a fast composite splitting algorithm

A Robust Algorithm for High-Resolution Dynamic Mri Based on the Partially Separable Functions Model

Linear Dynamic Sparse Modelling for Functional MR Imaging.

A Simulation Study of Dynamic MRI Based on Partially Separable Functions and Keyhole Techniques

PSF model simulation study using a cardiac phantom.

High-resolution imaging using partial separability of spatiotemporal signals with a novel data sampling scheme

Compressed Sensing-Based Simultaneous Recovery of Magnitude and Phase MR Images via Dual Trigonometric Sparsity

A Kernel Method for Higher Temporal Resolution Mri Using the Partial Separability (ps) Model

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

High resolution magnetic resonance thermometry based on the partial separable function model

MRI with phaseless encoding

High spatiotemporal resolution fMRI using partial separability model

Three-dimensional spatially resolved phase graph framework

Getting the phase consistent: The importance of phase description in balanced steady-state free precession MRI of multicompartment systems

High resolution dynamic cardiac MRI using partial separability of spatiotemporal signals with a novel sampling scheme

Accelerated partial separable model using dimension-reduced optimization technique for ultra-fast cardiac MRI

Accelerated phase contrast flow imaging with direct complex difference reconstruction.

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

Accelerating PS Model-Based Dynamic Cardiac MRI Using Compressed Sensing.