R: Средневзвешенное значение равно 0 для буферов, в которых нет данных

Я новичок в R, и мне трудно разобраться в этой проблеме, поэтому буду благодарен за любую помощь. У меня есть ежедневные растровые данные о затоплении, которые я использую для расчета средневзвешенного значения затопления в буферах, созданных вокруг определенных мест. Я очистил растры в QGIS и использую R для расчета средневзвешенного значения с помощью следующих строк:

library(sp)
library(raster)
library(ggplot2)
library(sf)

     # Load the raster file into R
      flooding_raster_data <- raster(raster_file)
      
      # Print information or perform further processing as needed
      print(paste("Loaded raster file:", raster_file))
      
      # Extract values from the raster that fall within the buffer polygons
      raster_values_in_buffer_weight <- extract(flooding_raster_data, buffer, fun = mean, na.rm = TRUE, weights=TRUE)
      
      # Create a new buffer object with the extracted floodshare values
      buffer_with_mean <- st_as_sf(buffer)
      buffer_with_mean$floodshare <- raster_values_in_buffer_weight
      
      # Write the sf object to a shapefile
      st_write(buffer_with_mean, output_file_path, append = FALSE)
      rm(flooding_raster_data, buffer_with_mean, raster_values_in_buffer_weight, raster_file)

Моя проблема заключается в том, что для растра, где нет данных (закодированных как отсутствие данных в QGIS), среднее значение буферов равно 0. Я проверил это по зональному среднему значению в QGIS, и те же самые буферы имеют зональное среднее значение NULL в QGIS. Я изо всех сил пытаюсь понять, почему это произойдет в R, поскольку я думаю, что R будет иметь среднее значение NA или пропустить, поэтому, если кто-то знает, почему, пожалуйста, дайте мне знать и как я могу это исправить. Опять же, я новичок в R, поэтому, возможно, я упустил что-то очень очевидное.

Я попытался сделать na.rm = FALSE, потому что подумал, что, возможно, проблема была в том, как я обрабатывал значения NA в коде R, но это просто сделало все взвешенное среднее значение NULL, поэтому я думаю, что не понимаю, как R обрабатывает недостающие данные. .

РЕДАКТИРОВАТЬ

Я пошел проверить документацию по извлечению R и взял их пример кода, но изменил значения растра на NA, и это дает средневзвешенное значение, равное 0, для значений извлечения в виде полигонов, если все значения растра равны NA. Почему он не выдает NA для средневзвешенного значения, хотя все остальное дает ответ как NA или NaN. Я добавил код репликации ниже.

###############################
# extract values with lines
###############################
r <- raster(ncol=36, nrow=18, vals=NA)
cds1 <- rbind(c(-50,0), c(0,60), c(40,5), c(15,-45), c(-10,-25))
cds2 <- rbind(c(80,20), c(140,60), c(160,0), c(140,-55))
lines <- spLines(cds1, cds2)

extract(r, lines)

###############################
# extract values with polygons
###############################
cds1 <- rbind(c(-180,-20), c(-160,5), c(-60, 0), c(-160,-60), c(-180,-20))
cds2 <- rbind(c(80,0), c(100,60), c(120,0), c(120,-55), c(80,0))
polys <- spPolygons(cds1, cds2)

v <- extract(r, polys)
# mean for each polygon
mean_polyg <- unlist(lapply(v, function(x) if (!is.null(x)) mean(x, na.rm=TRUE) else NA ))
print(mean_polyg)

# weighted mean
v <- extract(r, polys, weights=TRUE, fun=mean, na.rm = TRUE)
print(v)