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 13:44:27 浏览: 112
根据错误提示,第17行出现了意外的缩进。请检查该行以及前面的几行代码是否存在多余的空格或制表符。正确的代码应该是:
```python
# 创建空间连接对象
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)
```
阅读全文