r - 有没有办法将 MODIS 瓦片作为 R 中纬度/经度的函数？

Question

我正在寻找一种方法来自动返回 MODIS 瓦片（例如 h12 v4）作为经度和纬度的函数。

我希望 MODIS 包可以使用 getTile(***) 来完成。当我使用 getTile() 时，我看到有一个不错的 GUI，但我正在寻找一种自动方法作为 long/lat 的函数。

谢谢，杰夫

Answer 1

一个好建议！这目前包含在MODIS 'develop' 分支中，它肯定会成为未来软件包版本的一部分（从1.1.4开始）。在下一个 CRAN 版本之前，只需从 GitHub 安装当前的“开发”版本，看看这是否适合您：

# devtools::install_github("MatMatt/MODIS", ref = "develop")
library(MODIS)

## set point coordinates (taken from https://en.wikipedia.org/wiki/Frankfurt)
dat = data.frame(lon = 8.682222, lat = 50.110556, loc = "Frankfurt am Main")
pts = sf::st_as_sf(dat, coords = c("lon", "lat"), crs = 4326)

## get overlying tile
tls = getTile(pts)
tls@tile
# [1] "h18v03"

不过，一个重要的警告：目前，这只适用于sp（即Spatial*）和sf点输入。