ANALYSIS OF SPATIAL POINT PATTERNS BY KERNEL IDENTIFICATION

A major objective of plant ecology research is to determine the underlying processes responsible for the observed spatial distribution patterns of plant species. Plants can be approximated as points in space for this purpose, and thus, spatial point pattern analysis has become increasingly popular in ecological research. The basic piece of data for point pattern analysis is a point location of an ecological object in some study region. Therefore, point pattern analysis can only be performed if data can be collected. However, due to the lack of a convenient sampling method, a few previous studies have used point pattern analysis to examine the spatial patterns of grassland species. This is unfortunate because being able to explore point patterns in grassland systems has widespread implications for population dynamics, community‐level patterns, and ecological processes. In this study, we developed a new method to measure individual coordinates of species in grassland communities. This method records plant growing positions via digital picture samples that have been sub‐blocked within a geographical information system (GIS). Here, we tested out the new method by measuring the individual coordinates of Stipa grandis in grazed and ungrazed S. grandis communities in a temperate steppe ecosystem in China. Furthermore, we analyzed the pattern of S. grandis by using the pair correlation function g(r) with both a homogeneous Poisson process and a heterogeneous Poisson process. Our results showed that individuals of S. grandis were overdispersed according to the homogeneous Poisson process at 0–0.16 m in the ungrazed community, while they were clustered at 0.19 m according to the homogeneous and heterogeneous Poisson processes in the grazed community. These results suggest that competitive interactions dominated the ungrazed community, while facilitative interactions dominated the grazed community. In sum, we successfully executed a new sampling method, using digital photography and a geographical information system, to collect experimental data on the spatial point patterns for the populations in this grassland community.

There is considerable interest in cell biology in determining whether, and to what extent, the spatial arrangement of nuclear objects affects nuclear function. A common approach to address this issue involves analyzing a collection of images produced using some form of fluorescence microscopy. We assume that these images have been successfully pre-processed and a spatial point pattern representation of the objects of interest within the nuclear boundary is available. Typically in these scenarios, the number of objects per nucleus is low, which has consequences on the ability of standard analysis procedures to demonstrate the existence of spatial preference in the pattern. There are broadly two common approaches to look for structure in these spatial point patterns. First a spatial point pattern for each image is analyzed individually, or second a simple normalization is performed and the patterns are aggregated. In this paper we demonstrate using synthetic spatial point patterns drawn from predefined point processes how difficult it is to distinguish a pattern from complete spatial randomness using these techniques and hence how easy it is to miss interesting spatial preferences in the arrangement of nuclear objects. The impact of this problem is also illustrated on data related to the configuration of PML nuclear bodies in mammalian fibroblast cells.

This paper attempts to further the research by Odland and Ellis (1992) in applying event history methodology to the analysis of spatial point patterns (that is, event patterns). Its empirical focus is the event pattern derived from the adoption of an agricultural innovation, the Harvestore, in southern Ontario, Canada, from 1963 to 1986. Event history analysis involves the use of discrete‐state, continuous‐time stochastic models to investigate a temporal longitudinal record on discrete variables. Event history models are usually concerned with durations of time between events and the effects of intertemporal time dependencies on future event occurrences. As such, they are often referred to as duration models. Many of the methods used in event history analysis allow the use of other nonnegative interval measurements in place of standard temporal intervals to investigate a series of events. In particular, spatial intervals (or durations) of distances between events may also be accommodated by event history models. Our analysis extends the previous research of Odland and Ellis (1992) by using a wider range of parametric models to explore duration dependence, investigating the role of spatial censoring, and using a more extensive set of explanatory variables. In addition, simulation experiments and graphical tests are used to evaluate the empirical event pattern against one generated from Complete Spatial Randomness. Results indicate that the event pattern formed by the Harvestore adopter farms is clustered (that is, is described by positive duration dependency), the sales agent is a significant factor in the distribution of adopters, and that contrasting results are obtained from the analysis using censored data versus uncensored data.

There is currently a gap in theory for point patterns that lie on the surface of objects, with researchers focusing on patterns that lie in a Euclidean space, typically planar and spatial data. Methodology for planar and spatial data thus relies on Euclidean geometry and is therefore inappropriate for analysis of point patterns observed in non-Euclidean spaces. Recently, there has been extensions to the analysis of point patterns on a sphere, however, many other shapes are left unexplored. This is in part due to the challenge of defining the notion of stationarity for a point process existing on such a space due to the lack of rotational and translational isometries. Here, we construct functional summary statistics for Poisson processes defined on convex shapes in three dimensions. Using the Mapping Theorem, a Poisson process can be transformed from any convex shape to a Poisson process on the unit sphere which has rotational symmetries that allow for functional summary statistics to be constructed. We present the first and second order properties of such summary statistics and demonstrate how they can be used to construct a test statistics to determine whether an observed pattern exhibits complete spatial randomness or spatial preference on the original convex space. We compare this test statistic with one constructed from an analogue L-function for inhomogeneous point processes on the sphere. A study of the Type I and II errors of our test statistics are explored through simulations on ellipsoids of varying dimensions.

A review on anisotropy analysis of spatial point patterns

PDF HTML阅读 XML下载 导出引用 引用提醒 生态学中的点格局研究概况及其在国内的应用 DOI: 10.5846/stxb201607081399 作者: 作者单位: 作者简介: 通讯作者: 中图分类号: 基金项目: 国家自然基金项目（31370633）。 A review of point pattern analysis in ecology and its application in China Author: Affiliation: Fund Project: 摘要 | 图/表 | 访问统计 | 参考文献 | 相似文献 | 引证文献 | 资源附件 | 文章评论 摘要:点格局分析是研究生态学格局的工具之一，近年来在生态学中的应用越来越多。为深入了解点格局分析方法在国内的研究与应用情况，以所总结的研究进展、一般步骤和基本要点为背景，分析评述了1996-2015年期间以点格局为主题的国内中文核心期刊文献。结果表明，在国内生态学格局研究中，应用研究占据主导地位，研究对象广泛，包括以树木为主的乔、灌、草等不同生活型的植物，甚至包括景观；基础研究，包括概括性统计量、零模型与点过程模型等方面，以及专用软件工具包的开发等研究薄弱。在应用中存在一定问题，主要表现为：概括性统计量使用单一，且以Ripley的K-函数及其变形为主；零模型（或点过程模型）是科学问题的统计表达，但是有一半以上的研究未明确给出零模型。建议在未来应用研究中重视多种统计量的组合使用和原假设的建立，在探讨热带、亚热带森林等具有复杂空间结构系统的多样性格局时，考虑对象的不同世代和系统的不同垂直层次，并加强多变量或三维概括性统计量的开发、点格局分析方法与动态过程模型的结合研究等工作。 Abstract:Point pattern analysis, i.e., spatial point pattern analysis (SPPA) is a tool used in ecological pattern research and has recently become increasingly popular in ecology. To better understand SPPA research and its application in China, in this review we examined articles regarding SPPA published from 1996 to 2015 in China, which was based on a summary of the research progress, using steps and key points of SPPA. It was found that most articles focused on application research. Although there were extensive study objects, including forests, shrubs, grasslands, and landscapes, trees in forest communities were the main focus in all these application studies. Fundamental studies were found to be underdeveloped, such as studies regarding summary statistics and null models or point process models in SPPA, and the development of specific software packages for SPPA. Two important problems were found in SPPA applications: Firstly, although there are many alternatives, researchers preferred single summary statistics and only selected the K-function or its transformations; and secondly, although selecting null models as null hypotheses was a key step in SPPA, more than half of the research examined did not specify explicit null models. Therefore, it is recommended that different types of summary statistics are used to capture detailed spatial pattern structures, and the ecological questions or hypotheses should be clearly stated and correctly translated into statistical language, i.e., the null model and tests. Vertical layers and different generations should also be considered when studying systems with complex spatial structures, such as in tropical rainforests and subtropical evergreen broad-leaved forests. To enhance basic SPPA applications, these aspects discussed should be strengthened. 参考文献 相似文献 引证文献

Cox processes are natural models for point process phenomena that are environmentally driven, but much less natural for phenomena driven primarily by interactions amongst the points. The class of log-Gaussian Cox processes (LGCPs) has an elegant simplicity, and one of its attractive features is the tractability of the multivariate normal distribution carries over, to some extent, to the associated Cox process. In the statistical analysis of spatial point patterns, it is often assumed isotropy because of a simpler interpretation and ease of analysis. However, there are many cases in which the spatial structure depends on the direction. In this paper, we introduce new families of anisotropic spatial LGCPs that are useful to model spatial anisotropic point patterns that exhibit a degree of clustering. We propose classes of families consisting of elliptical and non-elliptical models. The former can be reduced to isotropic forms after some rotations, while the latter family goes beyond this property. We derive analytical forms for the covariance of the associated random field, and some second-order characteristics. A moment-based estimation procedure is followed to make inference on the parameters that control the degree of anisotropy. The estimation procedure is evaluated through a simulation study under a variety of scenarios and various degrees of anisotropy. Our methodology is illustrated on two real datasets of earthquakes in South America and the Mediterranean Europe.

A very common way of analyzing different and complicated plant behaviors is to use spatial point pattern analysis, which allows us to assess whether there is any structure present. To test the complete spatial randomness hypothesis, Diggle (1979) proposed a Monte Carlo test whose test statistic is the discrepancy between the estimated and the theoretical form of some summary function, such as the Ripley K-function. In this article, we improve this test by adding various weight functions and get more powerful tests if decreasing and increasing weight functions are used for processes with short and long, respectively, range of interaction.

To study the quantified indices for describing the distributional status of diseases in the spatial point pattern analysis, and provide the a statistic in disease prevention and control. G function, F function, J function and K function were summarized based on the inter-case distances from the view of spatial point pattern analysis. Through the introduction of the basic principles, these were used to analyze the data of acute schistosomiasis in the Guichi District, Chizhou City, Anhui province, with the study distances being from 0 to 3000 meters with 50-meter intervals. The findings were also validated by means of spatial moving scan window performed in SaTScan software. A total of 83 cases of acute schistosomiasis identified in Guichi District, and the point map showed that these cases were mainly distributed around the Yangtze and Qiupu rivers. The computational methods and characteristics of the four quantified indices were obtained. These acute schistosomiasis cases were also explored by using these indices, and the results showed that C and K functions were above 95% confidence interval. While, F and J functions were below 95% confidence interval. Ml these four indices showed that spatial clustering existed in the acute cases, which was consistent with the results of spatial moving scan window method. The latter method also found a most likely cluster, the coordinate of the circle center is (30.65 N, 117.44 E), radius is 2.69 km, and relative risk is 12.78 (DIR = 32. 80, P = 0. 0001). The quantified indices to describe the distributional status of diseases have not only solved the obstacle that spatial point pattern map which could only be analyzed qualitatively, but also supplied a theoretical foundation to deepen spatial clustering analysis. To study the quantified indices for describing the distributional status of

This study presents a massively parallel spatial computing approach that uses general-purpose graphics processing units (GPUs) to accelerate Ripley’s K function for univariate spatial point pattern analysis. Ripley’s K function is a representative spatial point pattern analysis approach that allows for quantitatively evaluating the spatial dispersion characteristics of point patterns. However, considerable computation is often required when analyzing large spatial data using Ripley’s K function. In this study, we developed a massively parallel approach of Ripley’s K function for accelerating spatial point pattern analysis. GPUs serve as a massively parallel platform that is built on many-core architecture for speeding up Ripley’s K function. Variable-grained domain decomposition and thread-level synchronization based on shared memory are parallel strategies designed to exploit concurrency in the spatial algorithm of Ripley’s K function for efficient parallelization. Experimental results demonstrate that substantial acceleration is obtained for Ripley’s K function parallelized within GPU environments.

SUMMARY We present a deviance residual for a spatial heterogeneous Poisson process. The residual is derived from an approximation of the normalising integral of the process. This yields a saturated model estimate of the local intensity at a point, which is the reciprocal of the Dirichlet tile area associated with the point. A number of examples of the use of such a residual in the analysis of spatial point patterns is given.

The wetland plants are very sensitive to hydrological regimes. In this study, the individual sizes of a widely distributed species (i.e., Deyeuxia angustifolia) at three typical marshes with different water table depths (i.e., wet meadow (WM) marsh; seasonal inundated (SI) marsh; perennial inundated (PI) marsh) were investigated in the Sanjiang Plain of Northeast China. Concurrently, three primary point pattern processes (homogeneous Poisson (HP) process, homogeneous Thomas (HT) process, and inhomogeneous Thomas (IT) process) were used to model spatial patterns in the distribution at 0-50cm scale for this tillering-cloning species. The plant height, diameter at breast height (d.b.h), internode number, branches number, and individual aboveground biomass of D. angustifolia decreased sharply with rising water level; however, its density and coverage increased first and then decreased as water level increases. The distribution of D. angustifolia totally diverged from the complete spatial randomness (CSR) model (i.e., HP process) suggesting strong aggregation at 0-50cm scale in all marshes, and aggregated intensity enhanced with increasing water level. Interestingly, the spatial distribution of D. angustifolia fits better with the nested double-cluster model (i.e., IT process) at all scales in WM and SI marshes, indicating that there is a series of clustered patterns under the slight flood stress. However, the spatial pattern fits well with the Poisson cluster model (i.e., HT process) at all scales in PI marsh, implying the small-scale clustering disappeared with the intensification of flooding stress. Our results highlight that the D. angustifolia population could adapt to flooding stress in a certain degree via individual miniaturization strategies and multi-aggregation mechanisms in the freshwater wetlands.

Associations between two spatial processes can be due to a real dependence between the two processes or to the dependence on common underlying variables. We propose to test the existence of a real dependence by use of local tests, leading to a global test of real dependence and a map of local interactions. We present first how classical interaction tests based on random rotations between completely observed processes such as those developed by Berman (Berman. Appl. Statist. (1986) 35, 54-62), can be integrated in local analyses. For this purpose, tests are first performed locally, and the distribution of their p-values is then compared to the corresponding value under the null hypothesis. A similar approach is proposed to test non-stationarity of a point pattern by using distance statistics popularized by Diggle (Diggle. Statistical Analysis of Spatial Point Patterns. (1983) Academic Press, New York). The problem of testing the interaction between a random field and a censoring area pattern process is discussed and an approach similar to the preceding ones is then proposed. The methods are mainly applied to agricultural examples but they can be applied to any microscopical images for which one wishes to analyse the spatial structure.

Single tree detection is important for monitoring tree plantations in order to estimate crop yield, which can be used for assessing food and biofuels. In the field of remote sensing, there have been studies regarding single tree extraction using Light Detection and Ranging (LiDAR) point cloud data. Canopy Height Model (CHM) can be simply derived from these LiDAR data where trees can be analysed using different physical characteristics such as tree crown width, and height. However, using these physical characteristics can’t be very useful for identifying which trees are within a plantation, a forest, or a single tree. This is the motivation for developing algorithms of automated differentiation of random trees from plantation trees. Object Based Image Analysis (OBIA) was used to acquire the locations of individual trees. Trees were separated from other land cover in the image using a CHM derived from LiDAR point cloud data. Using watershed algorithm in eCognition, trees with good canopy spacing was easily distinguished. A spatial point pattern analysis was done from the extracted locations of trees to differentiate regularly patterned points from randomly distributed points. Nearest neighbour statistics was done to provide measurement of the distribution of the extracted points. The Z-Score parameter from nearest neighbor was then added to the layers in eCognition as supplement to the available height derivatives such as mean, standard deviation, and different textural parameters. Supervised type of classification was used in order to classify the trees on the image. Training samples of each tree from the image were entered in Support Vector Machine (SVM) to classify the different trees. Two sets of validation were done to compare the results of using the object parameters only, and the object parameters and Z-Score. Addition of the Z-Score from the object parameters increased the obtained accuracy on both validation sets. Classification accuracy based on validation set 1 and set 2 increased by about 20% and 6%, respectively. Additional spatial pattern analysis methods and validation can further prove the use of such layers in improving classification of trees.

Spatial point pattern analysis applied to bubble nucleation in silicate melts