Empirical comparison of diffusion kurtosis imaging and diffusion basis spectrum imaging using the same acquisition in healthy young adults

As diffusion tensor imaging gains widespread use, many researchers have been motivated to go beyond the tensor model and fit more complex diffusion models, to gain a more complete description of white matter microstructure and associated pathology. Two such models are diffusion kurtosis imaging (DKI) and diffusion basis spectrum imaging (DBSI). It is not clear which DKI parameters are most closely related to DBSI parameters, so in the interest of enabling comparisons between DKI and DBSI studies, we conducted an empirical survey of the interrelation of these models in 12 healthy volunteers using the same diffusion acquisition. We found that mean kurtosis is positively associated with the DBSI fiber ratio and negatively associated with the hindered ratio. This was primarily driven by the radial component of kurtosis. The axial component of kurtosis was strongly and specifically correlated with the restricted ratio. The joint spatial distributions of DBSI and DKI parameters are tissue-dependent and stable across healthy individuals. Our contribution is a better understanding of the biological interpretability of the parameters generated by the two models in healthy individuals

Correcting for Motion between Acquisitions in Diffusion MR Imaging

The diffusion tensor (DT) and other diffusion models assume that each voxel corresponds to the same anatomical location in all the measurements. Movements and distortions violate this assumption and typically the images are realigned before model fitting. We propose a set of model-based methods to improve motion correction and avoid the errors that the traditional method introduces. The new methods are based on a three-step procedure to register DWI datasets, and use different reference images for DWIs with different gradient directions for registration, so the registrations take into account the contrast differences of measurements. Performance of the model-based registration techniques depends critically on outlier rejection. We develop new methods for fitting the diffusion tensor to diffusion MRI measurements in the presence of outliers by drawing on the RANSAC algorithm from computer vision. We compareone popularly used outlier rejection method RESTORE in the diffusion MRI literature with our new method. Then, we combine outlier rejection methods with model-based registration schemes, and compare the performance of motion correction with other methods. After aligning the dataset, we also update diffusion gradients for the registered datasets from both traditional and our methods, according to the transformations used in registrations. We develop and discuss a variety of registration evaluation methods using both synthetic and human-brain diffusion MRI datasets. Experiments demonstrate both quantitative and qualitative improvements using our new model-based methods

Deep learning‐based method for reducing residual motion effects in diffusion parameter estimation

PURPOSE: Conventional motion-correction techniques for diffusion MRI can introduce motion-level-dependent bias in derived metrics. To address this challenge, a deep learning-based technique was developed to minimize such residual motion effects. METHODS: The data-rejection approach was adopted in which motion-corrupted data are discarded before model-fitting. A deep learning-based parameter estimation algorithm, using a hierarchical convolutional neural network (H-CNN), was combined with motion assessment and corrupted volume rejection. The method was designed to overcome the limitations of existing methods of this kind that produce parameter estimations whose quality depends strongly on a proportion of the data discarded. Evaluation experiments were conducted for the estimation of diffusion kurtosis and diffusion-tensor-derived measures at both the individual and group levels. The performance was compared with the robust approach of iteratively reweighted linear least squares (IRLLS) after motion correction with and without outlier replacement. RESULTS: Compared with IRLLS, the H-CNN-based technique is minimally sensitive to motion effects. It was tested at severe motion levels when 70% to 90% of the data are rejected and when random motion is present. The technique had a stable performance independent of the numbers and schemes of data rejection. A further test on a data set from children with attention-deficit hyperactivity disorder shows the technique can potentially ameliorate spurious group-level difference caused by head motion. CONCLUSION: This method shows great potential for reducing residual motion effects in motion-corrupted diffusion-weighted-imaging data, bringing benefits that include reduced bias in derived metrics in individual scans and reduced motion-level-dependent bias in population studies employing diffusion MRI

Automated detection of brain abnormalities in neonatal hypoxia ischemic injury from MR images.

We compared the efficacy of three automated brain injury detection methods, namely symmetry-integrated region growing (SIRG), hierarchical region splitting (HRS) and modified watershed segmentation (MWS) in human and animal magnetic resonance imaging (MRI) datasets for the detection of hypoxic ischemic injuries (HIIs). Diffusion weighted imaging (DWI, 1.5T) data from neonatal arterial ischemic stroke (AIS) patients, as well as T2-weighted imaging (T2WI, 11.7T, 4.7T) at seven different time-points (1, 4, 7, 10, 17, 24 and 31 days post HII) in rat-pup model of hypoxic ischemic injury were used to assess the temporal efficacy of our computational approaches. Sensitivity, specificity, and similarity were used as performance metrics based on manual ('gold standard') injury detection to quantify comparisons. When compared to the manual gold standard, automated injury location results from SIRG performed the best in 62% of the data, while 29% for HRS and 9% for MWS. Injury severity detection revealed that SIRG performed the best in 67% cases while 33% for HRS. Prior information is required by HRS and MWS, but not by SIRG. However, SIRG is sensitive to parameter-tuning, while HRS and MWS are not. Among these methods, SIRG performs the best in detecting lesion volumes; HRS is the most robust, while MWS lags behind in both respects

PVR: Patch-to-Volume Reconstruction for Large Area Motion Correction of Fetal MRI

In this paper we present a novel method for the correction of motion artifacts that are present in fetal Magnetic Resonance Imaging (MRI) scans of the whole uterus. Contrary to current slice-to-volume registration (SVR) methods, requiring an inflexible anatomical enclosure of a single investigated organ, the proposed patch-to-volume reconstruction (PVR) approach is able to reconstruct a large field of view of non-rigidly deforming structures. It relaxes rigid motion assumptions by introducing a specific amount of redundant information that is exploited with parallelized patch-wise optimization, super-resolution, and automatic outlier rejection. We further describe and provide an efficient parallel implementation of PVR allowing its execution within reasonable time on commercially available graphics processing units (GPU), enabling its use in the clinical practice. We evaluate PVR's computational overhead compared to standard methods and observe improved reconstruction accuracy in presence of affine motion artifacts of approximately 30% compared to conventional SVR in synthetic experiments. Furthermore, we have evaluated our method qualitatively and quantitatively on real fetal MRI data subject to maternal breathing and sudden fetal movements. We evaluate peak-signal-to-noise ratio (PSNR), structural similarity index (SSIM), and cross correlation (CC) with respect to the originally acquired data and provide a method for visual inspection of reconstruction uncertainty. With these experiments we demonstrate successful application of PVR motion compensation to the whole uterus, the human fetus, and the human placenta.Comment: 10 pages, 13 figures, submitted to IEEE Transactions on Medical Imaging. v2: wadded funders acknowledgements to preprin

Fast and accurate Slicewise OutLIer Detection (SOLID) with informed model estimation for diffusion MRI data

The accurate characterization of the diffusion process in tissue using diffusion MRI is greatly challenged by the presence of artefacts. Subject motion causes not only spatial misalignments between diffusion weighted images, but often also slicewise signal intensity errors. Voxelwise robust model estimation is commonly used to exclude intensity errors as outliers. Slicewise outliers, however, become distributed over multiple adjacent slices after image registration and transformation. This challenges outlier detection with voxelwise procedures due to partial volume effects. Detecting the outlier slices before any transformations are applied to diffusion weighted images is therefore required. In this work, we present i) an automated tool coined SOLID for slicewise outlier detection prior to geometrical image transformation, and ii) a framework to naturally interpret data uncertainty information from SOLID and include it as such in model estimators. SOLID uses a straightforward intensity metric, is independent of the choice of the diffusion MRI model, and can handle datasets with a few or irregularly distributed gradient directions. The SOLID-informed estimation framework prevents the need to completely reject diffusion weighted images or individual voxel measurements by downweighting measurements with their degree of uncertainty, thereby supporting convergence and well-conditioning of iterative estimation algorithms. In comprehensive simulation experiments, SOLID detects outliers with a high sensitivity and specificity, and can achieve higher or at least similar sensitivity and specificity compared to other tools that are based on more complex and time-consuming procedures for the scenarios investigated. SOLID was further validated on data from 54 neonatal subjects which were visually inspected for outlier slices with the interactive tool developed as part of this study, showing its potential to quickly highlight problematic volumes and slices in large population studies. The informed model estimation framework was evaluated both in simulations and in vivo human data.Peer reviewe

From Diffusion to Tracts

Diffusion of water molecules within the brain tissue can be used to modulate the nuclear magnetic resonance signal that is used to form magnetic resonance images (MRI). As the signal itself can be noisy and its meaning challenging to interpret, mathematical models are generally fitted to these measurements to obtain the more accurate characterization of the brain microstructure. This, of course, requires that the mathematical model itself is sound in respect to the measurement setup. This dissertation focuses on the extensively used tensor models as they have been shown to unravel details of the physical diffusion phenomena along with various applications in the basic neuroscience, the clinical research, and even in the neurosurgery. One of the greatest challenges in the diffusion weighted MRI measurements is subject motion during the image acquisition as that can cause a complete loss of the measurement which is especially highlighted in ill or uncooperative patients studies. Due to the used acquisition technique, this loss extends to multiple measurements simultaneously resulting in an enormous gap in the sampling. Such gaps can be problematic for any model fitting, even for the currently available robust means developed to exclude outlier measurements from affecting the estimate. Hence in this dissertation, a tool coined as SOLID was developed to detect these outliers and to robustly process them during the tensor based model estimation. SOLID was implemented as a part of the widely used ExploreDTI toolbox to allow the rapid international distribution of the tool. Unfortunately, any reduction in the measurement sampling will lead to increasing error propagation during the model estimation. Mathematically this is detailed in terms of a condition number for the matrix inversion in the linear least squares fitting. Previously, the condition number has been used to optimize the diffusion weighted MRI acquisition gradient scheme but in this dissertation it was renovated into a novel quality control tool. The condition number of the matrix inversion that provides the model estimate can be calculated after the outliers are excluded to assess spatially and directionally varying error propagation to obviate any bias in subject or population studies. To motivate the importance of the robust methods and diffusion weighted MRI at large, neurocognitive studies with neonates’ visual abilities and bilinguals’ acquisition age of the second language were conducted as a part of this thesis. The findings in these studies indicated that premature birth affects the white matter structures across the brain whereas the age of acquisition of the second language affects only the speech related brain structures.Aivojen rakenteessa tapahtuvien muutosten mittaaminen on avainasemassa tutkittaessa esimerkiksi keskosena syntyneen lapsen kehitystä tai uusien taitojen, kuten kielten, oppimista. Ihmisaivojen tutkiminen on aiemmin rajoittunut aivojen toiminnan arviointiin aivosähkökäyrän ja neurokognitiivisten testien avulla. Viime vuosikymmenten kehitys magneettikuvaustekniikassa on tuonut mahdollisuuden tutkia kajoamattomasti myös aivojen rakennetta ja jopa seurata sen muutosta lapsen kasvaessa tai ihmisen oppiessa uusia taitoja. Yksi lupaavimmista aivojen tutkimusmenetelmistä on diffuusiopainotettu magneettikuvaus, jolle on löytynyt lukuisia käyttökohteita niin neurotieteessä, lääketieteellisissä tutkimuksissa kuin neurokirurgiassakin. Menetelmä perustuu vesimolekyylien lämpöliikkeen mittaamiseen aivoissa. Molekyylien liike on vapaata muun muassa valkean aineen rakenteiden myötäisesti, mutta lähes mahdotonta kohtisuoraan niiden lävitse. Jäljittämällä nämä reitit voidaan muodostaa tarkka malli aivojen rakenteesta. Mallin pohjalta on mahdollista laskea kuvaavia arvoja, jotka auttavat esimerkiksi määrittämään aivovaurion astetta. Diffuusiopainotetun magneettikuvauksen suurin haaste on menetelmän monimutkaisuus sekä mittauksen että analyysin osalta. Vain hyvin yksinkertaisissa tapauksissa asiantuntija voi arvioida suoraan diffuusiopainotetusta magneettikuvasta poikkeamia aivoissa. Yleensä käytetään matemaattisia menetelmiä kuvan tarkempaan analysointiin. Tällöin keskeistä on inversio-ongelman ratkaisu, missä potilaasta tehdyt mittaukset sovitetaan aivoja kuvaavaan matemaattiseen malliin. Sopivan mallin valinnalla on siis suuri vaikutus lopputuloksen hyödyllisyyteen. Diffuusiopainotettu magneettikuvaus on myös häiriöherkkä ja mittaukset sisältävät luonnostaan paljon kohinaa, jonka vaikutusta vähennetään tekemällä toistomittauksia. Toistomittaukset pidentävät kuvausaikaa, joka puolestaan voi olla haasteellinen potilaalle, koska potilaan pitää olla liikkumatta koko kuvauksen ajan. Potilaan pään pienikin liike voi johtaa huomattaviin mittavirheisiin, koska menetelmällä mitataan vesimolekyylien liikettä, jonka suuruus on vain kymmenien mikrometrien luokkaa. Tässä fysiikan väitöskirjassa keskityttiin diffuusiopainotetun magneettikuvauksen mallintamismenetelmien kehitystyöhön ja niiden käyttöönottoon Helsingin yliopistollisessa sairaalassa. Kehitimme kansainvälistä huomiota herättäneen SOLID-työkalun, jolla voidaan havaita sekä korjata potilaan liikkeestä aiheutuvia virheitä mittaustuloksissa. Tämän lisäksi esitimme laadunvalvonta menetelmän, jolla voidaan arvioida esimerkiksi potilaiden välisten mallinnustulosten vertailukelpoisuutta. Kehitettyjä menetelmiä testattiin ja sovellettiin kahdessa tutkimuksessa: Osoitimme, että vastasyntyneen lapsen kyky seurata katseellaan liikkuvaa kohdetta liittyy laaja-alaisiin muutoksiin aivojen valkean aineen rakenteessa. Lisäksi näytimme, että toisen kielen oppimisajankohta vaikuttaa aivojen puheentuottoon liittyvien aivorakenteiden muodostumiseen