有没有办法让MODIS瓦片在R中作为lat / long的函数？

Question

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

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

谢谢，杰夫

Answer 1

一个很好的建议！这目前包含在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点输入。