## Meta-QSM: An Image-Resolution-Arbitrary Network for QSM Reconstruction

Juan Liu ab and Kevin M. Koch abc a Center for Imaging Research, Medical College of Wisconsin, Milwaukee, WI, USA b Joint Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, WI, USA

c Radiology, Medical College of Wisconsin, Milwaukee, WI, USA

## ABSTRACT

Quantitative Susceptibility Mapping (QSM) can estimate the underlying tissue magnetic susceptibility and reveal pathology. Current deep-learning-based approaches to solve the QSM inverse problem are restricted on fixed image resolution. They trained a specific model for each image resolution which is inefficient in computing. In this work, we proposed a novel method called Meta-QSM to firstly solve QSM reconstruction of arbitrary image resolution with a single model. In Meta-QSM, weight prediction was used to predict the weights of kernels by taking the image resolution as input. The proposed method was evaluated on synthetic data and clinical data with comparison to existing QSM reconstruction methods. The experimental results showed the Meta-QSM can effectively reconstruct susceptibility maps with different image resolution using one neural network training.

Keywords: Quantitative Susceptibility Mapping, Deep Learning, Meta Learning

## 1. INTRODUCTION

Quantitative Susceptibility Mapping (QSM) can estimate the underlying tissue magnetic susceptibility to provide a novel image contrast and reveal pathology. QSM requires addressing a challenging post-processing problem: filtering of image phase estimates and inversion of the phase to susceptibility relationship. Both steps require solving ill-posed inverse computational problems, causing a wide variety of quantification errors, robustness limitations, and artifacts.

To date, all QSM methods rely on a dipolar convolution that relates source susceptibility to induced Larmor frequency offsets . 1, 2 While the forward relationship of this model (source to field) can be efficiently computed using Fast-Fourier-Transforms (FFT), a k-space singularity in the applied convolution filter results in an illconditioned relationship in the inverse model (field to source). Acquiring multiple orientations data to the magnetic field can sufficiently improve the conditioning of the inversion algorithm, serving as the empirical goldstandard for in vivo QSM assessment. 3 Single-orientation susceptibility maps computation is more challenging. Thresholding of the convolution operator 4-6 or use of more sophisticated regularization methods 7-10 are applied. However, these methods suffer from long computation time, regularization parameter tuning, conspicuity loss of fine details, and streaking artifacts, causing technical implementation challenges of QSM clinical translation and practice.

Current deep-learning-based methods such as QSMnet, 11 DeepQSM 12 have shown improved QSM reconstruction results compared with conventional methods. However, they trained a specific model for each image resolution, which requires train multiple models for different image resolution. If we train a specific model for each image resolution, it is impossible to store all these models and it is inefficient in computing. Thus, a method to solve QSM reconstruction of arbitrary image resolution is important for putting the deep-learning-based QSM approaches into practical use.

To solve these drawbacks, we proposed a method called Meta-QSM. It adopts the weight prediction method used in the meta-learning, to dynamically predict weights of the filters for each image resolution. Quantitative evaluation was conducted using synthetic test data with comparison to conventional QSM reconstruction methods. Besides, qualitative assessment was performed using two QSM datasets acquired with different image resolution.

## 2. METHODS

## 2.1 Training Data

We used one QSM dataset and data augmentation techniques to generate whole synthetic training data. The COSMOS dataset from 2016 ISMRM QSM challenge 13 was used, which was acquired with a fast 3-dimensional gradient-echo scans, 12 different head orientations, 1.06mm isotropic voxels on a 3T scanner. The COSMOS data was resampled to image resolution 0.75x0.75x3.0 mm . Elastic distortions are applied to geometrically transform 3 the susceptibility map. ? In addition, randomly sized geometric shapes, such as ellipsoids, spheres, cuboids, and cylinders with random susceptibility values are randomly placed on the augmented susceptibility maps to mimic the bleeding and calcifications, in order to increase the training data variability. Image contrast changes are applied to susceptibility maps as well to increases the training data variability in different image contrast. The image resolution of generated QSM data was randomly set to (0.5-1.0, 0.5-1.0, 2.0-4.0)mm 3 . The local fields are calculated using the well-defined forward dipole convolution relationship. The local fields and QSM images are cropped to image size 128x128x64 for training.

## (iii) Sagittal

Figure 1. One example of training data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-1acfcc97-661e-4fbe-ae92-2bcce8a3a988.jpg)

## 2.2 Neural Network Architecture

A 3D encoder-decoder convolutional neural network, based on a modified version of an established U-Net architecture, 14 was trained to perform pixel-wise susceptibility estimation. All vanilla convolution layers were replaced by weight-predict convolutional layers (WP-Conv), as shown in Fig.2. For WP-Conv, it consists of several fully connected layers and several activation layers. WP-Conv takes the image resolution as input, and the learn resolution specific weights and bias with the shape (inC, k, k, k, outC) and (outC) respectively. Here the inC is the number of channels of the input feature map. k is the kernel size, which is 3 in the paper. The outC is the number of channels of the output feature map.

(a) Neural Network Architecture

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-d9ec9ccc-e269-46fc-b581-c69d96a782fa.jpg)

(b) Weight-Prediction Convolutional Layer

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-7bf28644-7487-4bf1-85bd-81ef41d6e860.jpg)

Figure 2. Network structure of Meta-QSM. A 3D convolutional neural network was designed with 6 WP-Conv (kernel size 3x3x3, dilated rate 1x1x1), 3 WP-Conv (kernel size 3x3x3, dilated rate 2x2x2), 4 max pooling layer (pool size 2x2x2, stride size 2x2x2), 4 deconvolutional layers (stride size 2x2x2), 9 batch normalization layers, 4 feature concatenations, and 1 WP-Conv (kernel size 3x3x3, linear activation).

## 2.3 Training Details

L1 loss between the susceptibility map and label was utilized as a loss function for training. L2 regularization was added in all WP-Conv layers with regularization strength 0.01. 30,000 synthetic data were spilt into training dataset and validation dataset with ratio 9:1. The RMSprop optimizer was used in the deep learning training. The initial learning rate was set as 0.0001, with exponential decay at every 200 steps. Two NVIDIA tesla k40

graphics processing units (GPUs) were used for training with batch size 4. The neural network was trained and evaluated using Keras with Tensorflow as backend.

## 3. PERFORMANCE EVALUATION

## 3.1 Synthetic Data

100 synthetic data sets with matrix size 256x256x64 generated in similar fashion to the training data were used to evaluate the performance of Meta-QSM with comparison with Truncated K-Space Division (TKD) inversion, 4 iLSQR, 15 and Morphology Enabled Dipole Inversion (MEDI). 16 The results were evaluated against the groundtruth using root mean squared error (RMSE), high-frequency error norm (HFEN), and structural similarity (SSIM) index.

## 3.2 Comparison with Mesa-QSM

Wetrained two neural networks, denoted as Mesa-QSM, for image resolution 0.5x0.5x2.0mm 3 and 0.7639x0.7639x3.0mm . 3 Mesa-QSM replaces the WP-Conv with vanilla convolution layers in Meta-QSM. The neural networks was trained using the same hyperparameters.

After training, 100 synthetic data generated with image resolution 0.5x0.5x2.0mm 3 and 0.7639x0.7639x3.0mm 3 in in similar fashion to the training data were used to evaluate the performance of Meta-QSM and Mesa-QSM. The results were evaluated against the ground-truth using RMSE, HFEN, and SSIM index.

## 3.3 3T Research Data

MR imaging for a large local sports concussion study approved by the local institutional human research review board was performed on a 3T MRI scanner (GE Healthcare MR750) using a 32ch MRI head receive array using a commercially available susceptibility-weighted software application (SWAN, GE Healthcare). The data acquisition parameters were as follows: in-plane data matrix 320x256, field of view 24 cm, voxel size -0.5x0.5x2.0 mm , echo spacing - 7 ms, 4 echo times - [10.4, 3 17.4, 24.4, 31.4] ms, repetition time - 58.6 ms, autocalibrated parallel imaging factors - 3x1, acquisition time - 4 min.

Complex multi-echo images were reconstructed from raw k-space data using GE Orchestra. The brain masks were obtained using the SPM tool. 17 After background field removal using the Regularization-enabled Sophisticated Harmonic Artifact Reduction on Phase data (RESHARP) 18 method, susceptibility inversion was performed using the TKD, iLSQR, MEDI, Mesa-QSM, and Meta-QSM.

## 3.4 3T Clinical SWI Data

The clinical QSM data were acquired using susceptibility weighted angiography (SWAN, GE Healthcare), at a 3T MRI scanner (GE Healthcare MR750) with data acquisition parameters: in-plane data acquisition matrix 288x224, field of view 22 cm, slice thickness 3 mm, autocalibrated parallel imaging factors 1x2 or 1x3, number of slices 46-54, first echo time 12.6 ms, echo spacing 4.1 ms, number of echoes 7, flip angle 15 o , repetition time 39.7 ms, total scan time about 2 minutes.

The SWI images were processed by vendor reconstruction algorithms. The raw k-space data were saved for offline QSM processing. Multi-echo real and imaginary data were reconstructed from k-space data, with reconstruction matrix size 288x288, voxel size 0.76x0.76x3.0 mm 3 . The field map was obtained by the fitting of multi-echo phases. Brain masks were obtained using the SPM brain extraction tool. RESHARP were applied to remove the background field, with spherical kernel radius set as 6mm. Susceptibility inversion was performed using the TKD, iLSQR, MEDI, Mesa-QSM, and Meta-QSM.

In the above, for RESHARP, spherical kernel radius was set as 6mm to trade off the background removal performance and brain erosion. For TKD, the threshold was set to 0.20; for MEDI, the regularization factor was set to 1000. MEDI toolbox and STI toolbox publicly provided by their respective authors were used to calculate the QSM images use for comparison analyses in this report.

Table 1. Quantitative comparison Meta-QSM with conventional QSM reconstruction methods

|           | TKD               | iLSQR             | Meta-QSM          |
|-----------|-------------------|-------------------|-------------------|
| RMSE(%)   | 32 . 90 ± 1 . 55  | 52 . 94 ± 3 . 54  | 21 . 96 ± 3 . 84  |
| HFEN(%)   | 33 . 07 ± 1 . 65  | 51 . 11 ± 3 . 77  | 16 . 75 ± 3 . 84  |
| SSIM(0-1) | 0 . 960 ± 0 . 013 | 0 . 941 ± 0 . 023 | 0 . 977 ± 0 . 013 |

## 4. RESULTS

## 4.1 Synthetic Data

Table.1 illustrates error metrics using TKD, iLSQR, Meta-QSM from 100 synthetic data. The proposed method achieved the best score in RMSE, HFEN, and SSIM.

## 4.2 Comparison with Mesa-QSM

Table 2. Quantitative comparison Meta-QSM with Mesa-QSM on image resolution 0.5x0.5x2.0mm . 3

|           | Mesa-QSM          | Meta-QSM          |
|-----------|-------------------|-------------------|
| RMSE(%)   | 25 . 13 ± 1 . 75  | 28 . 33 ± 1 . 70  |
| HFEN(%)   | 17 . 07 ± 2 . 27  | 22 . 10 ± 2 . 11  |
| SSIM(0-1) | 0 . 968 ± 0 . 012 | 0 . 962 ± 0 . 014 |

Table 3. Quantitative comparison Meta-QSM with Mesa-QSM on image resolution 0.76x0.76x3.0mm . 3

|           | Mesa-QSM          | Meta-QSM          |
|-----------|-------------------|-------------------|
| RMSE(%)   | 25 . 34 ± 1 . 74  | 18 . 87 ± 2 . 00  |
| HFEN(%)   | 17 . 43 ± 2 . 38  | 14 . 04 ± 2 . 05  |
| SSIM(0-1) | 0 . 968 ± 0 . 012 | 0 . 984 ± 0 . 007 |

Table.2 and Table.3 shows the quantitative metrics of comparison Meta-QSM with two Mesa-QSM. For image resolution 0.76x0.76x3.0mm , Meta-QSM achieved better metric scores than Mesa-QSM. For image resolution 3 0.5x0.5x2.0mm , Mesa-QSM achieved slightly better metric scores than Meta-QSM. 3

## 4.3 3T Research Data

Fig.3 shows QSM images reconstructed by TKD, iLSQR, MEDI, Mesa-QSM, and Meta-QSM. Image blurring showed in TKD and MEDI results. Mesa-QSM and Meta-QSM can produce images with high image sharpness. Streaking artifacts is clearly observable in TKD and MEDI results. Meta-QSM showed a slight underestimated susceptibility values.

## 4.4 3T Clinical SWI Data

Fig.4 shows QSM images reconstructed by five methods and SWI images on a 26-year-old patient with dizziness and cognitive impairment after stroke. Meta-QSM can produce high quality QSM images.

Fig.5 shows QSM images and SWI images from a 28-year-old patient with left mesial temporal lesion and neurofibromatosis type 1. TKD, iLSQR, and MEDI results suffer from image blurring and shading artifacts around hemorrhage. Meta-QSM results show the best image quality.

Fig.6 shows QSM images and SWI images from a 36-year-old patient with surgical planning. Meta-QSM results can preserve the microstructures and suppress streaking artifacts.

Figure 3. Tissue fields (a) and QSM images (b-f) in a 3T MRI data with image resolution 0.5x0.5x2.0mm . 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-76871816-79d7-40a4-8823-1cf3e8e435be.jpg)

Figure 4. Tissue fields and QSM images (b-f), SWI images from a 26-year-old patient with dizziness and cognitive impairment after stroke.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-8796eb74-c205-48ac-bea7-a966fcefcd36.jpg)

## 5. DISCUSSION

Quantitative and qualitative analysis within this study has demonstrated that Meta-QSM can outperform conventional QSM reconstruction methods. In qualitative comparisons with TKD, iLSQR, and MEDI, the proposed method can outperform based on quantitative metrics evaluation. Qualitatively, Meta-QSM can not only preserve the tissue microstructure but also greatly suppress the streaking artifacts. Meta-QSM and Mesa-QSM show comparable performance in quantitative evaluation and qualitative visual assessment. For image resolution

Figure 5. Tissue fields and QSM images (b-f), SWI images from a 28-year-old patient with left mesial temporal lesion and neurofibromatosis type 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-899423c7-5e07-49de-be46-08b643dc4c6f.jpg)

Figure 6. Tissue fields and QSM images (b-f), SWI images from a 36-year-old patient with surgical planning.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00206v1-figure-633eba93-602a-48e4-bf6f-1edefe9b2e51.jpg)

0.5x0.5x2.0mm , Meta-QSM achieved slightly worse metric scores than Mesa-QSM, which may be due to the 3 training data with image resolution (0.5-1.0, 0.5-1.0, 2.0-4.0)mm 3 .

Compared QSMnet and DeepQSM, Meta-QSM has several advantages. First, Meta-QSM applied weight prediction to solve the fixed image resolution input. Second, Meta-QSM ultilized one in-vivo QSM dataset, data augmentation techniques and well-established physical model to generate whole training data. Third, Meta-QSM training data includes simulated bleeding/calcifications, which can greatly improve performance in pathological clinical cases.

There are several limitations of Meta-QSM. First, it is trained on synthetic data, which is due to the physical impossibility of developing a true experimental gold standard. The data distribution difference between the synthetic training data and collected test data can introduce inference errors. Second, Meta-QSM takes the local fields after background field removal to predict the susceptibility distribution. The background field removal errors can sequentially introduce susceptibility quantification errors.

## 6. CONCLUSION

We propose a novel QSM methods to solve the QSM reconstruction of arbitrary image resolution with a single model. The proposed Meta-Prediction Convolution could dynamically predict the weights of the filters. For each image resolution, the proposed Meta-Prediction Convolution generates a group of weights for the convolution layers. Thanks to the weight prediction, we can train a single model for QSM reconstruction of arbitrary image resolution.

## REFERENCES

- [1] Salomir, R., de Senneville, B. D., and Moonen, C. T., 'A fast calculation method for magnetic field inhomogeneity due to an arbitrary distribution of bulk susceptibility,' Concepts in Magnetic Resonance Part B 19 (1), 26-34 (2003).
- [2] Marques, J. and Bowtell, R., 'Application of a fourier-based method for rapid calculation of field inhomogeneity due to spatial variation of magnetic susceptibility,' Concepts in Magnetic Resonance Part B 25 (1), 65-78 (2005).
- [3] Liu, T., Spincemaille, P., De Rochefort, L., Kressler, B., and Wang, Y., 'Calculation of susceptibility through multiple orientation sampling (cosmos): a method for conditioning the inverse problem from measured magnetic field map to susceptibility source image in mri,' Magnetic Resonance in Medicine 61 (1), 196-204 (2009).
- [4] Shmueli, K., de Zwart, J. A., van Gelderen, P., Li, T.-Q., Dodd, S. J., and Duyn, J. H., 'Magnetic susceptibility mapping of brain tissue in vivo using mri phase data,' Magnetic Resonance in Medicine 62 (6), 1510-1522 (2009).
- [5] Wharton, S., Sch¨fer, A., and Bowtell, R., 'Susceptibility mapping in the human brain using threshold-based a k-space division,' Magnetic Resonance in Medicine 63 (5), 1292-1304 (2010).
- [6] Haacke, E., Tang, J., Neelavalli, J., and Cheng, Y., 'Susceptibility mapping as a means to visualize veins and quantify oxygen saturation,' Journal of Magnetic Resonance Imaging 32 (3), 663-676 (2010).
- [7] De Rochefort, L., Brown, R., Prince, M. R., and Wang, Y., 'Quantitative mr susceptibility mapping using piece-wise constant regularized inversion of the magnetic field,' Magnetic Resonance in Medicine 60 (4), 1003-1009 (2008).
- [8] de Rochefort, L., Liu, T., Kressler, B., Liu, J., Spincemaille, P., Lebon, V., Wu, J., and Wang, Y., 'Quantitative susceptibility map reconstruction from mr phase data using bayesian regularization: validation and application to brain imaging,' Magnetic Resonance in Medicine 63 (1), 194-206 (2010).
- [9] Liu, T., Liu, J., De Rochefort, L., Spincemaille, P., Khalidov, I., Ledoux, J. R., and Wang, Y., 'Morphology enabled dipole inversion (medi) from a single-angle acquisition: comparison with cosmos in human brain imaging,' Magnetic Resonance in Medicine 66 (3), 777-783 (2011).
- [10] Bilgic, B., Chatnuntawech, I., Fan, A. P., Setsompop, K., Cauley, S. F., Wald, L. L., and Adalsteinsson, E., 'Fast image reconstruction with l2-regularization,' Fast image reconstruction with L2regularization 40 (1), 181-191 (2014).
- [11] Yoon, J., Gong, E., Chatnuntawech, I., Bilgic, B., Lee, J., Jung, W., Ko, J., Jung, H., Setsompop, K., Zaharchuk, G., et al., 'Quantitative susceptibility mapping using deep neural network: Qsmnet,' NeuroImage (2018).
- [12] Bollmann, S., Rasmussen, K. G. B., Kristensen, M., Blendal, R. G., Østergaard, L. R., Plocharski, M., O'Brien, K., Langkammer, C., Janke, A., and Barth, M., 'Deepqsm-using deep learning to solve the dipole inversion for quantitative susceptibility mapping,' NeuroImage 195 , 373-383 (2019).
- [13] Langkammer, C., Schweser, F., Shmueli, K., Kames, C., Li, X., Guo, L., Milovic, C., Kim, J., Wei, H., Bredies, K., et al., 'Quantitative susceptibility mapping: report from the 2016 reconstruction challenge,' Magnetic Resonance in Medicine 79 (3), 1661-1673 (2018).
- [14] Ronneberger, O., Fischer, P., and Brox, T., 'U-net: Convolutional networks for biomedical image segmentation,' in [ International Conference on Medical Image Computing and Computer-assisted Intervention ], 234-241, Springer (2015).
- [15] Li, W., Wang, N., Yu, F., Han, H., Cao, W., Romero, R., Tantiwongkosi, B., Duong, T. Q., and Liu, C., 'A method for estimating and removing streaking artifacts in quantitative susceptibility mapping,' NeuroImage 108 , 111-122 (2015).
- [16] Liu, J., Liu, T., de Rochefort, L., Ledoux, J., Khalidov, I., Chen, W., Tsiouris, A. J., Wisnieff, C., Spincemaille, P., Prince, M. R., et al., 'Morphology enabled dipole inversion for quantitative susceptibility mapping using structural consistency between the magnitude image and the susceptibility map,' NeuroImage 59 (3), 2560-2568 (2012).
- [17] Brett, M., Anton, J.-L., Valabregue, R., Poline, J.-B., et al., 'Region of interest analysis using an spm toolbox,' in [ 8th International Conference on Functional Mapping of the Human Brain ], 16 (2), 497, Sendai (2002).
- [18] Sun, H. and Wilman, A. H., 'Background field removal using spherical mean value filtering and tikhonov regularization,' Magnetic Resonance in Medicine 71 (3), 1151-1157 (2014).
Juan Liu Kevin M. Koch
Meta-QSM: An Image-Resolution-Arbitrary Network for QSM Reconstruction
Quantitative Susceptibility Mapping (QSM) can estimate the underlying tissue
magnetic susceptibility and reveal pathology. Current deep-learning-based
approaches to solve the QSM inverse problem are restricted on fixed image
resolution. They trained a specific model for each image resolution which is
inefficient in computing. In this work, we proposed a novel method called
Meta-QSM to firstly solve QSM reconstruction of arbitrary image resolution with
a single model. In Meta-QSM, weight prediction was used to predict the weights
of kernels by taking the image resolution as input. The proposed method was
evaluated on synthetic data and clinical data with comparison to existing QSM
reconstruction methods. The experimental results showed the Meta-QSM can
effectively reconstruct susceptibility maps with different image resolution
using one neural network training.