2155

Optimisation of Poisson-disk Sampling Pattern for Highly Accelerated Femoral NCE-MRA
Hao Li1, Martin John Graves2, David John Lomas1, and Andrew Nicholas Priest2
1Department of Radiology, University of Cambridge, Cambridge, United Kingdom, 2Department of Radiology, Addenbrooke’s Hospital, Cambridge, United Kingdom

Synopsis

NCE-MRA techniques have a long acquisition time but can be accelerated by compressed sensing, parallel imaging and partial Fourier sampling using, for example, a Poisson-disk sampling pattern. In this study, we optimised the parameters in sampling pattern design with different acceleration factors for the 3D accelerated femoral fresh-blood-imaging sequence. The NCE-MRA data were reconstructed using weighted subtraction in k-space combined with phase correction. In a comparison using retrospective acceleration by subsampling of a full dataset, the optimised patterns outperformed the non-optimised patterns.

Introduction

NCE-MRA techniques such as fresh blood imaging (FBI)1 suffer from long acquisition times, which not only limits the clinical acceptance of these techniques but also renders them sensitive to artefacts from patient motion. Compressed sensing (CS), parallel imaging (PI) and partial Fourier (PF) sampling can be used to accelerate the acquisition by using, for example, a Poisson-disk sampling pattern (Fig. 1)2. However, the design of the optimal sampling pattern, and parameters still needs to be evaluated.
Conventionally the sampling pattern reflects the underlying signal density of k-space, which is most densely sampled in the centre of k-space3. Typically, a variable sampling density of the Poisson-disk distribution is preferred4. The density decay factor, determining the reduction of sampling density towards the periphery of k-space, needs to be optimised.
PF Sampling can be applied by covering only a fraction of k-space in both ky and kz dimensions (for 3D acquisitions) and using conjugate symmetry to fill the missing sections5. It remains unclear whether the introduction of PF can further improve acceleration performance or not in comparison with only using CS and PI, and what fraction of k-space in two dimensions should be filled in pattern design.
The fully sampled centre region of k-space is used for calibration in PI and phase estimation for PF homodyne detection. The shape and size of the calibration region can also influence the reconstruction quality and needs to be optimised.
This study investigates the optimal values for the above characteristics of sampling pattern design. Retrospective sub-sampling, i.e. acceleration, was performed on ten fully sampled datasets using different sampling patterns to determine the optimal ranges for parameters. The performance of the optimised patterns was then compared with typical non-optimised patterns.

Methods

Fully sampled femoral FBI datasets were acquired from ten healthy volunteers using a 1.5 T MR450 system (GE Healthcare, Waukesha, WI). Parameters included 224×224×80, ETL 80-90, FOV 40-44 cm, slice thickness 2 mm, TE 45ms, TR 2 or 3 heartbeats and total acquisition time 360 TRs. Retrospective simulated acceleration was performed using patterns with different parameters but the same total number of sampling points. Images were reconstructed using the k-space subtraction with phase and intensity correction method6.
Peak signal-to-noise ratio (PSNR) and structural similarity index measure (SSIM)7 were calculated to evaluate reconstruction accuracy, using fully sampled images as the reference. Contrast-to-noise ratio (CNR) of artery-to-background8 was used as a no-reference matrix to evaluate the arterial signal intensity level. Arterial and background signal was calculated by applying arterial and background tissue masks obtained from bright-blood images and subtracted images respectively. The standard deviation of the difference image (subtraction images of the reconstructed angiograms from the reference fully sampled angiograms) was used for noise estimation.
An automatic method was implemented to assess the sharpness of the vessel edges. Arterial regions were detected on the individual axial slices by the circular Hough transform9. The sharpness is then evaluated on the profiles of the arterial edges in two directions using a least-squares curve fit based on the sigmoid function10.

Results

Fig. 2 shows two examples of the parameter optimisation process. The density decay factor optimisation (A and B) shows that the four metrics all returned high values when the factor is between 0.08 and 0.14, indicating an optimal value range in this study. The PF sampling ratio optimisation (C and D) demonstrates that the range between 0.571 and 0.75 would be the optimal range for PF sampling ratio in the ky dimension. It can also be observed that a smaller PF ratio is more desirable for a larger AF (14x), so that a sufficient density of the sampling pattern can be maintained. Similarly, we determined the optimal value ranges of the other parameters under different AFs in this study (Table 1). Sampling patterns under different AFs were designed using these ranges of optimal values. As shown in Fig. 3, the performance of sampling patterns with AFs from 4 to 20 was evaluated on the ten fully sampled datasets. A typical series of non-optimised patterns with uniform density, square fully sampled centre regions and without PF were also evaluated for comparison. It can be observed that the optimised patterns improved the image quality in terms of signal loss, blurring and depiction of small branches. The corresponding quantitative evaluation results are summarised in Fig. 5. Optimised patterns outperformed non-optimised patterns in all four metrics.

Discussion and conclusion

In this study, we evaluated the influence of the parameters in Poisson-disk sampling pattern design and explored their optimal values. Optimised patterns outperformed non-optimised patterns in simulated acceleration. It was shown that compared to only using CS and PI, the implementation of PF in two directions (ky and kz) can further improve acceleration performance. Poisson-disk patterns using variable density with optimised k-space density decay factor performed better than uniformly distributed density patterns. It is possible that the optimal parameter values may vary for different vascular study regions, imaging protocols and reconstruction methods. Further work is planned to determine if these values are applicable generally or require specific technique- and region-based optimisation.

Acknowledgements

The authors acknowledge the support of the Addenbrooke’s Charitable Trust and the NIHR Cambridge Biomedical Research Centre. Hao Li acknowledges the China Scholarship Council and Cambridge Trust for fellowship support.

References

1. Miyazaki M, Sugiura S, Tateishi F, Wada H, Kassai Y, Abe H. Non-contrast-enhanced MR angiography using 3D ECG-synchronized half-Fourier fast spin echo. J. Magn. Reson. Imaging 2000;12:776–783.

2. Rapacchi S, Han F, Natsuaki Y, et al. High spatial and temporal resolution dynamic contrast-enhanced magnetic resonance angiography using compressed sensing with magnitude image subtraction. Magn. Reson. Med. 2014;71:1771–1783.

3. Hollingsworth KG. Reducing acquisition time in clinical MRI by data undersampling and compressed sensing reconstruction. Phys. Med. Biol. 2015;60:R297–R322.

4. Margosian PM, Seiberlich N, Takizawa M. Exploitation of Sparsity in MRI. 2013;2:289–306.

5. Xu Y, Mark Haacke E. Partial fourier imaging in multi-dimensions: A means to save a full factor of two in time. J. Magn. Reson. Imaging 2001;14:628–635.

6. Li H, Priest AN, Graves MJ, Lomas DJ. Highly Accelerated NCE-MRA Using Complex Subtraction with Intensity Correction: Improved Reconstruction Accuracy and Background Tissue Suppression. In: Proceedings of the 27th Annual Meeting of ISMRM, Montreal, Canada. ; 2019. p. 834.

7. Horé A, Ziou D. Image quality metrics: PSNR vs. SSIM. In: Proceedings - International Conference on Pattern Recognition. ; 2010. pp. 2366–2369.

8. Akasaka T, Fujimoto K, Yamamoto T, et al. Optimization of regularization parameters in compressed sensing of magnetic resonance angiography: Can statistical image metrics mimic radiologists’ perception? PLoS One 2016;11:1–14.

9. Atherton TJ, Kerbyson DJ. Size invariant circle detection. Image Vis. Comput. 1999;17:795–803.

10. Ahmad R, Ding Y, Simonetti OP. Edge sharpness assessment by parametric modeling: Application to magnetic resonance imaging. Concepts Magn. Reson. Part A 2015;44:138–149.

Figures

Fig. 1. A typical Poisson-disk sampling pattern for 3D accelerated imaging combining PI, CS and PF.

Fig. 2. Two examples of the parameter optimisation process. (a) and (b): Sampling patterns with varying density decay factor and corresponding quantitative measurements. (c) and (d): Sampling patterns with varying PF ratio in the ky dimension and corresponding quantitative measurements.

Table 1. The optimal value ranges and suggested values for the sampling pattern design (matrix size: 224×224×80)

Figure 4. Optimised (a) and non-optimised (b) sampling patterns and the corresponding example MIPs (right half) of reconstructed images (c and d) with different AFs. Blue arrowheads on (c) and (d) denote signal loss and blurring on images using non-optimised patterns. Yellow arrowheads denote improved delineation of small branches on images using optimised patterns.

Figure 5. Quantitative evaluation of acceleration using optimised and non-optimised sampling patterns with different AFs. The values are the mean and standard deviation calculated over the 10 volunteer datasets.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
2155