arcpy获取shp面积

使用ArcPy计算shp文件中面的面积，可以使用以下代码：

import arcpy

# 设置工作空间
arcpy.env.workspace = r"C:\data"

# 输入shp文件路径
in_shp = "my_shapefile.shp"

# 计算面积并输出到新的字段
arcpy.AddField_management(in_shp, "AREA", "DOUBLE")
expression = "!SHAPE.AREA@SQUAREKILOMETERS!"  # 可以根据需要修改为其他面积单位
arcpy.CalculateField_management(in_shp, "AREA", expression, "PYTHON_9.3")

# 输出结果
with arcpy.da.SearchCursor(in_shp, ["SHAPE@", "AREA"]) as cursor:
    for row in cursor:
        print("面积为 {} 平方千米。".format(row[1]))

其中，`AddField_management`函数用于添加一个新的双精度浮点型字段来存储面积；`CalculateField_management`函数用于计算面积并将结果输出到新的字段中；最后，使用`SearchCursor`函数遍历每个要素并打印其面积。

相关问题

# 导入 ArcPy 模块 import arcpy # 获取输入参数 input_layer = arcpy.GetParameterAsText(0) field_name = arcpy.GetParameterAsText(1) field_value = arcpy.GetParameterAsText(2) output_folder = arcpy.GetParameterAsText(3) # 定义筛选条件 where_clause = "{} = '{}'".format(arcpy.AddFieldDelimiters(input_layer, field_name), field_value) # 使用游标进行筛选并导出要素图层 with arcpy.da.SearchCursor(input_layer, ["OID@", "SHAPE@"], where_clause) as cursor: for row in cursor: output_name = "feature_{}.shp".format(row[0]) output_path = output_folder + "\\" + output_name arcpy.CopyFeatures_management(row[1], output_path)运行错误：Traceback (most recent call last): File "D:\多对一空间连接\批量删除小面积.py", line 13, in <module> NameError: name 'where_clause' is not defined 执行(删除小面)失败。请改正代码

看起来出错的原因是 `where_clause` 变量在定义之前被调用了。你可以将 `where_clause` 的定义放在获取输入参数之后，如下所示：

import arcpy

# 获取输入参数
input_layer = arcpy.GetParameterAsText(0)
field_name = arcpy.GetParameterAsText(1)
field_value = arcpy.GetParameterAsText(2)
output_folder = arcpy.GetParameterAsText(3)

# 定义筛选条件
where_clause = "{} = '{}'".format(arcpy.AddFieldDelimiters(input_layer, field_name), field_value)

# 使用游标进行筛选并导出要素图层
with arcpy.da.SearchCursor(input_layer, ["OID@", "SHAPE@"], where_clause) as cursor:
    for row in cursor:
        output_name = "feature_{}.shp".format(row[0])
        output_path = output_folder + "\\" + output_name
        arcpy.CopyFeatures_management(row[1], output_path)

这样应该就能够正常运行了。

import arcpy # 获取用户输入参数 input_feature = arcpy.GetParameterAsText(0) # 输入要素 query_expression = arcpy.GetParameterAsText(1) # 查询表达式 output_feature = arcpy.GetParameterAsText(2) # 输出要素 # 构建查询语句 query = "".format(arcpy.AddFieldDelimiters(input_feature, "面积"), query_expression) # 使用查询语句选择要素 arcpy.MakeFeatureLayer_management(input_feature, "lyr") arcpy.SelectLayerByAttribute_management("lyr", "NEW_SELECTION", query) # 将选择的要素导出到输出位置 arcpy.CopyFeatures_management("lyr", output_feature) # 清除选择 arcpy.SelectLayerByAttribute_management("lyr", "CLEAR_SELECTION")运行错误：Traceback (most recent call last): File "D:\多对一空间连接\批量删除小面积.py", line 15, in <module> NameError: name 'output_feature' is not defined 执行(删除小面)失败。请改正代码

在你的代码中，第9行定义了一个变量`output_feature`，但是在第17行使用该变量时出现了一个`NameError`错误，意味着该变量未定义或未被识别。这可能是因为你没有指定输出要素的路径和名称。

你可以尝试将第9行更改为以下代码：

output_feature = arcpy.GetParameterAsText(2) + "\\output.shp"

这将为输出要素添加一个默认名称`output.shp`，并将其保存到指定的输出路径中。然后在第17行使用该变量时，就不会出现`NameError`错误了。