Search

Although there exist some reference and stereotactic anatomical human brain atlases, there is no equivalent cerebrovascular atlas. A 3D reference atlas of the human cerebrovasculature that is interactive, stereotactic, very detailed, completely parcellated, fully labeled, extendable, dissectible, deformable, and user friendly, is needed in education, training, research, and clinical applications. We constructed a sophisticated and very detailed cerebrovascular atlas with 1325 vascular segments, the smallest of 80 μm in diameter. The atlas was created from multiple 3 and 7T scans by employing tools and methods developed in our lab. In contrast to a sort of "sampled" vascular geometry cited in the literature, we provide information about the "continuous" geometry of the vascular model measurable in 3D at every location and along any vascular segment for both hemispheres. This cerebrovascular atlas was validated taking into account domain knowledge from various fields including angiography, neurosurgery, neuroradiology, (clinical) neuroanatomy, and terminology. Validation was considered as a confirmation of the created vascular model to a typical (conventional) cerebrovascular system with nearly average dimensions. It was done in terms of vessel subdivision, origin, course, surrounding anatomy, patterns, length, diameter, and branches. This validation demonstrates that the constructed cerebrovascular model generally represents a normal, typical cerebrovasculature with its dimensions close to average. This 7T atlas along with the vascular model is a rich source of knowledge about the human cerebrovasculature. The atlas is applicable in education, research, and clinical applications; it has already been applied in deep brain stimulation., (Copyright © 2011 Elsevier Inc. All rights reserved.)

A knowledge-driven algorithm for a rapid, robust, accurate, and automatic extraction of the human cerebral ventricular system from MR neuroimages is proposed. Its novelty is in combination of neuroanatomy, radiological properties, and variability of the ventricular system with image processing techniques. The ventricular system is divided into six 3D regions: bodies and inferior horns of the lateral ventricles, third ventricle, and fourth ventricle. Within each ventricular region, a 2D region of interest (ROI) is defined based on anatomy and variability. Each ventricular region is further subdivided into subregions, and conditions detecting and preventing leakage into the extra-ventricular space are specified for each subregion. The algorithm extracts the ventricular system by (1) processing each ROI (to calculate its local statistics, determine local intensity ranges of cerebrospinal fluid and gray and white matters, set a seed point within the ROI, grow region directionally in 3D, check anti-leakage conditions, and correct growing if leakage occurred) and (2) connecting all unconnected regions grown by relaxing growing conditions. The algorithm was validated qualitatively on 68 and quantitatively on 38 MRI normal and pathological cases (30 clinical, 20 MGH Brain Repository, and 18 MNI BrainWeb data sets). It runs successfully for normal and pathological cases provided that the slice thickness is less than 3.0 mm in axial and less than 2.0 mm in coronal directions, and the data do not have a high inter-slice intensity variability. The algorithm also works satisfactorily in the presence of up to 9% noise and up to 40% RF inhomogeneity for the BrainWeb data. The running time is less than 5 s on a Pentium 4, 2.0 GHz PC. The best overlap metric between the results of a radiology expert and the algorithm is 0.9879 and the worst 0.9527; the mean and standard deviation of the overlap metric are 0.9723 and 0.01087, respectively.

A rapid algorithm for robust, accurate, and automatic extraction of the midsagittal plane (MSP) of the human cerebrum from normal and pathological neuroimages is proposed. The MSP is defined as a plane formed from the interhemispheric fissure line segments having the dominant orientation. The algorithm extracts the MSP in four steps: (1) determine suitable axial slices for processing, (2) localize the fissure line segments on them, (3) select inliers from the extracted fissure line segments through histogram-based outlier removal, and (4) calculate the equation of the MSP from the selected inliers. The fissure line segments are localized by minimizing the local symmetry index characterizing anatomical properties of images in the vicinity of the interhemispheric fissure. A two-stage angular and distance outlier removal is introduced to handle abnormalities. The algorithm has been validated quantitatively with 125 structural MRI and CT cases from 10 centers on three continents by studying its accuracy; tolerance to rotation, noise, asymmetry, and bias field; sensitivity to parameters; and performance. A statistical relationship between algorithm accuracy and the data's adherence to planarity is also determined. The algorithm extracts the MSP below 6 s on Pentium 4 (2.4 GHz) with the average angular and distance errors of (0.40 degrees; 0.63 mm) for normal and (0.59 degrees; 0.73 mm) for pathological cases. The robustness to noise, asymmetry, rotation, and bias field is achieved by extracting the MSP based on the dominant orientation and local symmetry index. A low computational cost results from applying simple operations capturing intrinsic anatomic features, constraining the searching space to the local vicinity of the interhemispheric fissure, and formulating a noniterative algorithm with a coarse and fine fixed-step searching. In comparison to the existing methods, our algorithm is much faster, performs accurately and robustly for a wide range of diversified data, and is fully automatic and thoroughly validated, which make it suitable for clinical applications.

The article introduces an atlas-assisted method and a tool called the Cerefy Neuroradiology Atlas (CNA), available over the Internet for neuroradiology and human brain mapping. The CNA contains an enhanced, extended, and fully segmented and labeled electronic version of the Talairach-Tournoux brain atlas, including parcelated gyri and Brodmann's areas. To our best knowledge, this is the first online, publicly available application with the Talairach-Tournoux atlas. The process of atlas-assisted neuroimage analysis is done in five steps: image data loading, Talairach landmark setting, atlas normalization, image data exploration and analysis, and result saving. Neuroimage analysis is supported by a near-real-time, atlas-to-data warping based on the Talairach transformation. The CNA runs on multiple platforms; is able to process simultaneously multiple anatomical and functional data sets; and provides functions for a rapid atlas-to-data registration, interactive structure labeling and annotating, and mensuration. It is also empowered with several unique features, including interactive atlas warping facilitating fine tuning of atlas-to-data fit, navigation on the triplanar formed by the image data and the atlas, multiple-images-in-one display with interactive atlas-anatomy-function blending, multiple label display, and saving of labeled and annotated image data. The CNA is useful for fast atlas-assisted analysis of neuroimage data sets. It increases accuracy and reduces time in localization analysis of activation regions; facilitates to communicate the information on the interpreted scans from the neuroradiologist to other clinicians and medical students; increases the neuroradiologist's confidence in terms of anatomy and spatial relationships; and serves as a user-friendly, public domain tool for neuroeducation. At present, more than 700 users from five continents have subscribed to the CNA.

The paper introduces an optimal algorithm for rapid calculation of a probabilistic functional atlas (PFA) of subcortical structures from data collected during functional neurosurgery procedures. The PFA is calculated based on combined intraoperative electrophysiology, pre- and intraoperative neuroimaging, and postoperative neurological verification. The algorithm converts the coordinates of the neurologically most effective contacts into probabilistic functional maps taking into account the geometry of a stimulating electrode. The PFA calculation comprises the reconstruction of the contact coordinates from two orthogonal projections, normalizing (warping) the contacts modeled as cylinders, voxelizing the contact models, calculating the atlas, and computing probability. In addition, an analytical representation of the PFA is formulated based on Gaussian modeling. The initial PFA has been calculated from the data collected during the treatment of 274 Parkinson's disease patients, most of them operated bilaterally (487 operated hemispheres). It contains the most popular stereotactic targets, the subthalamic nucleus, globus pallidus internus, and ventral intermedius nucleus. The key application of the algorithm is targeting in stereotactic and functional neurosurgery, and it also can be employed in human and animal brain research.