python空间数据处理：part3-分区统计以及表格整理

对栅格数据需要按照县域分区计算：

分区统计以及表格整理参考以下代码：
利用python快速进行批量区域统计(zonal statistics)

Python 分区统计 rasterstats zonal statistics - 简书

修改自己的参数后：

import glob
import numpy as np
import rasterio
import rasterstats
import pandas as pd
Input_folder = 'E:/seasonyield_predict/datapre/SH/vpd_growth_cliped_tif'
data_list = glob.glob(Input_folder + '/*.tif')
outpath='E:/seasonyield_predict/datapre/SH/vpd_growth_csv/'

for i in data_list:
    stats=rasterstats.zonal_stats('E:/seasonyield_predict/datapre/region/clip4326/db.shp',i)
    df_out=pd.DataFrame()
    df_out=pd.DataFrame.from_dict(data=stats)
    out=outpath+i.split('\\')[-1].split('.')[0] +'.csv'
    df_out.to_csv(out, header=True, index_label='fid', encoding='gbk')
    print(out)

 当一个tif文件中有多个波段，只想提取其中的某个波段栅格进行分析：

import os
import glob
import numpy as np
import rasterio
import rasterstats
import pandas as pd
from osgeo import gdal, ogr
from affine import Affine
Input_folder = 'E:/test'
data_list = glob.glob(Input_folder + '/*.tif')
outpath='E:/test/'

for i in data_list:
    ###查看有几个波段
    # td=gdal.Open(data_list[i])
    # # for i in range(1, td.RasterCount+1):
    # #     band = td.GetRasterBand(i)
    # #     band_name = band.GetDescription()
    # #     print("Band", i, "name:", band_name)
    with rasterio.open(i) as src:
        band = src.read(1)
        transform = src.transform
        affine_transform = Affine(transform.a, transform.b, transform.c, transform.d, transform.e, transform.f)
    stats=rasterstats.zonal_stats('E:/seasonyield_predict/datapre/region/clip4326/db.shp',band,affine=affine_transform)
    df_out=pd.DataFrame()
    df_out=pd.DataFrame.from_dict(data=stats)
    out=outpath+i.split('\\')[-1].split('.')[0] +'.csv'
    df_out.to_csv(out, header=True, index_label='fid', encoding='gbk')
    print(out)

【感慨，感谢愿意分享自己学习成果的博主们！同时也非常感谢身边的同学们，帮我解决了很多问题！感谢淋过雨的同学帮我撑伞，自己也会尽力帮别人撑伞！】