Dongwook Lee, Jaejun Yoo, Sungho Tak, Jong Chul Ye
8 min
Magnetic Resonance Imaging (MRI) is essential for diagnostics but constrained by long scan times. Accelerated acquisition techniques like parallel imaging and compressed sensing are used to shorten scans, but they are computationally expensive and sometimes struggle to fully remove coherent aliasing artifacts. The researchers investigate whether deep residual learning can provide faster, higher-quality image reconstruction by treating artifact removal as an image-restoration task.
The authors utilize a dual frame U-Net architecture, which uses a multi-scale approach to capture global artifact patterns. A novel aspect of this method is the separate training of magnitude and phase networks. This decomposition simplifies training and allows the magnitude network to act as a post-processing tool for images where only magnitude data is available. If phase information is present, the framework can perform an iterative k-space interpolation, using the learned network to replace handcrafted regularizers used in traditional compressed sensing algorithms. The networks were trained using a combination of public datasets (HCP) and in vivo clinical MR data.
The proposed model successfully removed strong, coherent aliasing artifacts from both single-coil and multi-channel MRI data, consistently outperforming traditional compressed sensing (ALOHA) and parallel imaging (GRAPPA) methods. A key advantage is the computational efficiency: once trained, the network reconstructs images in milliseconds, orders of magnitude faster than conventional iterative methods that take seconds or minutes to complete. The authors further show that pre-training the magnitude network on large datasets and using recursive, multi-scale network structures significantly enhances the final image quality.
This research demonstrates that deep learning can overcome the computational bottleneck inherent in current accelerated MRI techniques. By generating accurate results in near real-time, this framework could enable faster clinical workflows, reducing patient discomfort while maintaining the diagnostic quality required for medical imaging. The modular design, allowing for the separate handling of magnitude and phase, offers flexibility for different clinical data types.
OBJECTIVE: Accelerated magnetic resonance (MR) image acquisition with compressed sensing (CS) and parallel imaging is a powerful method to reduce MR imaging scan time. However, many reconstruction algorithms have high computational costs. To address this, we investigate deep residual learning networks to remove aliasing artifacts from artifact corrupted images. METHODS: The deep residual learning networks are composed of magnitude and phase networks that are separately trained. If both phase and magnitude information are available, the proposed algorithm can work as an iterative k-space interpolation algorithm using framelet representation. When only magnitude data are available, the proposed approach works as an image domain postprocessing algorithm. RESULTS: Even with strong coherent aliasing artifacts, the proposed network successfully learned and removed the aliasing artifacts, whereas current parallel and CS reconstruction methods were unable to remove these artifacts. CONCLUSION: Comparisons using single and multiple coil acquisition show that the proposed residual network provides good reconstruction results with orders of magnitude faster computational time than existing CS methods. SIGNIFICANCE: The proposed deep learning framework may have a great potential for accelerated MR reconstruction by generating accurate results immediately.
Alex: Splitting magnitude from phase sounds like a convenience. Why does it buy anything?
Sam: Two reasons. First, a lot of clinical work only uses magnitude images, often on a workstation with no access to raw data. A magnitude network can run there as pure post-processing. It also means public magnitude datasets become training data. They pre-trained on about twelve thousand slices from the Human Connectome Project.
Alex: And the second reason?
Sam: The magnitude result helps the phase network. Inside the brain, phase varies smoothly. Outside the head, the signal is basically zero, but the phase there is random noise, jumping all over its range. That noise made the phase network hard to train. So they run the magnitude network first, threshold its output to get a brain mask, and zero the phase outside it. A threshold of ten percent of the maximum brightness worked well.
Alex: A hand-picked threshold makes me nervous. How sensitive is the result to that choice?
Sam: They checked. Below about five percent, the phase network trained poorly. Between five and thirty percent, the phase error barely changed. So within a wide band, the choice doesn't matter much.
Alex: You said earlier the method can also work with raw data. How does that fit?
Sam: If you have both magnitude and phase, you have the full complex image. Then you can go back to k-space, keep the measurements you actually took, and replace only the missing ones with the network's estimate. Repeat that, with a damping step called Krasnoselskii-Mann iteration, which blends each new estimate with the last one to help it converge. They also retrained the networks on these nearly clean intermediate images, so the networks learn to handle both heavily aliased and almost artifact-free inputs.
Alex: So one pair of networks does two jobs: an image filter when you only have pictures, and the inner step of a k-space solver when you have raw data.
Sam: That's the design, yes.
Alex: Now the evidence. What was it tested on, and how big?
Sam: The real scans are brain images from a 3 Tesla scanner with four receiver coils. There were 81 slices from nine people in total. Seven people were used for training, one for validation, and one for testing. So the test set is twelve images from a single subject.
Alex: That's thin. One test subject makes me wonder how general any of this is.
Sam: I'd share that worry. The training data did mix spin echo and gradient echo scans at various timings, which the authors say made it more robust. But the testing is one subject, one acceleration factor, and undersampling simulated after the fact from fully sampled data.
Alex: Did splitting magnitude and phase cost anything, compared with just feeding the network real and imaginary parts?
Sam: They built exactly that baseline, same architecture, same training. The two came out about the same. The normalized error on magnitude dropped from about 0.24 for the aliased input to about 0.155 for both networks. The split version had a slightly better structural similarity score. So, no cost, and the practical advantages come for free.
Alex: And against GRAPPA and ALOHA, which is the comparison people will care about?
Sam: On single-coil data, ALOHA improved on the input but left visible aliasing, while the network removed it with the lowest error. On four-coil data, GRAPPA showed many errors, ALOHA was better but imperfect, and the network did best. In the parts of the paper I have, that head-to-head is shown as images and error maps, not a numeric table. Pre-training on the Connectome data helped modestly, and adding the k-space iteration gave the lowest error of all.
Alex: Put the speed in concrete terms.
Sam: The network needed under 41 milliseconds for a multi-coil magnitude image. With phase, since magnitude has to run first, about 61 milliseconds. The existing methods took from about thirty seconds to several minutes. Even fifteen rounds of the k-space iteration took around nine seconds. Training is the slow part: hours per dataset, done once.
Alex: What's the limitation that matters most?
Sam: Generalization. Because it works on images, the artifact pattern it learned depends on scan settings like field of view and matrix size. The authors say training on varied acquisitions is the remedy, but that wasn't tested broadly here. They also note some remaining blur. They attribute it to the squared-error loss, which favours averaged outputs, and suggest adversarial or perceptual losses as future work.
Alex: So who should sit down with the full paper, and where should they start?
Sam: People building MRI reconstruction, especially anyone stuck with magnitude-only data. Read the reconstruction flow section, where the masking and the iterative k-space scheme are laid out. Then look at the threshold experiment and the comparison against the real-and-imaginary baseline. Treat the quality comparisons as promising, given the single test subject.
Alex: And for everyone else, the one line to keep?
Sam: Split the image into brightness and phase, and one network becomes both a quick filter and the core of a k-space solver.
Alex: Fast scans, fast cleanup, and a test set that still needs to grow.