2621

LoopNet: A New Baseline Network for QSM Dipole Inversion
Chen Chen1, Yang Gao1, Min Li1, Zhuang Xiong2, Feng Liu2, and Hongfu Sun2
1School of Computer Science and Engineering, Central South University, Changsha, China, 2School of EECS, The University of Queensland, Brisbane, Australia

Synopsis

Keywords: Susceptibility/QSM, Quantitative Susceptibility mapping, LoopNet, Bidirectional U-net, QSM dipole inversion

Motivation: Most current deep learning (DL) QSM methods were developed based on U-net, whose performances might not be sufficiently good.

Goal(s): To proposed a new network baseline for deep learning QSM methods development.

Approach: We developed a LoopNet, by applying the proposed bidirectional loop and a self-tailored GHPA attention module into a Unet backbone, making better use of the latent information in deep networks.

Results: Simulated and in vivo experiments showed that the propoed LoopNet led to improved results than U-net.

Impact: This work introduces a novel deep neural network backbone, allowing researchers to develop innovative QSM methods easily by upgrading their original U-net to LoopNet, thanks to the plug-and-play design.

Introduction

QSM is a valuable technique that can calculate tissue susceptibility distribution from MRI phases. However, the QSM dipole inversion is a mathematically ill-posed inverse problem. Many deep learning (DL) algorithms [1-3] have been developed based on the conventional U-net [4] backbone. In this work, we proposed a novel baseline network, named LoopNet, for DLQSM reconstruction from the local field data. The LoopNet has more efficient information (hidden features) usage dataflow compared with traditional U-net. The major contribution is the introduction of a “reverse” concatenation connection from the U-net’s expanding path to the extracting path as well as a tailored Grouped multi-axis Hadamard Product Attention (GHPA) module, which is a memory-efficient to capture long-range information, which may be beneficial for QSM reconstruction from local field data. We conducted both simulated and in vivo experiments to validate the performances of the proposed LoopNet.

Method

LoopNet Architecture
As shown in Fig. 1(a), the proposed LoopNet was constructed based on a U-net backbone by adding a reverse concatenation connection from the expanding path of the U-net to the previous layers, thus forming a bidirectional loop (Fig. 1(b)). The GHPA is shown in Fig 1(c), which is designed to capture long-range information in the latent features. Different from conventional U-nets, LoopNet iterates the U-net body several times to make full use of the encoded latent maps, and the U-net extracting and expanding paths were shared for each iteration. The U-net backbone comprises 19 convolutional layers, 8GHPA modules, 18 batch normalization layers, 4 max-pooling, 4 unpooling layers, 4 forward concatenation layers, 4 reverse concatenation connections, and a final skip connection.
Training data preparation and network training
A total of 14400 cropped small patches (size: 483) were generated from 96 in vivo QSM subjects. The training inputs were the local field data simulated from the QSM label patches using the QSM forward model. All network parameters were initialized with random numbers drawn from a normal distribution with a mean of 0 and a standard deviation of 0.01. Subsequently, training for all networks was conducted for 100 epochs on two Nvidia Tesla A6000 GPUs using the Adam optimizer. The learning rate was set to 0.001, and mean-squared-error was employed as the loss function.
Validation on simulated and in vivo datasets
In this study, we compared our proposed LoopNet with several existing methods, including DL-based xQSM [1], QSMNet [2], LPCNN [3], and traditional iLSQR [5]. Simulated data generated based on a COSMOS label, and the in vivo data from two patients with intracranial hemorrhage acquired at 3T were adopted to investigate the performance of the proposed LoopNet.

Results

The effectiveness of the proposed bidirectional loop and the GHPA module was investigated using a simulated data, as shown in Fig. 2. As pointed out by the blue arrows, the proposed LoopNet successfully restored a small vein in the frontal white matter region, while the traditional U-net backbone failed to present it clearly. It is also found from the numerical metrics that the performances of the proposed LoopNet gradually increased with the iteration numbers.
Figure 3 compares the proposed LoopNet with multiple U-net based QSM methods on a simulated local field map. The proposed method resulted in the best numerical metrics with the highest PSNR and SSIM, and the lowest MSE. Besides, the LoopNet also achieved the minimum reconstruction error compared with previous DLQSM methods, as shown in the zoomed-in images.
Figure 4 compares different QSM reconstruction methods using two in vivo brain data from two patients with intracranial hemorrhage. All DLQSM methods led to similar QSM results, while according to the mIP images, xQSM and the proposed LoopNet seemed to locate the hemorrhage lesion more precisely, as pointed out by the blue arrows, maybe because they both introduced advanced network architectures (octave convolution and iterative architecture) to make full use of the latent information in the U-net.

Discussion and Conclusion

This paper introduces a novel deep neural network backbone, namely LoopNet, for solving QSM dipole inversion problem. Comparative studies were carried out on both simulated and in vivo data, and the results showed that the proposed LoopNet led to better QSM images with improved numerical metrics.

Acknowledgements

YG acknowledges support from the National Natural Science Foundation of China under Grant No. 62301616. HS acknowledges support from the Australian Research Council (DE210101297,DP230101628).

References

1. Gao, Y., Zhu, X., Moffat, B.A., Glarin, R., Wilman, A.H., Pike, G.B., Crozier, S., Liu, F., Sun, H., 2021. xQSM: quantitative susceptibility mapping with octave convolutional and noise-regularized neural networks. NMR Biomed. 34 .

2. J. Yoon, E. Gong, I. Chatnuntawech, B. Bilgic, J. Lee, W. Jung, J. Ko, H. Jung, K. Setsompop, G. Zaharchuk, E.Y. Kim, J. Pauly, J. Lee. Quantitative susceptibility mapping using deep neural network: QSMnet. Neuroimage. 2018.

3. Lai, K., M. Aggarwal, P.van, Xu.L,J,Sulam.(2020) Learned Proximal Networks for Quantitative Susceptibility Mapping. eess.IV.2008.05024.

4. Ronneberger O, Fischer P, Brox T. U-net: Convolutional networks for biomedical image segmentation. International Conference on Medical image computing and computer-assisted intervention. Springer, Cham, 2015: 234-241.

5. Li, W., Wang, N., Yu, F., Han, H., Cao, W., Romero, R., … Liu, C. (2015). A method for estimating and removing streaking artifacts in quantitative susceptibility mapping. NeuroImage.

Figures

Figure 1.The proposed LoopNet for QSM reconstruction, which is constructed upon (a) a conventional U-net backbone by incorporating (b) the proposed bidirectional Loop module and (c) the proposed Grouped multi-axis Hadamard Producet Attention (GHPA) module. Bidirectional skip connections are adopted to form the Loop architecture, while the GHPA module can enhance long-range information capture.

Figure 2. Ablation studies of the proposed bidirectional loop and GHAP modules. (a) The top two rows demonstrate the QSM images of different methods and the corresponding error maps, with MSE/PSNR/SSIM reported. (b) The bottom part compares the mIP maps of different QSM images. Blue arrows point out a vanishing vein in Unet reconstruction is successfully preserved by our methods, and the red arrows point to reconstruction error around Globus Pallidus.

Figure 3. Comparison of various QSM methods on a simulated human brain data at 3T. The reconstructions and corresponding error maps around the Golubus Pallidus region are demonstrated in the zoomed-in blue boxes. Numerical metrics including MSE, PSNR, and SSIM are also reported.

Figure 4. Comparison of the proposed QSM method with traditional iLSQR, and DL-based QSMnet, xQSM, and LPCNN on two in vivo brains scanned at 3T. The top rows (a) illustrate the corresponding reconstruction results of different methods, while the bottom part shows the corresponding mIP images. Blue arrows point to the hemorrhage lesion.

Proc. Intl. Soc. Mag. Reson. Med. 32 (2024)
2621
DOI: https://doi.org/10.58530/2024/2621