Unsupervised cell nucleus segmentation with active contours

The task of segmenting cell nuclei from cytoplasm in conventional Papanicolaou (Pap) stained cervical cell images is a classical image analysis problem which may prove to be crucial to the development of successful systems which automate the analysis of Pap smears for detection of cancer of the cervix. Although simple thresholding techniques will extract the nucleus in some cases, accurate unsupervised segmentation of very large image databases is elusive. Conventional active contour models as introduced by Kass, Witkin and Terzopoulos (1988) offer a number of advantages in this application, but suffer from the well-known drawbacks of initialisation and minimisation. Here we show that a Viterbi search-based dual active contour algorithm is able to overcome many of these problems and achieve over 99% accurate segmentation on a database of 20 130 Pap stained cell images. (C) 1998 Elsevier Science B.V. All rights reserved.

X-ray computed tomography to quantify tree rooting spatial distributions

Poor root development due to constraining soil conditions could be an important factor influencing health of urban trees. Therefore, there is a need for efficient techniques to analyze the spatial distribution of tree roots. An analytical procedure for describing tree rooting patterns from X-ray computed tomography (CT) data is described and illustrated. Large irregularly shaped specimens of undisturbed sandy soil were sampled from Various positions around the base of trees using field impregnation with epoxy resin, to stabilize the cohesionless soil. Cores approximately 200 mm in diameter by 500 mm in height were extracted from these specimens. These large core samples were scanned with a medical X-ray CT device, and contiguous images of soil slices (2 mm thick) were thus produced. X-ray CT images are regarded as regularly-spaced sections through the soil although they are not actual 2D sections but matrices of voxels similar to 0.5 mm x 0.5 mm x 2 mm. The images were used to generate the equivalent of horizontal root contact maps from which three-dimensional objects, assumed to be roots, were reconstructed. The resulting connected objects were used to derive indices of the spatial organization of roots, namely: root length distribution, root length density, root growth angle distribution, root spatial distribution, and branching intensity. The successive steps of the method, from sampling to generation of indices of tree root organization, are illustrated through a case study examining rooting patterns of valuable urban trees. (C) 1999 Elsevier Science B.V. All rights reserved.

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.

Mapping cortical change in Alzheimer's disease, brain development, and schizophrenia

This paper describes algorithms that can identify patterns of brain structure and function associated with Alzheimer's disease, schizophrenia, normal aging, and abnormal brain development based on imaging data collected in large human populations. Extraordinary information can be discovered with these techniques: dynamic brain maps reveal how the brain grows in childhood, how it changes in disease, and how it responds to medication. Genetic brain maps can reveal genetic influences on brain structure, shedding light on the nature-nurture debate, and the mechanisms underlying inherited neurobehavioral disorders. Recently, we created time-lapse movies of brain structure for a variety of diseases. These identify complex, shifting patterns of brain structural deficits, revealing where, and at what rate, the path of brain deterioration in illness deviates from normal. Statistical criteria can then identify situations in which these changes are abnormally accelerated, or when medication or other interventions slow them. In this paper, we focus on describing our approaches to map structural changes in the cortex. These methods have already been used to reveal the profile of brain anomalies in studies of dementia, epilepsy, depression, childhood and adult-onset schizophrenia, bipolar disorder, attention-deficit/ hyperactivity disorder, fetal alcohol syndrome, Tourette syndrome, Williams syndrome, and in methamphetamine abusers. Specifically, we describe an image analysis pipeline known as cortical pattern matching that helps compare and pool cortical data over time and across subjects. Statistics are then defined to identify brain structural differences between groups, including localized alterations in cortical thickness, gray matter density (GMD), and asymmetries in cortical organization. Subtle features, not seen in individual brain scans, often emerge when population-based brain data are averaged in this way. Illustrative examples are presented to show the profound effects of development and various diseases on the human cortex. Dynamically spreading waves of gray matter loss are tracked in dementia and schizophrenia, and these sequences are related to normally occurring changes in healthy subjects of various ages. (C) 2004 Published by Elsevier Inc.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

Modeling of subcutaneous absorption kinetics of infusion solutions in the elderly using technetium

Absorption kinetics of solutes given with the subcutaneous administration of fluids is ill-defined. The gamma emitter, technitium pertechnetate, enabled estimates of absorption rate to be estimated independently using two approaches. In the first approach, the counts remaining at the site were estimated by imaging above the subcutaneous administration site, whereas in the second approach, the plasma technetium concentration-time profiles were monitored up to 8 hr after technetium administration. Boluses of technetium pertechnetate were given both intravenously and subcutaneously on separate occasions with a multiple dosing regimen using three doses on each occasion. The disposition of technetium after iv administration was best described by biexponential kinetics with a V-ss of 0.30 +/- 0.11 L/kg and a clearance of 30.0 +/- 13.1 ml/min. The subcutaneous absorption kinetics was best described as a single exponential process with a half-life of 18.16 +/- 3.97 min by image analysis and a half-life of 11.58 +/- 2.48 min using plasma technetium time data. The bioavailability of technetium by the subcutaneous route was estimated to be 0.96 +/- 0.12. The absorption half-life showed no consistent change with the duration of the subcutaneous infusion. The amount remaining at the absorption site with time was similar when analyzed using image analysis, and plasma concentrations assuming multiexponential disposition kinetics and a first-order absorption process. Profiles of fraction remaining at the absorption sire generated by deconvolution analysis, image analysis, and assumption of a constant first-order absorption process were similar. Slowing of absorption from the subcutaneous administration site is apparent after the last bolus dose in three of the subjects and can De associated with the stopping of the infusion. In a fourth subject, the retention of technetium at the subcutaneous site is more consistent with accumulation of technetium near the absorption site as a result of systemic recirculation.

Cupped lesions of early onset dental erosion in young southeast Queensland adults

Background: Dental erosion manifests as cupped lesions on cusp apices and in fissures of teeth in patients from southeast Queensland referred with excessive tooth wear When found in young adults, these lesions may indicate early onset of active dental erosion. If the numbers and extent of cupped lesions increase with age, erosion may be a slow cumulative process. Methods: This cross-sectional study recorded the presence or absence and the relative sizes of cupped lesions from all cusps and occlusal fissures on premolar and permanent molar teeth from study models by image analysis. Type-specimens of cupped lesions were examined. Results: The Incidence by tooth reflected time in the mouth, post-tooth emergence. A linear increase in lesion number and size, with age, was found. However, cupped lesions occurred on mandibular first molar cusp apices as often, and attained greater extent, in adults under 27 years compared with older subjects. Conclusion: Marked differences were found between lesion number and size, between maxillary and mandibular molar sites that reflect differences in salivary protection against dental erosion. The significance of this study is that the mandibular first permanent molar indicates the age of onset and severity of dental erosion.

A nanometre-scale non-periodic structural variation in high temperature superconducting ceramics and the implications for properties

Non-periodic structural variation has been found in the high T-c cuprates, YBa2Cu3O7-x and Hg0.67Pb0.33Ba2Ca2Cu3O8+delta, by image analysis of high resolution transmission electron microscope (HRTEM) images. We use two methods for analysis of the HRTEM images. The first method is a means for measuring the bending of lattice fringes at twin planes. The second method is a low-pass filter technique which enhances information contained by diffuse-scattered electrons and reveals what appears to be an interference effect between domains of differing lattice parameter in the top and bottom of the thin foil. We believe that these methods of image analysis could be usefully applied to the many thousands of HRTEM images that have been collected by other workers in the high temperature superconductor field. This work provides direct structural evidence for phase separation in high T-c cuprates, and gives support to recent stripes models that have been proposed to explain various angle resolved photoelectron spectroscopy and nuclear magnetic resonance data. We believe that the structural variation is a response to an opening of an electronic solubility gap where holes are not uniformly distributed in the material but are confined to metallic stripes. Optimum doping may occur as a consequence of the diffuse boundaries between stripes which arise from spinodal decomposition. Theoretical ideas about the high T-c cuprates which treat the cuprates as homogeneous may need to be modified in order to take account of this type of structural variation.

Liquid distribution in wet granulation: dimensionless spray flux

This study investigates binder distribution in wet granulation and focuses on the nucleation zone, which is the area where the liquid binder and powder surface come into contact and form the initial nuclei. An equipment independent parameter, dimensionless spray flux Psi (a), is defined to characterise the most important process parameters in the nucleation process: solution flowrate, powder flux, and binder drop size. Ex-granulator experiments are used to study the relationship between dimensionless spray flux, process variables and the coverage of binder fluid on the powder surface. Lactose monohydrate powder on a variable speed riffler passed under a flat spray once only. Water and 7% HPC solution at two spray pressures were used as binders. Experiments with red dye and image analysis demonstrate that changes in dimensionless spray flux correlate with a measurable difference in powder surface coverage. Nucleation experiments show that spray flux controls the size and shape of the nuclei size distribution. At low Psi (a), the system operates in the drop controlled regime, where one drop forms one nucleus and the nuclei size distribution is narrow. At higher Psi (a), the powder surface cakes creating a broader size distribution. For controlled nucleation with the narrowest possible size distribution, it is recommended that the dimensionless spray flux be less than 0.1 to be in the drop-controlled regime. (C) 2001 Elsevier Science S.A. All rights reserved.

Mass balance characterisation of Al-7Si-Mg alloy microstructures as a function of solution treatment time

Mass balance calculations were performed to model the effect of solution treatment time on A356 and A357 alloy microstructures. Image analysis and electron probe microanalysis were used to characterise microstructures and confirm model predictions. In as-cast microstructures, up to 8 times more Mg is tied up in the pi-phase than in Mg2Si. The dissolution of pi is accompanied by a corresponding increase in the amount of beta-phase. This causes the rate of pi dissolution to be limited by the rate of beta formation. It is predicted that solution treatments of the order of tens of minutes at 540degreesC produce near-maximum T6 yield strengths, and that Mg contents in excess of 0.52 wt% have no advantage.

Cationic drug pharmacokinetics in diseased livers determined by fibrosis index, hepatic protein content, microsomal activity, and nature of drug

The disposition kinetics of six cationic drugs in perfused diseased and normal rat livers were determined by multiple indicator dilution and related to the drug physicochemical properties and liver histopathology. A carbon tetrachloride (CCl4)induced acute hepatocellular injury model had a higher fibrosis index (FI), determined by computer-assisted image analysis, than did an alcohol-induced chronic hepatocellular injury model. The alcohol-treated group had the highest hepatic alpha(1)- acid glycoprotein, microsomal protein (MP), and cytochrome P450 (P450) concentrations. Various pharmacokinetic parameters could be related to the octanol-water partition coefficient (log P-app) of the drug as a surrogate for plasma membrane partition coefficient and affinity for MP or P450, the dependence being lower in the CCl4-treated group and higher in the alcohol-treated group relative to controls. Stepwise regression analysis showed that hepatic extraction ratio, permeability-surface area product, tissue-binding constant, intrinsic clearance, partition ratio of influx (k(in)) and efflux rate constant (k(out)), and k(in)/k(out) were related to physicochemical properties of drug (log P-app or pK(a)) and liver histopathology (FI, MP, or P450). In addition, hepatocyte organelle ion trapping of cationic drugs was evident in all groups. It is concluded that fibrosis-inducing hepatic disease effects on cationic drug disposition in the liver may be predicted from drug properties and liver histopathology.

Scale-up of mixer granulators for effective liquid distribution

There is considerable anecdotal evidence from industry that poor wetting and liquid distribution can lead to broad granule size distributions in mixer granulators. Current scale-up scenarios lead to poor liquid distribution and a wider product size distribution. There are two issues to consider when scaling up: the size and nature of the spray zone and the powder flow patterns as a function of granulator scale. Short, nucleation-only experiments in a 25L PMA Fielder mixer using lactose powder with water and HPC solutions demonstrated the existence of different nucleation regimes depending on the spray flux Psi(a)-from drop-controlled nucleation to caking. In the drop-controlled regime at low Psi(a) values. each drop forms a single nucleus and the nuclei distribution is controlled by the spray droplet size distribution. As Psi(a) increases, the distribution broadens rapidly as the droplets overlap and coalesce in the spray zone. The results are in excellent agreement with previous experiments and confirm that for drop-controlled nucleation. Psi(a) should be less than 0.1. Granulator flow studies showed that there are two powder flow regimes-bumping and roping. The powder flow goes through a transition from bumping to roping as impeller speed is increased. The roping regime gives good bed turn over and stable flow patterns. This regime is recommended for good liquid distribution and nucleation. Powder surface velocities as a function of impeller speed were measured using high-speed video equipment and MetaMorph image analysis software, Powder surface velocities were 0.2 to 1 ms(-1)-an order of magnitude lower than the impeller tip speed. Assuming geometrically similar granulators, impeller speed should be set to maintain constant Froude number during scale-up rather than constant tip speed to ensure operation in the roping regime. (C) 2002 Published by Elsevier Science B.V.

Induction of metallothionein in human skin by routine exposure to sunlight: Evidence for a systemic response and enhanced induction at certain body sites

Expression of metallothionein, an antioxidant induced by a variety of stimuli including ultraviolet light, was quantitated by immunohistochemistry in the skin of males aged over 50 who had known short- and long-term exposures to sunlight. Skin punch biopsies were taken from two sites in each subject: the hand in all subjects and a range of other sites matched to patients with a previously excised primary melanoma. Metallothionein expression (strongest in the basal layers of the epidermis and primarily nuclear) was associated with both short- and long-term exposure to sunlight. A plateau of staining intensity was reached after 3 h sun exposure, within the previous 3 d before biopsy. Expression was also elevated in the nonexposed skin sites of subjects who had recent sun exposure, indicating a systemic response to exposure of remote sites. Using the skin of the hand to normalize responses to chronic exposure between individuals, the systemically modulated response to sunlight was significantly greater on the unexposed back than on other sites. The possibility of ultraviolet-induced cytokines selectively modifying the response of skin on a site-specific basis was investigated. The circulating leukocytes, but not lymphocytes, of two individuals exposed to 1 minimal erythema dose whole-body solar-simulated ultraviolet showed increased interleukin-6 mRNA 4 h after exposure. Interleukin-6 was not directly induced in these cell populations 4 h after ultraviolet A or ultraviolet B irradiation ex vivo . Leukocytes may therefore contribute to and amplify the systemic effects of ultraviolet-induced interleukin-6 and metallothionein expression.