python快速实现shp点提取栅格值 (两种方法对比)

大量快速值提取至点的方法

在遥感数据的处理中，常常需要实现的一个功能就是“提取指定shp点对应栅格图像的value值”。实现该功能可以使用ArcGIS中的提取工具-值提取至点。但遇到大量数据，众多实验区的数据时，手动操作着实会带来很大的麻烦。这里我个人改进总结了两种方法，可实现python快速提取shp点对应的栅格值方法。

rasterio+GeoPandas

先上图说明：

import matplotlib.pyplot as plt
import geopandas as gpd
import rasterio
from rasterio.plot import show
pointData=gpd.read_file(r'\ABBY.shp')		
raster=rasterio.open(r'\ABBY_CHM.TIF')		#CRS为UTM, EPSG:32610
pointData=pointData.to_crs(