Select to view content in your preferred language

Is it possible to count points in polygon and print out what each polygon contains

1882
7
06-18-2021 03:26 AM
Cookie010103
Emerging Contributor
Viewed 9 times
 
0

Sorry for the weird title but here's my sample table for further description:

I have a original point layer:


IDNameX(latitude)Y(longtitude)
ID1DogA1
ID2CatA1
ID3BirdB1
ID4CatA2
ID5CatB1
ID6DogA1

Those points are inside a grid layer, lets say 2x2 for instance, and the points are inside the grid. I'll try my best to demonstrate it in the follow table:


B
1. . . (Three points, which are ID1, ID2, ID6). . (ID3, ID5)
2. (ID4)(NULL)

Next, I use the count_points_in_polygon functino and create a new layer, which by default adds a column NUM_POINT.

Grid_ID NUM_POINT
A13
A22
B11
B20

But what I actually want is a table that shows which 'Name's are in each polygon, like the following table shows. 

Grid_IDNUM_POINT Name
A13Dog, Cat
A21Dog
B12Bird, Cat
B20-

 

Any solution or keywords is much appriciated, thank you!

7 Replies
JohannesLindner
MVP Frequent Contributor
import arcpy, os

points = "Point_layer_or_feature_class"
polygons = "Polygon_layer_or_feature_class"
out_table = "PathToOutputTable"

# Intersect points and polygons
in_features = [points, polygons]
intersect_points = "memory/intersect_points"  # this is for ArcGIS Pro
# intersect_points = "in_memory/intersect_points"  # this is for ArcMap
arcpy.analysis.Intersect(in_features, intersect_points)

# you now have a point feature class with the attributes of the original points and the attributes of the polygon the points are in.

# read the data into a list
fields = ["Grid_ID", "Name"]
point_data = [row for row in arcpy.da.SearchCursor(intersect_points , fields)]

# get the grid ids from the polygons
grid_ids = [row[0] for row in arcpy.da.SearchCursor(polygons, ["Grid_ID"])]
# if you don't care about polygons with 0 points in them, use this instead:
# grid_ids = list(set([p[0] for p in point_data]))

# create out_table
try:
    arcpy.management.Delete(out_table)
except:
    pass
arcpy.management.CreateTable(os.path.dirname(out_table), os.path.basename(out_table))
arcpy.management.AddField(out_table, "Grid_ID", "LONG")
arcpy.management.AddField(out_table, "Count", "LONG")
arcpy.management.AddField(out_table, "Name", "Text")

# for each grid_id, search point_data for all points with that grid_id, get the count and all names, write the result into out_table
with arcpy.da.InsertCursor(out_table, ["Grid_ID", "Count", "Name"]) as cursor:
    for grid_id in grid_ids:
        names_in_polygon = [p[1] for p in point_data if p[0] == grid_id]
        count = len(names_in_polygon)
        unique_names = ", ".join(set(names_in_polygon))
        cursor.insertRow([grid_id, count, unique_names])

Have a great day!
Johannes
Cookie010103
Emerging Contributor

I will try this soon, should this work in normal ArcGIS or Qgis? Anyways Thank you so much!

0 Kudos
JohannesLindner
MVP Frequent Contributor

This is for ArcGIS! Comment out line 9 or 10, depending on whether you work with ArcGIS Pro or ArcMap.

I'm very sure you could do something like that in QGIS, but I dont know how, as I don't work with QGIS.


Have a great day!
Johannes
Cookie010103
Emerging Contributor

I must say I'm an absolute beginner to python in GIS, mostly using the default functions for analysis. During the process several problems occurred.

1. Assuming the 'PathToOutputTalbe' is a variable where the table should be generated, So I created an empty file for it, but I'm not so sure about it.

 

2. It returns this error message, and generates almost a blank table as a result.

 

Runtime error 
Traceback (most recent call last):
  File "<string>", line 30, in <module>
SystemError: error return without exception set

 

 

3. Also, I might not be descriptive enough that the grid layer is actually a polygon, which should have its spacial value showed in the final result so I can visualize it on the map. In my case it's the 'left', 'right', 'top', 'bottom' column I assume.

I came up with the idea that can I simply generate 4 fields, but it didnt work as expected.

 

arcpy.management.AddField(out_table, "Left", "FLOAT")
arcpy.management.AddField(out_table, "Right", "FLOAT")
arcpy.management.AddField(out_table, "Top", "FLOAT")
arcpy.management.AddField(out_table, "Bottom", "FLOAT")

 

  At this point I might just post my original working table for your better understanding.  The desired result is that I can clink on a polygon, and read how many 'Names' are there in the selected polygon. Thank you and sorry for the dumb questions.

1.png

2.png

0 Kudos
Cookie010103
Emerging Contributor
import arcpy, os

points = "merge_simplified_new"
polygons = "Intersect_5x5_reprojected"
out_table = "C:\Users\lcdc7\Desktop\ArcGIS\Output"

in_features = [points, polygons]
intersect_points = "in_memory/intersect_points" 
arcpy.analysis.Intersect(in_features, intersect_points)

fields = ["Grid_ID", "Name"]
point_data = [row for row in arcpy.da.SearchCursor(intersect_points , fields)]
grid_ids = list(set([p[0] for p in point_data]))

try:
    arcpy.management.Delete(out_table)
except:
    pass
arcpy.management.CreateTable(os.path.dirname(out_table), os.path.basename(out_table))
arcpy.management.AddField(out_table, "Grid_ID", "LONG")
arcpy.management.AddField(out_table, "Count", "LONG")
arcpy.management.AddField(out_table, "Name", "Text")


with arcpy.da.InsertCursor(out_table, ["Grid_ID", "Count", "Name"]) as cursor:
    for grid_id in grid_ids:
        names_in_polygon = [p[1] for p in point_data if p[0] == grid_id]
        count = len(names_in_polygon)
        unique_names = ", ".join(set(names_in_polygon))
        cursor.insertRow([grid_id, count, unique_names])

p.s. This is the final code I use.

0 Kudos
DuncanHornby
MVP Notable Contributor

Absolutely NO CODE is required to answer this question just a little knowledge on the correct usage of the SPATIAL JOIN tool.

The trick is to set the Merge rule to JOIN on your Name field in the field mapping as shown below. Also adjust field length as required.

DuncanHornby_1-1624550463219.png

The output is a new FeatureClass grid with join count and concatenation of values from your Name field, see example below.

DuncanHornby_2-1624550505311.png

Any code is overkill and I suspect slower if your datasets get very large.

 

 

JohannesLindner
MVP Frequent Contributor

@Cookie010103: Sorry, kinda forgot about this thread. You should absolutely use Duncan's solution.

@DuncanHornby: This is great, thanks! Sometimes I jump to using code without checking out tools I don't know... TIL


Have a great day!
Johannes
0 Kudos