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

Cookie010103

Viewed 9 times

 

0

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

I have a original point layer:

ID	Name	X(latitude)	Y(longtitude)
ID1	Dog	A	1
ID2	Cat	A	1
ID3	Bird	B	1
ID4	Cat	A	2
ID5	Cat	B	1
ID6	Dog	A	1

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:

-	A	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
A1	3
A2	2
B1	1
B2	0

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

Grid_ID	NUM_POINT	Name
A1	3	Dog, Cat
A2	1	Dog
B1	2	Bird, Cat
B2	0	-

 

Any solution or keywords is much appriciated, thank you!

JohannesLindner

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])