An improved seeded region growing algorithm

Recently Adams and Bischof (1994) proposed a novel region growing algorithm for segmenting intensity images. The inputs to the algorithm are the intensity image and a set of seeds - individual points or connected components - that identify the individual regions to be segmented. The algorithm grows these seed regions until all of the image pixels have been assimilated. Unfortunately the algorithm is inherently dependent on the order of pixel processing. This means, for example, that raster order processing and anti-raster order processing do not, in general, lead to the same tessellation. In this paper we propose an improved seeded region growing algorithm that retains the advantages of the Adams and Bischof algorithm fast execution, robust segmentation, and no tuning parameters - but is pixel order independent. (C) 1997 Elsevier Science B.V.

Prediction of MHC class II-binding peptides using an evolutionary algorithm and artificial neural network

Motivation: Prediction methods for identifying binding peptides could minimize the number of peptides required to be synthesized and assayed, and thereby facilitate the identification of potential T-cell epitopes. We developed a bioinformatic method for the prediction of peptide binding to MHC class II molecules. Results: Experimental binding data and expert knowledge of anchor positions and binding motifs were combined with an evolutionary algorithm (EA) and an artificial neural network (ANN): binding data extraction --> peptide alignment --> ANN training and classification. This method, termed PERUN, was implemented for the prediction of peptides that bind to HLA-DR4(B1*0401). The respective positive predictive values of PERUN predictions of high-, moderate-, low- and zero-affinity binder-a were assessed as 0.8, 0.7, 0.5 and 0.8 by cross-validation, and 1.0, 0.8, 0.3 and 0.7 by experimental binding. This illustrates the synergy between experimentation and computer modeling, and its application to the identification of potential immunotheraaeutic peptides.

Relationships between blood pressure and heart rate in the saltwater crocodile Crocodylus porosus

The cardiac limb of the baroreflex loop was studied in the saltwater crocodile Crocodylus porosus, The classical pharmacological methodology using phenylephrine and sodium nitroprusside was used to trigger blood pressure changes, and the resulting alterations in heart rate were analysed quantitatively using a logistic function. Interindividual differences in resting heart rates and blood pressures were observed, but all seven animals displayed clear baroreflex responses. Atropine and sotalol greatly attenuated the response. A maximal baroreflex gain of 7.2 beats min(-1) kPa(-1) was found at a mean aortic pressure of 6.1 kPa, indicating the active role of the baroreflex in a wide pressure range encompassing hypotensive and hypertensive states. At the lowest mean aortic pressures (5.0 kPa), the synergistic role of the pulmonary-to-systemic shunt in buffering the blood pressure drop also contributes to blood pressure regulation, Pulse pressure showed a better correlation,vith heart rate and also a higher gain than mean aortic, systolic or diastolic pressures, and this is taken as an indicator of the existence of a differential control element working simultaneously with a linear proportional element.

Heart rate, blood lactate and kinematic data of elite colts (under-19) rugby union players during competition

Physiological and kinematic data were collected from elite under-19 rugby union players to provide a greater understanding of the physical demands of rugby union. Heart rate, blood lactate and time-motion analysis data were collected from 24 players (mean +/- s((x) over bar): body mass 88.7 +/- 9.9 kg, height 185 +/- 7 cm, age 18.4 +/- 0.5 years) during six competitive premiership fixtures. Six players were chosen at random from each of four groups: props and locks, back row forwards, inside backs, outside backs. Heart rate records were classified based on percent time spent in four zones (>95%, 85-95%, 75-84%, <75% HRmax). Blood lactate concentration was measured periodically throughout each match, with movements being classified as standing, walking, jogging, cruising, sprinting, utility, rucking/mauling and scrummaging. The heart rate data indicated that props and locks (58.4%) and back row forwards (56.2%) spent significantly more time in high exertion (85-95% HRmax) than inside backs (40.5%) and outside backs (33.9%) (P < 0.001). Inside backs (36.5%) and outside backs (38.5%) spent significantly more time in moderate exertion (75-84% HRmax) than props and locks (22.6%) and back row forwards (19.8%) (P < 0.05). Outside backs (20.1%) spent significantly more time in low exertion (< 75% HRmax) than props and locks (5.8%) and back row forwards (5.6%) (P < 0.05). Mean blood lactate concentration did not differ significantly between groups (range: 4.67 mmol.l(-1) for outside backs to 7.22 mmol.l(-1) for back row forwards; P < 0.05). The motion analysis data indicated that outside backs (5750 m) covered a significantly greater total distance than either props and locks or back row forwards (4400 and 4080 m, respectively; P < 0.05). Inside backs and outside backs covered significantly greater distances walking (1740 and 1780 m, respectively; P < 0.001), in utility movements (417 and 475 m, respectively; P < 0.001) and sprinting (208 and 340 m, respectively; P < 0.001) than either props and locks or back row forwards (walking: 1000 and 991 m; utility movements: 106 and 154 m; sprinting: 72 and 94 m, respectively). Outside backs covered a significantly greater distance sprinting than inside backs (208 and 340 m, respectively; P < 0.001). Forwards maintained a higher level of exertion than backs, due to more constant motion and a large involvement in static high-intensity activities. A mean blood lactate concentration of 4.8-7.2 mmol.l(-1) indicated a need for 'lactate tolerance' training to improve hydrogen ion buffering and facilitate removal following high-intensity efforts. Furthermore, the large distances (4.2-5.6 km) covered during, and intermittent nature of, match-play indicated a need for sound aerobic conditioning in all groups (particularly backs) to minimize fatigue and facilitate recovery between high-intensity efforts.

A consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements

To translate and transfer solution data between two totally different meshes (i.e. mesh 1 and mesh 2), a consistent point-searching algorithm for solution interpolation in unstructured meshes consisting of 4-node bilinear quadrilateral elements is presented in this paper. The proposed algorithm has the following significant advantages: (1) The use of a point-searching strategy allows a point in one mesh to be accurately related to an element (containing this point) in another mesh. Thus, to translate/transfer the solution of any particular point from mesh 2 td mesh 1, only one element in mesh 2 needs to be inversely mapped. This certainly minimizes the number of elements, to which the inverse mapping is applied. In this regard, the present algorithm is very effective and efficient. (2) Analytical solutions to the local co ordinates of any point in a four-node quadrilateral element, which are derived in a rigorous mathematical manner in the context of this paper, make it possible to carry out an inverse mapping process very effectively and efficiently. (3) The use of consistent interpolation enables the interpolated solution to be compatible with an original solution and, therefore guarantees the interpolated solution of extremely high accuracy. After the mathematical formulations of the algorithm are presented, the algorithm is tested and validated through a challenging problem. The related results from the test problem have demonstrated the generality, accuracy, effectiveness, efficiency and robustness of the proposed consistent point-searching algorithm. Copyright (C) 1999 John Wiley & Sons, Ltd.

The crack tip strain field of AISI 4340 - Part I - Measurement technique

This is the first paper in a study on the influence of the environment on the crack tip strain field for AISI 4340. A stressing stage for the environmental scanning electron microscope (ESEM) was constructed which was capable of applying loads up to 60 kN to fracture-mechanics samples. The measurement of the crack tip strain field required preparation (by electron lithography or chemical etching) of a system of reference points spaced at similar to 5 mu m intervals on the sample surface, loading the sample inside an electron microscope, image processing procedures to measure the displacement at each reference point and calculation of the strain field. Two algorithms to calculate strain were evaluated. Possible sources of errors were calculation errors due to the algorithm, errors inherent in the image processing procedure and errors due to the limited precision of the displacement measurements. Estimation of the contribution of each source of error was performed. The technique allows measurement of the crack tip strain field over an area of 50 x 40 mu m with a strain precision better than +/- 0.02 at distances larger than 5 mu m from the crack tip. (C) 1999 Kluwer Academic Publishers.

Increased ovulation rate in Merino ewes immunized against small synthetic peptide fragments of the inhibin alpha subunit

Four experiments were carried out in Merino ewes during a period of 4 years to determine the long-term effects of immunization against different synthetic peptides mimicking the amine terminal of the or subunit of porcine inhibin. Peptides were conjugated to human serum albumin and 100-200 mu g emulsified in Freund's complete adjuvant for the primary immunization. Usually two booster injections were given at monthly intervals with 50-100 mu g conjugated peptide using either incomplete Freund's adjuvant or Montanide : Marcel. In some experiments a further immunization was carried in the next year. Blood samples were taken 10 days after each immunization, during the luteal phase, for estimation of gonadotrophin concentrations and determination of inhibin antibody titres. One day after blood sampling cloprostenol was used to induce luteolysis and laparoscopy was performed in the subsequent oestrous cycle. Immunization of ewes with synthetic peptides 1-32, 1-26, 7-26 and 8-30 resulted in large increases in the ovulation rate (OR). An approximately two-fold increase in OR was observed following the first booster immunization with these peptides and a three- to five-fold increase after the second booster immunization. Immunization with these large peptides resulted in a sustained increase in OR for a period of at least 1 year after the second booster immunization. Of the shorter peptides, peptides 10-26 and 13-26 gave a reasonable ovulatory response, although it was more difficult to obtain a response with peptides 1-16, 8-22, 13-25, 8-19 and 10-19; peptides 7-13 and 1-6 gave no response (but were examined for one breeding season only). The smaller peptides led to lower inhibin antibody titres that were not necessarily associated with increased follicle-stimulating hormone (FSH) or OR. More intensive blood sampling in one experiment showed that following primary immunization against peptide 1-32 there was a transient increase in plasma FSH which did not lead to an increased OR. Moreover, a prolonged period of raised FSH after the first booster was significantly correlated with increased OR. In these animals antibody titres were only slightly increased after primary immunization, but after the first booster immunization higher titres were observed that were significantly correlated with trough FSH values and the subsequent OR. These results are interpreted as showing that (1) to obtain an increase in OR peptides 1-32, 1-26 and 7-26 are suitable as immunogens; (2) smaller peptides are less reliable, often require multiple injections, and the response may be delayed; and (3) an extended period of raised plasma FSH is needed to give a large ovulatory response.

Between-year fluctuations in schizophrenia birth rate and associations with the Southern Oscillation Index, cloud cover and sunshine

Our group have recently proposed that low prenatal vitamin D may be a risk-modifying factor for schizophrenia. Climate variability impacts on vitamin D levels in a population via fluctuations in the amount of available UV radiation. In order to explore this hypothesis, we examined fluctuations in the birthrates for people with schizophrenia born between 1920 and 1967 and three sets of variables strongly associated with UV radiation. These included: (a) the Southern Oscillation Index (SOI), a marker of El Nino which is the most prominent meteorological factor that influences Queensland weather: (b) measures of cloud cover and (c) measures of sunshine. Schizophrenia births were extracted from the Queensland Mental Health register and corrected for background population birth rates. Schizophrenia birth rates had several apparently non-random features in common with the SO1. The prominent SO1 fluctuation event that occurred between 1937 and 1943 is congruent with the most prominent fluctuation in schizophrenia birth rates. The relatively flat profile of SOI activity between 1927 and 1936 also corresponds to the flattest period in the schizophrenia time series. Both time series have prominent oscillations in the 3 ~, year range between 1946 and 1960. Significant associations between schizophrenia birth rates and measures of both sunshine and cloud cover were identified,and all three time series shared periodicity in the 3-4 year range. The analyses suggest that the risk of schizophrenia is higher for those born during times of increased cloud cover,reduced sunshine and positive SO1. These ecological analyses provide initial support for the vitamin D hypothesis, however alternative non-genetic candidate exposures also need to be considered. Other sites with year-to-year fluctuations in cloud cover and sunshine should examine patterns of association between these climate variables and schizophrenia birth rates. The Stanley Foundation supported this project.

Trial-of-antibiotic algorithm for the diagnosis of tuberculosis in a district hospital in a developing country with high HIV prevalence

OBJECTIVE: To evaluate a diagnostic algorithm for pulmonary tuberculosis based on smear microscopy and objective response to trial of antibiotics. SETTING: Adult medical wards, Hlabisa Hospital, South Africa, 1996-1997. METHODS: Adults with chronic chest symptoms and abnormal chest X-ray had sputum examined for Ziehl-Neelsen stained acid-fast bacilli by light microscopy. Those with negative smears were treated with amoxycillin for 5 days and assessed. Those who had not improved were treated with erythromycin for 5 days and reassessed. Response was compared with mycobacterial culture. RESULTS: Of 280 suspects who completed the diagnostic pathway, 160 (57%) had a positive smear, 46 (17%) responded to amoxycillin, 34 (12%) responded to erythromycin and 40 (14%) were treated as smear-negative tuberculosis. The sensitivity (89%) and specificity (84%) of the full algorithm for culture-positive tuberculosis were high. However, 11 patients (positive predictive value [PPV] 95%) were incorrectly diagnosed with tuberculosis, and 24 cases of tuberculosis (negative predictive value [NPV] 70%) were not identified. NPV improved to 75% when anaemia was included as a predictor. Algorithm performance was independent of human immunodeficiency virus status. CONCLUSION: Sputum smear microscopy plus trial of antibiotic algorithm among a selected group of tuberculosis suspects may increase diagnostic accuracy in district hospitals in developing countries.

The effect of solidification rate on the structure of magnesium-aluminium eutectic grains

A magnesium-aluminium alloy of eutectic composition was solidified under two different cooling conditions, producing a low and a high growth rate of the eutectic solid-liquid interface. The high growth rate specimen contained smaller eutectic grains and cells, with a smaller interphase spacing compared with the low growth rate specimen. The high growth rate specimen also contained some primary Mg17Al12 dendrites, suggesting that the coupled zone is skewed towards the Mg phase with increased undercooling, A lamellar eutectic morphology was observed in the low growth rate specimen, while the morphology was fibrous in the high growth rate specimen.

The interaction between the equity premium and the risk-free rate

The equity premium arises from the interaction between the atemporal risk premium for equity, the risk-free rate of intertemporal substitution and the impact of risk on the precautionary motive for saving. Depending on parameter values, the equity premium may either be increased or reduced by the presence of undiversifiable background risk. (C) 2000 Elsevier Science S.A. All rights reserved.

Exposure of rats to a high but not low dose of ethanol during early postnatal life increases the rate of loss of optic nerve axons and decreases the rate of myelination

Visual system abnormalities are commonly encountered in the fetal alcohol syndrome although the level of exposure at which they become manifest is uncertain. In this study we have examined the effects of either low (ETLD) or high dose (ETHD) ethanol, given between postnatal days 4-9, on the axons of the rat optic nerve. Rats were exposed to ethanol vapour in a special chamber for a period of 3 h per day during the treatment period. The blood alcohol concentration in the ETLD animals averaged similar to 171 mg/dl and in the ETHD animals similar to 430 mg/dl at the end of the treatment on any given day. Groups of 10 and 30-d-old mother-reared control (MRC), separation control (SC), ETLD and ETHD rats were anaesthetised with an intraperitoneal injection or ketamine and xylazine, and killed by intracardiac perfusion with phosphate-buffered glutaraldehyde. In the 10-d-old rat optic nerves there was a total of similar to 145000-165000 axons in MRC, SC and ETLD animals. About 4 % of these fibres were myelinated. The differences between these groups were not statistically significant. However, the 10-d-old ETHD animals had only about 75000 optic nerve axone (P < 0.05) of which about 2.8 % were myelinated. By 30 d of age there was a total of between 75000 90000 optic nerve axons, irrespective of the group examined. The proportion of axons which were myelinated at this age was still significantly lower (P < 0.001) in the ETHD animals (similar to 77 %) than in the other groups (about 98 %). It is concluded that the normal stages of development and maturation of the rat optic nerve axons, as assessed in this study, can be severely compromised by exposure to a relatively high (but not low) dose of ethanol between postnatal d 4 and 9.

Finite element modeling of fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins

In order to use the finite element method for solving fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins effectively and efficiently, we have presented, in this paper, the new concept and numerical algorithms to deal with the fundamental issues associated with the fluid-rock interaction problems. These fundamental issues are often overlooked by some purely numerical modelers. (1) Since the fluid-rock interaction problem involves heterogeneous chemical reactions between reactive aqueous chemical species in the pore-fluid and solid minerals in the rock masses, it is necessary to develop the new concept of the generalized concentration of a solid mineral, so that two types of reactive mass transport equations, namely, the conventional mass transport equation for the aqueous chemical species in the pore-fluid and the degenerated mass transport equation for the solid minerals in the rock mass, can be solved simultaneously in computation. (2) Since the reaction area between the pore-fluid and mineral surfaces is basically a function of the generalized concentration of the solid mineral, there is a definite need to appropriately consider the dependence of the dissolution rate of a dissolving mineral on its generalized concentration in the numerical analysis. (3) Considering the direct consequence of the porosity evolution with time in the transient analysis of fluid-rock interaction problems; we have proposed the term splitting algorithm and the concept of the equivalent source/sink terms in mass transport equations so that the problem of variable mesh Peclet number and Courant number has been successfully converted into the problem of constant mesh Peclet and Courant numbers. The numerical results from an application example have demonstrated the usefulness of the proposed concepts and the robustness of the proposed numerical algorithms in dealing with fluid-rock interaction problems in pore-fluid saturated hydrothermal/sedimentary basins. (C) 2001 Elsevier Science B.V. All rights reserved.

Minimum-order stable recursive filter design via the genetic algorithm approach

In this paper, the minimum-order stable recursive filter design problem is proposed and investigated. This problem is playing an important role in pipeline implementation sin signal processing. Here, the existence of a high-order stable recursive filter is proved theoretically, in which the upper bound for the highest order of stable filters is given. Then the minimum-order stable linear predictor is obtained via solving an optimization problem. In this paper, the popular genetic algorithm approach is adopted since it is a heuristic probabilistic optimization technique and has been widely used in engineering designs. Finally, an illustrative example is sued to show the effectiveness of the proposed algorithm.

Identifying rate-limiting nodes in large-scale cortical networks for visuospatial processing: an illustration using fMRI

With the advent of functional neuroimaging techniques, in particular functional magnetic resonance imaging (fMRI), we have gained greater insight into the neural correlates of visuospatial function. However, it may not always be easy to identify the cerebral regions most specifically associated with performance on a given task. One approach is to examine the quantitative relationships between regional activation and behavioral performance measures. In the present study, we investigated the functional neuroanatomy of two different visuospatial processing tasks, judgement of line orientation and mental rotation. Twenty-four normal participants were scanned with fMRI using blocked periodic designs for experimental task presentation. Accuracy and reaction time (RT) to each trial of both activation and baseline conditions in each experiment was recorded. Both experiments activated dorsal and ventral visual cortical areas as well as dorsolateral prefrontal cortex. More regionally specific associations with task performance were identified by estimating the association between (sinusoidal) power of functional response and mean RT to the activation condition; a permutation test based on spatial statistics was used for inference. There was significant behavioral-physiological association in right ventral extrastriate cortex for the line orientation task and in bilateral (predominantly right) superior parietal lobule for the mental rotation task. Comparable associations were not found between power of response and RT to the baseline conditions of the tasks. These data suggest that one region in a neurocognitive network may be most strongly associated with behavioral performance and this may be regarded as the computationally least efficient or rate-limiting node of the network.