Evaluation of different high-pass filter on the susceptibility in patients Parkinson’s disease and controls
Gerd Melkus1,2, Santanu Chakraborty1,2, and Fahad A Essbaiheen 1,2,3

1Medical Imaging, The Ottawa Hospital, Ottawa, ON, Canada, 2Radiology, University of Ottawa, Ottawa, ON, Canada, 3King Saud University, Riyadh, Saudi Arabia

Synopsis

Quantitative susceptibility mapping (QSM) was found as a useful method to evaluate neurodegenerative diseases such as Parkinson’s disease. For QSM reconstruction background field removed phase data is needed, but for retrospective studies only high-pass filtered data might be available. In this study we analyzed the influence of different high-pass filtered phase images on the susceptibility assessment for volunteers and Parkinson’s disease patients and compared the results to QSM estimation using background field removed phase data. With increasing high-pass filter strength consistently lower susceptibility results, but up to a certain filter strength differences in susceptibility can still be distinguished.

Purpose

Quantitative susceptibility mapping (QSM) is an emerging post-processing technique that uses gradient echo phase images to calculate the tissue magnetic susceptibility and it was found as a valuable method to characterize diseases such as Parkinson’s disease (PD) 1,2. It was also shown that Susceptibility Weighted Imaging (SWI) is useful to evaluate neurodegenerative diseases, but susceptibility quantification is hindered by the non-locality and orient dependency of high-pass filtered phase images 3,4. For retrospective studies, QSM might be obtained from available magnitude and vendor-provided high-pass filtered phase images of SWI. The purpose of this study is to estimate to what extent current QSM values may be affected when using vendor-provided high-pass filtered phase data. We analyzed the influence of different high-pass filtered phase images on the susceptibility assessment and compared the results to quantitative susceptibility estimation using a background field removed phase data.

Methods

All MRI experiments were performed at 3 Tesla (Siemens Trio, Siemens Medical, Erlangen, Germany) using a 32-channel head-coil. Flow compensated 3D multi gradient echo images were acquired on five volunteers and seven PD patients with the following parameters: TR = 44ms, TEs = 4.6, 11.2, 18.0, 24.8, 31.6, 38.4ms, α = 20°, FOV = 200 x 200mm2, matrix = 256 x 256, resolution = 0.75 x 0.75mm2, slice thickness = 1.5mm, number of slices = 80. Magnitude and phase data were stored on the scanner. Post processing of the data at TE = 18ms was performed in Matlab (Matlab 2014a, The MathWorks) in the following steps: (1) phase unwrapping, (2) high-pass filtering using different filter strengths (ranging from 0.01 to 0.3 of 3dB cutoff point, (Fig. 1)) and (3) susceptibility estimation using the iLSQR 5 method (an algorithm for sparse linear equations and sparse least squares) available from the STISuite toolbox 6. As reference, QSM post-processing was performed using a background field removal method (the V-SHARP algorithm 7) and susceptibility quantification applying iLSQR. Region of interest (ROIs) were draw in the red nucleus (RN) and the substantia nigra (SN) and the results for the two subject groups were compared for the different high-pass filtered images.

Results

Fig. 1 shows the eight different high-pass filter used in this analysis (from h=0.01 (the weakest) to h=0.3 (the strongest)). Fig. 2a shows the original phase, Fig. 2b the phase after phase unwrapping. The purple rectangle in Fig. 2c displays the SN/RN region of the brain for which the results of the remaining sub-figures of Fig. 2 are shown. Fig. 2d-g shows the phase images after different high-pass filter are applied. With increasing filter strength from h=0.025 to h=0.3 the phase is reduced for larger tissues (like SN and RN), the contours of the object remain. A similar trend can be seen for the calculated susceptibility using these high-pass filtered images (Fig. 2h-k). The results of the ROI analysis in SN and RN (Fig. 3a) can be seen in Fig. 3b and c, respectively. The calculated susceptibility values are in both tissues decreasing with increasing high-pass filter strength. While for the RN no significant difference was found between the PD and the control group, the SN showed a significant difference between both groups for susceptibility calculation using the background field removed phase and a high-pass filter up to h=0.05. Fig. 4 shows a comparison between a high-pass filtered phase image provided by the vendor (Fig. 4a) and using a high pass filter h=0.05 (Fig. 4b) on two different subjects. These images appear visually similar, though actual filter parameters used by the vendor are not known.

Discussion

This study evaluated the effect of different high-pass filtered phase images on the outcome of the susceptibility estimation for two different brain tissues (SN and RN) in PD and controls. The analysis showed, that the quantification using high-pass filtered susceptibility estimation resulted in consistently lower susceptibility values, but significant differences in the susceptibility could be detected up to a certain filter strength (here h=0.05 for the SN). For analyzing susceptibility differences in a retrospective study, where only the magnitude and the high-pass filter phase data is available, the high-pass filter strength needs to be evaluated and the effect of the filter needs to be considered.

Conclusion

The increasing strength of the high-pass filter results in consistently lower susceptibility values in SN and RN. Up to a certain filter strength (here h=0.05) differences in susceptibility can still be distinguished. However, as the vendor provided images visually appear close to our h=0.05 images, retrospective QSM reconstruction using vendor-provided SWI data will likely yield fair estimate of susceptibility.

Acknowledgements

No acknowledgement found.

References

1. Barbosa JHO, Santos AC, Tumas V, et al. Quantifying brain iron deposition in patients with Parkinson's disease using quantitative susceptibility mapping, R2 and R2. Magn Reson Imaging. 2015;33(5):559-65.

2. Lotfipour AK, Wharton S, Schwarz ST, et al. High resolution magnetic susceptibility mapping of the substantia nigra in Parkinson's disease. J Magn Reson Imaging. 2012;35(1):48-55.

3. Shmueli K, de Zwart JA, van Gelderen P, et al. Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data. Magnetic susceptibility mapping of brain tissue in vivo using MRI phase data. Magn Reson Med. 2009;62(6):1510-22.

4. Schäfer A, Wharton S, Gowland P, et al. Using magnetic field simulation to study susceptibility-related phase contrast in gradient echo MRI. Neuroimage. 2009;48(1):126-37.

5. Li W, Avram AV, Wu B, et al. Integrated Laplacian-based phase unwrapping and background phase removal for quantitative susceptibility mapping. NMR Biomed. 2014;27(2):219-27.

6. http://people.duke.edu/~cl160/index.html

7. Wu B, Li W, Guidon A, et al. Whole brain susceptibility mapping using compressed sensing. Magn Reson Med. 2012;67(1):137-47.

Figures

Figure 1. Eight different high pass filter used for the phase filtering. Ranging from 0.01 to 0.3 (h = 0.01, 0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.3) of the 3dB cutoff point.

Figure 2. Post-processing and susceptibility results of different high-pass filtered phase images for the SN/RN-region of a control subject. (a) Original phase image, (b) phase after phase unwrapping. (c) The purple rectangle indicates the location of images shown in (d-k). (d-g) Phase images after high-pass filtering. (h-k) Corresponding susceptibility maps.

Figure 3. (a) ROIs for the SN and RN shown as a blue and green overlay, respectively. (b) Susceptibility for SN for different phase filter. BFR = background field removal using V-SHARP. (c) Susceptibility for RN after high-pass phase filtering. * = significant difference (p<0.05) between PD and control group.


Figure 4. Exemplary comparison of high-pass filtered phase images from two different subjects. (a) Vendor-provided high-pass filtered image. (b) High-pass filtered image using the filter strength h = 0.05.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
3384