Latest Contributions by Cookie010103

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. ... View more

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 "", line 30, in 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. ... View more

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