Imaging the functioning human brain

One of the most exciting methodological advances for brain research field arises in functional brain imaging, which enables us to localize and characterize neural activity and biochemical events in the living human brain. Recently developed event-related functional MRI makes it possible to visualize the brain activity associated with cognitive processes with the temporal resolution of the hemodynamic response. In addition, the high sensitivity and selectivity of positron-emission tomography allow us to probe the neurochemical processes at the molecular level. Positron-emission tomography also has been applied to investigate the effects of therapeutic drugs as well as the effects of drugs of abuse.

Information processing in the human brain: magnetoencephalographic approach.

Rapid progress in effective methods to image brain functions has revolutionized neuroscience. It is now possible to study noninvasively in humans neural processes that were previously only accessible in experimental animals and in brain-injured patients. In this endeavor, positron emission tomography has been the leader, but the superconducting quantum interference device-based magnetoencephalography (MEG) is gaining a firm role, too. With the advent of instruments covering the whole scalp, MEG, typically with 5-mm spatial and 1-ms temporal resolution, allows neuroscientists to track cortical functions accurately in time and space. We present five representative examples of recent MEG studies in our laboratory that demonstrate the usefulness of whole-head magnetoencephalography in investigations of spatiotemporal dynamics of cortical signal processing.

Neural-network-based classification of cognitively normal, demented, Alzheimer disease and vascular dementia from single photon emission with computed tomography image data from brain.

Single photon emission with computed tomography (SPECT) hexamethylphenylethyleneamineoxime technetium-99 images were analyzed by an optimal interpolative neural network (OINN) algorithm to determine whether the network could discriminate among clinically diagnosed groups of elderly normal, Alzheimer disease (AD), and vascular dementia (VD) subjects. After initial image preprocessing and registration, image features were obtained that were representative of the mean regional tissue uptake. These features were extracted from a given image by averaging the intensities over various regions defined by suitable masks. After training, the network classified independent trials of patients whose clinical diagnoses conformed to published criteria for probable AD or probable/possible VD. For the SPECT data used in the current tests, the OINN agreement was 80 and 86% for probable AD and probable/possible VD, respectively. These results suggest that artificial neural network methods offer potential in diagnoses from brain images and possibly in other areas of scientific research where complex patterns of data may have scientifically meaningful groupings that are not easily identifiable by the researcher.

Convergence analysis and validation of low cost distance metrics for computational cost reduction of the Iterative Closest Point algorithm

The Iterative Closest Point algorithm (ICP) is commonly used in engineering applications to solve the rigid registration problem of partially overlapped point sets which are pre-aligned with a coarse estimate of their relative positions. This iterative algorithm is applied in many areas such as the medicine for volumetric reconstruction of tomography data, in robotics to reconstruct surfaces or scenes using range sensor information, in industrial systems for quality control of manufactured objects or even in biology to study the structure and folding of proteins. One of the algorithm’s main problems is its high computational complexity (quadratic in the number of points with the non-optimized original variant) in a context where high density point sets, acquired by high resolution scanners, are processed. Many variants have been proposed in the literature whose goal is the performance improvement either by reducing the number of points or the required iterations or even enhancing the complexity of the most expensive phase: the closest neighbor search. In spite of decreasing its complexity, some of the variants tend to have a negative impact on the final registration precision or the convergence domain thus limiting the possible application scenarios. The goal of this work is the improvement of the algorithm’s computational cost so that a wider range of computationally demanding problems from among the ones described before can be addressed. For that purpose, an experimental and mathematical convergence analysis and validation of point-to-point distance metrics has been performed taking into account those distances with lower computational cost than the Euclidean one, which is used as the de facto standard for the algorithm’s implementations in the literature. In that analysis, the functioning of the algorithm in diverse topological spaces, characterized by different metrics, has been studied to check the convergence, efficacy and cost of the method in order to determine the one which offers the best results. Given that the distance calculation represents a significant part of the whole set of computations performed by the algorithm, it is expected that any reduction of that operation affects significantly and positively the overall performance of the method. As a result, a performance improvement has been achieved by the application of those reduced cost metrics whose quality in terms of convergence and error has been analyzed and validated experimentally as comparable with respect to the Euclidean distance using a heterogeneous set of objects, scenarios and initial situations.

Ultrasonic and X-ray computed tomography characterization of progressive fracture damage in low-porous carbonate rocks

This paper studies the fracturing process in low-porous rocks during uniaxial compressive tests considering the original defects and the new mechanical cracks in the material. For this purpose, five different kinds of rocks have been chosen with carbonate mineralogy and low porosity (lower than 2%). The characterization of the fracture damage is carried out using three different techniques: ultrasounds, mercury porosimetry and X-ray computed tomography. The proposed methodology allows quantifying the evolution of the porous system as well as studying the location of new cracks in the rock samples. Intercrystalline porosity (the smallest pores with pore radius < 1 μm) shows a limited development during loading, disappearing rapidly from the porosimetry curves and it is directly related to the initial plastic behaviour in the stress–strain patterns. However, the biggest pores (corresponding to the cracks) suffer a continuous enlargement until the unstable propagation of fractures. The measured crack initiation stress varies between 0.25 σp and 0.50 σp for marbles and between 0.50 σp and 0.85 σp for micrite limestone. The unstable propagation of cracks is assumed to occur very close to the peak strength. Crack propagation through the sample is completely independent of pre-existing defects (porous bands, stylolites, fractures and veins). The ultrasonic response in the time-domain is less sensitive to the fracture damage than the frequency-domain. P-wave velocity increases during loading test until the beginning of the unstable crack propagation. This increase is higher for marbles (between 15% and 30% from initial vp values) and lower for micrite limestones (between 5% and 10%). When the mechanical cracks propagate unstably, the velocity stops to increase and decreases only when rock damage is very high. Frequency analysis of the ultrasonic signals shows clear changes during the loading process. The spectrum of treated waveforms shows two main frequency peaks centred at low (~ 20 kHz) and high (~ 35 kHz) values. When new fractures appear and grow the amplitude of the high-frequency peak decreases, while that of the low-frequency peak increases. Besides, a slight frequency shift is observed towards higher frequencies.

A Three-Dimensional Stereotaxic MRI Brain Atlas of the Cichlid Fish Oreochromis mossambicus

The African cichlid Oreochromis mossambicus (Mozambique tilapia) has been used as a model system in a wide range of behavioural and neurobiological studies. The increasing number of genetic tools available for this species, together with the emerging interest in its use for neurobiological studies, increased the need for an accurate hodological mapping of the tilapia brain to supplement the available histological data. The goal of our study was to elaborate a three-dimensional, high-resolution digital atlas using magnetic resonance imaging, supported by Nissl staining. Resulting images were viewed and analysed in all orientations (transverse, sagittal, and horizontal) and manually labelled to reveal structures in the olfactory bulb, telencephalon, diencephalon, optic tectum, and cerebellum. This high resolution tilapia brain atlas is expected to become a very useful tool for neuroscientists using this fish model and will certainly expand their use in future studies regarding the central nervous system.

Comment on “Low time resolution analysis of polar ice cores cannot detect impulsive nitrate events” by D.F. Smart et al.

Smart et al. (2014) suggested that the detection of nitrate spikes in polar ice cores from solar energetic particle (SEP) events could be achieved if an analytical system with sufficiently high resolution was used. Here we show that the spikes they associate with SEP events are not reliably recorded in cores from the same location, even when the resolution is clearly adequate. We explain the processes that limit the effective resolution of ice cores. Liquid conductivity data suggest that the observed spikes are associated with sodium or another nonacidic cation, making it likely that they result from deposition of sea salt or similar aerosol that has scavenged nitrate, rather than from a primary input of nitrate in the troposphere. We consider that there is no evidence at present to support the identification of any spikes in nitrate as representing SEP events. Although such events undoubtedly create nitrate in the atmosphere, we see no plausible route to using nitrate spikes to document the statistics of such events.

Quarterly technical progress and schedule report "high resolution, low cost solar cell contact development" : for the quarter ending 24 August 1979 /

"21 September 1979."

Widespread brain activity during an abdominal task markedly reduced after pain physiology education: fMRI evaluation of a single patient with chronic low back pain

The way people with chronic low back pain think about pain can affect the way they move. This case report concerns a patient with chronic disabling low back pain who underwent functional magnetic resonance imaging scans during performance of a voluntary trunk muscle task under three conditions: directly after training in the task and, after one week of practice, before and after a 2.5 hour pain physiology education session. Before education there was widespread brain activity during performance of the task, including activity in cortical regions known to be involved in pain, although the task was not painful. After education widespread activity was absent so that there was no brain activation outside of the primary somatosensory cortex. The results suggest that pain physiology education markedly altered brain activity during performance of the task. The data offer a possible mechanism for difficulty in acquisition of trunk muscle training in people with pain and suggest that the change in activity associated with education may reflect reduced threat value of the task.

The DVD animal model of schizophrenia: Effects of low maternal vitamin D on brain dopamine

Background: We have previously shown that the offspring of vitamin D3 depleted rats have enlarged ventricles and altered neurotrophin profiles (reduced NGF and GDNF). These findings enhance the biological plausibility that low prenatal vitamin D may be a risk factor for schizophrenia. Our recent behavioural studies have found that adult rats with developmental vitamin D deficiency (DVD) have a subtle increase in baseline locomotor activity and a heightened response to dopamine (DA) antagonists. The aim of this study was to investigate brain DA neurochemistry in the DVD model. Methods: We examined cerebrums and striatal tissue from neonates and a variety of brain tissues from the remaining littermates at adulthood. DA, DOPAC, HVA, serotonin and 5HIAA were analysed by HPLC. Single point comparisons for DA1, DA2 and NMDA receptors were also assessed in these tissues. Results: Significant increases in DA and HVA were found in brains from DVD deplete neonates (P=0.01). However, DA and its metabolites were not increased in either the neonate or adult striatum, however there was a trend towards increased DA and its metabolites in the accumbens (P=0.1). Receptor densities were unaffected by prenatal vitamin D levels. Conclusions: Although the effect of maternal diet appears to increase DA production and turnover in neonatal brain, this does not persist into adulthood. Thus other factors must underlie the increased locomotor activity noted in these animals. Future experiments will concentrate on monitoring accumbens and striatal DA release and turnover using microdialysis in pharmacologically challenged behavioural paradigms. References: Eyles D, Brown J; Mackay-Sim A, McGrath J, Feron F. (2003) Vitamin D3 and brain development. Neuroscience 118 (3) 641–653. Burne T, McGrath J, Eyles D, Mackay-Sim A. Behavioural characterization of vitamin D receptor knockout mice. (2005) Behavioural Brain Res: 157 299–308.

Femtosecond laser micro-inscription of optical coherence tomography resolution test artifacts

Optical coherence tomography (OCT) systems are becoming more commonly used in biomedical imaging and, to enable continued uptake, a reliable method of characterizing their performance and validating their operation is required. This paper outlines the use of femtosecond laser subsurface micro-inscription techniques to fabricate an OCT test artifact for validating the resolution performance of a commercial OCT system. The key advantage of this approach is that by utilizing the nonlinear absorption a three dimensional grid of highly localized point and line defects can be written in clear fused silica substrates.

Optical coherence tomography confirmed resolution of macular edema in a diabetic patient after hemodialysis

Purpose: The authors report the first case, to their knowledge, of resolution of diffuse macular edema after hemodialysis, which has been confirmed by optical coherence tomography. Methods: A 53-year-old white woman with type 2 diabetes developed worsening macular edema and was examined in the ophthalmology clinic and scheduled for macular grid laser photocoagulation. The laser, however, was deferred for 4 weeks because she had also developed end-stage renal failure and required hemodialysis. Results: When she was reviewed in the ophthalmology clinic 4 weeks later for laser therapy, it was found that her macular edema had resolved, vision had improved, and laser was unnecessary. Review at 6 months showed that the macular edema remains resolved. Conclusion: Ophthalmologists should be aware that hemodialysis may reduce macular edema in such patients. Copyright © by Ophthalmic Communications Society Inc.

Advanced Techniques for Image Quality Assessment of Modern X-ray Computed Tomography Systems

X-ray computed tomography (CT) imaging constitutes one of the most widely used diagnostic tools in radiology today with nearly 85 million CT examinations performed in the U.S in 2011. CT imparts a relatively high amount of radiation dose to the patient compared to other x-ray imaging modalities and as a result of this fact, coupled with its popularity, CT is currently the single largest source of medical radiation exposure to the U.S. population. For this reason, there is a critical need to optimize CT examinations such that the dose is minimized while the quality of the CT images is not degraded. This optimization can be difficult to achieve due to the relationship between dose and image quality. All things being held equal, reducing the dose degrades image quality and can impact the diagnostic value of the CT examination.

A recent push from the medical and scientific community towards using lower doses has spawned new dose reduction technologies such as automatic exposure control (i.e., tube current modulation) and iterative reconstruction algorithms. In theory, these technologies could allow for scanning at reduced doses while maintaining the image quality of the exam at an acceptable level. Therefore, there is a scientific need to establish the dose reduction potential of these new technologies in an objective and rigorous manner. Establishing these dose reduction potentials requires precise and clinically relevant metrics of CT image quality, as well as practical and efficient methodologies to measure such metrics on real CT systems. The currently established methodologies for assessing CT image quality are not appropriate to assess modern CT scanners that have implemented those aforementioned dose reduction technologies.

Thus the purpose of this doctoral project was to develop, assess, and implement new phantoms, image quality metrics, analysis techniques, and modeling tools that are appropriate for image quality assessment of modern clinical CT systems. The project developed image quality assessment methods in the context of three distinct paradigms, (a) uniform phantoms, (b) textured phantoms, and (c) clinical images.

The work in this dissertation used the “task-based” definition of image quality. That is, image quality was broadly defined as the effectiveness by which an image can be used for its intended task. Under this definition, any assessment of image quality requires three components: (1) A well defined imaging task (e.g., detection of subtle lesions), (2) an “observer” to perform the task (e.g., a radiologists or a detection algorithm), and (3) a way to measure the observer’s performance in completing the task at hand (e.g., detection sensitivity/specificity).

First, this task-based image quality paradigm was implemented using a novel multi-sized phantom platform (with uniform background) developed specifically to assess modern CT systems (Mercury Phantom, v3.0, Duke University). A comprehensive evaluation was performed on a state-of-the-art CT system (SOMATOM Definition Force, Siemens Healthcare) in terms of noise, resolution, and detectability as a function of patient size, dose, tube energy (i.e., kVp), automatic exposure control, and reconstruction algorithm (i.e., Filtered Back-Projection– FPB vs Advanced Modeled Iterative Reconstruction– ADMIRE). A mathematical observer model (i.e., computer detection algorithm) was implemented and used as the basis of image quality comparisons. It was found that image quality increased with increasing dose and decreasing phantom size. The CT system exhibited nonlinear noise and resolution properties, especially at very low-doses, large phantom sizes, and for low-contrast objects. Objective image quality metrics generally increased with increasing dose and ADMIRE strength, and with decreasing phantom size. The ADMIRE algorithm could offer comparable image quality at reduced doses or improved image quality at the same dose (increase in detectability index by up to 163% depending on iterative strength). The use of automatic exposure control resulted in more consistent image quality with changing phantom size.

Based on those results, the dose reduction potential of ADMIRE was further assessed specifically for the task of detecting small (<=6 mm) low-contrast (<=20 HU) lesions. A new low-contrast detectability phantom (with uniform background) was designed and fabricated using a multi-material 3D printer. The phantom was imaged at multiple dose levels and images were reconstructed with FBP and ADMIRE. Human perception experiments were performed to measure the detection accuracy from FBP and ADMIRE images. It was found that ADMIRE had equivalent performance to FBP at 56% less dose.

Using the same image data as the previous study, a number of different mathematical observer models were implemented to assess which models would result in image quality metrics that best correlated with human detection performance. The models included naïve simple metrics of image quality such as contrast-to-noise ratio (CNR) and more sophisticated observer models such as the non-prewhitening matched filter observer model family and the channelized Hotelling observer model family. It was found that non-prewhitening matched filter observers and the channelized Hotelling observers both correlated strongly with human performance. Conversely, CNR was found to not correlate strongly with human performance, especially when comparing different reconstruction algorithms.

The uniform background phantoms used in the previous studies provided a good first-order approximation of image quality. However, due to their simplicity and due to the complexity of iterative reconstruction algorithms, it is possible that such phantoms are not fully adequate to assess the clinical impact of iterative algorithms because patient images obviously do not have smooth uniform backgrounds. To test this hypothesis, two textured phantoms (classified as gross texture and fine texture) and a uniform phantom of similar size were built and imaged on a SOMATOM Flash scanner (Siemens Healthcare). Images were reconstructed using FBP and a Sinogram Affirmed Iterative Reconstruction (SAFIRE). Using an image subtraction technique, quantum noise was measured in all images of each phantom. It was found that in FBP, the noise was independent of the background (textured vs uniform). However, for SAFIRE, noise increased by up to 44% in the textured phantoms compared to the uniform phantom. As a result, the noise reduction from SAFIRE was found to be up to 66% in the uniform phantom but as low as 29% in the textured phantoms. Based on this result, it clear that further investigation was needed into to understand the impact that background texture has on image quality when iterative reconstruction algorithms are used.

To further investigate this phenomenon with more realistic textures, two anthropomorphic textured phantoms were designed to mimic lung vasculature and fatty soft tissue texture. The phantoms (along with a corresponding uniform phantom) were fabricated with a multi-material 3D printer and imaged on the SOMATOM Flash scanner. Scans were repeated a total of 50 times in order to get ensemble statistics of the noise. A novel method of estimating the noise power spectrum (NPS) from irregularly shaped ROIs was developed. It was found that SAFIRE images had highly locally non-stationary noise patterns with pixels near edges having higher noise than pixels in more uniform regions. Compared to FBP, SAFIRE images had 60% less noise on average in uniform regions for edge pixels, noise was between 20% higher and 40% lower. The noise texture (i.e., NPS) was also highly dependent on the background texture for SAFIRE. Therefore, it was concluded that quantum noise properties in the uniform phantoms are not representative of those in patients for iterative reconstruction algorithms and texture should be considered when assessing image quality of iterative algorithms.

The move beyond just assessing noise properties in textured phantoms towards assessing detectability, a series of new phantoms were designed specifically to measure low-contrast detectability in the presence of background texture. The textures used were optimized to match the texture in the liver regions actual patient CT images using a genetic algorithm. The so called “Clustured Lumpy Background” texture synthesis framework was used to generate the modeled texture. Three textured phantoms and a corresponding uniform phantom were fabricated with a multi-material 3D printer and imaged on the SOMATOM Flash scanner. Images were reconstructed with FBP and SAFIRE and analyzed using a multi-slice channelized Hotelling observer to measure detectability and the dose reduction potential of SAFIRE based on the uniform and textured phantoms. It was found that at the same dose, the improvement in detectability from SAFIRE (compared to FBP) was higher when measured in a uniform phantom compared to textured phantoms.

The final trajectory of this project aimed at developing methods to mathematically model lesions, as a means to help assess image quality directly from patient images. The mathematical modeling framework is first presented. The models describe a lesion’s morphology in terms of size, shape, contrast, and edge profile as an analytical equation. The models can be voxelized and inserted into patient images to create so-called “hybrid” images. These hybrid images can then be used to assess detectability or estimability with the advantage that the ground truth of the lesion morphology and location is known exactly. Based on this framework, a series of liver lesions, lung nodules, and kidney stones were modeled based on images of real lesions. The lesion models were virtually inserted into patient images to create a database of hybrid images to go along with the original database of real lesion images. ROI images from each database were assessed by radiologists in a blinded fashion to determine the realism of the hybrid images. It was found that the radiologists could not readily distinguish between real and virtual lesion images (area under the ROC curve was 0.55). This study provided evidence that the proposed mathematical lesion modeling framework could produce reasonably realistic lesion images.

Based on that result, two studies were conducted which demonstrated the utility of the lesion models. The first study used the modeling framework as a measurement tool to determine how dose and reconstruction algorithm affected the quantitative analysis of liver lesions, lung nodules, and renal stones in terms of their size, shape, attenuation, edge profile, and texture features. The same database of real lesion images used in the previous study was used for this study. That database contained images of the same patient at 2 dose levels (50% and 100%) along with 3 reconstruction algorithms from a GE 750HD CT system (GE Healthcare). The algorithms in question were FBP, Adaptive Statistical Iterative Reconstruction (ASiR), and Model-Based Iterative Reconstruction (MBIR). A total of 23 quantitative features were extracted from the lesions under each condition. It was found that both dose and reconstruction algorithm had a statistically significant effect on the feature measurements. In particular, radiation dose affected five, three, and four of the 23 features (related to lesion size, conspicuity, and pixel-value distribution) for liver lesions, lung nodules, and renal stones, respectively. MBIR significantly affected 9, 11, and 15 of the 23 features (including size, attenuation, and texture features) for liver lesions, lung nodules, and renal stones, respectively. Lesion texture was not significantly affected by radiation dose.

The second study demonstrating the utility of the lesion modeling framework focused on assessing detectability of very low-contrast liver lesions in abdominal imaging. Specifically, detectability was assessed as a function of dose and reconstruction algorithm. As part of a parallel clinical trial, images from 21 patients were collected at 6 dose levels per patient on a SOMATOM Flash scanner. Subtle liver lesion models (contrast = -15 HU) were inserted into the raw projection data from the patient scans. The projections were then reconstructed with FBP and SAFIRE (strength 5). Also, lesion-less images were reconstructed. Noise, contrast, CNR, and detectability index of an observer model (non-prewhitening matched filter) were assessed. It was found that SAFIRE reduced noise by 52%, reduced contrast by 12%, increased CNR by 87%. and increased detectability index by 65% compared to FBP. Further, a 2AFC human perception experiment was performed to assess the dose reduction potential of SAFIRE, which was found to be 22% compared to the standard of care dose.

In conclusion, this dissertation provides to the scientific community a series of new methodologies, phantoms, analysis techniques, and modeling tools that can be used to rigorously assess image quality from modern CT systems. Specifically, methods to properly evaluate iterative reconstruction have been developed and are expected to aid in the safe clinical implementation of dose reduction technologies.