python语言基于geopandas和rasterio包进行矢量提取栅格值操作

导包

import geopandas as gpd
from osgeo import gdal
import rasterio as rio
import rasterio.mask
import pandas as pd
import numpy as np
from tqdm import tqdm
import os

读取数据

shpdatafile=r'C:\Users\Administrator\Desktop\数据采样降低空间自相关\100m/100m_yhln.shp'
pocl_h_path = r'D:\000-毕业论文改\提取特征\点云特征\按多边形裁剪_根据地面点归一化_IntensityMetrics'

shpdata=gpd.read_file(shpdatafile)

栅格数据存储如下，所有的数据都在一个文件夹里：

矢量提取栅格主代码

基本逻辑是使用for循环遍历矢量数据内的每一个矢量，使用rasterio中的mask.mask函数通过每一个矢量提取对应栅格的像元均值。最终生成结果为csv文件，如果想要生成矢量文件可以使用geopandas的merge函数将csv文件合并在原矢量上。

ee = []
val_n