Есть ли способ получить плитки MODIS как функцию lat/long в R?

Question

Я ищу способ автоматического возврата плитки MODIS (например, h12 v4) в зависимости от долготы и широты.

Я надеялся, что пакет MODIS сможет сделать это, используя getTile(***). Я вижу хороший графический интерфейс, когда я использую getTile(), но я ищу автоматический метод как функцию long / lat.

спасибо Джефф

Answer 1

Хорошее предложение! В настоящее время он включен в ветку MODIS "разработка" и определенно будет частью будущих версий пакетов (начиная с 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 точечные входы.