1. 裁剪栅格(cropping)
- 统一投影坐标系
vector = st_transform(vector,projection(raster))
- 设定矢量数据范围外栅格数据的值为
NA
raster_masked= mask(raster,as(vector,'Spatial')
2. 栅格数据提取
- 点
vector$feature=raster::extract(raster,as(vector,'Spatial')
vector$buffer_feature=raster::extract(raster,