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(