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

Question

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

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

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

有人知道这在 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)

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

Overlap (or convert) a raster in shapefile with R software

tiff

r

converters

r-raster