A Stereoscopic Volume Rendered Brain Atlas

The anatomy of the human brain is organized as a complex arrangement of interrelated structures in three dimensional space. To facilitate the understanding of both structure and function, we have created a volume rendered brain atlas (VRBA) with an intuitive interface that allows real-time stereoscopic rendering of brain anatomy. The VRBA incorporates 2-dimensional and 3-dimensional texture mapping to display segmented brain anatomy co-registered with a T1 MRI. The interface allows the user to remove and add any of the 62 brain structures, as well as control the display of the MRI dataset. The atlas also contains brief verbal and written descriptions of the different anatomical regions to correlate structure with function. A variety of stereoscopic projection methods are supported by the VRBA and provide an abstract, yet simple, way of visualizing brain anatomy and 3-dimensional relationships between different nuclei.

Advanced luminance control and black offset correction for multi-projector display systems

In order to display a homogeneous image using multiple projectors, differences in the projected intensities must be compensated. In this paper, we present novel approaches to combine and extend existing techniques for edge blending and luminance harmonization to achieve a detailed luminance control. Furthermore, we apply techniques for improving the contrast ratio of multi-segmented displays also to the black offset correction. We also present a simple scheme to involve the displayed context in the correction process to dynamically improve the contrast in brighter images. In addition, we present a metric to evaluate the different methods and their influence on the visual quality.

Monitoring the Time Course of Perception without Awareness: A Comparison of Mirror Masked Words and Nonwords

Mirror masked words are embedded into a context that makes them appear as senseless patterns or as strings of unfamiliar letters. Thus, mirror masked words can be shown for several hundreds of milliseconds without being recognised as words. We sought to further investigate effects of nonsconscious reading by monitoring event-related brain potentials (ERPs) while participants observed mirror masked letter strings. ERPs were recorded while participants observed mirror masked words and nonwords. Data of 15 participants was segmented into periods of quasi-stable field topography (microstates). Microstates for masked words and nonwords were compared using randomization tests, statistical parametric scalp maps and Low Resolution Electromagnetic Tomography (LORETA). ERPs to masked words and nonwords showed significant topographic differences between 136 and 256 ms, indicating that stimuli were nonconsciously discriminated. A LORETA model localised sources of activation discriminating between masked words and nonwords in left operculum, the right superior parietal lobe and right superior temporal gyrus indicating higher current density for nonwords than for words in these areas. ERPs of mirror masked stimuli can indicate unconscious discrimination even in cases where behavioural priming is unreliable. This approach might be useful for investigating differences in early, nonconscious stages of word perception.

Prestimulus EEG microstates influence visual event-related potential microstates in field maps with 47 channels

The influence of the immediate prestimulus EEG microstate (sub-second epoch of stable topography/map landscape) on the map landscape of visually evoked 47-channel event-related potential (ERP) microstates was examined using the frequent, non-target stimuli of a cognitive paradigm (12 volunteers). For the two most frequent prestimulus microstate classes (oriented left anterior-right posterior and right anterior-left posterior), ERP map series were selectively averaged. The post-stimulus ERP grand average map series was segmented into microstates; 10 were found. The centroid locations of positive and negative map areas were extracted as landscape descriptors. Significant differences (MANOVAs and t-tests) between the two prestimulus classes were found in four of the ten ERP microstates. The relative orientation of the two ERP microstate classes was the same as prestimulus in some ERP microstates, but reversed in others. — Thus, brain electric microstates at stimulus arrival influence the landscapes of the post-stimulus ERP maps and therefore, information processing; prestimulus microstate effects differed for different post-stimulus ERP microstates.

Brain electric microstates and momentary conscious mind states as building blocks of spontaneous thinking: I. Visual imagery and abstract thoughts

Prompted reports of recall of spontaneous, conscious experiences were collected in a no-input, no-task, no-response paradigm (30 random prompts to each of 13 healthy volunteers). The mentation reports were classified into visual imagery and abstract thought. Spontaneous 19-channel brain electric activity (EEG) was continuously recorded, viewed as series of momentary spatial distributions (maps) of the brain electric field and segmented into microstates, i.e. into time segments characterized by quasi-stable landscapes of potential distribution maps which showed varying durations in the sub-second range. Microstate segmentation used a data-driven strategy. Different microstates, i.e. different brain electric landscapes must have been generated by activity of different neural assemblies and therefore are hypothesized to constitute different functions. The two types of reported experiences were associated with significantly different microstates (mean duration 121 ms) immediately preceding the prompts; these microstates showed, across subjects, for abstract thought (compared to visual imagery) a shift of the electric gravity center to the left and a clockwise rotation of the field axis. Contrariwise, the microstates 2 s before the prompt did not differ between the two types of experiences. The results support the hypothesis that different microstates of the brain as recognized in its electric field implement different conscious, reportable mind states, i.e. different classes (types) of thoughts (mentations); thus, the microstates might be candidates for the `atoms of thought'.

Quantification of Popliteal Artery Deformation During Leg Flexion in Subjects With Peripheral Artery Disease: A Pilot Study

Purpose: To quantify the in vivo deformations of the popliteal artery during leg flexion in subjects with clinically relevant peripheral artery disease (PAD). Methods: Five patients (4 men; mean age 69 years, range 56–79) with varying calcification levels of the popliteal artery undergoing endovascular revascularization underwent 3-dimensional (3D) rotational angiography. Image acquisition was performed with the leg straight and with a flexion of 70°/20° in the knee/hip joints. The arterial centerline and the corresponding branches in both positions were segmented to create 3D reconstructions of the arterial trees. Axial deformation, twisting, and curvatures were quantified. Furthermore, the relationships between the calcification levels and the deformations were investigated. Results: An average shortening of 5.9%±2.5% and twist rate of 3.8±2.2°/cm in the popliteal artery were observed. Maximal curvatures in the straight and flexed positions were 0.12±0.04 cm−1 and 0.24±0.09 cm−1, respectively. As the severity of calcification increased, the maximal curvature in the straight position increased from 0.08 to 0.17 cm−1, while an increase from 0.17 to 0.39 cm−1 was observed for the flexed position. Axial elongations and arterial twisting were not affected by the calcification levels. Conclusion: The popliteal artery of patients with symptomatic PAD is exposed to significant deformations during flexion of the knee joint. The severity of calcification directly affects curvature, but not arterial length or twisting angles. This pilot study also showed the ability of rotational angiography to quantify the 3D deformations of the popliteal artery in patients with various levels of calcification.

In vitro accuracy evaluation of image-guided robot system for direct cochlear access

HYPOTHESIS A previously developed image-guided robot system can safely drill a tunnel from the lateral mastoid surface, through the facial recess, to the middle ear, as a viable alternative to conventional mastoidectomy for cochlear electrode insertion. BACKGROUND Direct cochlear access (DCA) provides a minimally invasive tunnel from the lateral surface of the mastoid through the facial recess to the middle ear for cochlear electrode insertion. A safe and effective tunnel drilled through the narrow facial recess requires a highly accurate image-guided surgical system. Previous attempts have relied on patient-specific templates and robotic systems to guide drilling tools. In this study, we report on improvements made to an image-guided surgical robot system developed specifically for this purpose and the resulting accuracy achieved in vitro. MATERIALS AND METHODS The proposed image-guided robotic DCA procedure was carried out bilaterally on 4 whole head cadaver specimens. Specimens were implanted with titanium fiducial markers and imaged with cone-beam CT. A preoperative plan was created using a custom software package wherein relevant anatomical structures of the facial recess were segmented, and a drill trajectory targeting the round window was defined. Patient-to-image registration was performed with the custom robot system to reference the preoperative plan, and the DCA tunnel was drilled in 3 stages with progressively longer drill bits. The position of the drilled tunnel was defined as a line fitted to a point cloud of the segmented tunnel using principle component analysis (PCA function in MatLab). The accuracy of the DCA was then assessed by coregistering preoperative and postoperative image data and measuring the deviation of the drilled tunnel from the plan. The final step of electrode insertion was also performed through the DCA tunnel after manual removal of the promontory through the external auditory canal. RESULTS Drilling error was defined as the lateral deviation of the tool in the plane perpendicular to the drill axis (excluding depth error). Errors of 0.08 ± 0.05 mm and 0.15 ± 0.08 mm were measured on the lateral mastoid surface and at the target on the round window, respectively (n =8). Full electrode insertion was possible for 7 cases. In 1 case, the electrode was partially inserted with 1 contact pair external to the cochlea. CONCLUSION The purpose-built robot system was able to perform a safe and reliable DCA for cochlear implantation. The workflow implemented in this study mimics the envisioned clinical procedure showing the feasibility of future clinical implementation.

Increasing the dynamic range for the analysis of boron in PGAA

Prompt gamma activation analysis (PGAA) is especially sensitive for elements with high neutron-capture cross sections, like boron, which can be detected down to a level of ng/g. However, if it is a major component, the high count rate from its signal will distort the spectra, making the evaluation difficult. A lead attenuator was introduced in front of the HPGe-detector to reduce low-energy gamma radiation and specifically the boron gamma rays reaching the detector, whose thickness was found to be optimal at 10 mm. Detection efficiencies with and without the lead attenuator were compared, and it was shown that the dynamic range of the PGAA technique was significantly increased. The method was verified with the analyses of stoichiometric compounds: TiB2, NiB, PVC, Alborex, and Alborite.

Topographic bone thickness maps for Bonebridge implantations.

Bonebridge™ (BB) implantation relies on optimal anchoring of the bone-conduction implant in the temporal bone. Preoperative position planning has to account for the available bone thickness minimizing unwanted interference with underlying anatomical structures. This study describes the first clinical experience with a planning method based on topographic bone thickness maps (TBTM) for presigmoid BB implantations. The temporal bone was segmented enabling three-dimensional surface generation. Distances between the external and internal surface were color encoded and mapped to a TBTM. Suitable implant positions were planned with reference to the TBTM. Surgery was performed according to the standard procedure (n = 7). Computation of the TBTM and consecutive implant position planning took 70 min on average for a trained technician. Surgical time for implantations under passive TBTM image guidance was 60 min, on average. The sigmoid sinus (n = 5) and dura mater (n = 1) were exposed, as predicted with the TBTM. Feasibility of the TBTM method was shown for standard presigmoid BB implantations. The projection of three-dimensional bone thickness information into a single topographic map provides the surgeon with an intuitive display of the anatomical situation prior to implantation. Nevertheless, TBTM generation time has to be significantly reduced to simplify integration in clinical routine.

Insular and caudate lesions release abnormal yawning in stroke patients.

Abnormal yawning is an underappreciated phenomenon in patients with ischemic stroke. We aimed at identifying frequently affected core regions in the supratentorial brain of stroke patients with abnormal yawning and contributing to the anatomical network concept of yawning control. Ten patients with acute anterior circulation stroke and ≥3 yawns/15 min without obvious cause were analyzed. The NIH stroke scale (NIHSS), Glasgow Coma Scale (GCS), symptom onset, period with abnormal yawning, blood oxygen saturation, glucose, body temperature, blood pressure, heart rate, and modified Rankin scale (mRS) were assessed for all patients. MRI lesion maps were segmented on diffusion-weighted images, spatially normalized, and the extent of overlap between the different stroke patterns was determined. Correlations between the period with abnormal yawning and the apparent diffusion coefficient (ADC) in the overlapping regions, total stroke volume, NIHSS and mRS were performed. Periods in which patients presented with episodes of abnormal yawning lasted on average for 58 h. Average GCS, NIHSS, and mRS scores were 12.6, 11.6, and 3.5, respectively. Clinical parameters were within normal limits. Ischemic brain lesions overlapped in nine out of ten patients: in seven patients in the insula and in seven in the caudate nucleus. The decrease of the ADC within the lesions correlated with the period with abnormal yawing (r = -0.76, Bonferroni-corrected p = 0.02). The stroke lesion intensity of the common overlapping regions in the insula and the caudate nucleus correlates with the period with abnormal yawning. The insula might be the long sought-after brain region for serotonin-mediated yawning.

Automatic quantitation of multiple sclerosis lesions on MR images

A two-pronged approach for the automatic quantitation of multiple sclerosis (MS) lesions on magnetic resonance (MR) images has been developed. This method includes the design and use of a pulse sequence for improved lesion-to-tissue contrast (LTC) and seeks to identify and minimize the sources of false lesion classifications in segmented images. The new pulse sequence, referred to as AFFIRMATIVE (Attenuation of Fluid by Fast Inversion Recovery with MAgnetization Transfer Imaging with Variable Echoes), improves the LTC, relative to spin-echo images, by combining Fluid-Attenuated Inversion Recovery (FLAIR) and Magnetization Transfer Contrast (MTC). In addition to acquiring fast FLAIR/MTC images, the AFFIRMATIVE sequence simultaneously acquires fast spin-echo (FSE) images for spatial registration of images, which is necessary for accurate lesion quantitation. Flow has been found to be a primary source of false lesion classifications. Therefore, an imaging protocol and reconstruction methods are developed to generate "flow images" which depict both coherent (vascular) and incoherent (CSF) flow. An automatic technique is designed for the removal of extra-meningeal tissues, since these are known to be sources of false lesion classifications. A retrospective, three-dimensional (3D) registration algorithm is implemented to correct for patient movement which may have occurred between AFFIRMATIVE and flow imaging scans. Following application of these pre-processing steps, images are segmented into white matter, gray matter, cerebrospinal fluid, and MS lesions based on AFFIRMATIVE and flow images using an automatic algorithm. All algorithms are seamlessly integrated into a single MR image analysis software package. Lesion quantitation has been performed on images from 15 patient volunteers. The total processing time is less than two hours per patient on a SPARCstation 20. The automated nature of this approach should provide an objective means of monitoring the progression, stabilization, and/or regression of MS lesions in large-scale, multi-center clinical trials. ^

Thermal Conditions at the Bed of the Laurentide Ice Sheet in Maine During Deglaciation: Implications for Esker Formation

The University of Maine Ice Sheet Model was used to study basal conditions during retreat of the Laurentide ice sheet in Maine. Within 150 km of the margin, basal melt rates average similar to 5 mm a(-1) during retreat. They decline over the next 100km, so areas of frozen bed develop in northern Maine during retreat. By integrating the melt rate over the drainage area typically subtended by an esker, we obtained a discharge at the margin of similar to 1.2 m(3) s(-1). While such a discharge could have moved the material in the Katahdin esker, it was likely too low to build the esker in the time available. Additional water from the glacier surface was required. Temperature gradients in the basal ice increase rapidly with distance from the margin. By conducting upward into the ice all of the additional viscous heat produced by any perturbation that increases the depth of flow in a flat conduit in a distributed drainage system, these gradients inhibit the formation of sharply arched conduits in which an esker can form. This may explain why eskers commonly seem to form near the margin and are typically segmented, with later segments overlapping onto earlier ones.

Segmentation and Recognition of Multi-Food Meal Images for Carbohydrate Counting

In this paper, we propose novel methodologies for the automatic segmentation and recognition of multi-food images. The proposed methods implement the first modules of a carbohydrate counting and insulin advisory system for type 1 diabetic patients. Initially the plate is segmented using pyramidal mean-shift filtering and a region growing algorithm. Then each of the resulted segments is described by both color and texture features and classified by a support vector machine into one of six different major food classes. Finally, a modified version of the Huang and Dom evaluation index was proposed, addressing the particular needs of the food segmentation problem. The experimental results prove the effectiveness of the proposed method achieving a segmentation accuracy of 88.5% and recognition rate equal to 87%

Landmark Detection for Fusion of Fundus and MRI Towards a Patient-Specific Multi-Modal Eye Model

Ophthalmologists typically acquire different image modalities to diagnose eye pathologies. They comprise e.g., Fundus photography, Optical Coherence Tomography (OCT), Computed Tomography (CT) and Magnetic Resonance Imaging (MRI). Yet, these images are often complementary and do express the same pathologies in a different way. Some pathologies are only visible in a particular modality. Thus, it is beneficial for the ophthalmologist to have these modalities fused into a single patient-specific model. The presented article’s goal is a fusion of Fundus photography with segmented MRI volumes. This adds information to MRI which was not visible before like vessels and the macula. This article’s contributions include automatic detection of the optic disc, the fovea, the optic axis and an automatic segmentation of the vitreous humor of the eye.

3D morphometry using automated aortic segmentation in native MR angiography: an alternative to contrast enhanced MRA?

INTRODUCTION Native-MR angiography (N-MRA) is considered an imaging alternative to contrast enhanced MR angiography (CE-MRA) for patients with renal insufficiency. Lower intraluminal contrast in N-MRA often leads to failure of the segmentation process in commercial algorithms. This study introduces an in-house 3D model-based segmentation approach used to compare both sequences by automatic 3D lumen segmentation, allowing for evaluation of differences of aortic lumen diameters as well as differences in length comparing both acquisition techniques at every possible location. METHODS AND MATERIALS Sixteen healthy volunteers underwent 1.5-T-MR Angiography (MRA). For each volunteer, two different MR sequences were performed, CE-MRA: gradient echo Turbo FLASH sequence and N-MRA: respiratory-and-cardiac-gated, T2-weighted 3D SSFP. Datasets were segmented using a 3D model-based ellipse-fitting approach with a single seed point placed manually above the celiac trunk. The segmented volumes were manually cropped from left subclavian artery to celiac trunk to avoid error due to side branches. Diameters, volumes and centerline length were computed for intraindividual comparison. For statistical analysis the Wilcoxon-Signed-Ranked-Test was used. RESULTS Average centerline length obtained based on N-MRA was 239.0±23.4 mm compared to 238.6±23.5 mm for CE-MRA without significant difference (P=0.877). Average maximum diameter obtained based on N-MRA was 25.7±3.3 mm compared to 24.1±3.2 mm for CE-MRA (P<0.001). In agreement with the difference in diameters, volumes obtained based on N-MRA (100.1±35.4 cm(3)) were consistently and significantly larger compared to CE-MRA (89.2±30.0 cm(3)) (P<0.001). CONCLUSIONS 3D morphometry shows highly similar centerline lengths for N-MRA and CE-MRA, but systematically higher diameters and volumes for N-MRA.