This paper presents an efficient algorithm for large deformation diffeomorphic metric mapping (LDDMM) with geodesic shooting for image registration. We introduce a novel finite dimensional Fourier representation of diffeomorphic deformations based on the key fact that the high frequency components of a diffeomorphism remain stationary throughout the integration process when computing the deformation associated with smooth velocity fields. We show that manipulating high dimensional diffeomorphisms can be carried out entirely in the bandlimited space by integrating the nonstationary low frequency components of the displacement field. This insight substantially reduces the computational cost of the registration problem. Experimental results show that our method is significantly faster than the state-of-the-art diffeomorphic image registration methods while producing equally accurate alignment. We demonstrate our algorithm in two different applications of image registration: neuroimaging and in-utero imaging.
The neural correlates of spaceflight-induced sensorimotor impairments are unknown. Head down-tilt bed rest (HDBR) serves as a microgravity analog because it mimics the headward fluid shift and axial body unloading of spaceflight. We investigated focal brain white matter (WM) changes and fluid shifts during 70 days of 6° HDBR in 16 subjects who were assessed pre (2x), during (3x), and post-HDBR (2x). Changes over time were compared to those in control subjects (n = 12) assessed four times over 90 days. Diffusion MRI was used to assess WM microstructure and fluid shifts. Free-Water Imaging was used to quantify distribution of intracranial extracellular free water (FW). Additionally, we tested whether WM and FW changes correlated with changes in functional mobility and balance measures. HDBR resulted in FW increases in fronto-temporal regions and decreases in posterior-parietal regions that largely recovered by two weeks post-HDBR. WM microstructure was unaffected by HDBR. FW decreases in the post-central gyrus and precuneus correlated negatively with balance changes. We previously reported that gray matter increases in these regions were associated with less HDBR-induced balance impairment, suggesting adaptive structural neuroplasticity. Future studies are warranted to determine causality and underlying mechanisms.
We propose a non-parametric approach for characterizing heterogeneous diseases in large-scale studies. We target diseases where multiple types of pathology present simultaneously in each subject and a more severe disease manifests as a higher level of tissue destruction. For each subject, we model theof local image descriptors as samples generated by an unknown subject-specific probability density. Instead of approximating the probability density via a parametric family, we propose to side step the parametric inference by directly estimating the divergence between subject densities. Our method maps the collection of local image descriptors to a signaturethat is used to predict a clinical measurement. We are able to interpret the prediction of the clinical variable in the population and individual levels by carefully studying the divergences. We illustrate an application this method on simulated data as well as on a large-scale lung CT study of Chronic Obstructive Pulmonary Disease (COPD). Our approach outperforms classical methods on both simulated and COPD data and demonstrates the state-of-the-art prediction on an important physiologic measure of airflow (the forced respiratory volume in one second, FEV1).
Diffusion MRI tractography is increasingly used in pre-operative neurosurgical planning to visualize critical fiber tracts. However, a major challenge for conventional tractography, especially in patients with brain tumors, is tracing fiber tracts that are affected by vasogenic edema, which increases water content in the tissue and lowers diffusion anisotropy. One strategy for improving fiber tracking is to use a tractography method that is more sensitive than the traditional single-tensor streamline tractography. We performed experiments to assess the performance of two-tensor unscented Kalman filter (UKF) tractography in edema. UKF tractography fits a diffusion model to the data during fiber tracking, taking advantage of prior information from the previous step along the fiber. We studied UKF performance in a synthetic diffusion MRI digital phantom with simulated edema and in retrospective data from two neurosurgical patients with edema affecting the arcuate fasciculus and corticospinal tracts. We compared the performance of several tractography methods including traditional streamline, UKF single-tensor, and UKF two-tensor. To provide practical guidance on how the UKF method could be employed, we evaluated the impact of using various seed regions both inside and outside the edematous regions, as well as the impact of parameter settings on the tractography sensitivity. We quantified the sensitivity of different methods by measuring the percentage of the patient-specific fMRI activation that was reached by the tractography. We expected that diffusion anisotropy threshold parameters, as well as the inclusion of a free water model, would significantly influence the reconstruction of edematous WM fiber tracts, because edema increases water content in the tissue and lowers anisotropy. Contrary to our initial expectations, varying the fractional anisotropy threshold and including a free water model did not affect the UKF two-tensor tractography output appreciably in these two patient datasets. The most effective parameter for increasing tracking sensitivity was the generalized anisotropy (GA) threshold, which increased the length of tracked fibers when reduced to 0.075. In addition, the most effective seeding strategy was seeding in the whole brain or in a large region outside of the edema. Overall, the main contribution of this study is to provide insight into how UKF tractography can work, using a two-tensor model, to begin to address the challenge of fiber tract reconstruction in edematous regions near brain tumors.
Adrian V Dalca, Katherine L Bouman, William T. Freeman, Natalia S Rost, Mert R Sabuncu, and Polina Golland. 6/2017. “Population Based Image Imputation.” Inf Process Med Imaging, 10265, Pp. 659-671.Abstract
We present an algorithm for creating high resolution anatomically plausible images consistent with acquired clinical brain MRI scans with large inter-slice spacing. Although large databases of clinical images contain a wealth of information, medical acquisition constraints result in sparse scans that miss much of the anatomy. These characteristics often render computational analysis impractical as standard processing algorithms tend to fail when applied to such images. Highly specialized or application-specific algorithms that explicitly handle sparse slice spacing do not generalize well across problem domains. In contrast, our goal is to enable application of existing algorithms that were originally developed for high resolution research scans to significantly undersampled scans. We introduce a model that captures fine-scale anatomical similarity across subjects in clinical image collections and use it to fill in the missing data in scans with large slice spacing. Our experimental results demonstrate that the proposed method outperforms current upsampling methods and promises to facilitate subsequent analysis not previously possible with scans of this quality.
Diffusion imaging is critical for detecting acute brain injury. However, normal apparent diffusion coefficient (ADC) maps change rapidly in early childhood, making abnormality detection difficult. In this article, we explored clinical PACS and electronic healthcare records (EHR) to create age-specific ADC atlases for clinical radiology reference. Using the EHR and three rounds of multiexpert reviews, we found ADC maps from 201 children 0-6 years of age scanned between 2006 and 2013 who had brain MRIs with no reported abnormalities and normal clinical evaluations 2+ years later. These images were grouped in 10 age bins, densely sampling the first 1 year of life (5 bins, including neonates and 4 quarters) and representing the 1-6 year age range (an age bin per year). Unbiased group-wise registration was used to construct ADC atlases for 10 age bins. We used the atlases to quantify (a) cross-sectional normative ADC variations; (b) spatiotemporal heterogeneous ADC changes; and (c) spatiotemporal heterogeneous volumetric changes. The quantified age-specific whole-brain and region-wise ADC values were compared to those from age-matched individual subjects in our study and in multiple existing independent studies. The significance of this study is that we have shown that clinically acquired images can be used to construct normative age-specific atlases. These first of their kind age-specific normative ADC atlases quantitatively characterize changes of myelination-related water diffusion in the first 6 years of life. The quantified voxel-wise spatiotemporal ADC variations provide standard references to assist radiologists toward more objective interpretation of abnormalities in clinical images. Our atlases are available at https://www.nitrc.org/projects/mgh_adcatlases.
We perform a review of the literature in the field of white matter tractography for neurosurgical planning, focusing on those works where tractography was correlated with clinical information such as patient outcome, clinical functional testing, or electro-cortical stimulation. We organize the review by anatomical location in the brain and by surgical procedure, including both supratentorial and infratentorial pathologies, and excluding spinal cord applications. Where possible, we discuss implications of tractography for clinical care, as well as clinically relevant technical considerations regarding the tractography methods. We find that tractography is a valuable tool in variable situations in modern neurosurgery. Our survey of recent reports demonstrates multiple potentially successful applications of white matter tractography in neurosurgery, with progress towards overcoming clinical challenges of standardization and interpretation.
The glymphatic pathway is a system which facilitates continuous cerebrospinal fluid (CSF) and interstitial fluid (ISF) exchange and plays a key role in removing waste products from the rodent brain. Dysfunction of the glymphatic pathway may be implicated in the pathophysiology of Alzheimer's disease. Intriguingly, the glymphatic system is most active during deep wave sleep general anesthesia. By using paramagnetic tracers administered into CSF of rodents, we previously showed the utility of MRI in characterizing a macroscopic whole brain view of glymphatic transport but we have yet to define and visualize the specific flow patterns. Here we have applied an alternative mathematical analysis approach to a dynamic time series of MRI images acquired every 4min over ∼3h in anesthetized rats, following administration of a small molecular weight paramagnetic tracer into the CSF reservoir of the cisterna magna. We use Optimal Mass Transport (OMT) to model the glymphatic flow vector field, and then analyze the flow to find the network of CSF-ISF flow channels. We use 3D visualization computational tools to visualize the OMT defined network of CSF-ISF flow channels in relation to anatomical and vascular key landmarks from the live rodent brain. The resulting OMT model of the glymphatic transport network agrees largely with the current understanding of the glymphatic transport patterns defined by dynamic contrast-enhanced MRI revealing key CSF transport pathways along the ventral surface of the brain with a trajectory towards the pineal gland, cerebellum, hypothalamus and olfactory bulb. In addition, the OMT analysis also revealed some interesting previously unnoticed behaviors regarding CSF transport involving parenchymal streamlines moving from ventral reservoirs towards the surface of the brain, olfactory bulb and large central veins.
In order to bridge microscopic molecular motion with macroscopic diffusion MR signal in complex structures, we propose a general stochastic model for molecular motion in a magnetic field. The Fokker-Planck equation of this model governs the probability density function describing the diffusion-magnetization propagator. From the propagator we derive a generalized version of the Bloch-Torrey equation and the relation to the random phase approach. This derivation does not require assumptions such as a spatially constant diffusion coefficient, or ad hoc selection of a propagator. In particular, the boundary conditions that implicitly incorporate the microstructure into the diffusion MR signal can now be included explicitly through a spatially varying diffusion coefficient. While our generalization is reduced to the conventional Bloch-Torrey equation for piecewise constant diffusion coefficients, it also predicts scenarios in which an additional term to the equation is required to fully describe the MR signal.
The hypothesis that brain pathways form 2D sheet-like structures layered in 3D as "pages of a book" has been a topic of debate in the recent literature. This hypothesis was mainly supported by a qualitative evaluation of "path neighborhoods" reconstructed with diffusion MRI (dMRI) tractography. Notwithstanding the potentially important implications of the sheet structure hypothesis for our understanding of brain structure and development, it is still considered controversial by many for lack of quantitative analysis. A means to quantify sheet structure is therefore necessary to reliably investigate its occurrence in the brain. Previous work has proposed the Lie bracket as a quantitative indicator of sheet structure, which could be computed by reconstructing path neighborhoods from the peak orientations of dMRI orientation density functions. Robust estimation of the Lie bracket, however, is challenging due to high noise levels and missing peak orientations. We propose a novel method to estimate the Lie bracket that does not involve the reconstruction of path neighborhoods with tractography. This method requires the computation of derivatives of the fiber peak orientations, for which we adopt an approach called normalized convolution. With simulations and experimental data we show that the new approach is more robust with respect to missing peaks and noise. We also demonstrate that the method is able to quantify to what extent sheet structure is supported for dMRI data of different species, acquired with different scanners, diffusion weightings, dMRI sampling schemes, and spatial resolutions. The proposed method can also be used with directional data derived from other techniques than dMRI, which will facilitate further validation of the existence of sheet structure.
Implant placement has been widely used in various kinds of surgery. However, accurate intraoperative drilling performance is essential to avoid injury to adjacent structures. Although some commercially-available surgical navigation systems have been approved for clinical applications, these systems are expensive and the source code is not available to researchers. 3D Slicer is a free, open source software platform for the research community of computer-aided surgery. In this study, a loadable module based on Slicer has been developed and validated to support surgical navigation. This research module allows reliable calibration of the surgical drill, point-based registration and surface matching registration, so that the position and orientation of the surgical drill can be tracked and displayed on the computer screen in real time, aiming at reducing risks. In accuracy verification experiments, the mean target registration error (TRE) for point-based and surface-based registration were 0.31±0.06mm and 1.01±0.06mm respectively, which should meet clinical requirements. Both phantom and cadaver experiments demonstrated the feasibility of our surgical navigation software module.
We consider diffusion within pores with general shapes in the presence of spatially linear magnetic field profiles. The evolution of local magnetization of the spin bearing particles can be described by the Bloch-Torrey equation. We study the diffusive process in the eigenbasis of the non-Hermitian Bloch-Torrey operator. It is possible to find expressions for some special temporal gradient waveforms employed to sensitize the nuclear magnetic resonance (NMR) signal to diffusion. For more general gradient waveforms, we derive an efficient numerical solution by introducing a novel matrix formalism. Compared to previous methods, this new approach requires a fewer number of eigenfunctions to achieve the same accuracy. This shows that these basis functions are better suited to the problem studied. The new framework could provide new important insights into the fundamentals of diffusion sensitization, which could further the development of the field of NMR.
The basal ganglia is part of a complex system of neuronal circuits that play a key role in the integration and execution of motor, cognitive and emotional function in the human brain. Parkinson's disease is a progressive neurological disorder of the motor circuit characterized by tremor, rigidity, and slowness of movement. Deep brain stimulation (DBS) of the subthalamic nucleus and the globus pallidus pars interna provides an efficient treatment to reduce symptoms and levodopa-induced side effects in Parkinson's disease patients. While the underlying mechanism of action of DBS is still unknown, the potential modulation of white matter tracts connecting the surgical targets has become an active area of research. With the introduction of advanced diffusion MRI acquisition sequences and sophisticated post-processing techniques, the architecture of the human brain white matter can be explored in vivo. The goal of this study is to investigate the white matter connectivity between the subthalamic nucleus and the globus pallidus. Two multi-fiber tractography methods were used to reconstruct pallido-subthalamic, subthalamo-pallidal and pyramidal fibers in five healthy subjects datasets of the Human Connectome Project. The anatomical accuracy of the tracts was assessed by four judges with expertise in neuroanatomy, functional neurosurgery, and diffusion MRI. The variability among subjects was evaluated based on the fractional anisotropy and mean diffusivity of the tracts. Both multi-fiber approaches enabled the detection of complex fiber architecture in the basal ganglia. The qualitative evaluation by experts showed that the identified tracts were in agreement with the expected anatomy. Tract-derived measurements demonstrated relatively low variability among subjects. False-negative tracts demonstrated the current limitations of both methods for clinical decision-making. Multi-fiber tractography methods combined with state-of-the-art diffusion MRI data have the potential to help identify white matter tracts connecting DBS targets in functional neurosurgery intervention.
INTRODUCTION: Huntington's disease (HD) is a genetic neurodegenerative disorder that primarily affects striatal neurons. Striatal volume loss is present years before clinical diagnosis; however, white matter degradation may also occur prior to diagnosis. Diffusion-weighted imaging (DWI) can measure microstructural changes associated with degeneration that precede macrostructural changes. DWI derived measures enhance understanding of degeneration in prodromal HD (pre-HD). METHODS: As part of the PREDICT-HD study, N = 191 pre-HD individuals and 70 healthy controls underwent two or more (baseline and 1-5 year follow-up) DWI, with n = 649 total sessions. Images were processed using cutting-edge DWI analysis methods for large multicenter studies. Diffusion tensor imaging (DTI) metrics were computed in selected tracts connecting the primary motor, primary somato-sensory, and premotor areas of the cortex with the subcortical caudate and putamen. Pre-HD participants were divided into three CAG-Age Product (CAP) score groups reflecting clinical diagnosis probability (low, medium, or high probabilities). Baseline and longitudinal group differences were examined using linear mixed models. RESULTS: Cross-sectional and longitudinal differences in DTI measures were present in all three CAP groups compared with controls. The high CAP group was most affected. CONCLUSIONS: This is the largest longitudinal DWI study of pre-HD to date. Findings showed DTI differences, consistent with white matter degeneration, were present up to a decade before predicted HD diagnosis. Our findings indicate a unique role for disrupted connectivity between the premotor area and the putamen, which may be closely tied to the onset of motor symptoms in HD.
The Open Anatomy Browser (OABrowser) is an open source, web-based, zero-installation anatomy atlas viewer based on current web browser technologies and evolving anatomy atlas interoperability standards. OABrowser displays three-dimensional anatomical models, image cross-sections of labeled structures and source radiological imaging, and a text-based hierarchy of structures. The viewer includes novel collaborative tools: users can save bookmarks of atlas views for later access and exchange those bookmarks with other users, and dynamic shared views allow groups of users can participate in a collaborative interactive atlas viewing session. We have published several anatomy atlases (an MRI-derived brain atlas and atlases of other parts of the anatomy) to demonstrate OABrowser's functionality. The atlas source data, processing tools, and the source for OABrowser are freely available through GitHub and are distributed under a liberal open source license.
PURPOSE: To investigate the heating of EEG electrodes during magnetic resonance imaging (MRI) scans and to better understand the underlying physical mechanisms with a focus on the antenna effect. MATERIALS AND METHODS: Gold cup and conductive plastic electrodes were placed on small watermelons with fiberoptic probes used to measure electrode temperature changes during a variety of 1.5T and 3T MRI scans. A subset of these experiments was repeated on a healthy human volunteer. RESULTS: The differences between gold and plastic electrodes did not appear to be practically significant. For both electrode types, we observed heating below 4°C for straight wires whose lengths were multiples of ½ the radiofrequency (RF) wavelength and stronger heating (over 15°C) for wire lengths that were odd multiples of ¼ RF wavelength, consistent with the antenna effect. CONCLUSIONS: The antenna effect, which has received little attention so far in the context of EEG-MRI safety, can play as significant a role as the loop effect (from electromagnetic induction) in the heating of EEG electrodes, and therefore wire lengths that are odd multiples of ¼ RF wavelength should be avoided. These results have important implications for the design of EEG electrodes and MRI studies as they help to minimize the risk to patients undergoing MRI with EEG electrodes in place.
Inferring the microstructure of complex media from the diffusive motion of molecules is a challenging problem in diffusion physics. In this paper, we introduce a novel representation of diffusion MRI (dMRI) signal from tissue with spatially-varying diffusivity using a diffusion disturbance function. This disturbance function contains information about the (intra-voxel) spatial fluctuations in diffusivity due to restrictions, hindrances and tissue heterogeneity of the underlying tissue substrate. We derive the short- and long-range disturbance coefficients from this disturbance function to characterize the tissue structure and organization. Moreover, we provide an exact relation between the disturbance coefficients and the time-varying moments of the diffusion propagator, as well as their relation to specific tissue microstructural information such as the intra-axonal volume fraction and the apparent axon radius. The proposed approach is quite general and can model dMRI signal for any type of gradient sequence (rectangular, oscillating, etc.) without using the Gaussian phase approximation. The relevance of the proposed PICASO model is explored using Monte-Carlo simulations and in-vivo dMRI data. The results show that the estimated disturbance coefficients can distinguish different types of microstructural organization of axons.
In schizophrenia, abnormalities in structural connectivity between brain regions known to contain mirror neurons and their relationship to negative symptoms related to a domain of social cognition are not well understood. Diffusion tensor imaging (DTI) scans were acquired in 16 patients with first episode schizophrenia and 16 matched healthy controls. FA and Trace of the tracts interconnecting regions known to be rich in mirror neurons, i.e., anterior cingulate cortex (ACC), inferior parietal lobe (IPL) and premotor cortex (PMC) were evaluated. A significant group effect for Trace was observed in IPL-PMC white matter fiber tract (F (1, 28) = 7.13, p = .012), as well as in the PMC-ACC white matter fiber tract (F (1, 28) = 4.64, p = .040). There were no group differences in FA. In addition, patients with schizophrenia showed a significant positive correlation between the Trace of the left IPL-PMC white matter fiber tract, and the Ability to Feel Intimacy and Closeness score (rho = .57, p = 0.034), and a negative correlation between the Trace of the left PMC-ACC and the Relationships with Friends and Peers score (rho = remove -.54, p = 0.049). We have demonstrated disrupted white mater microstructure within the white matter tracts subserving brain regions containing mirror neurons. We further showed that such structural disruptions might impact negative symptoms and, more specifically, contribute to the inability to feel intimacy (a measure conceptually related to theory of mind) in first episode schizophrenia. Further studies are needed to understand the potential of our results for diagnosis, prognosis and therapeutic interventions.
PURPOSE: Diffusion encoding with asymmetric gradient waveforms is appealing because the asymmetry provides superior efficiency. However, concomitant gradients may cause a residual gradient moment at the end of the waveform, which can cause significant signal error and image artifacts. The purpose of this study was to develop an asymmetric waveform designs for tensor-valued diffusion encoding that is not sensitive to concomitant gradients. METHODS: The "Maxwell index" was proposed as a scalar invariant to capture the effect of concomitant gradients. Optimization of "Maxwell-compensated" waveforms was performed in which this index was constrained. Resulting waveforms were compared to waveforms from literature, in terms of the measured and predicted impact of concomitant gradients, by numerical analysis as well as experiments in a phantom and in a healthy human brain. RESULTS: Maxwell-compensated waveforms with Maxwell indices below 100 (mT/m) ms showed negligible signal bias in both numerical analysis and experiments. By contrast, several waveforms from literature showed gross signal bias under the same conditions, leading to a signal bias that was large enough to markedly affect parameter maps. Experimental results were accurately predicted by theory. CONCLUSION: Constraining the Maxwell index in the optimization of asymmetric gradient waveforms yields efficient diffusion encoding that negates the effects of concomitant fields while enabling arbitrary shapes of the b-tensor. This waveform design is especially useful in combination with strong gradients, long encoding times, thick slices, simultaneous multi-slice acquisition, and large FOVs.
Computational biomechanics of the brain for neurosurgery is an emerging area of research recently gaining in importance and practical applications. This review paper presents the contributions of the Intelligent Systems for Medicine Laboratory and its collaborators to this field, discussing the modeling approaches adopted and the methods developed for obtaining the numerical solutions. We adopt a physics-based modeling approach and describe the brain deformation in mechanical terms (such as displacements, strains, and stresses), which can be computed using a biomechanical model, by solving a continuum mechanics problem. We present our modeling approaches related to geometry creation, boundary conditions, loading, and material properties. From the point of view of solution methods, we advocate the use of fully nonlinear modeling approaches, capable of capturing very large deformations and nonlinear material behavior. We discuss finite element and meshless domain discretization, the use of the total Lagrangian formulation of continuum mechanics, and explicit time integration for solving both time-accurate and steady-state problems. We present the methods developed for handling contacts and for warping 3D medical images using the results of our simulations. We present two examples to showcase these methods: brain shift estimation for image registration and brain deformation computation for neuronavigation in epilepsy treatment.
Jean-Jacques Lemaire, Antonio De Salles, Guillaume Coll, Youssef El Ouadih, Rémi Chaix, Jérôme Coste, Franck Durif, Nikos Makris, and Ron Kikinis. 8/2019. “MRI Atlas of the Human Deep Brain.” Front Neurol, 10, Pp. 851.Abstract
Mastering detailed anatomy of the human deep brain in clinical neurosciences is challenging. Although numerous pioneering works have gathered a large dataset of structural and topographic information, it is still difficult to transfer this knowledge into practice, even with advanced magnetic resonance imaging techniques. Thus, classical histological atlases continue to be used to identify structures for stereotactic targeting in functional neurosurgery. Physicians mainly use these atlases as a template co-registered with the patient's brain. However, it is possible to directly identify stereotactic targets on MRI scans, enabling personalized targeting. In order to help clinicians directly identify deep brain structures relevant to present and future medical applications, we built a volumetric MRI atlas of the deep brain (MDBA) on a large scale (infra millimetric). Twelve hypothalamic, 39 subthalamic, 36 telencephalic, and 32 thalamic structures were identified, contoured, and labeled. Nineteen coronal, 18 axial, and 15 sagittal MRI plates were created. Although primarily designed for direct labeling, the anatomic space was also subdivided in twelfths of AC-PC distance, leading to proportional scaling in the coronal, axial, and sagittal planes. This extensive work is now available to clinicians and neuroscientists, offering another representation of the human deep brain ([https://hal.archives-ouvertes.fr/] [hal-02116633]). The atlas may also be used by computer scientists who are interested in deciphering the topography of this complex region.
Intraoperative tissue deformation, known as brain shift, decreases the benefit of using preoperative images to guide neurosurgery. Non-rigid registration of preoperative magnetic resonance (MR) to intraoperative ultrasound (US) has been proposed as a means to compensate for brain shift. We focus on the initial registration from MR to predurotomy US. We present a method that builds on previous work to address the need for accuracy and generality of MR-iUS registration algorithms in multi-site clinical data. To improve accuracy of registration, we use high-dimensional texture attributes instead of image intensities and propose to replace the standard difference-based attribute matching with correlation-based attribute matching. We also present a strategy that deals explicitly with the large field-of-view mismatch between MR and iUS images. We optimize key parameters across independent MR-iUS brain tumor datasets acquired at three different institutions, with a total of 43 tumor patients and 758 corresponding landmarks to validate the registration algorithm. Despite differences in imaging protocols, patient demographics and landmark distributions, our algorithm was able to reduce landmark errors prior to registration in three data sets (5.37 ± 4.27, 4.18 ± 1.97 and 6.18 ± 3.38 mm, respectively) to a consistently low level (2.28 ± 0.71, 2.08 ± 0.37 and 2.24 ± 0.78 mm, respectively). Our algorithm is compared to 15 other algorithms that have been previously tested on MR-iUS registration and it is competitive with the state-of-the-art on multiple datasets. We show that our algorithm has one of the lowest errors in all datasets (accuracy), and this is achieved while sticking to a fixed set of parameters for multi-site data (generality). In contrast, other algorithms/tools of similar performance need per-dataset parameter tuning (high accuracy but lower generality), and those that stick to fixed parameters have larger errors or inconsistent performance (generality but not the top accuracy). We further characterized landmark errors according to brain regions and tumor types, a topic so far missing in the literature. We found that landmark errors were higher in high-grade than low-grade glioma patients, and higher in tumor regions than in other brain regions.
PURPOSE: In image-guided surgery for glioma removal, neurosurgeons usually plan the resection on images acquired before surgery and use them for guidance during the subsequent intervention. However, after the surgical procedure has begun, the preplanning images become unreliable due to the brain shift phenomenon, caused by modifications of anatomical structures and imprecisions in the neuronavigation system. To obtain an updated view of the resection cavity, a solution is to collect intraoperative data, which can be additionally acquired at different stages of the procedure in order to provide a better understanding of the resection. A spatial mapping between structures identified in subsequent acquisitions would be beneficial. We propose here a fully automated segmentation-based registration method to register ultrasound (US) volumes acquired at multiple stages of neurosurgery. METHODS: We chose to segment sulci and falx cerebri in US volumes, which remain visible during resection. To automatically segment these elements, first we trained a convolutional neural network on manually annotated structures in volumes acquired before the opening of the dura mater and then we applied it to segment corresponding structures in different surgical phases. Finally, the obtained masks are used to register US volumes acquired at multiple resection stages. RESULTS: Our method reduces the mean target registration error (mTRE) between volumes acquired before the opening of the dura mater and during resection from 3.49 mm (± 1.55 mm) to 1.36 mm (± 0.61 mm). Moreover, the mTRE between volumes acquired before opening the dura mater and at the end of the resection is reduced from 3.54 mm (± 1.75 mm) to 2.05 mm (± 1.12 mm). CONCLUSION: The segmented structures demonstrated to be good candidates to register US volumes acquired at different neurosurgical phases. Therefore, our solution can compensate brain shift in neurosurgical procedures involving intraoperative US data.