R 范围不重叠

r

由森林火灾数据组成，土地覆盖是包含栖息地指标的

modis

数据，而

bcr

是地理区域。在实现此代码时

fire <- bcr %>% 
     st_buffer(dist = 10000) %>% 
     st_transform(crs = projection(r)) %>% 
     crop(r, .) %>% 
     projectRaster(crs = projection(landcover))

#error comes from crop(r, .)

您可以从此处下载必要的文件：

库（sf）
图书馆（光栅）
图书馆（ExactExtractor）
图书馆（绿色）
图书馆（tidyverse）
#解决命名空间冲突
选择%
设置名称（土地覆盖）

library(MODIS)

hdf <- getHdf(HdfName = 'MCD64A1.A2000306.h12v11.006.2017012010432.hdf',forceDownload=T)

gdalUtils::get_subdatasets(hdf)

sds <- gdalUtils::get_subdatasets(hdf)

r <- raster(sds[1])

library(sf)
library(raster)
library(exactextractr)
library(viridis)
library(tidyverse)
# resolve namespace conflicts
select <- dplyr::select
map <- purrr::map
projection <- raster::projection

bcr <- read_sf("data/gis-data.gpkg", "bcr") %>% filter(bcr_code == 27) %>% st_transform(crs = paste("+proj=sinu +lon_0=0 +x_0=0 +y_0=0",

landcover <- list.files("insert_folder_name_here", "^modis_mcd12q1_umd", 
                        full.names = TRUE) %>% 
             stack()
# label layers with year
landcover <- names(landcover) %>% 
             str_extract("(?<=modis_mcd12q1_umd_)[0-9]{4}") %>% 
             paste0("y", .) %>% 
             setNames(landcover, .)