Create Thiessen Polygons (Analysis)
Creates Thiessen polygons from point input features.
Each Thiessen polygon contains only a single point input feature. Any location within a Thiessen polygon is closer to its associated point than to any other point input feature.
This tool is used to divide the area covered by the point input features into Thiessen or proximal zones. These zones represent full areas where any location within the zone is closer to its associated input point than to any other input point.
- Where S is a set of points in coordinate or Euclidean space (x,y), for any point p in that space, there is one point of S closest to p, except where point p is equidistant to two or more points of S.
- A single proximal polygon (Voronoi cell) is defined by all points p closest to a single point in S, that is, the total area in which all points p are closer to a given point in S than to any other point in S.
Thiessen proximal polygons are constructed as follows:
- All points are triangulated into a triangulated irregular network (TIN) that meets the Delaunay criterion.
- The perpendicular bisectors for each triangle edge are generated, forming the edges of the Thiessen polygons. The location at which the bisectors intersect determine the locations of the Thiessen polygon vertices.
The outside boundary of the output Thiessen polygon feature class is the extent of the point input features plus an additional 10%. If the Extent environment is set to a specific extent window, this tool tool will use the environment setting to set its outside boundary.
The theoretical background for creating Thiessen polygons is as follows:
This tool may produce unexpected results with data in a geographic coordinate system since the Delaunay triangulation method used by the tool works best with data in a projected coordinate system.
The point input features from which Thiessen polygons will be generated.
The output feature class containing the Thiessen polygons that are generated from the point input features.
Determines which attributes from the point input features will be transferred to the output feature class.
The following Python window script demonstrates how to use the CreateThiessenPolygons tool in immediate mode.
import arcpy arcpy.env.workspace = "C:/data/data.gdb" arcpy.CreateThiessenPolygons_analysis("schools", "c:/output/output.gdb/thiessen1", "ALL")
The following stand-alone script demonstrates how to use the CreateThiessenPolygons function.
# Name: CreateThiessenPolygons_Example2.py # Description: Create symmetrical difference between input and update features # Author: ESRI # Import system modules import arcpy from arcpy import env # Set environment settings env.workspace = "C:/data/data.gdb" # Set local variables inFeatures = "schools" outFeatureClass = "c:/output/output.gdb/thiessen1" outFields = "ALL" # Execute CreateThiessenPolygons arcpy.CreateThiessenPolygons_analysis(inFeatures, outFeatureClass, outFields)