У меня есть виртуальная растровая мозаика (.vrt) с четырьмя полосами, показывающими 4 типа землепользования.У меня также есть шейп-файл с 520 точками, который я хочу нарисовать в 2-километровом буфере и извлечь общую площадь каждого типа землепользования на буфер.
До сих пор я использовал следующий код:
ek_sample<- raster("G:/Heather_tiff_files/use_S/raster_mosaic_3006.vrt")
sample_vec<- st_read("G:/Heather_tiff_files/bat_final_6_sweref.shp")
ex <- extract(ek_sample, sample_vec, fun=sum ,buffer = 2000, df=TRUE)
Это просто создает кадр данных с общим количеством пикселей, но не числом на полосу.Как я могу написать, что я хочу сумму на группу?
Большое спасибо,
Хизер.