最新贡献:Cookie010103

import arcpy, os points = "merge_simplified_new" polyons = " intersect_5x5_reprojection " out_table = "C:\Users\lcdc7\Desktop\ArcGIS\Output" in_features = [points, polyons] intersect_points = "in_memory/intersect_points" arcpy.analysis。Intersect(in_features, intersect_points) fields = ["Grid_ID"， "Name"] point_data = [arcpy.da中的行与行。搜索Cursor(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. ．.．查看更多

我必须说，我绝对是GIS中的python初学者，主要使用默认函数进行分析。在这个过程中出现了几个问题。1.假设“PathToOutputTalbe”是一个变量，表应该生成，所以我为它创建了一个空文件，但我不确定。2.它返回这个错误消息，结果生成一个几乎空白的表。运行时错误Traceback(最近一次调用):File ""， line 30, in SystemError: error return without exception set 3。此外，我可能没有足够的描述网格层实际上是一个多边形，它应该在最终结果中显示其空间值，以便我可以在地图上可视化它。在我的例子中，我假设是“左”，“右”，“上”，“下”列。我想到了这个想法，我可以简单地生成4个字段，但它没有像预期的那样工作。 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. ．.．查看更多

我将尝试这个很快，这应该在正常的ArcGIS或Qgis工作吗?无论如何，非常感谢! ．.．查看更多