使用 R 軟件在 shapefile 中重疊（或轉換）柵格

Question

我有一個用於拉丁美洲的 .TIFF 文件（柵格）和一個用於巴西市政當局的 .shp 文件（shapefile），我想將兩者重疊，但我想總結柵格擁有的所有信息。

就我而言，我有一個關於光污染的 TIFF 文件。 柵格中的每個點代表一個光信息，范圍為 0-63。 shapefile 只有幾何圖形。

最后，我希望市政當局的每條信息都返回所有柵格點的總和（或平均值）。

有人知道這在 R 軟件中是可能的嗎？

非常感謝

我無法舉例，因為 TIFF 文件太大，但我發送了圖像

我想像這樣覆蓋一個 map ，但尺寸更小

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)