1

我正在尝试下载和可视化肯尼亚（以及其他类似国家）的 NDVI 数据。

为此，我检索了肯尼亚使用的边界框，osmdata并在从MODIStsp

library(osmdata)
library(MODIStsp)

kenya_bb <- osmdata::getbb("Kenya")

MODIStsp(gui             = FALSE,
         out_folder      = "Data",
         out_folder_mod  = "Data",
         selprod         = "Vegetation_Indexes_16Days_1Km (M*D13A2)",
         bbox            =  kenya_bb,
         bandsel         = "NDVI",
         user            = "mstp_test" ,
         password        = "MSTP_test_01",
         start_date      = "2020.06.01", 
         end_date        = "2020.06.15", 
         verbose         = FALSE)

然后我只是将下载的数据可视化，但似乎没有裁剪出肯尼亚

library(raster)
library(here)
library(ggplot2)

NDVI_raster <- raster(here::here("Data/VI_16Days_1Km_v6/NDVI/MOD13A2_NDVI_2020_161.tif"))

NDVI_df <- as.data.frame(NDVI_raster, xy = TRUE, na.rm = TRUE)
rownames(NDVI_df) <- c()

ggplot(data = NDVI_df,
       aes(x=x,y=y)) +
  geom_raster(aes(fill = MOD13A2_NDVI_2020_161))

任何关于如何裁剪单个国家的想法都值得赞赏

4

1 回答 1

0

这种行为的原因是参数的默认值spatmeth是"tiles"，因此使用参数start_x、end_x和（start_y以及end_y，如果没有明确指定，它们的默认值）指定的值。

从功能文档：

bbox numeric(4) 在 out_proj 坐标系中输出边界框 (xmin, ymin, xmax, ymax)。如果 spatmeth == "tiles" 则忽略，默认值：NULL

spatmeth字符 ["tiles" | "bbox" | "file"]，表示如何检索处理范围。如果是“tiles”，则使用指定的tiles（start_x ....）。如果是“文件”，则从 spafile 中指定的空间文件中检索范围。如果为“bbox”，则使用指定的边界框，默认：“tiles”

要使用kenya_bb范围，您必须：

  1. 添加spatmeth = "tiles";
  2. 正确bbox = kenya_bb到bbox = as.vector(kenya_bb)

此外，最近修复了一个影响此案例的错误。您还应该从 GitHub 重新安装包：

remotes::install_github("ropensci/MODIStsp")
于 2021-04-14T14:47:38.377 回答