Spatial Statistics II

# Spatial Statistics II

Télécharger la présentation

## Spatial Statistics II

- - - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - - -
##### Presentation Transcript

1. Spatial Statistics II RESM 575 Spring 2011 Lecture 11

2. What we will use spatial stats to find • How features are distributed • Where are the clusters • What is the pattern created by the features • What are the relationships between sets of features or values

3. Notes on Nearest Neighbor • Use if there is a direct interaction between features • Good for data analyzed along a line (ie plants or wildlife observations collected along a transect) • Stat is calculated using the distance between features

4. . . . . . . . . . . . . . . . . . . . . . . . . Counting the number of features within defined distances • Known as the K-function, or Ripley’s K-function • Here you specify a distance interval and use GIS to calculate the average number of neighboring features within the distance of each feature • If the ave # of features found at a distance is greater than the ave concentration of features throughout the study area, then the distribution is considered clustered at that distance K is the distance bands

5. Notes on K-function • Good for if you are interested in how the pattern changes at different scales of analysis • Includes all distances to neighbors within a given distance • Emergency calls, bird nests

6. K function in GIS • Finds the distance from each point to every other point • Then, for each point, counts up the number of surrounding points within a given distance

7. Interpreting the K stat

8. Paper on class website for today.

9. Factors influencing the K function • Points near the edge of a study area are likely to have less neighboring points • At larger distances its more likely that fewer neighboring points will be found than actually exist

10. Spatial Statistics III

11. Geostatistical Analyst of ArcGIS 10 • For advanced surface modeling • Extension of ArcGIS 10 • Tools for creating a statistically valid surface

12. Loading the Geostatistical Extension Customize - > Extensions Customize - > Tool bars -> Geostatistical Analyst

13. Further reading • Armstrong, M. 1998. Basic Linear Geostatistics. Springer, Berlin. • Chiles, J. and Delfiner, P. 1999. Geostatistics. Modeling Spatial Uncertainty. John Wiley and Sons, New York. • Cressie, N. 1988. Spatial prediction and ordinary kriging. Mathematical Geology 20:405-421. (Erratum, Mathematical Geology 21: 493-494) • Cressie, N. 1990. The origins of kriging. Mathematical Geology 22:239-252. • Isaaks, E.H. and Srivastrava, R.M. 1989. An introduction to Applied Geostatistics. Oxford University Press, New York. • Johnston, Kevin, Jay M. Ver Hoef, Konstantin Krivoruchko, and Neil Lucas. Using ArcGis Geostatistical Analyst, 2001. Environmental System Research Institute, Redlands, CA. • Shaw, Gareth and Dennis Wheeler.  Statistical Techniques in Geographical Analysis, 1994. David Fulton Publishers, London.

14. Part A. Background on Interpolation Techniques Deterministic methods Geostatistical methods Some important principles

15. Interpolating a surface • Generate the most accurate surface • Sample point data as input • Characterize the error and variability of the predicted surface

16. Interpolation techniques • Deterministic • Use mathematical functions for interpolation • IDW, global and local polynomial, radial basis • Geostatistical • Relies on both statistical and mathematical methods • Can be used to assess the uncertainty of the predictions NOTE: Both rely on similarity of nearby points to create the surface

17. Deterministic techniques • Inverse distance weighted • Global polynomial • Local polynomial • Radial bias functions

18. Inverse distance weighted • Reasonably accurate if the points are evenly distributed and the surface characteristics do not change across the landscape • Values of closer points are weighted more heavily than those further away

19. IDW (in the Geostatistical Analyst)

20. Global polynomial • Identify and model local structures and surface trends • Fit a plane between the sample points One bend = 2nd order Two bends = 3rd order Etc… Plane = first order

21. GP

22. GP

23. Local polynomial • Fitting many smaller overlapping planes

24. LP

25. Radial basis • Captures global trends and picks up local variation (bending and stretching of surface to match all the measured values)

26. RB

27. Geostatistical methods • Based on statistical methods not just mathematical • Include spatial autocorrelation • Provide a measure of certainty or accuracy • Kriging • Cokriging

28. Principals of Geostat Methods • Unlike the deterministic methods, geostatistics assumes that all values are a result of a random process with dependence • What does this mean?

29. Ex • Flip three coins and determine if H or T • The fourth coin will not be flipped; it will be laid down based on what the 2nd and 3rd are • Rule to lay the 4th: • if the 2nd and 3rd are tails, the fourth is the opposite of the first, if not then the 4th is same as first

30. How does this relate to predicting locations in an interpolation? • In coin ex, dependence rules were given • In reality, dependence rules are not known • In geostats, there are two key tasks • To uncover the dependence rules • To make predictions KEY: the predictions come from knowing the dependency rules!

31. Principles of Geostat Methods • Besides random process with dependence… • Stationarity • Mean stationarity • mean is constant between samples and is independent of location • Second order stationarity for covariance • covariance is the same between any two points that are at the same distance and direction apart no matter which points you choose • Intrinsic stationarity for semivariograms • variance of the difference is the same between any two points that are at the same distance and direction apart no matter which two points you choose

32. Kriging • In geostats, there are two key tasks • To uncover the dependence rules • To make predictions Semivariogram and covariance functions Interpolate areas

33. Kriging • Similar to IDW (weights surrounding values to derive a prediction) • Different in that it incorporates the spatial arrangement among the measured points (must calculate spatial autocorrelation)

34. Kriging

35. Cokriging • Uses information on several variable types • Requires much more estimation (autocorrelation for each variable and cross-correlations)

36. Kriging process • Calculate the empirical semivariogram • Fit a model • Make a prediction

37. Empirical semivariogram • Tests for spatial autocorrelation (things closer are more alike) spatial modeling, structural analysis or variography Combinations of the points low on both the x and y axis have more autocorrelation Increasing dissimilarity Increasing distance

38. Fit a Model • Defining a line (weighted least squares) that provides the best fit through the points in the empirical semivariogram cloud • Line is considered a model quantifying the spatial autocorrelation in a model

39. Make a prediction • From the kriging weights for the measured values, you can calculate a prediction for the location with the unknown value.

40. Part B. The Geostatistical Process Explore the data Fit a model Perform diagnostics Compare the models