使用Python中的点从栅格中提取值

Question

是否有一个程序包，我可以在其中使用一组栅格文件并根据我在shapefile中的点从栅格文件中提取值。

我在Jupyter笔记本中加载了shapefile和栅格，但是没有找到可以使用上载点从上载栅格中提取数据的功能

Answer 1

您可以在crop程序包中使用raster函数来创建一个新的栅格对象，该对象限于多边形的大小。

这是raster包的帮助示例。

library(sp)
library(raster)
r <- raster(nrow=45, ncol=90)
x <- runif(ncell(r)) * 360 - 180
y <- runif(ncell(r)) * 180 - 90
xy <- cbind(x, y)
r0 <- rasterize(xy, r)
b <- as(extent(0, 20, 10, 30), 'SpatialPolygons')
crs(b) <- crs(r)
rb <- crop(r0, b)

如果形状很多，虽然返回列表而不是空间对象，但是over函数可能会有所帮助。