HYDI-DSI revisited: Constrained non-parametric EAP imaging without q-space re-gridding

Producción CientíficaHybrid Diffusion Imaging (HYDI) was one of the first attempts to use multi-shell samplings of the q-space to infer diffusion properties beyond Diffusion Tensor Imaging (DTI) or High Angular Resolution Diffusion Imaging (HARDI). HYDI was intended as a flexible protocol embedding both DTI (for lower -values) and HARDI (for higher -values) processing, as well as Diffusion Spectrum Imaging (DSI) when the entire data set was exploited. In the latter case, the spherical sampling of the q-space is re-gridded by interpolation to a Cartesian lattice whose extent covers the range of acquired b-values, hence being acquisition-dependent. The Discrete Fourier Transform (DFT) is afterwards used to compute the corresponding Cartesian sampling of the Ensemble Average Propagator (EAP) in an entirely non-parametric way. From this lattice, diffusion markers such as the Return To Origin Probability (RTOP) or the Mean Squared Displacement (MSD) can be numerically estimated. We aim at re-formulating this scheme by means of a Fourier Transform encoding matrix that eliminates the need for q-space re-gridding at the same time it preserves the non-parametric nature of HYDI-DSI. The encoding matrix is adaptively designed at each voxel according to the underlying DTI approximation, so that an optimal sampling of the EAP can be pursued without being conditioned by the particular acquisition protocol. The estimation of the EAP is afterwards carried out as a regularized Quadratic Programming (QP) problem, which allows to impose positivity constraints that cannot be trivially embedded within the conventional HYDI-DSI. We demonstrate that the definition of the encoding matrix in the adaptive space allows to analytically (as opposed to numerically) compute several popular descriptors of diffusion with the unique source of error being the cropping of high frequency harmonics in the Fourier analysis of the attenuation signal. They include not only RTOP and MSD, but also Return to Axis/Plane Probabilities (RTAP/RTPP), which are defined in terms of specific spatial directions and are not available with the former HYDI-DSI. We report extensive experiments that suggest the benefits of our proposal in terms of accuracy, robustness and computational efficiency, especially when only standard, non-dedicated q-space samplings are available.Ministerio de Ciencia e Innovación (PID2021-124407NB-I00 and TED2021-130758B-I00)Ministry of Science and Higher Education (Poland) (PPN/BEK/ 2019/1/00421

Efficient estimation of propagator anisotropy and non‐Gaussianity in multishell diffusion MRI with micro‐structure adaptive convolution kernels and dual Fourier integral transforms

Producción CientíficaPurpose:We seek to reformulate the so-called Propagator Anisotropy (PA) andNon-Gaussianity (NG), originally conceived for the Mean Apparent Propagatordiffusion MRI (MAP-MRI), to the Micro-Structure adaptive convolution ker-nels and dual Fourier Integral Transforms (MiSFIT). These measures describerelevant normalized features of the Ensemble Average Propagator (EAP).Theory and Methods:First, the indices, which are defined as the EAP’sdissimilarity from an isotropic (PA) or a Gaussian (NG) one, are analyticallyreformulated within the MiSFIT framework. Then a comparison between theresulting maps is drawn by means of a visual analysis, a quantitative assess-ment via numerical simulations, a test-retest study across the MICRA dataset (6subjects scanned five times) and, finally, a computational time evaluation.Results:Findings illustrate the visual similarity between the indices computedwith either technique. Evaluation against synthetic ground truth data, however,demonstrates MiSFIT’s improved accuracy. In addition, the test–retest studyreveals MiSFIT’s higher degree of reliability in most of white matter regions.Finally, the computational time evaluation shows MiSFIT’s time reduction upto two orders of magnitude.Conclusions:Despite being a direct development on the MAP-MRI represen-tation, the PA and the NG can be reliably and efficiently computed withinMiSFIT’s framework. This, together with the previous findings in the originalMiSFIT’s article, could mean the difference that definitely qualifies diffusionMRI to be incorporated into regular clinical settings.Ministerio de Educación, Junta de Castilla y León y Fondo Social Europeo, (Grant/Award Number: OrdenEDU/1100/2017 12/12)Ministerio de Ciencia e Innovación, Grant/AwardNumbers: (RTI2018-094569-B-I00),(PID2021-124407NB-I00)Ministry of Science and Higher Education of Poland,(Grant/Award Number:692/STYP/13/2018)Narodowa Agencja Wymiany Akademickiej, (Grant/AwardNumber: PPN/BEK/2019/1/00421

Spatially-variant noise filtering in magnetic resonance imaging : a consensus-based approach

In order to accelerate the acquisition process in multiple-coil Magnetic Resonance scanners, parallel techniques were developed. These techniques reduce the acquisition time via a sub-sampling of the k-space and a reconstruction process. From a signal and noise perspective, the use of a acceleration techniques modify the structure of the noise within the image. In the most common algorithms, like SENSE, the final magnitude image after the reconstruction is known to follow a Rician distribution for each pixel, just like single coil systems. However, the noise is spatially non-stationary, i.e. the variance of noise becomes x-dependent. This effect can also be found in magnitude images due to other processing inside the scanner. In this work we propose a method to adapt well-known noise filtering techniques initially designed to deal with stationary noise to the case of spatially variant Rician noise. The method copes with inaccurate estimates of variant noise patterns in the image, showing its robustness in realistic cases. The method employs a consensus strategy in conjunction with a set of aggregation functions and a penalty function. Multiple possible outputs are generated for each pixel assuming different unknown input parameters. The consensus approach merges them into a unique filtered image. As a filtering technique, we have selected the Linear Minimum Mean Square Error (LMMSE) estimator for Rician data, which has been used to test our methodology due to its simplicity and robustness. Results with synthetic and in vivo data confirm the good behavior of our approach

Anisotropic Diffusion Filter with Memory based on Speckle Statistics for Ultrasound Images

Ultrasound imaging exhibits considerable difficulties for medical visual inspection and for the 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 work, we propose an anisotropic diffusion filter with a probabilistic-driven memory mechanism to overcome the over-filtering problem by following a tissue selective philosophy. Specifically, 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

Apparent propagator anisotropy from single-shell diffusion MRI acquisitions

Purpose The apparent propagator anisotropy (APA) is a new diffusion MRI metric that, while drawing on the benefits of the ensemble averaged propagator anisotropy (PA) compared to the fractional anisotropy (FA), can be estimated from single‐shell data. Theory and Methods Computation of the full PA requires acquisition of large datasets with many diffusion directions and different b‐values, and results in extremely long processing times. This has hindered adoption of the PA by the community, despite evidence that it provides meaningful information beyond the FA. Calculation of the complete propagator can be avoided under the hypothesis that a similar sensitivity/specificity may be achieved from apparent measurements at a given shell. Assuming that diffusion anisotropy (DiA) is nondependent on the b‐value, a closed‐form expression using information from one single shell (ie, b‐value) is reported. Results Publicly available databases with healthy and diseased subjects are used to compare the APA against other anisotropy measures. The structural information provided by the APA correlates with that provided by the PA for healthy subjects, while it also reveals statistically relevant differences in white matter regions for two pathologies, with a higher reliability than the FA. Additionally, APA has a computational complexity similar to the FA, with processing‐times several orders of magnitude below the PA. Conclusions The APA can extract more relevant white matter information than the FA, without any additional demands on data acquisition. This makes APA an attractive option for adoption into existing diffusion MRI analysis pipelines

Robust Estimation of the Apparent Diffusion Coefficient Invariant to Acquisition Noise and Physiological Motion

Purpose: In this work we have proposed a methodology for the estimation of the apparent diffusion coeffcient in the body from multiple breath hold diffusion weighted images, which is robust to two preeminent confounding factors: noise and motion during acquisition. Methods: We have extended a method for the joint groupwise multimodal registration and apparent diffusion coefficient estimation, previously proposed by the authors, in order to correct the bias that arises from the non-Gaussianity of the data and the registration procedure. Results: Results show that the proposed methodology provides a statistically signi ficant improvement both in robustness for displacement elds calculation and in terms of accuracy for the apparent diffusion coefficient estimation as compared with traditional sequential approaches. Reproducibility has also been measured on real data in terms of the distribution of apparent diffusion coefficient differences obtained from different b-values subsets. Conclusions: Our proposal has shown to be able to effectively correct the estimation bias by introducing additional computationally light procedures to the original method, thus providing robust apparent diffusion coefficient maps in the liver and allowing an accurate and reproducible analysis of the tissue

Micro-structure diffusion scalar measures from reduced MRI acquisitions

In diffusion MRI, the Ensemble Average diffusion Propagator (EAP) provides relevant microstructural information and meaningful descriptive maps of the white matter previously obscured by traditional techniques like the Diffusion Tensor. The direct estimation of the EAP, however, requires a dense sampling of the Cartesian q-space. Due to the huge amount of samples needed for an accurate reconstruction, more efficient alternative techniques have been proposed in the last decade. Even so, all of them imply acquiring a large number of diffusion gradients with different b-values. In order to use the EAP in practical studies, scalar measures must be directly derived, being the most common the return-to-origin probability (RTOP) and the return-to-plane and return-to-axis probabilities (RTPP, RTAP). In this work, we propose the so-called “Apparent Measures Using Reduced Acquisitions” (AMURA) to drastically reduce the number of samples needed for the estimation of diffusion properties. AMURA avoids the calculation of the whole EAP by assuming the diffusion anisotropy is roughly independent from the radial direction. With such an assumption, and as opposed to common multi-shell procedures based on iterative optimization, we achieve closed-form expressions for the measures using information from one single shell. This way, the new methodology remains compatible with standard acquisition protocols commonly used for HARDI (based on just one b-value). We report extensive results showing the potential of AMURA to reveal microstructural properties of the tissues compared to state of the art EAP estimators, and is well above that of Diffusion Tensor techniques. At the same time, the closed forms provided for RTOP, RTPP, and RTAP-like magnitudes make AMURA both computationally efficient and robust

Influence of ultrasound speckle tracking strategies for motion and strain estimation

Speckle Tracking is one of the most prominent techniques used to estimate the regional movement of the heart based on ultrasound acquisitions. Many different approaches have been proposed, proving their suitability to obtain quantitative and qualitative information regarding myocardial deformation, motion and function assessment. New proposals to improve the basic algorithm usually focus on one of these three steps: (1) the similarity measure between images and the speckle model; (2) the transformation model, i.e. the type of motion considered between images; (3) the optimization strategies, such as the use of different optimization techniques in the transformation step or the inclusion of structural information. While many contributions have shown their good performance independently, it is not always clear how they perform when integrated in a whole pipeline. Every step will have a degree of influence over the following and hence over the final result. Thus, a Speckle Tracking pipeline must be analyzed as a whole when developing novel methods, since improvements in a particular step might be undermined by the choices taken in further steps. This work presents two main contributions: (1) We provide a complete analysis of the influence of the different steps in a Speckle Tracking pipeline over the motion and strain estimation accuracy. (2) The study proposes a methodology for the analysis of Speckle Tracking systems specifically designed to provide an easy and systematic way to include other strategies. We close the analysis with some conclusions and recommendations that can be used as an orientation of the degree of influence of the models for speckle, the transformation models, interpolation schemes and optimization strategies over the estimation of motion features. They can be further use to evaluate and design new strategy into a Speckle Tracking system