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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Impact of MR Acquisition Parameters on DTI Scalar Indexes: A Tractography Based Approach.

Acquisition parameters play a crucial role in Diffusion Tensor Imaging (DTI), as they have a major impact on the values of scalar measures such as Fractional Anisotropy (FA) or Mean Diffusivity (MD) that are usually the focus of clinical studies based on white matter analysis. This paper presents an analysis on the impact of the variation of several acquisition parameters on these scalar measures with a novel double focus. First, a tractography-based approach is employed, motivated by the significant number of clinical studies that are carried out using this technique. Second, the consequences of simultaneous changes in multiple parameters are analyzed: number of gradient directions, b-value and voxel resolution. Results indicate that the FA is most affected by changes in the number of gradients and voxel resolution, while MD is specially influenced by variations in the b-value. Even if the choice of a tractography algorithm has an effect on the numerical values of the final scalar measures, the evolution of these measures when acquisition parameters are modified is parallel.

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.

Gamma mixture classifier for plaque detection in intravascular ultrasonic images.

Carotid and coronary vascular incidents are mostly caused by vulnerable plaques. Detection and characterization of vulnerable plaques are important for early disease diagnosis and treatment. For this purpose, the echomorphology and composition have been studied. Several distributions have been used to describe ultrasonic data depending on tissues, acquisition conditions, and equipment. Among them, the Rayleigh distribution is a one-parameter model used to describe the raw envelope RF ultrasound signal for its simplicity, whereas the Nakagami distribution (a generalization of the Rayleigh distribution) is the two-parameter model which is commonly accepted. However, it fails to describe B-mode images or Cartesian interpolated or subsampled RF images because linear filtering changes the statistics of the signal. In this work, a gamma mixture model (GMM) is proposed to describe the subsampled/interpolated RF images and it is shown that the parameters and coefficients of the mixture are useful descriptors of speckle pattern for different types of plaque tissues. This new model outperforms recently proposed probabilistic and textural methods with respect to plaque description and characterization of echogenic contents. Classification results provide an overall accuracy of 86.56% for four classes and 95.16% for three classes. These results evidence the classifier usefulness for plaque characterization. Additionally, the classifier provides probability maps according to each tissue type, which can be displayed for inspecting local tissue composition, or used for automatic filtering and segmentation.

Multi-Stencil Streamline Fast Marching: A General 3-D Framework to Determine Myocardial Thickness and Transmurality in Late Enhancement Images.

We propose a fully 3-D methodology for the computation of myocardial nonviable tissue transmurality in contrast enhanced magnetic resonance images. The outcome is a continuous map defined within the myocardium where not only current state-of-the-art measures of transmurality can be calculated, but also information on the location of nonviable tissue is preserved. The computation is done by means of a partial differential equation framework we have called multi-stencil streamline fast marching. Using it, the myocardial and scarred tissue thickness is simultaneously computed. Experimental results show that the proposed 3-D method allows for the computation of transmurality in myocardial regions where current 2-D methods are not able to as conceived, and it also provides more robust and accurate results in situations where the assumptions on which current 2-D methods are based-i.e., there is a visible endocardial contour and its corresponding epicardial points lie on the same slice-, are not met.

A stochastic modelling framework for the reconstruction of cardiovascular signals.

This paper presents a common stochastic modelling framework for physiological signals which allows patient simulation following a synthesis-by-analysis approach. Within this framework, we propose a general model-based methodology able to reconstruct missing or artifacted signal intervals in cardiovascular monitoring applications. The proposed model consists of independent stages which provide high flexibility to incorporate signals of different nature in terms of shape, cross-correlation and variability. The reconstruction methodology is based on model sampling and selection based on a wide range of boundary conditions, which include prior information. Results on real data show how the proposed methodology fits the particular approaches presented so far for electrocardiogram (ECG) reconstruction and how a simple extension within the framework can significantly improve their performance.

Stochastic modeling of the PPG signal: a synthesis-by-analysis approach with applications.

In this paper, we propose a stochastic model of photoplethysmographic signals that is able to synthesize an arbitrary number of other statistically equivalent signals to the one under analysis. To that end, we first preprocess the pulse signal to normalize and time-align pulses. In a second stage, we design a single-pulse model, which consists of ten parameters. In the third stage, the time evolution of this ten-parameter vector is approximated by means of two autoregressive moving average models, one for the trend and one for the residue; this model is applied after a decorrelation step which let us to process each vector component in parallel. The experiments carried out show that the model we here propose is able to maintain the main features of the original signal; this is accomplished by means of both a linear spectral analysis and also by comparing two measures obtained from a nonlinear analysis. Finally, we explore the capability of the model to: 1) track physical activity; 2) obtain statistics of clinical parameters by model sampling; and 3) recover corrupted or missing signal epochs by synthesis.