我有一个虚拟的栅格镶嵌（.vrt），其中有四个波段显示4种土地利用类型。我也有shapepoint文件，该文件有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)
这只是产生一个具有像素总数而不是每个波段数目的数据帧。我该如何编码每个频段我想要的金额？
非常感谢，
希瑟。