MOCCA: A Fast Algorithm for Parallel MRI Reconstruction Using Model Based Coil Calibration

Gerlind Plonka,Yannick Riebe
2024-10-29
Abstract:We propose a new fast algorithm for simultaneous recovery of the coil sensitivities and of the magnetization image from incomplete Fourier measurements in parallel MRI. Our approach is based on a parameter model for the coil sensitivities using bivariate trigonometric polynomials of small degree. The derived MOCCA algorithm has low computational complexity of $O(N_c N^2 \log N)$ for $N \times N$ images and $N_c$ coils and achieves very good performance for incomplete MRI data. We present a complete mathematical analysis of the proposed reconstruction method. Further, we show that MOCCA achieves similarly good reconstruction results as ESPIRiT with a considerably smaller numerical effort which is due to the employed parameter model. Our numerical examples show that MOCCA can outperform several other reconstruction methods.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in parallel magnetic resonance imaging (parallel MRI), how to simultaneously recover coil sensitivities and magnetization images from incomplete Fourier measurement data. Specifically, the author proposes a new fast algorithm - MOCCA (Model - based Coil Calibration Algorithm) to solve the following problems: 1. **Simultaneous recovery of coil sensitivities and magnetization images**: Traditional parallel MRI reconstruction methods usually require pre - estimating coil sensitivities, but these sensitivities are often unknown in practical applications and have a complex relationship with magnetization images. Therefore, how to simultaneously recover coil sensitivities and magnetization images from incomplete k - space data is a challenge. 2. **Reduction of computational complexity**: Existing methods such as ESPIRiT can achieve good reconstruction results, but have high computational complexity, especially when dealing with high - resolution images. Therefore, developing an algorithm with lower computational complexity is crucial for improving reconstruction speed and efficiency. 3. **Multi - channel blind deconvolution problem**: Parallel MRI reconstruction can be regarded as a multi - channel blind deconvolution problem because the relationship between coil sensitivities and magnetization images is unknown, and there is also a correlation between the sensitivities of different coils. How to effectively solve this problem is a key challenge. ### Main contributions of the MOCCA algorithm - **Model assumption based on bivariate trigonometric polynomials**: The MOCCA algorithm assumes that coil sensitivities can be represented by low - order bivariate trigonometric polynomials, so that the sensitivities can be determined by a small number of parameters. This not only reduces the computational complexity but also improves the reconstruction accuracy. - **Low computational complexity**: The computational complexity of the MOCCA algorithm is \(O(N_c N^2 \log N)\), where \(N_c\) is the number of coils and \(N\) is the side length of the image. Compared with other methods, the MOCCA algorithm significantly reduces the amount of computation while maintaining good reconstruction results. - **No need to solve eigenvalue problems pixel - by - pixel**: Unlike methods such as ESPIRiT, the MOCCA algorithm does not need to solve eigenvalue problems for each pixel. Instead, it recovers all sensitivity parameters by calculating an eigenvector of a medium - sized MOCCA matrix. This avoids problems such as phase correction and further simplifies the calculation process. ### Summary The MOCCA algorithm successfully solves multiple challenges in parallel MRI reconstruction, including the simultaneous recovery of coil sensitivities and magnetization images, the reduction of computational complexity, and the multi - channel blind deconvolution problem, by introducing new mathematical models and optimizing calculation methods. Experimental results show that the MOCCA algorithm can be comparable to or even better than existing advanced methods in performance, while having lower computational complexity.