import arcpy# 定义输入图层和输出图层input_layer = arcpy.GetParameterAsText(0)output_layer = arcpy.GetParameterAsText(1)# 定义合并条件merge_field = arcpy.GetParameterAsText(2) # 要合并的字段名merge_gap = arcpy.GetParameterAsText(3) # 相邻图斑面积差距# 定义面积筛选条件selection_field = arcpy.GetParameterAsText(4) # 用于筛选的字段名selection_value = arcpy.GetParameterAsText(5) # 筛选条件selection_area = float(arcpy.GetParameterAsText(6)) # 面积筛选阈值# 进行面积筛选arcpy.MakeFeatureLayer_management(input_layer, "layer", "{}='{}' AND SHAPE_AREA > {}".format(selection_field, selection_value, selection_area))# 寻找相邻图斑arcpy.PolygonNeighbors_analysis("layer", "neighbors", ["FID"])# 定义合并列表merge_list = []# 进行合并with arcpy.da.SearchCursor("neighbors", ['src_FID', 'nbr_FID', 'src_{}'.format(merge_field)]) as cursor: for row in cursor: if row[0] < row[1]: src_geom = arcpy.da.SearchCursor("layer", ["SHAPE@"], "FID={}".format(row[0])).next()[0] nbr_geom = arcpy.da.SearchCursor("layer", ["SHAPE@"], "FID={}".format(row[1])).next()[0] if src_geom.area < nbr_geom.area: src_geom, nbr_geom = nbr_geom, src_geom if src_geom.area - nbr_geom.area > merge_gap: merge_list.append(row[0]) merge_list.append(row[1])arcpy.management.Dissolve("layer", output_layer, "{}".format(merge_field), None, "MULTI_PART", "DISSOLVE_LINES")# 删除被合并的图斑with arcpy.da.UpdateCursor(output_layer, ['FID']) as cursor: for row in cursor: if row[0] in merge_list: cursor.deleteRow()运行上面代码报错Traceback (most recent call last): File "D:\实验YY\新建文件夹 (2)\批量合并小图斑.py", line 13, in <module> File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\analysis.py", line 984, in PolygonNeighbors raise e ExecuteError: 执行失败。参数无效。 ERROR 000732: 输入要素: 数据集 layer 不存在或不受支持 执行(PolygonNeighbors)失败。 执行(批量合并小图斑)失败。请改正代码

时间: 2023-06-27 21:04:13 浏览: 50
根据错误提示,看起来是找不到名为 "layer" 的图层。可能是因为输入的图层名称不正确或图层不存在。请检查代码中的输入和输出参数名称是否正确,并确保输入的图层存在于数据中。此外,您还可以尝试在代码中添加一些调试语句来确定出错的位置,例如在“定义输入图层和输出图层”后添加一个 print 语句以确保正确获取了输入和输出参数。
相关问题

# 导入需要使用的模块 import arcpy # 设置输入参数 input_features = arcpy.GetParameterAsText(0) # 输入要素类 output_features = arcpy.GetParameterAsText(1) # 输出要素类 minimum_area = arcpy.GetParameterAsText(2) # 面积阈值,单位为平方米 # 将面积阈值转换为浮点型 try: minimum_area = float(minimum_area) except ValueError: arcpy.AddError("面积阈值必须为数字!") sys.exit() # 创建要素图层并添加要素 arcpy.MakeFeatureLayer_management(input_features, "input_layer") arcpy.SelectLayerByAttribute_management("input_layer", "NEW_SELECTION", "Shape_Area < {}".format(minimum_area)) arcpy.CopyFeatures_management("input_layer", output_features) # 输出结果 arcpy.SetParameterAsText(3, output_features)运行错误:IndentationError: unexpected indent (批量删除小面积.py, line 8) 执行(删除小面)失败。请改正代码

这个错误通常是由于代码的缩进不正确导致的。可以检查一下第8行的缩进是否正确。 正确的代码应该像这样: ```python # 导入需要使用的模块 import arcpy # 设置输入参数 input_features = arcpy.GetParameterAsText(0) # 输入要素类 output_features = arcpy.GetParameterAsText(1) # 输出要素类 minimum_area = arcpy.GetParameterAsText(2) # 面积阈值,单位为平方米 # 将面积阈值转换为浮点型 try: minimum_area = float(minimum_area) except ValueError: arcpy.AddError("面积阈值必须为数字!") sys.exit() # 创建要素图层并添加要素 arcpy.MakeFeatureLayer_management(input_features, "input_layer") arcpy.SelectLayerByAttribute_management("input_layer", "NEW_SELECTION", "Shape_Area < {}".format(minimum_area)) arcpy.CopyFeatures_management("input_layer", output_features) # 输出结果 arcpy.SetParameterAsText(3, output_features) ```

import arcpy import sys # 设置输入参数 input_features = arcpy.GetParameterAsText(0) # 输入要素类 output_features = arcpy.GetParameterAsText(1) # 输出要素类 minimum_area = arcpy.GetParameterAsText(2) # 面积阈值,单位为平方米 # 将面积阈值转换为浮点型 try: minimum_area = float(minimum_area) except ValueError: arcpy.AddError("面积阈值必须为数字!") sys.exit() # 创建要素图层并添加要素 arcpy.MakeFeatureLayer_management(input_features, "input_layer") arcpy.SelectLayerByAttribute_management("input_layer", "NEW_SELECTION", "Shape_Area < {}".format(minimum_area)) arcpy.CopyFeatures_management("input_layer", output_features) # 输出结果 arcpy.SetParameterAsText(3, output_features)运行错误:IndentationError: unexpected indent (批量删除小面积.py, line 8) 执行(删除小面)失败。请改正代码

这个错误是由于代码的缩进不正确导致的,建议按照以下方式修改代码: ```python import arcpy import sys # 设置输入参数 input_features = arcpy.GetParameterAsText(0) # 输入要素类 output_features = arcpy.GetParameterAsText(1) # 输出要素类 minimum_area = arcpy.GetParameterAsText(2) # 面积阈值,单位为平方米 # 将面积阈值转换为浮点型 try: minimum_area = float(minimum_area) except ValueError: arcpy.AddError("面积阈值必须为数字!") sys.exit() # 创建要素图层并添加要素 arcpy.MakeFeatureLayer_management(input_features, "input_layer") arcpy.SelectLayerByAttribute_management("input_layer", "NEW_SELECTION", "Shape_Area < {}".format(minimum_area)) arcpy.CopyFeatures_management("input_layer", output_features) # 输出结果 arcpy.SetParameterAsText(3, output_features) ``` 这样就可以避免缩进错误了。

相关推荐

运行import arcpy # 定义输入图层和输出图层 input_layer = arcpy.GetParameterAsText(0) output_layer = arcpy.GetParameterAsText(1) # 定义合并条件 merge_field = arcpy.GetParameterAsText(2) # 要合并的字段名 merge_gap = arcpy.GetParameterAsText(3) # 相邻图斑面积差距 # 定义面积筛选条件 selection_field = arcpy.GetParameterAsText(4) # 用于筛选的字段名 selection_value = arcpy.GetParameterAsText(5) # 筛选条件 selection_area = float(arcpy.GetParameterAsText(6)) # 面积筛选阈值 # 进行面积筛选 arcpy.MakeFeatureLayer_management(input_layer, "layer", "{}='{}' AND SHAPE_AREA > {}".format(selection_field, selection_value, selection_area)) # 寻找相邻图斑 arcpy.PolygonNeighbors_analysis("layer", "neighbors", ["FID"]) # 定义合并列表 merge_list = [] # 进行合并 with arcpy.da.SearchCursor("neighbors", ['src_FID', 'nbr_FID', 'src_{}'.format(merge_field)]) as cursor: for row in cursor: if row[0] < row[1]: src_geom = arcpy.da.SearchCursor("layer", ["SHAPE@"], "FID={}".format(row[0])).next()[0] nbr_geom = arcpy.da.SearchCursor("layer", ["SHAPE@"], "FID={}".format(row[1])).next()[0] if src_geom.area < nbr_geom.area: src_geom, nbr_geom = nbr_geom, src_geom if src_geom.area - nbr_geom.area > merge_gap: merge_list.append(row[0]) merge_list.append(row[1]) arcpy.management.Dissolve("layer", output_layer, "{}".format(merge_field), None, "MULTI_PART", "DISSOLVE_LINES") # 删除被合并的图斑 with arcpy.da.UpdateCursor(output_layer, ['FID']) as cursor: for row in cursor: if row[0] in merge_list: cursor.deleteRow()时报错Traceback (most recent call last): File "D:\实验YY\新建文件夹 (2)\批量合并小图斑.py", line 13, in <module> File "c:\program files (x86)\arcgis\desktop10.2\arcpy\arcpy\analysis.py", line 984, in PolygonNeighbors raise e ExecuteError: 执行失败。参数无效。 ERROR 000732: 输入要素: 数据集 layer 不存在或不受支持 执行(PolygonNeighbors)失败。 执行(批量合并小图斑)失败。请改正代码

import sys sys.tracebacklimit = 0 import os os.environ['PYTHONUNBUFFERED'] = '1'import arcpy # 获取参数 input_features = arcpy.GetParameterAsText(0) join_field = arcpy.GetParameterAsText(1) target_feature = arcpy.GetParameterAsText(2) target_field = arcpy.GetParameterAsText(3) area_threshold = arcpy.GetParameterAsText(4) # 创建空间连接 join_result = arcpy.SpatialJoin_analysis(input_features, target_feature, "in_memory/spatial_join", "JOIN_ONE_TO_ONE", "KEEP_ALL", "", "INTERSECT") # 使用MakeFeatureLayer创建要素图层,并使用AddFieldDelimiters处理字段名称 join_layer = arcpy.management.MakeFeatureLayer(join_result, "join_layer").getOutput(0) join_field_name = arcpy.AddFieldDelimiters(join_layer, join_field) # 使用SelectLayerByAttribute选择重叠面积大于阈值的要素 arcpy.management.SelectLayerByAttribute(join_layer, "NEW_SELECTION", "Shape_Area > " + str(area_threshold)) # 使用SummaryStatistics工具进行面积求和 summary_table = arcpy.Statistics_analysis(join_layer, "in_memory/summary_table", [["Shape_Area", "SUM"]], join_field_name) # 使用TableToNumPyArray将结果转换为字典 sum_dict = {} with arcpy.da.TableToNumPyArray(summary_table, [join_field, "SUM_Shape_Area"]) as arr: for row in arr: sum_dict[row[0]] = row[1] # 使用UpdateCursor更新目标要素类的目标字段 with arcpy.da.UpdateCursor(target_feature, [target_field, join_field], sql_clause=(None, "ORDER BY OBJECTID")) as cursor: for row in cursor: join_value = row[1] if join_value in sum_dict: area_sum = sum_dict[join_value] row[0] = area_sum cursor.updateRow(row) # 导出结果 output_feature = arcpy.GetParameterAsText(5) arcpy.CopyFeatures_management(target_feature, output_feature) # 删除游标对象和要素图层对象 del cursor, join_layer运行错误SyntaxError: invalid syntax (空间连接.py, line 4) 执行(ccc)失败。请改正代码

最新推荐

recommend-type

pre_o_1csdn63m9a1bs0e1rr51niuu33e.a

pre_o_1csdn63m9a1bs0e1rr51niuu33e.a
recommend-type

matlab建立计算力学课程的笔记和文件.zip

matlab建立计算力学课程的笔记和文件.zip
recommend-type

FT-Prog-v3.12.38.643-FTD USB 工作模式设定及eprom读写

FT_Prog_v3.12.38.643--FTD USB 工作模式设定及eprom读写
recommend-type

matlab基于RRT和人工势场法混合算法的路径规划.zip

matlab基于RRT和人工势场法混合算法的路径规划.zip
recommend-type

matlab基于matlab的两步定位软件定义接收机的开源GNSS直接位置估计插件模块.zip

matlab基于matlab的两步定位软件定义接收机的开源GNSS直接位置估计插件模块.zip
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

SPDK_NVMF_DISCOVERY_NQN是什么 有什么作用

SPDK_NVMF_DISCOVERY_NQN 是 SPDK (Storage Performance Development Kit) 中用于查询 NVMf (Non-Volatile Memory express over Fabrics) 存储设备名称的协议。NVMf 是一种基于网络的存储协议,可用于连接远程非易失性内存存储器。 SPDK_NVMF_DISCOVERY_NQN 的作用是让存储应用程序能够通过 SPDK 查询 NVMf 存储设备的名称,以便能够访问这些存储设备。通过查询 NVMf 存储设备名称,存储应用程序可以获取必要的信息,例如存储设备的IP地址、端口号、名称等,以便能
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。