写一个arcgis10.2.2工具箱python2.7脚本,用来:国有农用地使用权宗地分类面积调查表批量生成
时间: 2024-02-06 22:12:59 浏览: 69
以下是一个简单的arcgis10.2.2工具箱python2.7脚本,用于批量生成国有农用地使用权宗地分类面积调查表。该脚本假设您已经有一个包含所有必需字段的要素类,如宗地代码、宗地面积、土地用途等。
```python
# 导入模块
import arcpy
# 设置工作空间和要素类
arcpy.env.workspace = "C:/data"
fc = "land_parcels.shp"
# 定义要素类字段
fields = ["Parcel_ID", "Parcel_Area", "Land_Use"]
# 创建游标,依次访问所有要素
with arcpy.da.SearchCursor(fc, fields) as cursor:
for row in cursor:
# 提取字段信息
parcel_id = row[0]
parcel_area = row[1]
land_use = row[2]
# 计算面积和分类信息
# 这里的计算公式需要根据实际需求进行修改
cultivated_area = parcel_area * 0.8
forest_area = parcel_area * 0.2
other_area = parcel_area - cultivated_area - forest_area
classification = "农用地"
if land_use == "林地":
classification = "林地"
# 创建新的调查表
output_file = "survey_table_{}.xls".format(parcel_id)
with open(output_file, 'w') as f:
# 写入表头
f.write("宗地代码\t宗地面积\t分类\t农用地面积\t林地面积\t其他面积\n")
# 写入数据
f.write("{}\t{}\t{}\t{}\t{}\t{}\n".format(parcel_id, parcel_area, classification, cultivated_area, forest_area, other_area))
arcpy.AddMessage("Created survey table for parcel {}".format(parcel_id))
```
这个脚本将遍历所有要素,并为每个要素创建一个新的调查表。每个调查表将包含宗地代码、宗地面积、土地分类和各类面积信息。您可以根据实际需求修改计算公式和输出格式。
阅读全文