Comparison of electric field strength and spatial distribution of electroconvulsive therapy and magnetic seizure therapy in a realistic human head model.

BACKGROUND: This study examines the strength and spatial distribution of the electric field induced in the brain by electroconvulsive therapy (ECT) and magnetic seizure therapy (MST). METHODS: The electric field induced by standard (bilateral, right unilateral, and bifrontal) and experimental (focal electrically administered seizure therapy and frontomedial) ECT electrode configurations as well as a circular MST coil configuration was simulated in an anatomically realistic finite element model of the human head. Maps of the electric field strength relative to an estimated neural activation threshold were used to evaluate the stimulation strength and focality in specific brain regions of interest for these ECT and MST paradigms and various stimulus current amplitudes. RESULTS: The standard ECT configurations and current amplitude of 800-900mA produced the strongest overall stimulation with median of 1.8-2.9 times neural activation threshold and more than 94% of the brain volume stimulated at suprathreshold level. All standard ECT electrode placements exposed the hippocampi to suprathreshold electric field, although there were differences across modalities with bilateral and right unilateral producing respectively the strongest and weakest hippocampal stimulation. MST stimulation is up to 9 times weaker compared to conventional ECT, resulting in direct activation of only 21% of the brain. Reducing the stimulus current amplitude can make ECT as focal as MST. CONCLUSIONS: The relative differences in electric field strength may be a contributing factor for the cognitive sparing observed with right unilateral compared to bilateral ECT, and MST compared to right unilateral ECT. These simulations could help understand the mechanisms of seizure therapies and develop interventions with superior risk/benefit ratio.

In vivo snapshot hyperspectral image analysis of age-related macular degeneration.

Drusen, the hallmark lesions of age related macular degeneration (AMD), are biochemically heterogeneous and the identification of their biochemical distribution is key to the understanding of AMD. Yet the challenges are to develop imaging technology and analytics, which respect the physical generation of the hyperspectral signal in the presence of noise, artifacts, and multiple mixed sources while maximally exploiting the full data dimensionality to uncover clinically relevant spectral signatures. This paper reports on the statistical analysis of hyperspectral signatures of drusen and anatomical regions of interest using snapshot hyperspectral imaging and non-negative matrix factorization (NMF). We propose physical meaningful priors as initialization schemes to NMF for finding low-rank decompositions that capture the underlying physiology of drusen and the macular pigment. Preliminary results show that snapshot hyperspectral imaging in combination with NMF is able to detect biochemically meaningful components of drusen and the macular pigment. To our knowledge, this is the first reported demonstration in vivo of the separate absorbance peaks for lutein and zeaxanthin in macular pigment.

Using perturbation theory to compute the morphological similarity of diffusion tensors.

Computing the morphological similarity of diffusion tensors (DTs) at neighboring voxels within a DT image, or at corresponding locations across different DT images, is a fundamental and ubiquitous operation in the postprocessing of DT images. The morphological similarity of DTs typically has been computed using either the principal directions (PDs) of DTs (i.e., the direction along which water molecules diffuse preferentially) or their tensor elements. Although comparing PDs allows the similarity of one morphological feature of DTs to be visualized directly in eigenspace, this method takes into account only a single eigenvector, and it is therefore sensitive to the presence of noise in the images that can introduce error intothe estimation of that vector. Although comparing tensor elements, rather than PDs, is comparatively more robust to the effects of noise, the individual elements of a given tensor do not directly reflect the diffusion properties of water molecules. We propose a measure for computing the morphological similarity of DTs that uses both their eigenvalues and eigenvectors, and that also accounts for the noise levels present in DT images. Our measure presupposes that DTs in a homogeneous region within or across DT images are random perturbations of one another in the presence of noise. The similarity values that are computed using our method are smooth (in the sense that small changes in eigenvalues and eigenvectors cause only small changes in similarity), and they are symmetric when differences in eigenvalues and eigenvectors are also symmetric. In addition, our method does not presuppose that the corresponding eigenvectors across two DTs have been identified accurately, an assumption that is problematic in the presence of noise. Because we compute the similarity between DTs using their eigenspace components, our similarity measure relates directly to both the magnitude and the direction of the diffusion of water molecules. The favorable performance characteristics of our measure offer the prospect of substantially improving additional postprocessing operations that are commonly performed on DTI datasets, such as image segmentation, fiber tracking, noise filtering, and spatial normalization.

LV volume quantification via spatiotemporal analysis of real-time 3-D echocardiography.

This paper presents a method of four-dimensional (4-D) (3-D + Time) space-frequency analysis for directional denoising and enhancement of real-time three-dimensional (RT3D) ultrasound and quantitative measures in diagnostic cardiac ultrasound. Expansion of echocardiographic volumes is performed with complex exponential wavelet-like basis functions called brushlets. These functions offer good localization in time and frequency and decompose a signal into distinct patterns of oriented harmonics, which are invariant to intensity and contrast range. Deformable-model segmentation is carried out on denoised data after thresholding of transform coefficients. This process attenuates speckle noise while preserving cardiac structure location. The superiority of 4-D over 3-D analysis for decorrelating additive white noise and multiplicative speckle noise on a 4-D phantom volume expanding in time is demonstrated. Quantitative validation, computed for contours and volumes, is performed on in vitro balloon phantoms. Clinical applications of this spaciotemporal analysis tool are reported for six patient cases providing measures of left ventricular volumes and ejection fraction.

Improving PET-based physiological quantification through methods of wavelet denoising.

The goal of this study was to evaluate methods of multidimensional wavelet denoising on restoring the fidelity of biological signals hidden within dynamic positron emission tomography (PET) images. A reduction of noise within pixels, between adjacent regions, and time-serial frames was achieved via redundant multiscale representations. In analyzing dynamic PET data of healthy volunteers, a multiscale method improved the estimate-to-error ratio of flows fivefold without loss of detail. This technique also maintained accuracy of flow estimates in comparison with the "gold standard," using dynamic PET with O15-water. In addition, in studies of coronary disease patients, flow patterns were preserved and infarcted regions were well differentiated from normal regions. The results show that a wavelet-based noise-suppression method produced reliable approximations of salient underlying signals and led to an accurate quantification of myocardial perfusion. The described protocol can be generalized to other temporal biomedical imaging modalities including functional magnetic resonance imaging and ultrasound.

Use of wavelet transforms in analysis of time-activity data from cardiac PET.

UNLABELLED: Because of its intrinsic quantitative properties, PET permits measurement of myocardial perfusion and metabolism in absolute terms (i.e., mL/g/min). However, quantification has been limited by errors produced in image acquisition, selection of regions of interest, and data analysis. The goal of this study was to evaluate a newly developed, novel, wavelet-based noise-reduction approach that can objectively extract biologic signals hidden within dynamic PET data. METHODS: Quantification of myocardial perfusion using dynamic PET imaging with 82Rb, H2(15)O, and 13NH3 was selected to evaluate the effects of the wavelet-based noise-reduction protocol. Dynamic PET data were fitted to appropriate mathematic models before and after wavelet-based noise reduction to get flow estimates. Time-activity curves, precision, accuracy, and differentiating capacity derived from the wavelet protocol were compared with those obtained from unmodified data processing. A total of 84 human studies was analyzed, including 43 at rest (18 82Rb scans, 18 H2(15)O scans, and 7 13NH3 scans) and 41 after coronary hyperemia with dipyridamole (17 82Rb scans, 17 H2(15)O scans, and 7 13NH3 scans). RESULTS: For every tracer tested under all conditions, the wavelet method improved the shape of blood and tissue time-activity curves, increased estimate-to-error ratios, and maintained fidelity of flow in regions as small as 0.85 cm3. It also improved the accuracy of flow estimates derived from 82Rb to the level of that achieved with H2(15)O, which was not affected markedly by the wavelet process. In studies of patients with coronary disease, regional heterogeneity of myocardial perfusion was preserved and flow estimates in infarcted regions were differentiated more easily from normal regions. CONCLUSION: The wavelet-based noise-reduction method effectively and objectively extracted tracer time-activity curves from data with low signal-to-noise ratios and improved the accuracy and precision of measurements with all tracer techniques studied. The approach should be generalizable to other image modalities such as functional MRI and CT and, therefore, improve the ability to quantify dynamic physiologic processes.

Quantification of myocardial perfusion in human subjects using 82Rb and wavelet-based noise reduction.

UNLABELLED: Quantification of myocardial perfusion with 82Rb has been difficult to achieve because of the low signal-to-noise ratio of the dynamic data curves. This study evaluated the accuracy of flow estimates after the application of a novel multidimensional wavelet-based noise-reduction protocol. METHODS: Myocardial perfusion was estimated using 82Rb and a two-compartment model from dynamic PET scans on 11 healthy volunteers at rest and after hyperemic stress with dipyridamole. Midventricular planes were divided into eight regions of interest, and a wavelet transform protocol was applied to images and time-activity curves. Flow estimates without and with the wavelet approach were compared with those obtained using H2(15)O. RESULTS: Over a wide flow range (0.45-2.75 mL/g/min), flow achieved with the wavelet approach correlated extremely closely with values obtained with H2(15)O (y = 1.03 x -0.12; n = 23 studies, r = 0.94, P < 0.001). If the wavelet noise-reduction technique was not used, the correlation was less strong (y = 1.11 x + 0.24; n = 23 studies, r = 0.79, P < 0.001). In addition, the wavelet approach reduced the regional variation from 75% to 12% and from 62% to 11% (P < 0.001 for each comparison) for resting and stress studies, respectively. CONCLUSION: The use of a wavelet protocol allows near-optimal noise reduction, markedly enhances the physiologic flow signal within the PET images, and enables accurate measurement of myocardial perfusion with 82Rb in human subjects over a wide range of flows.

Wavelets in temporal and spatial processing of biomedical images.

We review some of the most recent advances in the area of wavelet applications in medical imaging. We first review key concepts in the processing of medical images with wavelet transforms and multiscale analysis, including time-frequency tiling, overcomplete representations, higher dimensional bases, symmetry, boundary effects, translational invariance, orientation selectivity, and best-basis selection. We next describe some applications in magnetic resonance imaging, including activation detection and denoising of functional magnetic resonance imaging and encoding schemes. We then present an overview in the area of ultrasound, including computational anatomy with three-dimensional cardiac ultrasound. Next, wavelets in tomography are reviewed, including their relationship to the radon transform and applications in position emission tomography imaging. Finally, wavelet applications in digital mammography are reviewed, including computer-assisted diagnostic systems that support the detection and classification of small masses and methods of contrast enhancement.

Speckle reduction and contrast enhancement of echocardiograms via multiscale nonlinear processing.

This paper presents an algorithm for speckle reduction and contrast enhancement of echocardiographic images. Within a framework of multiscale wavelet analysis, we apply wavelet shrinkage techniques to eliminate noise while preserving the sharpness of salient features. In addition, nonlinear processing of feature energy is carried out to enhance contrast within local structures and along object boundaries. We show that the algorithm is capable of not only reducing speckle, but also enhancing features of diagnostic importance, such as myocardial walls in two-dimensional echocardiograms obtained from the parasternal short-axis view. Shrinkage of wavelet coefficients via soft thresholding within finer levels of scale is carried out on coefficients of logarithmically transformed echocardiograms. Enhancement of echocardiographic features is accomplished via nonlinear stretching followed by hard thresholding of wavelet coefficients within selected (midrange) spatial-frequency levels of analysis. We formulate the denoising and enhancement problem, introduce a class of dyadic wavelets, and describe our implementation of a dyadic wavelet transform. Our approach for speckle reduction and contrast enhancement was shown to be less affected by pseudo-Gibbs phenomena. We show experimentally that this technique produced superior results both qualitatively and quantitatively when compared to results obtained from existing denoising methods alone. A study using a database of clinical echocardiographic images suggests that such denoising and enhancement may improve the overall consistency of expert observers to manually defined borders.

A wavelet-based metric for visual texture discrimination with applications in evolutionary ecology.

Much work on natural and sexual selection is concerned with the conspicuousness of visual patterns (textures) on animal and plant surfaces. Previous attempts by evolutionary biologists to quantify apparency of such textures have involved subjective estimates of conspicuousness or statistical analyses based on transect samples. We present a method based on wavelet analysis that avoids subjectivity and that uses more of the information in image textures than transects do. Like the human visual system for texture discrimination, and probably like that of other vertebrates, this method is based on localized analysis of orientation and frequency components of the patterns composing visual textures. As examples of the metric's utility, we present analyses of crypsis for tigers, zebras, and peppered moth morphs.

Mammographic feature enhancement by multiscale analysis.

Introduces a novel approach for accomplishing mammographic feature analysis by overcomplete multiresolution representations. The authors show that efficient representations may be identified within a continuum of scale-space and used to enhance features of importance to mammography. Methods of contrast enhancement are described based on three overcomplete multiscale representations: 1) the dyadic wavelet transform (separable), 2) the phi-transform (nonseparable, nonorthogonal), and 3) the hexagonal wavelet transform (nonseparable). Multiscale edges identified within distinct levels of transform space provide local support for image enhancement. Mammograms are reconstructed from wavelet coefficients modified at one or more levels by local and global nonlinear operators. In each case, edges and gain parameters are identified adaptively by a measure of energy within each level of scale-space. The authors show quantitatively that transform coefficients, modified by adaptive nonlinear operators, can make more obvious unseen or barely seen features of mammography without requiring additional radiation. The authors' results are compared with traditional image enhancement techniques by measuring the local contrast of known mammographic features. They demonstrate that features extracted from multiresolution representations can provide an adaptive mechanism for accomplishing local contrast enhancement. By improving the visualization of breast pathology, one can improve chances of early detection while requiring less time to evaluate mammograms for most patients.

Visual textures, machine vision and animal camouflage.

The psychophysics of visual texture perception and texture discrimination have been investigated extensively during the past 30 years. Humans have been the main study subjects, but some research on texture perception has involved other species, and there is good reason to think that the most general results from humans apply to other vertebrates as well. Psychophysicists have suggested that some of their findings on human vision reflect adaptive 'tricks' for countering prey camouflage, but this possibility has not been widely communicated to evolutionary biologists. We review the psychophysicists' main conclusions on texture discrimination, and list additional questions that their results raise when animal coats are considered as visual textures. We also suggest ways in which advances in computer vision can be combined with psychophysics to provide new perspectives on the function of animal coat patterns.