1639

An integrated model-based framework for the correction of signal pile-up and translational offsets in prostate diffusion MRI
Muhammad Usman1, Lebina Kakkar2, Karin Shmueli3, Simon Arridge1, and David Atkinson2

1Centre for Medical Image Computing, Department of Computer Science, University College London, London, United Kingdom, 2Centre for Medical Imaging, Division of Medicine, University College London, London, United Kingdom, London, United Kingdom, 3Medical Physics and Biomedical Engineering, University College London, London, United Kingdom, London, United Kingdom

Synopsis

Prostate diffusion EPI scans suffer from geometric distortions, signal pile-up and signal drop-out due to differences in susceptibility values at interface between prostate and rectal-air. In this work, an integrated model based framework is proposed that can correct for signal pile-up in regions of severe distortions and can compensate for any translational offsets that may exist between different scans. In-vivo validation of the proposed method is done in patients.

Introduction

Prostate diffusion EPI scans suffer from geometric distortions, signal pile-up and signal drop-out due to differences in susceptibility values at interface between prostate and rectal-air. Techniques have been proposed that can correct for geometric distortions using an additional B0 map that measures frequency offset from resonance at different spatial locations. In regions of severe signal pile-up, these techniques are prone to failure in recovering the underlying true signal. Furthermore, due to drifts and inaccuracies in the measurement of the center frequency [1], the reconstructed images for the EPI scans can be shifted in the phase encoding direction. These shifts are not present for the B0 scan that has a higher bandwidth, thus making use of a B0 map difficult. In this work, using a set of EPI data acquired in blip-up and blip-down phase encoding directions, an integrated model based framework is proposed that can correct for signal pile-up in regions of severe distortions and compensate for any translational offsets between different scans. Validation of proposed method is done in two patients.

Method

In presence of translational offset along the phase encoding direction, the distorted k-space Y can be related to ideal image f via following model [2]:

$$\mathbf{Y}(k,l)=\sum_{n=0}^{N-1}\sum_{m=0}^{M-1} \mathbf{f}(m,n)e^{-i2\pi(\frac{mk}{M}+\frac{(n-Δn_0)l}{N})}e^{-i2\pi(\triangle \mathbf{B_0} (m,n)t(k,l))} $$

$$\mathbf{Y}=\mathbf{Ef}\space\space\space Eq(1)$$

where E is the encoding operator that contains translational offset and off-resonance related effects, m, n are image indices and k, l are k-space indices, t(k,l) is the sample time for location (k,l) in k-space, ΔB0 is the B0 map, Δn0 is the translational offset whose direction varies with direction of phase encoding (+ve for blip up and –ve for blip down scans).

The proposed method has following two steps:

1) Estimating the translational offset: The translational offset Δn0 is estimated as the parameter that maximizes the mutual information (MI) similarity measure between preliminary model based reconstructions for blip-up and blip-down EPI data by solving the following unconstrained problem:

$$ \widehat{Δn_0}=arg_{Δn_0}max(MI(\mathbf{E_{1}}^{H}\mathbf{Y_{1}}, \mathbf{E_{2}}^{H}\mathbf{Y_{2}}))\space\space\space Eq(2)$$

where E1 and E2 are encoding operators for blip-up and blip-down data that differ in times of the k-space samples, Y1 and Y2 are the acquired blip-up and blip-down k-space data for b-value of 0 .

2) Model based Reconstruction: For each b-value, the EPI data acquired in blip-up and blip-down scans can be combined into a single formulation by setting Y=[Y1 Y2]T and E=[E1 E2]T in Eq (1). Model based reconstruction is done (Fig. 1) from combined k-space data Y in a manner analogous to iterative SENSE reconstruction [3] that solves Hermitian-symmetric form of Eq (1) based on conjugate gradient method [4] .

Experiments

Prostate scanning was performed on 3T scanner (Achieva, Philips Healthcare) equipped with 32-channel cardiac receive coil array. Multi-slice EPI data in both blip-up and blip-down phase encoding directions were acquired in two patients with parameters: resolution: 2 x 2 x 4 mm3, FOV: 180-200 x 180-200 x 55-80 mm3, partial Fourier acquisition, phase encoding direction AP/PA (for blip-up/down), b-values of 0, 500 and 1000 s/mm2. B0 map was calculated from a separate dual-echo gradient echo scan with resolution: 2x2x2 mm3, FOV: 200 x 200 x 90 mm3, echo time difference between two echoes=2.3msec, scan time: 1-1.2 min. The translational offset was estimated using Eq (2) followed by model based reconstruction. The proposed model based reconstruction from both blip-up and blip-down data was compared against uncorrected blip-up and blip-down data, correction using blip-up data only and blip-down data only. A modified T2W scan with fat saturation and resolution similar to the DWI was acquired for comparison.

Results and Discussion

For patient 1, the values of optimal translational offset (in pixels) as a function of slice number is shown in Fig.2a together with convergence plot as function of iteration number in Fig. 2b. The median optimal translational offset was in the range of 3-4 pixels across all the slices. In Fig. 3 and 4, the model based reconstructions from both blip-up and blip-down data are shown for b-values of 0 and 1000 s/mm2, respectively. Results for patient 2 are summarized in Fig. 5. In the region around the prostate rectal-air interface, the variation in B0 was around 120 Hz that corresponds to shift of 7-8 pixels at bandwidth/pixel of ~15.80 Hz in phase encoding direction. This resulted in pile-up in regions where B0 gradient magnitude was high and its direction was opposite to that of phase encoding gradient (see regions pointed out by red arrows in Fig. 4 and Fig.5). The proposed method corrected for translational offset as well as pile-up artefacts that cannot be corrected using only one-directional phase encoding. In future works, we will investigate the effects due to phase changes that may occur between diffusion weighted blip-up and blip-down EPI scans.

Acknowledgements

This work is supported by Cancer Research UK (grant A21099) and the NIHR Biomedical Research Unit at UCLH. We acknowledge research support provided by Philips Healthcare.

References

[1] Andersson et al, NeuroImag 2003

[2] Munger et al, TMI 2000

[3] Stoer et al, Num Math, 1990

[4] Pruessman et al, MRM 2001

Figures

Fig 1: Implementation of iterative model based reconstruction method using data from both blip up and blip down EPI scans. In each iteration of the conjugate gradient (CG )algorithm, we iterate back and forth between the raw k-space data (Y1,Y2) and image f via encoding operators E1 and E2 that contain the translational offset and B0 distortion information.

Fig.2: Invivo Results a) translational offset ( measured in number of pixels) estimated for a patient as a function of slice number, b) Convergence of minimization algorithm as a function of iteration number

Fig 3: Prostate patient 1 results for EPI data (b-value of 0): a) B0 map (in Hz) estimated from dual echo gradient echo scan, b and c) Uncorrected images reconstructed from blip-up and blip down EPI scans, red arrows indicate the regions of pile up, d and e) Corrected images using blip up only and blip down only, translational offsets are corrected but areas of pile up are not corrected. f) proposed method: both areas of pile up and stretching are corrected with close correspondence to modified T2W reference scan in (g)

Fig 4: Prostate patient 1 results for diffusion EPI data (b-value of 1000): a and b) Uncorrected images reconstructed from blip-up and blip down EPI scans, red arrows indicate the regions of pile up, c and d) corrected images using blip up only and blip down only, e) proposed method, f) modified T2W reference scan. Red arrows indicate the regions of pile up which were corrected with proposed method

Fig 5: Prostate patient 2 results for EPI data (b-value of 0 and 500): Uncorrected images reconstructed from blip-up and blip down EPI scans are shown in first and second columns, respectively. The proposed method reconstructions are shown in third column, Standard T2W scan as reference is shown in fourth column. Red arrows indicate the region of pile up that was corrected with the proposed technique.

Proc. Intl. Soc. Mag. Reson. Med. 26 (2018)
1639