有没有办法在 R 中将 MODIS 切片作为 lat/long 的函数？
IS there a way to get MODIS tiles as a function of lat/long in R?
我正在寻找一种方法来自动 return MODIS 瓦片（例如 h12 v4）作为经度和纬度的函数。
我希望 MODIS 包可以使用 getTile(***) 来完成。当我使用 getTile() 时，我看到有一个很好的 GUI，但我正在寻找一种自动方法作为 long/lat 的函数。
谢谢，
杰夫
好建议！这目前包含在 MODIS 'develop' 分支中，它肯定会成为未来包版本的一部分（从 1.1.4 ).在下一个 CRAN 版本之前，只需从 GitHub 安装当前的 'develop' 版本，看看这是否适合你：
# 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点输入。
我正在寻找一种方法来自动 return MODIS 瓦片（例如 h12 v4）作为经度和纬度的函数。
我希望 MODIS 包可以使用 getTile(***) 来完成。当我使用 getTile() 时，我看到有一个很好的 GUI，但我正在寻找一种自动方法作为 long/lat 的函数。
谢谢， 杰夫
好建议！这目前包含在 MODIS 'develop' 分支中，它肯定会成为未来包版本的一部分（从 1.1.4 ).在下一个 CRAN 版本之前，只需从 GitHub 安装当前的 'develop' 版本，看看这是否适合你：
# 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点输入。