8 research outputs found
Inferring Species Richness and Turnover by Statistical Multiresolution Texture Analysis of Satellite Imagery
<div><h3>Background</h3><p>The quantification of species-richness and species-turnover is essential to effective monitoring of ecosystems. Wetland ecosystems are particularly in need of such monitoring due to their sensitivity to rainfall, water management and other external factors that affect hydrology, soil, and species patterns. A key challenge for environmental scientists is determining the linkage between natural and human stressors, and the effect of that linkage at the species level in space and time. We propose pixel intensity based Shannon entropy for estimating species-richness, and introduce a method based on statistical wavelet multiresolution texture analysis to quantitatively assess interseasonal and interannual species turnover.</p> <h3>Methodology/Principal Findings</h3><p>We model satellite images of regions of interest as textures. We define a texture in an image as a spatial domain where the variations in pixel intensity across the image are both stochastic and multiscale. To compare two textures quantitatively, we first obtain a multiresolution wavelet decomposition of each. Either an appropriate probability density function (pdf) model for the coefficients at each subband is selected, and its parameters estimated, or, a non-parametric approach using histograms is adopted. We choose the former, where the wavelet coefficients of the multiresolution decomposition at each subband are modeled as samples from the generalized Gaussian pdf. We then obtain the joint pdf for the coefficients for all subbands, assuming independence across subbands; an approximation that simplifies the computational burden significantly without sacrificing the ability to statistically distinguish textures. We measure the difference between two textures' representative pdf's via the Kullback-Leibler divergence (KL). Species turnover, or diversity, is estimated using both this KL divergence and the difference in Shannon entropy. Additionally, we predict species richness, or diversity, based on the Shannon entropy of pixel intensity.To test our approach, we specifically use the green band of Landsat images for a water conservation area in the Florida Everglades. We validate our predictions against data of species occurrences for a twenty-eight years long period for both wet and dry seasons. Our method correctly predicts 73% of species richness. For species turnover, the newly proposed KL divergence prediction performance is near 100% accurate. This represents a significant improvement over the more conventional Shannon entropy difference, which provides 85% accuracy. Furthermore, we find that changes in soil and water patterns, as measured by fluctuations of the Shannon entropy for the red and blue bands respectively, are positively correlated with changes in vegetation. The fluctuations are smaller in the wet season when compared to the dry season.</p> <h3>Conclusions/Significance</h3><p>Texture-based statistical multiresolution image analysis is a promising method for quantifying interseasonal differences and, consequently, the degree to which vegetation, soil, and water patterns vary. The proposed automated method for quantifying species richness and turnover can also provide analysis at higher spatial and temporal resolution than is currently obtainable from expensive monitoring campaigns, thus enabling more prompt, more cost effective inference and decision making support regarding anomalous variations in biodiversity. Additionally, a matrix-based visualization of the statistical multiresolution analysis is presented to facilitate both insight and quick recognition of anomalous data.</p> </div
Estimated and measured local species-richness.
<p>The local species-richness ( diversity) and the estimated local species-richness (i.e. the Shannon entropy of the green band, using <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone.0046616.e106" target="_blank">Equation 5</a>) are reported in plot (a) and (b) from 1984 to 2011 respectively. In plot (c) the functional relationship between the Shannon entropy of the green band and the local species-richness for WCA 1 is reported. The inset, which reports the normalized entropy, shows the ability of the Shannon entropy to capture at least 70% of the measured local species-richness. This percentage is 80% and 77% in the dry and wet season respectively. The dashed grey curves are the 95% confidence interval of the linear regression exponent. Variabilities of measured exponents are found by bootstrapping over points and deriving slopes by the linear and the Jackknife models <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone.0046616-Warton1" target="_blank">[92]</a>.</p
Remote-sensed images for the Arthur R. Marshall Loxahatchee National Wildlife Refuge (WCA 1) during the dry season for the period 1987β2011.
<p>The first three years (1984β1986) images are not represented. The representative region in which the texture analysis is performed is delineated in red for each image. The red regions are characterized by a cloud cover lower than 20%. The green regions identify where the data of species are available. <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone.0046616.s001" target="_blank">Figure S1</a> reports the images for the wet season.</p
Graphical explanation of the analysis performed for the WCA 1.
<p>Landsat images (from 1984 to 2011) are acquired for the dry and the wet season. The example is reported for the years 1997 and 1998. The changes in vegetation composition are analyzed using and diversity among seasons and among years. For the interseasonal analysis the time-scale is on average six months between seasons of the same year, while for the interannual analysis the time-scale is about a year between the same season of different years. and diversity from data are compared to the estimates of the Shannon entropy (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone.0046616.e106" target="_blank">Equation 5</a>) and of the KL divergence (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone.0046616.e083" target="_blank">Equation 4</a>) for the green band of the images respectively.</p
Interannual entropy of WCA 1 Landsat images and rainfall in the period 1984β2011 considered.
<p>(a, b) Shannon entropy of the representative regions of WCA 1 in the dry and in the wet season respectively (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone-0046616-g001" target="_blank">Figure 1</a>, and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0046616#pone.0046616.s001" target="_blank">Figure S1</a>) for the red, green, and blue bands. The Shannon entropy for the green bad is proportional to the diversity of plant species. (c) Average annual rainfall (in mm) in the dry (red line) and wet season (blue line).</p
Cells in different regions of colonies respond differently to DNA damage.
<p>(AβC) In order to investigate regional heterogeneity, colonies of differentiated and NCS-treated cells were computationally divided into windows, which were then classified according to the density of cells contained. 25 colonies were divided into 164 windows, containing a total of 13,133 cells. As demonstrated, high cell density regions (purple) tend to have low Oct4 (B), and reduced induction of cPARP (C) upon DNA damage, than low (blue) and mixed (red) density regions. (D) Oct4 distribution of the cells in the above classified subregions. Black line shows the distribution of all cells. (E) Proportion of Oct4 positive cells in each subpopulation. (F) proportion of cPARP positive cells.</p
Large-scale mapping of single-cell heterogeneity in Oct4 mRNA expression.
<p>(A) Quantification of Oct4 mRNA level with smFISH in a hES cell. Top: raw image; Bottom: segmented FISH spots. (B) Segmentation of Oct4 spots in hES cells with nuclear boundaries in red, cell boundaries in yellow, and colonies boundaries in green. (<b>C</b>) Method for setting smFISH spot detection thresholds to ensure spot detection consistency across adjacent fields in a large region. (D) Reconstruction of Oct4 mRNA levels across 6 contiguous fields of view using the spot matching procedure detailed in the main text. Each circle represents an individual cell and circle color denotes the expression level of Oct4. Overall, Oct4 levels diminish through differentiation, although small clusters of cells with elevated Oct4 levels can be observed after 2 days of differentiation.</p
Spatial enrichment of G1 and G2 phase cells in colonies.
<p>(A) Illustration of the distance transform applied to cells in all colonies. Outside of the boundaries of the colony, the distance from edge is β0β, and within the border of the colony, the distance ranges between 0 and the maximum radius of the colony. (B) Frequencies of cells as a function of distance from edge for the G1 (red), S (magenta), and G2 (blue) and M (green) subpopulations. Error bars indicate the 95% confidence interval based on 1000 bootstrap samples. Abscissa tick marks indicate the distance intervals from which the population and bootstrap frequencies were calculated. Each distance interval contains approximately 1500 cells. The data points are located halfway in between each interval. The frequencies at each point add up to 1. Of note is the statistically significant peak in the G2/M subpopulation at 25 microns. (C) A cellular map illustrating the periphery (blue), Cell Layer 1 (green), Cell Layer 2 (red), and Cell Layer 3 (blue) subpopulations. (D) Distributions of cell cycle phases for each cell layer across many colonies. There is a significant enrichment of G2 cells in Cell Layer 1 relative to Cell Layers 2 and 3.</p