Highly accelerated Point-Spread Function mapping based on Finite Rate of Innovation for EPI distortion correction.

Hybrid MR-PET scans enable acquisition of both types of images within a session. Despite consistent subject positioning, image alignment can still be challenging. Functional MR images rely on echo-planar imaging (EPI) and present geometric distortions due to static B0 field inhomogeneities. Direct B0 [1] and Point Spread Function (PSF) mapping [2] (Figure 1) have been proposed for distortion correction. The PSF method is more robust [3], but acquisition times are long even with previous acceleration approaches [4-7].

Figure 1

The PSF mapping sequence is repeated while varying the amplitude of the prewinder phase encode gradient (in red).

We used the Finite Rate of Innovation (FRI) framework [8, 9] to detect the PSF peak position to sub-pixel precision using as few ks samples as possible.

Images were acquired on a 3.0T Philips Achieva: 2.5×2.5×4.0 mm3, matrix 96×95, 95 ks steps with under-sampling retrospectively performed. PSF peak position was estimated using both the fully-sampled (zero-filled by a factor of 1000) and highly-undersampled data.

To determine the PSF peak location, pattern matching was performed. The signal measured at each spatial location was compared to a predicted signal pattern accounting for the ks sampling scheme and the search progressively refined up to the intended precision. EPI images were undistorted as in [10].

Figure 2 shows example PSF peak shift maps (relative to expected undistorted positions). Figure 3A shows the original EPI slice matching Figure 2, and Figures 3B-E show the corrected images using each displacement map. Comparison with the GE image (Figure 3F) confirms accurate geometrical corrections.

Figure 2

Maps of the PSF shifts (in pixels) estimated with: A0 2; B) 3; C) 4; D) all 95 ks

Figure 3

EPI image (with corresponding outer contour in yellow); Undistorted EPI images estimating the displacement field from: B) 2; C) 3; D) 4 and E) all 95 ks samples and F) GE image (corresponding outer contour in yellow propagated to all undistorted EPI images)

Using the proposed approach the position of the PSF peak can be estimated from a very small number of samples. In the future distortion map estimation could easily be incorporated into standard preparation phases. Making distortion correction of EPI images more practical would facilitate combining functional PET and MR information as well as structural connectivity information from diffusion-weighted MR images.