Multi-Device Parallel MRI Reconstruction: Efficient Partitioning for Undersampled 5D Cardiac CINE.

Cardiac CINE, a form of dynamic cardiac MRI, is indispensable in the diagnosis and treatment of heart conditions, offering detailed visualization essential for the early detection of cardiac diseases. As the demand for higher-resolution images increases, so does the volume of data requiring processing, presenting significant computational challenges that can impede the efficiency of diagnostic imaging. Our research presents an approach that takes advantage of the computational power of multiple Graphics Processing Units (GPUs) to address these challenges. GPUs are devices capable of performing large volumes of computations in a short period, and have significantly improved the cardiac MRI reconstruction process, allowing images to be produced faster. The innovation of our work resides in utilizing a multi-device system capable of processing the substantial data volumes demanded by high-resolution, five-dimensional cardiac MRI. This system surpasses the memory capacity limitations of single GPUs by partitioning large datasets into smaller, manageable segments for parallel processing, thereby preserving image integrity and accelerating reconstruction times. Utilizing OpenCL technology, our system offers adaptability and cross-platform functionality, ensuring wider applicability. The proposed multi-device approach offers an advancement in medical imaging, accelerating the reconstruction process and facilitating faster and more effective cardiac health assessment.

Elastic AlignedSENSE for Dynamic MR Reconstruction: A Proof of Concept in Cardiac Cine.

Numerous methods in the extensive literature on magnetic resonance imaging (MRI) reconstruction exploit temporal redundancy to accelerate cardiac cine. Some of them include motion compensation, which involves high computational costs and long runtimes. In this work, we proposed a method-elastic alignedSENSE (EAS)-for the direct reconstruction of a motion-free image plus a set of nonrigid deformations to reconstruct a 2D cardiac sequence. The feasibility of the proposed approach was tested in 2D Cartesian and golden radial multi-coil breath-hold cardiac cine acquisitions. The proposed approach was compared against parallel imaging compressed sense (sPICS) and group-wise motion corrected compressed sense (GWCS) reconstructions. EAS provides better results on objective measures with considerable less runtime when an acceleration factor is higher than 10×. Subjective assessment of an expert, however, invited proposing the combination of EAS and GWCS as a preferable alternative to GWCS or EAS in isolation.

Jacobian weighted temporal total variation for motion compensated compressed sensing reconstruction of dynamic MRI.

PURPOSE: To eliminate the need of spatial intraframe regularization in a recently reported dynamic MRI compressed-sensing-based reconstruction method with motion compensation and to increase its performance. THEORY AND METHODS: We propose a new regularization metric based on the introduction of a spatial weighting measure given by the Jacobian of the estimated deformations. It shows convenient discretization properties and, as a byproduct, it also provides a theoretical support to a result reported by others based on an intuitive design. The method has been applied to the reconstruction of both short and long axis views of the heart of four healthy volunteers. Quantitative image quality metrics as well as straightforward visual assessment are reported. RESULTS: Short and long axis reconstructions of cardiac cine MRI sequences have shown superior results than previously reported methods both in terms of quantitative metrics and of visual assessment. Fine details are better preserved due to the lack of additional intraframe regularization, with no significant image artifacts even for an acceleration factor of 12. CONCLUSIONS: The proposed Jacobian Weighted temporal Total Variation results in better reconstructions of highly undersampled cardiac cine MRI than previously proposed methods and sets a theoretical ground for forward and backward predictors used elsewhere. Magn Reson Med 77:1208-1215, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Nonrigid groupwise registration for motion estimation and compensation in compressed sensing reconstruction of breath-hold cardiac cine MRI.

PURPOSE: Compressed sensing methods with motion estimation and compensation techniques have been proposed for the reconstruction of accelerated dynamic MRI. However, artifacts that naturally arise in compressed sensing reconstruction procedures hinder the estimation of motion from reconstructed images, especially at high acceleration factors. This work introduces a robust groupwise nonrigid motion estimation technique applied to the compressed sensing reconstruction of dynamic cardiac cine MRI sequences. THEORY AND METHODS: A spatio-temporal regularized, groupwise, nonrigid registration method based on a B-splines deformation model and a least squares metric is used to estimate and to compensate the movement of the heart in breath-hold cine acquisitions and to obtain a quasistatic sequence with highly sparse representation in temporally transformed domains. RESULTS: Short axis in vivo datasets are used for validation, both original multicoil as well as DICOM data. Fully sampled data were retrospectively undersampled with various acceleration factors and reconstructions were compared with the two well-known methods k-t FOCUSS and MASTeR. The proposed method achieves higher signal to error ratio and structure similarity index for medium to high acceleration factors. CONCLUSIONS: Reconstruction methods based on groupwise registration show higher quality reconstructions for cardiac cine images than the pairwise counterparts tested.

Computing thickness of irregularly-shaped thin walls using a locally semi-implicit scheme with extrapolation to solve the Laplace equation: Application to the right ventricle.

Cardiac Magnetic Resonance (CMR) Imaging is currently considered the gold standard imaging modality in cardiology. However, it is accompanied by a tradeoff between spatial resolution and acquisition time. Providing accurate measures of thin walls relative to the image resolution may prove challenging. One such anatomical structure is the cardiac right ventricle. Methods for measuring thickness of wall-like anatomical structures often rely on the Laplace equation to provide point-to-point correspondences between both boundaries. This work presents limex, a novel method to solve the Laplace equation using ghost nodes and providing extrapolated values, which is tested on three different datasets: a mathematical phantom, a set of biventricular segmentations from CMR images of ten pigs and the database used at the RV Segmentation Challenge held at MICCAI'12. Thickness measurements using the proposed methodology are more accurate than state-of-the-art methods, especially with the coarsest image resolutions, yielding mean L1 norms of the error between 43.28% and 86.52% lower than the second-best methods on the different test datasets. It is also computationally affordable. Limex has outperformed other state-of-the-art methods in classifying RV myocardial segments by their thickness.

Fast 4D elastic group-wise image registration. Convolutional interpolation revisited.

BACKGROUND AND OBJECTIVE: This paper proposes a new and highly efficient implementation of 3D+t groupwise registration based on the free-form deformation paradigm. METHODS: Deformation is posed as a cascade of 1D convolutions, achieving great reduction in execution time for evaluation of transformations and gradients. RESULTS: The proposed method has been applied to 4D cardiac MRI and 4D thoracic CT monomodal datasets. Results show an average runtime reduction above 90%, both in CPU and GPU executions, compared with the classical tensor product formulation. CONCLUSIONS: Our implementation, although fully developed for the metric sum of squared differences, can be extended to other metrics and its adaptation to multiresolution strategies is straightforward. Therefore, it can be extremely useful to speed up image registration procedures in different applications where high dimensional data are involved.

Time-efficient three-dimensional transmural scar assessment provides relevant substrate characterization for ventricular tachycardia features and long-term recurrences in ischemic cardiomyopathy.

Delayed gadolinium-enhanced cardiac magnetic resonance (LGE-CMR) imaging requires novel and time-efficient approaches to characterize the myocardial substrate associated with ventricular arrhythmia in patients with ischemic cardiomyopathy. Using a translational approach in pigs and patients with established myocardial infarction, we tested and validated a novel 3D methodology to assess ventricular scar using custom transmural criteria and a semiautomatic approach to obtain transmural scar maps in ventricular models reconstructed from both 3D-acquired and 3D-upsampled-2D-acquired LGE-CMR images. The results showed that 3D-upsampled models from 2D LGE-CMR images provided a time-efficient alternative to 3D-acquired sequences to assess the myocardial substrate associated with ischemic cardiomyopathy. Scar assessment from 2D-LGE-CMR sequences using 3D-upsampled models was superior to conventional 2D assessment to identify scar sizes associated with the cycle length of spontaneous ventricular tachycardia episodes and long-term ventricular tachycardia recurrences after catheter ablation. This novel methodology may represent an efficient approach in clinical practice after manual or automatic segmentation of myocardial borders in a small number of conventional 2D LGE-CMR slices and automatic scar detection.

Analysis of the pyramidal tract in tumor patients using diffusion tensor imaging.

In this work, we propose to use fiber tracking in order to analyze and quantify the state of the pyramidal tracts in patients affected by tumors. We introduce a framework that includes an automatic method to obtain the fibers involved in the pyramidal tract of any subject, in order to compare robustly fiber bundles affected by tumors with healthy fiber tracts from control subjects and also to quantify the relative state of degeneration between the fiber tracts in the two hemispheres of the same patient. The comparative analyses proposed in our methodology are based on a new set of measures on the pyramidal tract, which take into account intrinsic properties of the fibers involved in the bundle as well as the similarity with the pyramidal tract of a standard healthy subject, modeled as the average of a set of controls. In order to perform better comparison studies and to take into account more information of the whole bundle, a mapping technique is used to represent the fiber tracts in 2D. Here, we show a set of experiments using 5 tumor patients and 10 control subjects, including pre- and post-operative studies in patients that have been treated with partial or total tumor resection. The results obtained indicate the usefulness of the method showing good overall performance. A reproducibility study using several acquisitions of the same patient is also presented to validate the techniques employed.

Reconstruction techniques for cardiac cine MRI.

The present survey describes the state-of-the-art techniques for dynamic cardiac magnetic resonance image reconstruction. Additionally, clinical relevance, main challenges, and future trends of this image modality are outlined. Thus, this paper aims to provide a general vision about cine MRI as the standard procedure in functional evaluation of the heart, focusing on technical methodologies.

OpenCLIPER: An OpenCL-Based C++ Framework for Overhead-Reduced Medical Image Processing and Reconstruction on Heterogeneous Devices.

Medical image processing is often limited by the computational cost of the involved algorithms. Whereas dedicated computing devices (GPUs in particular) exist and do provide significant efficiency boosts, they have an extra cost of use in terms of housekeeping tasks (device selection and initialization, data streaming, synchronization with the CPU, and others), which may hinder developers from using them. This paper describes an OpenCL-based framework that is capable of handling dedicated computing devices seamlessly and that allows the developer to concentrate on image processing tasks. The framework handles automatically device discovery and initialization, data transfers to and from the device and the file system and kernel loading and compiling. Data structures need to be defined only once independently of the computing device; code is unique, consequently, for every device, including the host CPU. Pinned memory/buffer mapping is used to achieve maximum performance in data transfers. Code fragments included in the paper show how the computing device is almost immediately and effortlessly available to the users algorithms, so they can focus on productive work. Code required for device selection and initialization, data loading and streaming and kernel compilation is minimal and systematic. Algorithms can be thought of as mathematical operators (called processes), with input, output and parameters, and they may be chained one after another easily and efficiently. Also for efficiency, processes can have their initialization work split from their core workload, so process chains and loops do not incur in performance penalties. Algorithm code is independent of the device type targeted.

Space-time variant weighted regularization in compressed sensing cardiac cine MRI.

PURPOSE: To analyze the impact on image quality and motion fidelity of a motion-weighted space-time variant regularization term in compressed sensing cardiac cine MRI. METHODS: k-t SPARSE-SENSE with temporal total variation (tTV) is used as the base reconstruction algorithm. Motion in the dynamic image is estimated by means of a robust registration technique for non-rigid motion. The resulting deformation fields are used to leverage the regularization term. The results are compared with standard k-t SPARSE-SENSE with tTV regularization as well as with an improved version of this algorithm that makes use of tTV and temporal Fast Fourier Transform regularization in x-f domain. RESULTS: The proposed method with space-time variant regularization provides higher motion fidelity and image quality than the two previously reported methods. Difference images between undersampled reconstruction and fully sampled reference images show less systematic errors with the proposed approach. CONCLUSIONS: Usage of a space-time variant regularization offers reconstructions with better image quality than the state of the art approaches used for comparison.

Vortical features for myocardial rotation assessment in hypertrophic cardiomyopathy using cardiac tagged magnetic resonance.

Left ventricular rotational motion is a feature of normal and diseased cardiac function. However, classical torsion and twist measures rely on the definition of a rotational axis which may not exist. This paper reviews global and local rotation descriptors of myocardial motion and introduces new curl-based (vortical) features built from tensorial magnitudes, intended to provide better comprehension about fibrotic tissue characteristics mechanical properties. Fifty-six cardiomyopathy patients and twenty-two healthy volunteers have been studied using tagged magnetic resonance by means of harmonic phase analysis. Rotation descriptors are built, with no assumption about a regular geometrical model, from different approaches. The extracted vortical features have been tested by means of a sequential cardiomyopathy classification procedure; they have proven useful for the regional characterization of the left ventricular function by showing great separability not only between pathologic and healthy patients but also, and specifically, between heterogeneous phenotypes within cardiomyopathies.

A Second Order Multi-Stencil Fast Marching Method with a Non-Constant Local Cost Model.

The Fast Marching method is widely employed in several fields of image processing. Some years ago a Multi-Stencil version (MSFM) was introduced to improve its accuracy by solving the equation for a set of stencils and choosing the best solution at each considered node. The following work proposes a modified numerical scheme for MSFM to take into account the variation of the local cost, which has proven to be second order. The influence of the stencil set choice on the algorithm outcome with respect to stencil orthogonality and axis swapping is also explored, where stencils are taken from neighborhoods of varying radius. The experimental results show that the proposed schemes improve the accuracy of their original counterparts, and that the use of permutation-invariant stencil sets provides robustness against shifted vector coordinates in the stencil set.

Sequential anisotropic Wiener filtering applied to 3D MRI data.

We present three different sequential Wiener filters, namely, isotropic, orientation and anisotropic. The first one is similar to the classical Wiener filter in the sense that it uses an isotropic neighborhood to estimate its parameters. Here we present a sequential version of it. The orientation Wiener filter uses oriented neighborhoods to estimate the structure orientation present at each voxel, giving rise to a modified estimator of the parameters. Finally, the anisotropic Wiener filter combines both approaches adaptively so that the appropriate approach is locally selected. Several synthetic experiments are presented showing the performance of the filters with respect to their parameters. A mean square error analysis is performed using a publicly available magnetic resonance imaging (MRI) brain phantom and a comparison with other filtering approaches is carried out. In addition, results from filtering real MRI data are presented.

Weaning from mechanical ventilation: a retrospective analysis leading to a multimodal perspective.

Practitioners' decision for mechanical aid discontinuation is a challenging task that involves a complete knowledge of a great number of clinical parameters, as well as its evolution in time. Recently, an increasing interest on respiratory pattern variability as an extubation readiness indicator has appeared. Reliable assessment of this variability involves a set of signal processing and pattern recognition techniques. This paper presents a suitability analysis of different methods used for breathing pattern complexity assessment. The contribution of this analysis is threefold: 1) to serve as a review of the state of the art on the so-called weaning problem from a signal processing point of view; 2) to provide insight into the applied processing techniques and how they fit into the problem; 3) to propose additional methods and further processing in order to improve breathing pattern regularity assessment and weaning readiness decision. Results on experimental data show that sample entropy outperforms other complexity assessment methods and that multidimensional classification does improve weaning prediction. However, the obtained performance may be objectionable for real clinical practice, a fact that paves the way for a multimodal signal processing framework, including additional high-quality signals and more reliable statistical methods.

Multi-oriented windowed harmonic phase reconstruction for robust cardiac strain imaging.

The purpose of this paper is to develop a method for direct estimation of the cardiac strain tensor by extending the harmonic phase reconstruction on tagged magnetic resonance images to obtain more precise and robust measurements. The extension relies on the reconstruction of the local phase of the image by means of the windowed Fourier transform and the acquisition of an overdetermined set of stripe orientations in order to avoid the phase interferences from structures outside the myocardium and the instabilities arising from the application of a gradient operator. Results have shown that increasing the number of acquired orientations provides a significant improvement in the reproducibility of the strain measurements and that the acquisition of an extended set of orientations also improves the reproducibility when compared with acquiring repeated samples from a smaller set of orientations. Additionally, biases in local phase estimation when using the original harmonic phase formulation are greatly diminished by the one here proposed. The ideas here presented allow the design of new methods for motion sensitive magnetic resonance imaging, which could simultaneously improve the resolution, robustness and accuracy of motion estimates.

An approach for contour detection of human kidneys from ultrasound images using Markov random fields and active contours.

In this paper, a novel method for the boundary detection of human kidneys from three dimensional (3D) ultrasound (US) is proposed. The inherent difficulty of interpretation of such images, even by a trained expert, makes the problem unsuitable for classical methods. The method here proposed finds the kidney contours in each slice. It is a probabilistic Bayesian method. The prior defines a Markov field of deformations and imposes the restriction of contour smoothness. The likelihood function imposes a probabilistic behavior to the data, conditioned to the contour position. This second function, which is also Markov, uses an empirical model of distribution of the echographical data and a function of the gradient of the data. The model finally includes, as a volumetric extension of the prior, a term that forces smoothness along the depth coordinate. The experiments that have been carried out on echographies from real patients validate the model here proposed. A sensitivity analysis of the model parameters has also been carried out.

On low order moments of the homodyned-K distribution.

Fractional low order moments have been reported as beneficial for sampling computations using the K distribution. However, it has been recently pointed out that this it not the case for the homodyned-K distribution for a tissue discrimination problem. In this paper we show that such an statement is not fully justified. To that end, we follow a standard pattern recognition procedure both to determine class separability measures and to classify data with several classifiers. We conclude that the optimum order of the moments is intimately linked to the specific statistical properties of the tissues to be discriminated. Some ideas on how to choose the optimum order are discussed.

The EM Method in a Probabilistic Wavelet-Based MRI Denoising.

Human body heat emission and others external causes can interfere in magnetic resonance image acquisition and produce noise. In this kind of images, the noise, when no signal is present, is Rayleigh distributed and its wavelet coefficients can be approximately modeled by a Gaussian distribution. Noiseless magnetic resonance images can be modeled by a Laplacian distribution in the wavelet domain. This paper proposes a new magnetic resonance image denoising method to solve this fact. This method performs shrinkage of wavelet coefficients based on the conditioned probability of being noise or detail. The parameters involved in this filtering approach are calculated by means of the expectation maximization (EM) method, which avoids the need to use an estimator of noise variance. The efficiency of the proposed filter is studied and compared with other important filtering techniques, such as Nowak's, Donoho-Johnstone's, Awate-Whitaker's, and nonlocal means filters, in different 2D and 3D images.

Anisotropic diffusion filter with memory based on speckle statistics for ultrasound images.

Ultrasound (US) imaging exhibits considerable difficulties for medical visual inspection and for development of automatic analysis methods due to speckle, which negatively affects the perception of tissue boundaries and the performance of automatic segmentation methods. With the aim of alleviating the effect of speckle, many filtering techniques are usually considered as a preprocessing step prior to automatic analysis methods or visual inspection. Most of the state-of-the-art filters try to reduce the speckle effect without considering its relevance for the characterization of tissue nature. However, the speckle phenomenon is the inherent response of echo signals in tissues and can provide important features for clinical purposes. This loss of information is even magnified due to the iterative process of some speckle filters, e.g., diffusion filters, which tend to produce over-filtering because of the progressive loss of relevant information for diagnostic purposes during the diffusion process. In this paper, we propose an anisotropic diffusion filter with a probabilistic-driven memory mechanism to overcome the over-filtering problem by following a tissue selective philosophy. In particular, we formulate the memory mechanism as a delay differential equation for the diffusion tensor whose behavior depends on the statistics of the tissues, by accelerating the diffusion process in meaningless regions and including the memory effect in regions where relevant details should be preserved. Results both in synthetic and real US images support the inclusion of the probabilistic memory mechanism for maintaining clinical relevant structures, which are removed by the state-of-the-art filters.