RESUMEN
Combining multi-site data can strengthen and uncover trends, but is a task that is marred by the influence of site-specific covariates that can bias the data and, therefore, any downstream analyses. Post-hoc multi-site correction methods exist but have strong assumptions that often do not hold in real-world scenarios. Algorithms should be designed in a way that can account for site-specific effects, such as those that arise from sequence parameter choices, and in instances where generalisation fails, should be able to identify such a failure by means of explicit uncertainty modelling. This body of work showcases such an algorithm that can become robust to the physics of acquisition in the context of segmentation tasks while simultaneously modelling uncertainty. We demonstrate that our method not only generalises to complete holdout datasets, preserving segmentation quality but does so while also accounting for site-specific sequence choices, which also allows it to perform as a harmonisation tool.
Asunto(s)
Imagen por Resonancia Magnética , Neuroimagen , Humanos , Incertidumbre , Imagen por Resonancia Magnética/métodos , Algoritmos , Encéfalo/diagnóstico por imagen , Procesamiento de Imagen Asistido por Computador/métodosRESUMEN
Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.
Asunto(s)
Axones , Imagen de Difusión por Resonancia Magnética , Ratas , Animales , Imagen de Difusión por Resonancia Magnética/métodos , Simulación por Computador , Procesamiento de Imagen Asistido por Computador/métodosRESUMEN
BACKGROUND: The incidence of major surgery is on the rise globally, and more than 20% of patients are readmitted to hospital following discharge from hospital. During their hospital stay, patients are monitored for early detection of clinical deterioration, which includes regularly measuring physiological parameters such as blood pressure, heart rate, respiratory rate, temperature, and pulse oximetry. This monitoring ceases upon hospital discharge, as patients are deemed clinically stable. Monitoring after discharge is relevant to detect adverse events occurring in the home setting and can be made possible through the development of digital technologies and mobile networks. Smartwatches and other technological devices allow patients to self-measure physiological parameters in the home setting, and Bluetooth connectivity can facilitate the automatic collection and transfer of this data to a secure server with minimal input from the patient. OBJECTIVE: This paper presents the protocol for the DREAMPath (Domiciliary Recovery After Medicalization Pathway) study, which aims to measure compliance with a multidevice remote monitoring kit after discharge from hospital following major surgery. METHODS: DREAMPath is a single-center, prospective, observational, cohort study, comprising 30 patients undergoing major intracavity surgery. The primary outcome is to assess patient compliance with wearable and interactive smart technology in the first 30 days following discharge from hospital after major surgery. Secondary outcomes will explore the relation between unplanned health care events and physiological data collected in the study, as well as to explore a similar relationship with daily patient-reported outcome measures (Quality of Recovery-15 score). Secondary outcomes will be analyzed using appropriate regression methods. Cardiopulmonary exercise testing data will also be collected to assess correlations with wearable device data. RESULTS: Recruitment was halted due to COVID-19 restrictions and will progress once research staff are back from redeployment. We expect that the study will be completed in the first quarter of 2022. CONCLUSIONS: Digital health solutions have been recently made possible due to technological advances, but urgency in rollout has been expedited due to COVID-19. The DREAMPath study will inform readers about the feasibility of remote monitoring for a patient group that is at an increased risk of acute deterioration. TRIAL REGISTRATION: ISRCTN Registry ISRCTN62293620; https://www.isrctn.com/ISRCTN62293620. INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID): DERR1-10.2196/30638.
RESUMEN
Fiber tractography is widely used to non-invasively map white-matter bundles in vivo using diffusion-weighted magnetic resonance imaging (dMRI). As it is the case for all scientific methods, proper validation is a key prerequisite for the successful application of fiber tractography, be it in the area of basic neuroscience or in a clinical setting. It is well-known that the indirect estimation of the fiber tracts from the local diffusion signal is highly ambiguous and extremely challenging. Furthermore, the validation of fiber tractography methods is hampered by the lack of a real ground truth, which is caused by the extremely complex brain microstructure that is not directly observable non-invasively and that is the basis of the huge network of long-range fiber connections in the brain that are the actual target of fiber tractography methods. As a substitute for in vivo data with a real ground truth that could be used for validation, a widely and successfully employed approach is the use of synthetic phantoms. In this work, we are providing an overview of the state-of-the-art in the area of physical and digital phantoms, answering the following guiding questions: "What are dMRI phantoms and what are they good for?", "What would the ideal phantom for validation fiber tractography look like?" and "What phantoms, phantom datasets and tools used for their creation are available to the research community?". We will further discuss the limitations and opportunities that come with the use of dMRI phantoms, and what future direction this field of research might take.
Asunto(s)
Imagen de Difusión Tensora/métodos , Fantasmas de Imagen , Sustancia Blanca/diagnóstico por imagen , Artefactos , Humanos , Procesamiento de Imagen Asistido por ComputadorRESUMEN
The intra-axonal water exchange time (τi), a parameter associated with axonal permeability, could be an important biomarker for understanding and treating demyelinating pathologies such as Multiple Sclerosis. Diffusion-Weighted MRI (DW-MRI) is sensitive to changes in permeability; however, the parameter has so far remained elusive due to the lack of general biophysical models that incorporate it. Machine learning based computational models can potentially be used to estimate such parameters. Recently, for the first time, a theoretical framework using a random forest (RF) regressor suggests that this is a promising new approach for permeability estimation. In this study, we adopt such an approach and for the first time experimentally investigate it for demyelinating pathologies through direct comparison with histology. We construct a computational model using Monte Carlo simulations and an RF regressor in order to learn a mapping between features derived from DW-MRI signals and ground truth microstructure parameters. We test our model in simulations, and find strong correlations between the predicted and ground truth parameters (intra-axonal volume fraction f: R2 =0.99, τi: R2 =0.84, intrinsic diffusivity d: R2 =0.99). We then apply the model in-vivo, on a controlled cuprizone (CPZ) mouse model of demyelination, comparing the results from two cohorts of mice, CPZ (N=8) and healthy age-matched wild-type (WT, N=8). We find that the RF model estimates sensible microstructure parameters for both groups, matching values found in literature. Furthermore, we perform histology for both groups using electron microscopy (EM), measuring the thickness of the myelin sheath as a surrogate for exchange time. Histology results show that our RF model estimates are very strongly correlated with the EM measurements (ρ = 0.98 for f, ρ = 0.82 for τi). Finally, we find a statistically significant decrease in τi in all three regions of the corpus callosum (splenium/genu/body) of the CPZ cohort (<τi>=310ms/330ms/350ms) compared to the WT group (<τi>=370ms/370ms/380ms). This is in line with our expectations that τi is lower in regions where the myelin sheath is damaged, as axonal membranes become more permeable. Overall, these results demonstrate, for the first time experimentally and in vivo, that a computational model learned from simulations can reliably estimate microstructure parameters, including the axonal permeability .
Asunto(s)
Axones/patología , Cuerpo Calloso/patología , Enfermedades Desmielinizantes/diagnóstico por imagen , Aprendizaje Automático , Sustancia Blanca/diagnóstico por imagen , Animales , Axones/metabolismo , Axones/ultraestructura , Simulación por Computador , Cuerpo Calloso/ultraestructura , Cuprizona/toxicidad , Enfermedades Desmielinizantes/inducido químicamente , Enfermedades Desmielinizantes/patología , Imagen de Difusión por Resonancia Magnética , Modelos Animales de Enfermedad , Procesamiento de Imagen Asistido por Computador , Ratones , Microscopía Electrónica , Inhibidores de la Monoaminooxidasa/toxicidad , Método de Montecarlo , Permeabilidad , Sustancia Blanca/patologíaRESUMEN
PURPOSE: Time-dependence is a key feature of the diffusion-weighted (DW) signal, knowledge of which informs biophysical modelling. Here, we study time-dependence in the human spinal cord, as its axonal structure is specific and different from the brain. METHODS: We run Monte Carlo simulations using a synthetic model of spinal cord white matter (WM) (large axons), and of brain WM (smaller axons). Furthermore, we study clinically feasible multi-shell DW scans of the cervical spinal cord (b = 0; b = 711 s mm-2 ; b = 2855 s mm-2 ), obtained using three diffusion times (Δ of 29, 52 and 76 ms) from three volunteers. RESULTS: Both intra-/extra-axonal perpendicular diffusivities and kurtosis excess show time-dependence in our synthetic spinal cord model. This time-dependence is reflected mostly in the intra-axonal perpendicular DW signal, which also exhibits strong decay, unlike our brain model. Time-dependence of the total DW signal appears detectable in the presence of noise in our synthetic spinal cord model, but not in the brain. In WM in vivo, we observe time-dependent macroscopic and microscopic diffusivities and diffusion kurtosis, NODDI and two-compartment SMT metrics. Accounting for large axon calibers improves fitting of multi-compartment models to a minor extent. CONCLUSIONS: Time-dependence of clinically viable DW MRI metrics can be detected in vivo in spinal cord WM, thus providing new opportunities for the non-invasive estimation of microstructural properties. The time-dependence of the perpendicular DW signal may feature strong intra-axonal contributions due to large spinal axon caliber. Hence, a popular model known as "stick" (zero-radius cylinder) may be sub-optimal to describe signals from the largest spinal axons.
Asunto(s)
Axones/patología , Imagen de Difusión por Resonancia Magnética , Médula Espinal/diagnóstico por imagen , Adulto , Algoritmos , Encéfalo/diagnóstico por imagen , Simulación por Computador , Femenino , Voluntarios Sanos , Humanos , Procesamiento de Imagen Asistido por Computador/métodos , Masculino , Método de Montecarlo , Factores de TiempoRESUMEN
Microscopic diffusion anisotropy (µA) has been recently gaining increasing attention for its ability to decouple the average compartment anisotropy from orientation dispersion. Advanced diffusion MRI sequences, such as double diffusion encoding (DDE) and double oscillating diffusion encoding (DODE) have been used for mapping µA, usually using measurements from a single b shell. However, the accuracy of µA estimation vis-à-vis different b-values was not assessed. Moreover, the time-dependence of this metric, which could offer additional insights into tissue microstructure, has not been studied so far. Here, we investigate both these concepts using theory, simulation, and experiments performed at 16.4T in the mouse brain, ex-vivo. In the first part, simulations and experimental results show that the conventional estimation of microscopic anisotropy from the difference of D(O)DE sequences with parallel and orthogonal gradient directions yields values that highly depend on the choice of b-value. To mitigate this undesirable bias, we propose a multi-shell approach that harnesses a polynomial fit of the signal difference up to third order terms in b-value. In simulations, this approach yields more accurate µA metrics, which are similar to the ground-truth values. The second part of this work uses the proposed multi-shell method to estimate the time/frequency dependence of µA. The data shows either an increase or no change in µA with frequency depending on the region of interest, both in white and gray matter. When comparing the experimental results with simulations, it emerges that simple geometric models such as infinite cylinders with either negligible or finite radii cannot replicate the measured trend, and more complex models, which, for example, incorporate structure along the fibre direction are required. Thus, measuring the time dependence of microscopic anisotropy can provide valuable information for characterizing tissue microstructure.
Asunto(s)
Encéfalo/anatomía & histología , Imagen de Difusión por Resonancia Magnética/métodos , Procesamiento de Imagen Asistido por Computador/métodos , Animales , Anisotropía , Encéfalo/fisiología , Ratones , Ratones Endogámicos C57BLRESUMEN
Quality control (QC) is a fundamental component of any study. Diffusion MRI has unique challenges that make manual QC particularly difficult, including a greater number of artefacts than other MR modalities and a greater volume of data. The gold standard is manual inspection of the data, but this process is time-consuming and subjective. Recently supervised learning approaches based on convolutional neural networks have been shown to be competitive with manual inspection. A drawback of these approaches is they still require a manually labelled dataset for training, which is itself time-consuming to produce and still introduces an element of subjectivity. In this work we demonstrate the need for manual labelling can be greatly reduced by training on simulated data, and using a small amount of labelled data for a final calibration step. We demonstrate its potential for the detection of severe movement artefacts, and compare performance to a classifier trained on manually-labelled real data.
Asunto(s)
Artefactos , Mapeo Encefálico/métodos , Procesamiento de Imagen Asistido por Computador/métodos , Control de Calidad , Aprendizaje Automático Supervisado , Mapeo Encefálico/normas , Conectoma/métodos , Imagen de Difusión por Resonancia Magnética/métodos , Imagen de Difusión por Resonancia Magnética/normas , Femenino , Humanos , Procesamiento de Imagen Asistido por Computador/normas , Recién Nacido , Masculino , Aprendizaje Automático Supervisado/normasRESUMEN
Mapping axon diameters within the central and peripheral nervous system could play an important role in our understanding of nerve pathways, and help diagnose and monitor an array of neurological disorders. Numerous diffusion MRI methods have been proposed for imaging axon diameters, most of which use conventional single diffusion encoding (SDE) spin echo sequences. However, a growing number of studies show that oscillating gradient spin echo (OGSE) sequences can provide additional advantages over conventional SDE sequences. Recent theoretical results suggest that this is especially the case in realistic scenarios, such as when fibres have unknown or dispersed orientation. In the present study, we adopt the ActiveAx approach to experimentally investigate the extent of these advantages by comparing the performances of SDE and trapezoidal OGSE in viable nerve tissue. We optimise SDE and OGSE ActiveAx protocols for a rat peripheral nerve tissue and test their performance using Monte Carlo simulations and a 800 mT/m gradient strength pre-clinical imaging experiment. The imaging experiment uses excised sciatic nerve from a rat's leg placed in a MRI compatible viable isolated tissue (VIT) maintenance chamber, which keeps the tissue in a viable physiological state that preserves the structural complexity of the nerve and enables lengthy scan times. We compare model estimates to histology, which we perform on the nerve post scanning. Optimisation produces a three-shell SDE and OGSE ActiveAx protocol, with the OGSE protocol consisting of one SDE sequence and two low-frequency oscillating gradient waveform sequences. Both simulation and imaging results show that the OGSE ActiveAx estimates of the axon diameter index have a higher accuracy and a higher precision compared to those from SDE. Histology estimates of the axon diameter index in our nerve tissue samples are 4-5.8 µm and these are excellently matched with the OGSE estimates 4.2-6.5 µm, while SDE overestimates at 5.2-8 µm for the same sample. We found OGSE estimates to be more precise with on average a 0.5 µm standard deviation compared to the SDE estimates which have a 2 µm standard deviation. When testing the robustness of the estimates when the number of the diffusion gradient directions reduces, we found that both OGSE and SDE estimates are affected, however OGSE is more robust to these changes than the SDE. Overall, these results suggest, quantitatively and in in vivo conditions, that low-frequency OGSE sequences may provide improved accuracy of axon diameter mapping compared to standard SDE sequences.
Asunto(s)
Axones , Imagen por Resonancia Magnética/métodos , Neuroimagen/métodos , Nervio Ciático/diagnóstico por imagen , Animales , Simulación por Computador , Imagen de Difusión por Resonancia Magnética/métodos , Imagen de Difusión por Resonancia Magnética/normas , Imagen por Resonancia Magnética/normas , Método de Montecarlo , Neuroimagen/normas , Ratas , Ratas Sprague-Dawley , Sensibilidad y EspecificidadRESUMEN
Because of their low bandwidth in the phase-encode (PE) direction, the susceptibility-induced off-resonance field causes distortions in echo planar imaging (EPI) images. It is therefore crucial to correct for susceptibility-induced distortions when performing diffusion studies using EPI. The susceptibility-induced field is caused by the object (head) disrupting the field and it is typically assumed that it remains constant within a framework defined by the object, (i.e. it follows the object as it moves in the scanner). However, this is only approximately true. When a non-spherical object rotates around an axis other than that parallel with the magnetic flux (the z-axis) it changes the way it disrupts the field, leading to different distortions. Hence, if using a single field to correct for distortions there will be residual distortions in the volumes where the object orientation is substantially different to that when the field was measured. In this paper we present a post-processing method for estimating the field as it changes with motion during the course of an experiment. It only requires a single measured field and knowledge of the orientation of the subject when that field was acquired. The volume-to-volume changes of the field as a consequence of subject movement are estimated directly from the diffusion data without the need for any additional or special acquisitions. It uses a generative model that predicts how each volume would look predicated on field change and inverts that model to yield an estimate of the field changes. It has been validated on both simulations and experimental data. The results show that we are able to track the field with high accuracy and that we are able to correct the data for the adverse effects of the changing field.
Asunto(s)
Artefactos , Mapeo Encefálico/métodos , Imagen Eco-Planar/métodos , Procesamiento de Imagen Asistido por Computador/métodos , Algoritmos , Movimientos de la Cabeza , Humanos , Movimiento (Física)RESUMEN
In this paper we evaluate the three main methods for correcting the susceptibility-induced artefact in diffusion-weighted magnetic-resonance (DW-MR) data, and assess how correction is affected by the susceptibility field's interaction with motion. The susceptibility artefact adversely impacts analysis performed on the data and is typically corrected in post-processing. Correction strategies involve either registration to a structural image, the application of an acquired field-map or the use of additional images acquired with different phase-encoding. Unfortunately, the choice of which method to use is made difficult by the absence of any systematic comparisons of them. In this work we quantitatively evaluate these methods, by extending and employing a recently proposed framework that allows for the simulation of realistic DW-MR datasets with artefacts. Our analysis separately evaluates the ability for methods to correct for geometric distortions and to recover lost information in regions of signal compression. In terms of geometric distortions, we find that registration-based methods offer the poorest correction. Field-mapping techniques are better, but are influenced by noise and partial volume effects, whilst multiple phase-encode methods performed best. We use our simulations to validate a popular surrogate metric of correction quality, the comparison of corrected data acquired with AP and LR phase-encoding, and apply this surrogate to real datasets. Furthermore, we demonstrate that failing to account for the interaction of the susceptibility field with head movement leads to increased errors when analysing DW-MR data. None of the commonly used post-processing methods account for this interaction, and we suggest this may be a valuable area for future methods development.
Asunto(s)
Imagen de Difusión por Resonancia Magnética/métodos , Artefactos , Encéfalo/diagnóstico por imagen , Humanos , Interpretación de Imagen Asistida por Computador , Movimiento (Física)RESUMEN
Most motion correction methods work by aligning a set of volumes together, or to a volume that represents a reference location. These are based on an implicit assumption that the subject remains motionless during the several seconds it takes to acquire all slices in a volume, and that any movement occurs in the brief moment between acquiring the last slice of one volume and the first slice of the next. This is clearly an approximation that can be more or less good depending on how long it takes to acquire one volume and in how rapidly the subject moves. In this paper we present a method that increases the temporal resolution of the motion correction by modelling movement as a piecewise continous function over time. This intra-volume movement correction is implemented within a previously presented framework that simultaneously estimates distortions, movement and movement-induced signal dropout. We validate the method on highly realistic simulated data containing all of these effects. It is demonstrated that we can estimate the true movement with high accuracy, and that scalar parameters derived from the data, such as fractional anisotropy, are estimated with greater fidelity when data has been corrected for intra-volume movement. Importantly, we also show that the difference in fidelity between data affected by different amounts of movement is much reduced when taking intra-volume movement into account. Additional validation was performed on data from a healthy volunteer scanned when lying still and when performing deliberate movements. We show an increased correspondence between the "still" and the "movement" data when the latter is corrected for intra-volume movement. Finally we demonstrate a big reduction in the telltale signs of intra-volume movement in data acquired on elderly subjects.
Asunto(s)
Mapeo Encefálico/métodos , Encéfalo/anatomía & histología , Imagen de Difusión por Resonancia Magnética/métodos , Aumento de la Imagen/métodos , Artefactos , Simulación por Computador , Humanos , Movimiento , Reproducibilidad de los ResultadosRESUMEN
Diffusion MRI has been proposed as a non-invasive technique for axonal diameter mapping. However, accurate estimation of small diameters requires strong gradients, which is a challenge for the transition of the technique from preclinical to clinical MRI scanners, since these have weaker gradients. In this work, we develop a framework to estimate the lower bound for accurate diameter estimation, which we refer to as the resolution limit. We analyse only the contribution from the intra-axonal space and assume that axons can be represented by impermeable cylinders. To address the growing interest in using techniques for diffusion encoding that go beyond the conventional single diffusion encoding (SDE) sequence, we present a generalised analysis capable of predicting the resolution limit regardless of the gradient waveform. Using this framework, waveforms were optimised to minimise the resolution limit. The results show that, for parallel cylinders, the SDE experiment is optimal in terms of yielding the lowest possible resolution limit. In the presence of orientation dispersion, diffusion encoding sequences with square-wave oscillating gradients were optimal. The resolution limit for standard clinical MRI scanners (maximum gradient strength 60-80 mT/m) was found to be between 4 and 8 µm, depending on the noise levels and the level of orientation dispersion. For scanners with a maximum gradient strength of 300 mT/m, the limit was reduced to between 2 and 5 µm.
Asunto(s)
Axones/ultraestructura , Imagen de Difusión Tensora/métodos , Interpretación de Imagen Asistida por Computador/métodos , Modelos Neurológicos , Sustancia Blanca/citología , Sustancia Blanca/diagnóstico por imagen , Algoritmos , Anisotropía , Simulación por Computador , Humanos , Aumento de la Imagen/métodos , Campos Magnéticos , Modelos Anatómicos , Dosis de Radiación , Reproducibilidad de los Resultados , Sensibilidad y Especificidad , Relación Señal-RuidoRESUMEN
Some microstructure parameters, such as permeability, remain elusive because mathematical models that express their relationship to the MR signal accurately are intractable. Here, we propose to use computational models learned from simulations to estimate these parameters. We demonstrate the approach in an example which estimates water residence time in brain white matter. The residence time τi of water inside axons is a potentially important biomarker for white matter pathologies of the human central nervous system, as myelin damage is hypothesised to affect axonal permeability, and thus τi. We construct a computational model using Monte Carlo simulations and machine learning (specifically here a random forest regressor) in order to learn a mapping between features derived from diffusion weighted MR signals and ground truth microstructure parameters, including τi. We test our numerical model using simulated and in vivo human brain data. Simulation results show that estimated parameters have strong correlations with the ground truth parameters (R2={0.88,0.95,0.82,0.99}) for volume fraction, residence time, axon radius and diffusivity respectively), and provide a marked improvement over the most widely used Kärger model (R2={0.75,0.60,0.11,0.99}). The trained model also estimates sensible microstructure parameters from in vivo human brain data acquired from healthy controls, matching values found in literature, and provides better reproducibility than the Kärger model on both the voxel and ROI level. Finally, we acquire data from two Multiple Sclerosis (MS) patients and compare to the values in healthy subjects. We find that in the splenium of corpus callosum (CC-S) the estimate of the residence time is 0.57±0.05s for the healthy subjects, while in the MS patient with a lesion in CC-S it is 0.33±0.12s in the normal appearing white matter (NAWM) and 0.19±0.11s in the lesion. In the corticospinal tracts (CST) the estimate of the residence time is 0.52±0.09s for the healthy subjects, while in the MS patient with a lesion in CST it is 0.56±0.05s in the NAWM and 0.13±0.09s in the lesion. These results agree with our expectations that the residence time in lesions would be lower than in NAWM because the loss of myelin should increase permeability. Overall, we find parameter estimates in the two MS patients consistent with expectations from the pathology of MS lesions demonstrating the clinical potential of this new technique.
Asunto(s)
Encéfalo/diagnóstico por imagen , Simulación por Computador , Aprendizaje Automático , Modelos Teóricos , Sustancia Blanca/diagnóstico por imagen , Adulto , Encéfalo/patología , Imagen Eco-Planar , Femenino , Humanos , Interpretación de Imagen Asistida por Computador/métodos , Masculino , Persona de Mediana Edad , Método de Montecarlo , Esclerosis Múltiple/diagnóstico por imagen , Esclerosis Múltiple/patología , Permeabilidad , Sustancia Blanca/patología , Adulto JovenRESUMEN
Data in recordings obtained from ambulatory patients using wearable sensors are often corrupted by motion artefact and are, in general, noisier than the data obtained from the nonmobile patients. Identifying and ignoring erroneous measurements from these data is very important, if wearable sensors are to be incorporated into clinical practice. In this paper, we propose a novel Signal Quality Index, intended to assess whether reliable heart rates can be obtained from a single channel of ECG collected from ambulatory patients, using wearable sensors. The proposed system is based on wavelet entropy measurements of the heart rate variability signal. The system was trained and tested on expert-labeled data from a particular wearable sensor and was also tested on labeled data from a different sensor. The sensitivities and specificities achieved were 94% and 98%, respectively, on data from the same sensor as the training set, and 91% and 97%, respectively, on data from a different sensor, indicating the potential of the system to generalize across different sensors. Because the system relies on a single channel of ECG, it has the potential for inclusion in applications using wearable sensors and in the most basic clinical environments.
Asunto(s)
Electrocardiografía Ambulatoria/métodos , Frecuencia Cardíaca/fisiología , Procesamiento de Señales Asistido por Computador , Algoritmos , Entropía , Humanos , Aprendizaje AutomáticoRESUMEN
PURPOSE: To introduce a novel diffusion pulse sequence, namely double oscillating diffusion encoding (DODE), and to investigate whether it adds sensitivity to microscopic diffusion anisotropy (µA) compared to the well-established double diffusion encoding (DDE) methodology. METHODS: We simulate measurements from DODE and DDE sequences for different types of microstructures exhibiting restricted diffusion. First, we compare the effect of varying pulse sequence parameters on the DODE and DDE signal. Then, we analyse the sensitivity of the two sequences to the microstructural parameters (pore diameter and length) which determine µA. Finally, we investigate specificity of measurements to particular substrate configurations. RESULTS: Simulations show that DODE sequences exhibit similar signal dependence on the relative angle between the two gradients as DDE sequences, however, the effect of varying the mixing time is less pronounced. The sensitivity analysis shows that in substrates with elongated pores and various orientations, DODE sequences increase the sensitivity to pore diameter, while DDE sequences are more sensitive to pore length. Moreover, DDE and DODE sequence parameters can be tailored to enhance/suppress the signal from a particular range of substrates. CONCLUSIONS: A combination of DODE and DDE sequences maximize sensitivity to µA, compared to using just the DDE method. Magn Reson Med 78:550-564, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
Asunto(s)
Imagen de Difusión por Resonancia Magnética/métodos , Procesamiento de Imagen Asistido por Computador/métodos , Procesamiento de Señales Asistido por Computador , Anisotropía , Técnicas Citológicas , Microscopía , Modelos TeóricosRESUMEN
Non-invasive estimation of cell size and shape is a key challenge in diffusion MRI. This article presents a model-based approach that provides independent estimates of pore size and eccentricity from diffusion MRI data. The technique uses a geometric model of finite cylinders with gamma-distributed radii to represent pores of various sizes and elongations. We consider both macroscopically isotropic substrates and substrates of semi-coherently oriented anisotropic pores and we use Monte Carlo simulations to generate synthetic data. We compare the sensitivity of single and double diffusion encoding (SDE and DDE) sequences to the size distribution and eccentricity, and further analyse different protocols of DDE sequences with parallel and/or perpendicular pairs of gradients. We show that explicitly accounting for size distribution is necessary for accurate microstructural parameter estimates, and a model that assumes a single size yields biased eccentricity values. We also find that SDE sequences support estimates, although DDE sequences with mixed parallel and perpendicular gradients enhance accuracy. In the case of macroscopically anisotropic substrates, this model-based approach can be extended to a rotationally invariant framework to provide features of pore shape (specifically eccentricity) in the presence of size distribution and orientation dispersion.
Asunto(s)
Simulación por Computador , Imagen de Difusión por Resonancia Magnética/métodos , Microscopía , Modelos Biológicos , AnisotropíaRESUMEN
PURPOSE: To identify optimal pulsed gradient spin-echo (PGSE) and oscillating gradient spin-echo (OGSE) sequence settings for maximizing sensitivity to axon diameter in idealized and practical conditions. METHODS: Simulations on a simple two-compartment white matter model (with nonpermeable cylinders) are used to investigate a wide space of clinically plausible PGSE and OGSE sequence parameters with trapezoidal diffusion gradient waveforms. Signal sensitivity is measured as a derivative of the signal with respect to axon diameter. Models of parallel and dispersed fibers are investigated separately to represent idealized and practical conditions. RESULTS: Simulations show that, for the simple case of gradients perfectly perpendicular to straight parallel fibers, PGSE always gives maximum sensitivity. However, in real-world scenarios where fibers have unknown and dispersed orientation, low-frequency OGSE provides higher sensitivity. Maximum sensitivity results show that on current clinical scanners (Gmax = 60 mT/m, signal to noise ratio (SNR) = 20) axon diameters below 6 µm are indistinguishable from zero. Scanners with stronger gradient systems such as the Massachusetts General Hospital (MGH) Connectom scanner (Gmax = 300 mT/m) can extend this sensitivity limit down to 2-3 µm, probing a much greater proportion of the underlying axon diameter distribution. CONCLUSION: Low-frequency OGSE provides additional sensitivity to PGSE in practical situations. OGSE is particularly advantageous for systems with high performance gradients.
Asunto(s)
Axones/ultraestructura , Imagen de Difusión por Resonancia Magnética/métodos , Sustancia Blanca/anatomía & histología , Simulación por Computador , Humanos , Procesamiento de Imagen Asistido por Computador/métodos , Fantasmas de Imagen , Sensibilidad y EspecificidadRESUMEN
Stejskal and Tanner's ingenious pulsed field gradient design from 1965 has made diffusion NMR and MRI the mainstay of most studies seeking to resolve microstructural information in porous systems in general and biological systems in particular. Methods extending beyond Stejskal and Tanner's design, such as double diffusion encoding (DDE) NMR and MRI, may provide novel quantifiable metrics that are less easily inferred from conventional diffusion acquisitions. Despite the growing interest on the topic, the terminology for the pulse sequences, their parameters, and the metrics that can be derived from them remains inconsistent and disparate among groups active in DDE. Here, we present a consensus of those groups on terminology for DDE sequences and associated concepts. Furthermore, the regimes in which DDE metrics appear to provide microstructural information that cannot be achieved using more conventional counterparts (in a model-free fashion) are elucidated. We highlight in particular DDE's potential for determining microscopic diffusion anisotropy and microscopic fractional anisotropy, which offer metrics of microscopic features independent of orientation dispersion and thus provide information complementary to the standard, macroscopic, fractional anisotropy conventionally obtained by diffusion MR. Finally, we discuss future vistas and perspectives for DDE.
Asunto(s)
Imagen por Resonancia Magnética/clasificación , Imagen por Resonancia Magnética/normas , Espectroscopía de Resonancia Magnética/clasificación , Espectroscopía de Resonancia Magnética/normas , Procesamiento de Señales Asistido por Computador , Terminología como Asunto , Guías como AsuntoRESUMEN
In this paper we demonstrate a simulation framework that enables the direct and quantitative comparison of post-processing methods for diffusion weighted magnetic resonance (DW-MR) images. DW-MR datasets are employed in a range of techniques that enable estimates of local microstructure and global connectivity in the brain. These techniques require full alignment of images across the dataset, but this is rarely the case. Artefacts such as eddy-current (EC) distortion and motion lead to misalignment between images, which compromise the quality of the microstructural measures obtained from them. Numerous methods and software packages exist to correct these artefacts, some of which have become de-facto standards, but none have been subject to rigorous validation. In the literature, improved alignment is assessed using either qualitative visual measures or quantitative surrogate metrics. Here we introduce a simulation framework that allows for the direct, quantitative assessment of techniques, enabling objective comparisons of existing and future methods. DW-MR datasets are generated using a process that is based on the physics of MRI acquisition, which allows for the salient features of the images and their artefacts to be reproduced. We apply this framework in three ways. Firstly we assess the most commonly used method for artefact correction, FSL's eddy_correct, and compare it to a recently proposed alternative, eddy. We demonstrate quantitatively that using eddy_correct leads to significant errors in the corrected data, whilst eddy is able to provide much improved correction. Secondly we investigate the datasets required to achieve good correction with eddy, by looking at the minimum number of directions required and comparing the recommended full-sphere acquisitions to equivalent half-sphere protocols. Finally, we investigate the impact of correction quality by examining the fits from microstructure models to real and simulated data.