Dispersing geocoded points within a Polygon

1947
1
11-27-2013 02:49 AM
MTP
by
New Contributor II
Hi,

I have a problem with dispersing geocoded points (randomely/)equally within a polygon. More specific I have geocoded adresses using zip codes and would like to disperse them within a polygon. I know that it is possible to disperse markers but I would like to disperse the points. Dispersing them randomely would be not a problem as it is i.e. possible with Geowizard.

Therefore, i have searched the forums and have found some solutions for Avenue but in my ArcGIS 10 this obviously is not working. Based on the post here: http://arcpy.wordpress.com/2013/06/07/disperse-overlapping-points/#comment-220 I have tried to develop a script. But it is not working.

Can anyone advice what could be changed to make it work? Is it necessary to have the x, y already in the tael definded as columns?

Thank you very much!

import random
import arcpy

arcpy.env.workspace = �??D:\�?�\database.gdb�?�

in_points = �??D:\�?�\�?�\geodatabase.gdb\Geocoding_Results.shp�?�
polygon = �??D:\�?�\geodatabase.gdb\zipcode.shp�?�

def point_in_poly(poly, x, y):
�??�?�"Returns if the point is inside the polygon.

Parameters:
poly: arcpy.Polygon() geometry
x: x coordinate (float)
y: y coordinate (float)

�??�?�"
pg = arcpy.PointGeometry(arcpy.Point(x, y), poly.spatialReference)
return poly.contains(pg)

def disperse_points(in_points, polygon):
�??�?�"Randomly disperse points inside a polygon.

Parameters:
in_points: Point feature class/layer (with or without selection)
polygon: arcpy.Polygon() geometry

�??�?�"

lenx = polygon.extent.width
leny = polygon.extent.height

with arcpy.da.UpdateCursor(in_points, "shape@xy") as points:
for p in points:
x = (random.random() * lenx) + polygon.extent.XMin
y = (random.random() * leny) + polygon.extent.YMin
inside = point_in_poly(polygon, x, y)
while not inside:
x = (random.random() * lenx) + polygon.extent.XMin
y = (random.random() * leny) + polygon.extent.YMin
inside = point_in_poly(polygon, x, y)
points.updateRow([(x, y)])
Tags (2)
0 Kudos
1 Reply
MTP
by
New Contributor II
To visualize what I would like to achieve:

The geocoding result (with only one ZIP/normally I have around 500):
[ATTACH=CONFIG]29401[/ATTACH]

And the result I would like to receive:

[ATTACH=CONFIG]29402[/ATTACH]
0 Kudos