Запрос изображений Landsat

Question

Можно ли каким-либо образом автоматически отправить запрос в USGS по дате для получения изображений Landsat 8? На самом деле, я делаю это вручную, подавая массовый запрос в USGD через веб-интерфейс Earthexplorer. Кроме того, есть ли код R, который я могу использовать для извлечения информации из изображений Landsat 8? (Я знаю, что пакет Landsat уже есть)

Answer 1

Вы проверили API ESPA от USGS? Он позволяет вам запрашивать метаданные, заказы на сборку и т. Д. В репозитории есть демо-версия ноутбука Jupyter, которая даст вам хорошее представление о возможностях.

Вас также может заинтересовать API службы инвентаризации USGS/EROS для получения метаданных о приобретениях для загрузки.

Answer 2

Я не знаю ни кода, ни пакета для загрузки изображений Landsat. Но чтобы получить информацию о них, я использую функцию read.metadata внутри пакета plotKML.

Вот пример, где я читаю азимут в.xml данных

library(plotKML)
meta<-read.metadata(paste(DataPath,"LC82330852015268LGN00.xml",sep=""))
azi=as.numeric(meta$value[9])

Надеюсь, это полезно для вас

Answer 3

Вы должны проверить rLandsat пакет, https://github.com/socialcopsdev/rLandsat, который предоставляет простой способ запрашивать изображения LandSat в R, используя API-интерфейсы ESPA и sat-api.

Пример из пакета,

# get all the product IDs for India, alternatively can define path and row
result = landsat_search(min_date = "2018-01-01", max_date = "2018-01-16", country = "India")

# inputting espa creds
espa_creds("yourusername", "yourpassword")

# getting available products
prods = espa_products(result$product_id)
prods = prods$master

# placing an espa order
result_order = espa_order(result$product_id, product = c("sr","sr_ndvi"),
                          projection = "lonlat",
                          order_note = "All India Jan 2018")
order_id = result_order$order_details$orderid

# getting order status
durl = espa_status(order_id = order_id, getSize = TRUE)
downurl = durl$order_details

# download; after the order is complete
landsat_download(download_url = downurl$product_dload_url, dest_file = getwd())