1,721,158 research outputs found
Autocalibrating and calibrationless parallel magnetic resonance imaging as a bilinear inverse problem
No full textModern reconstruction methods for magnetic resonance imaging (MRI) exploit the spatially varying sensitivity profiles of receive-coil arrays as additional source of information. This allows to reduce the number of time-consuming Fourier-encoding steps by undersampling. The receive sensitivities are a priori unknown and influenced by geometry and electric properties of the (moving) subject. For optimal results, they need to be estimated jointly with the image from the same undersampled measurement data. Formulated as an inverse problem, this leads to a bilinear reconstruction problem related to multi-channel blind deconvolution. In this work, we will discuss some recently developed approaches for the solution of this problem
Bayesian MRI reconstruction with joint uncertainty estimation using diffusion models
No full textPurpose
We introduce a framework that enables efficient sampling from learned probability distributions for MRI reconstruction.
Method
Samples are drawn from the posterior distribution given the measured k-space using the Markov chain Monte Carlo (MCMC) method, different from conventional deep learning-based MRI reconstruction techniques. In addition to the maximum a posteriori estimate for the image, which can be obtained by maximizing the log-likelihood indirectly or directly, the minimum mean square error estimate and uncertainty maps can also be computed from those drawn samples. The data-driven Markov chains are constructed with the score-based generative model learned from a given image database and are independent of the forward operator that is used to model the k-space measurement.
Results
We numerically investigate the framework from these perspectives: (1) the interpretation of the uncertainty of the image reconstructed from undersampled k-space; (2) the effect of the number of noise scales used to train the generative models; (3) using a burn-in phase in MCMC sampling to reduce computation; (4) the comparison to conventional
-wavelet regularized reconstruction; (5) the transferability of learned information; and (6) the comparison to fastMRI challenge.
Conclusion
A framework is described that connects the diffusion process and advanced generative models with Markov chains. We demonstrate its flexibility in terms of contrasts and sampling patterns using advanced generative priors and the benefits of also quantifying the uncertainty for every pixel
Method and device for reconstructing a sequence of magnetic resonance images
No full textA method for reconstructing a sequence of magnetic resonance (MR) images of an object under investigation, includes the steps of (a) providing a series of sets of image raw data including an image content of the MR images to be reconstructed, the image raw data being collected with the use of at least one radiofrequency receiver coil of a magnetic resonance imaging (MRI) device, wherein each set of image raw data includes a plurality of data samples being generated with a gradient-echo sequence, in particular a FLASH sequence, that spatially encodes an MRI signal received with the at least one radiofrequency receiver coil using a non-Cartesian k-space trajectory, each set of image raw data includes a set of homogeneously distributed lines in k-space with equivalent spatial frequency content, the lines of each set of image raw data cross the center of k-space and cover a continuous range of spatial frequencies, and the positions of the lines of each set of image raw data differ in successive sets of image raw data, and (b) subjecting the sets of image raw data to a regularized nonlinear inverse reconstruction process to provide the sequence of MR images, wherein each of the MR images is created by a simultaneous estimation of a sensitivity of the at least one receiver coil and the image content and in dependency on a difference between a current estimation of the sensitivity of the at least one receiver coil and the image content and a preceding estimation of the sensitivity of the at least one receiver coil and the image content
Accelerating Non-Cartesian MRI Reconstruction Convergence Using k-Space Preconditioning
No full textWe propose a k-space preconditioning formulation for accelerating the convergence of iterative Magnetic Resonance Imaging (MRI) reconstructions from non-uniformly sampled k-space data. Existing methods either use sampling density compensations which sacrifice reconstruction accuracy, or circulant preconditioners which increase per-iteration computation. Our approach overcomes both shortcomings. Concretely, we show that viewing the reconstruction problem in the dual formulation allows us to precondition in k-space using density-compensation-like operations. Using the primal-dual hybrid gradient method, the proposed preconditioning method does not have inner loops and are competitive in accelerating convergence compared to existing algorithms. We derive l2 -optimized preconditioners, and demonstrate through experiments that the proposed method converges in about ten iterations in practice
Method and device for reconstructing a sequence of magnetic resonance images
No full textA method for reconstructing a sequence of magnetic resonance (MR) images of an object under investigation, includes the steps of (a) providing a series of sets of image raw data including an image content of the MR images to be reconstructed, the image raw data being collected with the use of at least one radiofrequency receiver coil of a magnetic resonance imaging (MRI) device, wherein each set of image raw data includes a plurality of data samples being generated with a gradient-echo sequence, in particular a FLASH sequence, that spatially encodes an MRI signal received with the at least one radiofrequency receiver coil using a non-Cartesian k-space trajectory, each set of image raw data includes a set of homogeneously distributed lines in k-space with equivalent spatial frequency content, the lines of each set of image raw data cross the center of k-space and cover a continuous range of spatial frequencies, and the positions of the lines of each set of image raw data differ in successive sets of image raw data, and (b) subjecting the sets of image raw data to a regularized nonlinear inverse reconstruction process to provide the sequence of MR images, wherein each of the MR images is created by a simultaneous estimation of a sensitivity of the at least one receiver coil and the image content and in dependency on a difference between a current estimation of the sensitivity of the at least one receiver coil and the image content and a preceding estimation of the sensitivity of the at least one receiver coil and the image content
A Multi-GPU Programming Library for Real-Time Applications
No full textWe present MGPU, a C++ programming library targeted at single-node multi-GPU systems. Such systems combine disproportionate floating point performance with high data locality and are thus well suited to implement real-time algorithms. We describe the library design, programming interface and implementation details in light of this specific problem domain. The core concepts of this work are a novel kind of container abstraction and MPI-like communication methods for intra-system communication. We further demonstrate how MGPU is used as a framework for porting existing GPU libraries to multi-device architectures. Putting our library to the test, we accelerate an iterative non-linear image reconstruction algorithm for real-time magnetic resonance imaging using multiple GPUs. We achieve a speed-up of about 1.7 using 2 GPUs and reach a final speed-up of 2.1 with 4 GPUs. These promising results lead us to conclude that multi-GPU systems are a viable solution for real-time MRI reconstruction as well as signal-processing applications in general
Physics-based Reconstruction Methods for Magnetic Resonance Imaging
No full textConventional Magnetic Resonance Imaging (MRI) is hampered by long scan times and only qualitative image contrasts that prohibit a direct comparison between different systems. To address these limitations, model-based reconstructions explicitly model the physical laws that govern the MRI signal generation. By formulating image reconstruction as an inverse problem, quantitative maps of the underlying physical parameters can then be extracted directly from efficiently acquired k-space signals without intermediate image reconstruction -- addressing both shortcomings of conventional MRI at the same time. This review will discuss basic concepts of model-based reconstructions and report about our experience in developing several model-based methods over the last decade using selected examples that are provided complete with data and code
- …
