import arcpy # 设置工具箱参数 input_features = arcpy.GetParameterAsText(0) # 输入要素图层 join_features = arcpy.GetParameterAsText(1) # 连接要素图层 output_features = arcpy.GetParameterAsText(2) # 输出要素图层 join_fields = arcpy.GetParameterAsText(3).split(';') # 连接要素图层中需要赋值的字段,用“;”隔开 area_threshold = arcpy.GetParameterAsText(4) # 面积阈值,只有面积大于该值的要素才会被连接 # 创建空间连接对象 spatial_join = arcpy.SpatialJoin_analysis(target_features=input_features, join_features=join_features, out_feature_class=output_features, join_operation="JOIN_ONE_TO_MANY", join_type="KEEP_ALL", match_option="INTERSECT") # 根据面积阈值进行筛选 if len(area_threshold) > 0: arcpy.SelectLayerByAttribute_management(in_layer_or_view=spatial_join, where_clause="Shape_Area >= {}".format(area_threshold)) # 设置字段映射 field_mapping = arcpy.FieldMappings() for field in join_fields: field_map = arcpy.FieldMap() field_map.addInputField(spatial_join, field) output_field = field_map.outputField output_field.name = "{}_{}".format(join_features, field) field_map.outputField = output_field field_mapping.addFieldMap(field_map) # 对要素进行赋值 arcpy.FeatureClassToFeatureClass_conversion(spatial_join, output_features, "", field_mapping)运行错误:IndentationError: unexpected indent (空间连接.py, line 17) 执行(空间连接多对一)失败。请改正代码
时间: 2023-12-06 18:44:31 浏览: 86
使用arcpy获取简单的矢量数据信息.zip_arcPy 查询_arcpy 数据读入_arcpy 读取数据_arcpy矢量拼接
这个错误通常是由于代码缩进不正确导致的。在你的代码中,第17行有一个额外的缩进,需要将其删除。请使用以下代码:
```python
import arcpy
# 设置工具箱参数
input_features = arcpy.GetParameterAsText(0) # 输入要素图层
join_features = arcpy.GetParameterAsText(1) # 连接要素图层
output_features = arcpy.GetParameterAsText(2) # 输出要素图层
join_fields = arcpy.GetParameterAsText(3).split(';') # 连接要素图层中需要赋值的字段,用“;”隔开
area_threshold = arcpy.GetParameterAsText(4) # 面积阈值,只有面积大于该值的要素才会被连接
# 创建空间连接对象
spatial_join = arcpy.SpatialJoin_analysis(target_features=input_features, join_features=join_features, out_feature_class=output_features, join_operation="JOIN_ONE_TO_MANY", join_type="KEEP_ALL", match_option="INTERSECT")
# 根据面积阈值进行筛选
if len(area_threshold) > 0:
arcpy.SelectLayerByAttribute_management(in_layer_or_view=spatial_join, where_clause="Shape_Area >= {}".format(area_threshold))
# 设置字段映射
field_mapping = arcpy.FieldMappings()
for field in join_fields:
field_map = arcpy.FieldMap()
field_map.addInputField(spatial_join, field)
output_field = field_map.outputField
output_field.name = "{}_{}".format(join_features, field)
field_map.outputField = output_field
field_mapping.addFieldMap(field_map)
# 对要素进行赋值
arcpy.FeatureClassToFeatureClass_conversion(spatial_join, output_features, "", field_mapping)
```
阅读全文