Search

The article presents a study that examines the role of automatic neural control on carbon dioxide reactivity and the dynamic pressure autoregulation in the frequency range between 0.005 and 0.40 hertz. The study analyzed the experimental data from healthy persons with the use of a nonlinear multivariate modeling approach. Results show that mean artevial blood pressure (MABP) and the mean cerebral blood flow velocity (MCBFV) decrease during ganglion blockade.

One of the fundamental principles of cortical brain regions, including the hippocampus, is that information is represented in the ensemble firing of populations of neurons, i.e., spatio-temporal patterns of electrophysiological activity. The hippocampus has long been known to be responsible for the formation of declarative, or fact-based, memories. Damage to the hippocampus disrupts the propagation of spatio-temporal patterns of activity through hippocampal internal circuitry, resulting in a severe anterograde amnesia. Developing a neural prosthesis for the damaged hippocampus requires restoring this multiple-input, multiple-output transformation of spatio-temporal patterns of activity. Because the, mechanisms underlying synaptic transmission and generation of electrical activity in neurons are inherently nonlinear, any such prosthesis must be based on a nonlinear multiple-input, multiple-output model. In this paper, we have formulated the transformational process of multi-site propagation of spike activity between two subregions of the hippocampus (CA3 and CA1) as the identification of a multiple-input, multiple-output (MIMO) system, and proposed that it can be decomposed into a series of multiple-input, single-output (MISO) systems. Each MISO system is modeled as a physiologically plausible structure that consists of 1) linear/nonlinear feedforward Volterra kernels modeling synaptic transmission and dendritic integration, 2) a linear feedback Volterra kernel modeling spike-triggered after-potentials, 3) a threshold for spike generation, 4) a summation process for somatic integration, and 5) a noise term representing intrinsic neuronal noise and the contributions of unobserved inputs. Input and output spike trains were recorded from hippocampal CA3 and CA1 regions of rats performing a spatial delayed-nonmatch-to-sample memory task that requires normal hippocampal function. Kernels were expanded with Laguerre basis functions and estimated using a maximum-likelihood method. Complexity of the feedforward kernel was progressively increased to capture higher-order system nonlinear dynamics. Results showed higher prediction accuracies as kernel complexity increased. Self-kernels describe the nonlinearities within each input. Cross-kernels capture the nonlinear interaction between inputs. Second- and third-order nonlinear models were found to successfully predict the CA1 output spike distribution based on CA3 input spike trains. First-order, linear models were shown to be insufficient. [ABSTRACT FROM AUTHOR]

Reports on the development of a hippocampal-cortical neural prostheses by researchers from the U.S. Use of the novel medical apparatus to restore lost cognitive function; Safety and efficacy of replacing the CA3 region of the hippocampal slice with the biomimetic device; Characterization and mathematical modeling of the nonlinear dynamic properties of the region.

Inversion layer ultrasonic transducers have been investigated recently as an interesting approach in wideband transducer design. In this paper we present an analytical model of multilayer ultrasonic transducers with an inversion layer. Our analysis of the wave propagation problem of an inversion layer transducer includes a functional decomposition of the electrical input impedance. It becomes clear from this decomposition that an inversion layer transducer can be modeled as three elements in series connection, i.e., a clamped capacitance, a classical motional impedance, and a coupled motional impedance. The first two elements make up the classical model of a single element transducer. The coupled motional impedance describes the coupled interaction between the regular and the inverted piezoelectric sublayers, and thus reflects the effect of an inversion layer. We present examples which show that inversion layer transducers are advantageous in achieving such useful features as dual-frequency operation mode as used in harmonic imaging or broadband performance desired in most ultrasonic applications. [ABSTRACT FROM AUTHOR]

Half-thickness inversion layer high-frequency ultrasonic transducers were fabricated using lithium niobate (LiNbO3) single crystal plate. The transducers developed for this study used a 36° rotated Y-cut LiNbO3 thin plate with an active element thickness of 115 μm. The designed center frequency was in the range of 30 to 60 MHz. Half-thickness inversion layer was formed after the sample was annealed at a high temperature, and it is shown that the inversion layer thickness can be controlled by the temperature. Silver powder/epoxy composite and parylene were used as acoustic matching layers. A lossy silver epoxy was used as the backing material. Using an analytical method, the electrical impedance for different inversion layer ratios was determined. The measured resonant frequency was consistent with the modeled data. Even-order higher frequency broadband ultrasonic transducers with a center frequency at 60 MHz was obtained using half-thickness inversion layer of LiNbO3 single crystal. [ABSTRACT FROM AUTHOR]

The effect of spontaneous beat-to-beat mean arterial blood pressure fluctuations and breath-to-breath end-tidal CO2 fluctuations on beat-to-beat cerebral blood flow velocity variations is studied using the Laguerre-Volterra network methodology for multiple-input nonlinear systems. The observations made from experimental measurements from ten healthy human subjects reveal that, whereas pressure fluctuations explain most of the high-frequency blood flow velocity variations (above 0.04 Hz), end-tidal CO2 fluctuations as well as nonlinear interactions between pressure and CO2 have a considerable effect in the lower frequencies (below 0.04 Hz). They also indicate that cerebral autoregulation is strongly nonlinear and dynamic (frequency-dependent). Nonlinearities are mainly active in the low-frequency range (below 0.04 Hz) and are more prominent in the dynamics of the end-tidal CO2-blood flow velocity relationship. Significant nonstationarities are also revealed by the obtained models, with greater variability evident for the effects of CO2 on blood flow velocity dynamics. [ABSTRACT FROM AUTHOR]

Proposes the use of a class of feedforward artificial neural networks with polynomial activation functions for practical modeling of highorder Volterra systems. Relationship between Volterra models and feedforward multilayer neural networks; Structure of discrete-time Volterra models; Examples of Volterra systems modeling.

Presents information on the comparison of the characteristics in renal autoregulation in normotensive Sprague-Dawey and spontaneously hypertensive rats, relating to the application of biomedical engineering and applied biophysics. Information on mathematical methods; Other detailed information.

Assesses the use of feedforward artificial neural networks (FANN) to model the nonlinear dynamics of the renal autoregulatory mechanisms in rats and compares their performance in terms of model prediction with the Volterra models. Details on the existence of nonlinear interactions between the myogenic and tubuloglomerular feedback (TGF) mechanism which regulate renal blood flow in response to an increase in blood pressure; Findings of studies conducted.

Objective: To compare the novel model‐based hemodynamic physiomarker of Dynamic Vasomotor Reactivity (DVR) with biomarkers based on Diffusion Tensor Imaging (DTI) and some widely used neurocognitive scores in terms of their ability to delineate patients with amnestic Mild Cognitive Impairment (MCI) from age‐matched cognitively normal controls. Materials & Methods: The model‐based DVR and MRI‐based DTI markers were obtained from 36 patients with amnestic MCI and 16 age‐matched controls without cognitive impairment, for whom widely used neurocognitive scores were available. These markers and scores were subsequently compared in terms of statistical delineation between patients and controls. Results: It was found that statistically significant delineation between MCI patients and controls was comparable for DVR or DTI markers (p < 0.01). The performance of both types of markers was consistent with the scores of some (but not all) widely used neurocognitive tests. Conclusion: Since DTI offers a measure of cerebral white matter integrity, the results suggest that the model‐based hemodynamic marker of DVR may correlate with cognitive impairment due to white matter lesions. This finding is consistent with the hypothesis that dysregulation of cerebral microcirculation may be an early cause of cognitive impairment, which has been recently corroborated by several studies. [ABSTRACT FROM AUTHOR]

Heart rate variability is a robust biomarker of emotional well-being, consistent with the shared brain networks regulating emotion regulation and heart rate. While high heart rate oscillatory activity clearly indicates healthy regulatory brain systems, can increasing this oscillatory activity also enhance brain function? To test this possibility, we randomly assigned 106 young adult participants to one of two 5-week interventions involving daily biofeedback that either increased heart rate oscillations (Osc+ condition) or had little effect on heart rate oscillations (Osc− condition) and examined effects on brain activity during rest and during regulating emotion. While there were no significant changes in the right amygdala-medial prefrontal cortex (MPFC) functional connectivity (our primary outcome), the Osc+ intervention increased left amygdala-MPFC functional connectivity and functional connectivity in emotion-related resting-state networks during rest. It also increased down-regulation of activity in somatosensory brain regions during an emotion regulation task. The Osc− intervention did not have these effects. In this healthy cohort, the two conditions did not differentially affect anxiety, depression, or mood. These findings indicate that modulating heart rate oscillatory activity changes emotion network coordination in the brain. [ABSTRACT FROM AUTHOR]

Cerebral Autoregulation (CA) is an important physiological mechanism stabilizing cerebral blood flow (CBF) in response to changes in cerebral perfusion pressure (CPP). By maintaining an adequate, relatively constant supply of blood flow, CA plays a critical role in brain function. Quantifying CA under different physiological and pathological states is crucial for understanding its implications. This knowledge may serve as a foundation for informed clinical decision-making, particularly in cases where CA may become impaired. The quantification of CA functionality typically involves constructing models that capture the relationship between CPP (or arterial blood pressure) and experimental measures of CBF. Besides describing normal CA function, these models provide a means to detect possible deviations from the latter. In this context, a recent white paper from the Cerebrovascular Research Network focused on Transfer Function Analysis (TFA), which obtains frequency domain estimates of dynamic CA. In the present paper, we consider the use of time-domain techniques as an alternative approach. Due to their increased flexibility, time-domain methods enable the mitigation of measurement/physiological noise and the incorporation of nonlinearities and time variations in CA dynamics. Here, we provide practical recommendations and guidelines to support researchers and clinicians in effectively utilizing these techniques to study CA. [ABSTRACT FROM AUTHOR]

Abstract: Developing a neural prosthesis for the damaged hippocampus requires restoring the transformation of population neural activities performed by the hippocampal circuitry. To bypass a damaged region, output spike trains need to be predicted from the input spike trains and then reinstated through stimulation. We formulate a multiple-input, multiple-output (MIMO) nonlinear dynamic model for the input–output transformation of spike trains. In this approach, a MIMO model comprises a series of physiologically-plausible multiple-input, single-output (MISO) neuron models that consist of five components each: (1) feedforward Volterra kernels transforming the input spike trains into the synaptic potential, (2) a feedback kernel transforming the output spikes into the spike-triggered after-potential, (3) a noise term capturing the system uncertainty, (4) an adder generating the pre-threshold potential, and (5) a threshold function generating output spikes. It is shown that this model is equivalent to a generalized linear model with a probit link function. To reduce model complexity and avoid overfitting, statistical model selection and cross-validation methods are employed to choose the significant inputs and interactions between inputs. The model is applied successfully to the hippocampal CA3–CA1 population dynamics. Such a model can serve as a computational basis for the development of hippocampal prostheses. [Copyright &y& Elsevier]

This paper presents the results of a computational study that compares simulated compartmental (differential equation) and Volterra models of the dynamic effects of insulin on blood glucose concentration in humans. In the first approach, we employ the widely accepted "minimal model" and an augmented form of it, which incorporates the effect of insulin secretion by the pancreas, in order to represent the actual closed-loop operating conditions of the system, and in the second modeling approach, we employ the general class of Volterra-type models that are estimated from input-output data. We demonstrate both the equivalence between the two approaches analytically and the feasibility of obtaining accurate Volterra models from insulin--glucose data generated from the compartmental models. The results corroborate the proposition that it may be preferable to obtain data-driven (i.e., inductive) models in a more general and realistic operating context, without resorting to the restrictive prior assumptions and simplifications regarding model structure and/or experimental protocols (e.g., glucose tolerance tests) that are necessary for the compartmental models proposed previously. These prior assumptions may lead to results that are improperly constrained or biased by preconceived (and possibly erroneous) notions--a risk that is avoided when we let the data guide the inductive selection of the appropriate model within the general class of Volterra-type models, as our simulation results suggest. [ABSTRACT FROM AUTHOR]

We tested the influence of blood pressure variability on the reproducibility of dynamic cerebral autoregulation (DCA) estimates. Data were analyzed from the 2nd CARNet bootstrap initiative, where mean arterial blood pressure (MABP), cerebral blood flow velocity (CBFV) and end tidal CO2 were measured twice in 75 healthy subjects. DCA was analyzed by 14 different centers with a variety of different analysis methods. Intraclass Correlation (ICC) values increased significantly when subjects with low power spectral density MABP (PSD-MABP) values were removed from the analysis for all gain, phase and autoregulation index (ARI) parameters. Gain in the low frequency band (LF) had the highest ICC, followed by phase LF and gain in the very low frequency band. No significant differences were found between analysis methods for gain parameters, but for phase and ARI parameters, significant differences between the analysis methods were found. Alternatively, the Spearman-Brown prediction formula indicated that prolongation of the measurement duration up to 35 minutes may be needed to achieve good reproducibility for some DCA parameters. We conclude that poor DCA reproducibility (ICC<0.4) can improve to good (ICC > 0.6) values when cases with low PSD-MABP are removed, and probably also when measurement duration is increased. [ABSTRACT FROM AUTHOR]

This letter proposes a novel method, multi-input, multi-output neuronal mode network (MIMO-NMN), for modeling encoding dynamics and functional connectivity in neural ensembles such as the hippocampus. Compared with conventional approaches such as the Volterra-Wiener model, linear-nonlinear-cascade (LNC) model, and generalized linear model (GLM), the NMN has several advantages in terms of estimation accuracy, model interpretation, and functional connectivity analysis. We point out the limitations of current neural spike modeling methods, especially the estimation biases caused by the imbalanced class problem when the number of zeros is significantly larger than ones in the spike data. We use synthetic data to test the performance of NMN with a comparison of the traditional methods, and the results indicate the NMN approach could reduce the imbalanced class problem and achieve better predictions. Subsequently, we apply the MIMO-NMN method to analyze data from the human hippocampus. The results indicate that the MIMO-NMN method is a promising approach to modeling neural dynamics and analyzing functional connectivity of multi-neuronal data. [ABSTRACT FROM AUTHOR]

Neurostimulation is a promising therapy for abating epileptic seizures. However, it is extremely difficult to identify optimal stimulation patterns experimentally. In this study, human recordings are used to develop a functional 24 neuron network statistical model of hippocampal connectivity and dynamics. Spontaneous seizure-like activity is induced in silico in this reconstructed neuronal network. The network is then used as a testbed to design and validate a wide range of neurostimulation patterns. Commonly used periodic trains were not able to permanently abate seizures at any frequency. A simulated annealing global optimization algorithm was then used to identify an optimal stimulation pattern, which successfully abated 92% of seizures. Finally, in a fully responsive, or closed-loop, neurostimulation paradigm, the optimal stimulation successfully prevented the network from entering the seizure state. We propose that the framework presented here for algorithmically identifying patient-specific neurostimulation patterns can greatly increase the efficacy of neurostimulation devices for seizures. [ABSTRACT FROM AUTHOR]

This letter examines the results of input-output (nonparametric) modeling based on the analysis of data generated by a mechanism-based (parametric) model of CA3-CA1 neuronal connections in the hippocampus. The motivation is to obtain biological insight into the interpretation of such input-output (Volterra-equivalent) models estimated from synthetic data. The insights obtained may be subsequently used to interpretat input-output models extracted from actual experimental data. Specifically, we found that a simplified parametric model may serve as a useful tool to study the signal transformations in the hippocampal CA3-CA1 regions. Input-output modeling of model-based synthetic data show that GABAergic interneurons are responsible for regulating neuronal excitation, controlling the precision of spike timing, and maintaining network oscillations, in a manner consistent with previous studies. The input-output model obtained from real data exhibits intriguing similarities with its synthetic-data counterpart, demonstrating the importance of a dynamic resonance in the system/model response around 2 Hz to 3 Hz. Using the input-output model from real data as a guide, we may be able to amend the parametric model by incorporating more mechanisms in order to yield better-matching input-output model. The approach we present can also be applied to the study of other neural systems and pathways. [ABSTRACT FROM AUTHOR]

Much of the research on cannabinoids (CBs) has focused on their effects at the molecular and synaptic level. However, the effects of CBs on the dynamics of neural circuits remains poorly understood. This study aims to disentangle the effects of CBs on the functional dynamics of the hippocampal Schaffer collateral synapse by using data-driven nonparametric modeling. Multi-unit activity was recorded from rats doing an working memory task in control sessions and under the influence of exogenously administered tetrahydrocannabinol (THC), the primary CB found in marijuana. It was found that THC left firing rate unaltered and only slightly reduced theta oscillations. Multivariate autoregressive models, estimated from spontaneous spiking activity, were then used to describe the dynamical transformation from CA3 to CA1. They revealed that THC served to functionally isolate CA1 from CA3 by reducing feedforward excitation and theta information flow. The functional isolation was compensated by increased feedback excitation within CA1, thus leading to unaltered firing rates. Finally, both of these effects were shown to be correlated with memory impairments in the working memory task. By elucidating the circuit mechanisms of CBs, these results help close the gap in knowledge between the cellular and behavioral effects of CBs. [ABSTRACT FROM AUTHOR]

Objective: As an extension to our study comparing a putative compartmental and data-based model of linear dynamic cerebral autoregulation (CA) and CO2-vasomotor reactivity (VR), we study the CA–VR process in a nonlinear context. Methods: We use the concept of principal dynamic modes (PDM) in order to obtain a compact and more easily interpretable input–output model. This in silico study permits the use of input data with a dynamic range large enough to simulate the classic homeostatic CA and VR curves using a putative structural model of the regulatory control of the cerebral circulation. The PDM model obtained using theoretical and experimental data are compared. Results: It was found that the PDM model was able to reflect accurately both the simulated static CA and VR curves in the associated nonlinear functions (ANFs). Similar to experimental observations, the PDM model essentially separates the pressure-flow relationship into a linear component with fast dynamics and nonlinear components with slow dynamics. In addition, we found good qualitative agreement between the PDMs representing the dynamic theoretical and experimental CO2-flow relationship. Conclusion: Under the modeling assumption and in light of other experimental findings, we hypothesize that PDMs obtained from experimental data correspond with passive fluid dynamical and active regulatory mechanisms. Significance: Both hypothesis-based and data-based modeling approaches can be combined to offer some insight into the physiological basis of PDM model obtained from human experimental data. The PDM modeling approach potentially offers a practical way to quantify the status of specific regulatory mechanisms in the CA–VR process. [ABSTRACT FROM AUTHOR]

A nonparametric model of smooth muscle tension response to electrical stimulation was estimated using the Laguerre expansion technique of nonlinear system kernel estimation. The experimental data consisted of force responses of smooth muscle to energy-matched alternating single pulse and burst current stimuli. The burst stimuli led to at least a 10-fold increase in peak force in smooth muscle from Mytilus edulis, despite the constant energy constraint. A linear model did not fit the data. However, a second-order model fit the data accurately, so the higher-order models were not required to fit the data. Results showed that smooth muscle force response is not linearly related to the stimulation power. [ABSTRACT FROM AUTHOR]

Very productive collaborative investigations characterized how multineuron hippocampal ensembles recorded in nonhuman primates (NHPs) encode short-term memory necessary for successful performance in a delayed match to sample (DMS) task and utilized that information to devise a unique nonlinear multi-input multi-output (MIMO) memory prosthesis device to enhance short-term memory in real-time during task performance. Investigations have characterized how the hippocampus in primate brain encodes information in a multi-item, rule-controlled, delayed match to sample (DMS) task. The MIMO model was applied via closed loop feedback micro-current stimulation during the task via conformal electrode arrays and enhanced performance of the complex memory requirements. These findings clearly indicate detection of a means by which the hippocampus encodes information and transmits this information to other brain regions involved in memory processing. By employing the nonlinear dynamic multi-input/multi-output (MIMO) model, developed and adapted to hippocampal neural ensemble firing patterns derived from simultaneous recorded multi-neuron CA1 and CA3 activity, it was possible to extract information encoded in the Sample phase of DMS trials that was necessary for successful performance in the subsequent Match phase of the task. The extension of this MIMO model to online delivery of electrical stimulation patterns to the same recording loci that exhibited successful CA1 firing in the DMS Sample Phase provided the means to increase task performance on a trial-by-trial basis. Increased utility of the MIMO model as a memory prosthesis was exhibited by the demonstration of cumulative increases in DMS task performance with repeated MIMO stimulation over many sessions. These results, reported below in this article, provide the necessary demonstrations to further the feasibility of the MIMO model as a memory prosthesis to recover and/or enhance encoding of cognitive information in humans with memory disruptions resulting from brain injury, disease or aging. [ABSTRACT FROM AUTHOR]

A modeling approach to characterize the nonlinear dynamic transformations of the dentate gyrus of the hippocampus is presented and experimentally validated. The dentate gyrus is the first region of the hippocampus which receives and integrates sensory information via the perforant path. The perforant path is composed of two distinct pathways: 1) the lateral path and 2) the medial perforant path. The proposed approach examines and captures the short-term dynamic characteristics of these two pathways using a nonparametric, third-order Poisson-Volterra model. The nonlinear characteristics of the two pathways are represented by Poisson-Volterra kernels, which are quantitative descriptors of the nonlinear dynamic transformations. The kernels were computed with experimental data from in vitro hippocampal slices. The electrophysiological activity was measured with custom-made multielectrode arrays, which allowed selective stimulation with random impulse trains and simultaneous recordings of extracellular field potential activity. The results demonstrate that this mathematically rigorous approach is suitable for the multipathway complexity of the hippocampus and yields interpretable models that have excellent predictive capabilities. The resulting models not only accurately predict previously reported electrophysiological descriptors, such as paired pulses, but more important, can be used to predict the electrophysiological activity of dentate granule cells to arbitrary stimulation patterns at the perforant path. [ABSTRACT FROM AUTHOR]

The encoding of mechanical stimuli into action potentials in two types of spider mechanoreceptor neurons is modeled by use of the principal dynamic modes (PDM) methodology. The PDM model is equivalent to the general Wiener–Bose model and consists of a minimum set of linear dynamic filters (PDMs), followed by a multivariate static nonlinearity and a threshold function. The PDMs are obtained by performing eigen-decomposition of a matrix constructed using the first-order and second-order Volterra kernels of the system, which are estimated by means of the Laguerre expansion technique, utilizing measurements of pseudorandom mechanical stimulation (input signal) and the resulting action potentials (output signal). The static nonlinearity, which can be viewed as a measure of the probability of action potential firing as a function of the PDM output values, is computed as the locus of points of the latter that correspond to output action potentials. The performance of the model is assessed by computing receiver operating characteristic (ROC) curves, akin to the ones used in decision theory and quantified by computing the area under the ROC curve. Three PDMs are revealed by the analysis. The first PDM exhibits a high-pass characteristic, illustrating the importance of the velocity of slit displacement in the generation of action potentials at the mechanoreceptor output, while the second and third PDMs exhibit band-pass and low-pass characteristics, respectively. The corresponding three-input nonlinearity exhibits asymmetric behavior with respect to its arguments, suggesting directional dependence of the mechanoreceptor response on the mechanical stimulation and the PDM outputs, in agreement to our findings from a previous study (Ann Biomed Eng 27:391–402, 1999). Differences between the Type A and B neurons are observed in the zeroth-order Volterra kernels (related to the average firing), as well as in the magnitudes of the second and third PDMs that perform band-pass and low-pass processing of the input signal, respectively. [ABSTRACT FROM AUTHOR]

This paper presents a novel segmentation methodology for automated classification and differentiation of soft tissues using multiband data obtained with the newly developed system of high-resolution ultrasonic transmission tomography (HUTT) for imaging biological organs. This methodology extends and combines two existing approaches: the L-level set active contour (AC) segmentation approach and the agglomerative hierarchical k-means approach for unsupervised Clustering (UC). To prevent the trapping of the current iterative minimization AC algorithm in a local minimum, we introduce a multiresolution approach that applies the level set functions at successively increasing resolutions of the image data. The resulting AC clusters are subsequently rearranged by the UC algorithm that seeks the optimal set of clusters yielding the minimum within-cluster distances in the feature space. The presented results from Monte Carlo simulations and experimental animal-tissue data demonstrate that the proposed methodology outperforms other existing methods without depending on heuristic parameters and provides a reliable means for soft tissue differentiation in HUTT images. [ABSTRACT FROM AUTHOR]

In this paper, we are interested in soft tissue differentiation by multiband images obtained from the High-Resolution Ultrasonic Transmission Tomography (HUTT) system using a spectral target detection method based on constrained energy minimization (CEM). We have developed a new tissue differentiation method (called "CEM filter bank") consisting of multiple CEM filters specially designed for detecting multiple types of tissues. Statistical inference on the output of the CEM filter bank is used to make a decision based on the maximum statistical significance rather than the magnitude of each CEM filter output. We test and validate this method through three-dimensional interphantom/intraphantom soft tissue classification where target profiles obtained from an arbitrary single slice are used for differentiation over multiple other tomographic slices. The performance of the proposed classifier is assessed using receiver operating characteristic analysis. We also apply our method to classify tiny structures inside a bovine kidney and sheep kidneys. Using the proposed method we can detect physical objects and biological tissues such as styrofoam balls, chicken tissue, calyces, and vessel-duct successfully. [ABSTRACT FROM AUTHOR]