用Arcpy批量对raster文件进行区统计

最新推荐文章于 2023-05-16 14:13:26 发布

小怪兽慢慢走

最新推荐文章于 2023-05-16 14:13:26 发布

阅读量1.6k

点赞数 2

分类专栏： Arcpy 文章标签： Arcpy zonal statistics 批量

本文链接：https://blog.csdn.net/darcy_vb/article/details/85919008

版权

Arcpy 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

import arcpy
from arcpy import env
from arcpy.sa import *
import arcpy.da as da


arcpy.CheckOutExtension("spatial")
env.overwriteOutput=1

workingDir="D:\\lyj\\code\\"
env.workspace=workingDir+"globe_temp_raster_result"


in_zone_data=workingDir+"country_188\\"+"country_188.shp"
zone_field="COUNTRY"
in_value_rasters=arcpy.ListRasters()


# 读入新的10年数据时，修改以下代码，startYear=1901 startYear=2001...
startYear=2001

for in_value_raster in in_value_rasters:
    print in_value_raster
    # out_table=env.workspace+"\\dbf\\" + in_value_raster + ".dbf"   # 这样方便，但是cal_1901.txt.tif.dbf这样的文件名python会报错，想了半天才发现是这里错了
    out_table=env.workspace+"\\dbf\\"+ str(startYear) +".dbf"   # 1901.dbf,1902.dbf这样的文件名不会报错

    outZSaT=ZonalStatisticsAsTable(in_zone_data,zone_field,in_value_raster,out_table,"DATA","SUM")

    txtfile=env.workspace+"\\txt\\"+str(startYear)+".txt"
    with open(txtfile,"w") as f:                         # 将print的内容写入.txt文件
        fields=['COUNTRY','SUM']                         # 指定要读入的字段
        with arcpy.da.SearchCursor(out_table, fields) as rows:
            for row in rows:
                print('{0},{1}'.format(row[0],row[1]))
                f.write('{0},{1}'.format(row[0],row[1]))
                f.write('\n')                            # 写完一行后换行
    f.close()

    startYear += 1