Create Spatially Balanced Points (Geostatisical Analyst)

Zusammenfassung

Generates a set of sample points based on inclusion probabilities resulting in a spatially balanced, i.e. maximized, and thus more efficient sample design.

How Create Spatially Balanced Points works

Verwendung

• The input probability raster contains probability values between 0 and 1. The higher the value, the more likely that the cell will be included in the sample design.
• All values in the study area should have inclusion probabilities >= 0, while all areas outside the study area should have Null values.
• The cell size of the inclusion probability raster determines the finest resolution at which samples will be generated (samples are located at the center of selected cells).
• When point, line, or polygon features are converted to raster (to obtain the input probability raster), the following should be considered:
  • The cell size (resolution) should be fine enough to distinguish all the important features in the population. To accomplish this, the cell size can be set to less than half the minimum distance between features.
  • For line and polygon features, the cell size should be set so that features (like meandering streams) are adequately represented in the resulting raster.
  • The precision with which sample locations can be located in the field should also be considered: if locations are to be found using a GPS with a positional accuracy of 10m, then the cell size should be of 10m.
  • The size of the inclusion probability raster should be considered: if it is very large (small cell size or many cells), then processing time will be longer.
• To avoid artifacts in the output, it is recommended that the number of sample locations be less than 1% of the number of cells in the inclusion probability raster.
• Only the Mersenne Twister random number generator (specified in Random_number_generator) is supported in this tool.
• The random number Seed value can be changed in the Environment Settings.
  • If a random number seed value of 0 is used (the default value), then each time the tool is run, a different set of random numbers will be used and a different set of sample locations will be generated.
  • If the random number seed is set to a fixed number other than 0, then the tool will produce the same set of sample locations each time it is run (until the seed value is changed). Setting the random number seed value to a fixed number > 0 is useful when you want to add to or substract from a set of sample locations by rerunning the tool using a different number of sample points each time.

Syntax

CreateSpatiallyBalancedPoints_ga (in_probability_raster, number_output_points, out_feature_class)

Parameter	Erläuterung	Datentyp
in_probability_raster	This raster defines the inclusion probabilities for each location in the area of interest. The location values range from 0 (low inclusion probability) to 1 (high inclusion probability).	Raster Layer; Mosaic Layer
number_output_points	Specify how many sample locations to generate.	Long
out_feature_class	The output feature class contains the selected sample locations and their inclusion probabilities.	Feature Class

Codebeispiel

CreateSpatiallyBalancedPoints example 1 (Python window)

Create a set of spatially balanced points based on an input inclusion probability raster.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.CreateSpatiallyBalancedPoints_ga("ca_prob", "10", "C:/gapyexamples/output/csbp")

CreateSpatiallyBalancedPoints example 2 (stand-alone script)

Create a set of spatially balanced points based on an input inclusion probability raster.

# Name: CreateSpatiallyBalancedPoints_Example_02.py
# Description: This tool generates a set of sample points based on inclusion
#   probabilities. The resulting sample design is spatially balanced, meaning
#   that the spatial independence between samples is maximized, making the 
#   design more efficient than sampling the study area at random.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"

# Set local variables
inProb = "ca_prob"
numberPoints = 10
outPoints = "C:/gapyexamples/output/csbp"

# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")

# Execute CreateSpatiallyBalancedPoints
arcpy.CreateSpatiallyBalancedPoints_ga(inProb, numberPoints, outPoints)

Umgebung

Output Coordinate System, Extent, Random number generator

Lizenzierungsinformationen

ArcGIS for Desktop Basic: Erfordert Geostatistical Analyst

ArcGIS for Desktop Standard: Erfordert Geostatistical Analyst

ArcGIS for Desktop Advanced: Erfordert Geostatistical Analyst