Наложение (или преобразование) растра в шейп-файл с помощью программного обеспечения R

Question

У меня есть файл .TIFF (растр) для Латинской Америки и файл .shp (шейп-файл) для бразильских муниципалитетов, и я хочу перекрыть оба, но я хочу суммировать всю информацию, имеющуюся в растре.

Вмой случай, у меня есть файл TIFF о загрязнении света. Каждая точка в растре представляет одну световую информацию в диапазоне от 0 до 63. Шейп-файл имеет только геометрию.

В конце я хочу, чтобы каждая информация в муниципалитете возвращала сумму (или среднее значение) всех растровых точек.

Кто-то знает, что это возможно в программном обеспечении R?

Большое спасибо

Я не могу привести пример, потому что файл TIFF слишком большой, но я отправляю изображение

Я хочу наложить карту как это , но с меньшими размерами

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)