把tif的值映射到shp中

问题描述

假如目前有一个（多个）tif文件和一个shp文件，想要把tif中每个像素的值集成到shp文件的新字段中。如果栅格和像素是一一对应的，问题将会变得非常简单：直接把每个像素的值映射到每个栅格中的新字段里即可。如果栅格和像素不是一一对应的，这里我们需要用到zonal_stats，把所有涉及到的像素值求平均值再映射到栅格里去。接下来是一个简单的示例，假设我们要把不同日期的NDVI值映射到小栅格里面去。

代码

import geopandas as gpd
import rasterio
from rasterstats import zonal_stats
import os

shapefile_path = 'your.shp'