使用 R 软件在 shapefile 中重叠（或转换）栅格

Question

我有一个用于拉丁美洲的 .TIFF 文件（栅格）和一个用于巴西市政当局的 .shp 文件（shapefile），我想将两者重叠，但我想总结一个栅格所拥有的所有信息。

就我而言，我有一个关于光污染的 TIFF 文件。栅格中的每个点代表一个光信息，范围为 0-63。 shapefile 只有几何图形。

最后，我希望市政当局的每个信息都返回所有栅格点的总和（或平均值）。

有人知道这在 R 软件中是可能的吗？

非常感谢

我无法举例，因为 TIFF 文件太大，但我发送了 image

我想覆盖像this one 这样的地图，但尺寸更小

Answer 1

试试这个：

require(raster)

rs=raster("pathToTiff.tif")
shp=shapefile("pathToShapefile.shp")

##to get the sum of all raster values within in each administrative area:
##assuming that they are in the same CRS, if not in same CRS run this first - shp=spTransform(shp,crs(rs))
extract(rs,shp,fun=sum)