0

有什么方法可以按日期自动向 USGS 发出请求以获取 Landsat 8 图像?实际上,我是通过 Earthexplorer Web 界面向 USGD 提交批量请求来手动执行此操作的。此外,是否有可用于从 Landsat 8 图像中提取信息的 R 代码?(我知道已经有 Landsat 包)

4

3 回答 3

0

你检查过 USGS的ESPA API吗?它允许您请求元数据、构建订单等。他们的 repo 中有一个jupyter notebook 演示,可以让您很好地了解这些功能。

您可能还对USGS/EROS Inventory Service API感兴趣,该 API用于获取有关获取的元数据以供下载。

于 2018-06-15T21:30:34.583 回答
0

我不知道下载 landsat 图像的任何代码或包。但是为了获取有关它们的信息,我使用了 plotKML 包中的 read.metadata 函数。

这是一个示例,我在 .xml 数据中读取方位角

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

希望这对你有用

于 2016-02-17T18:53:35.550 回答
0

您应该查看rLandsathttps://github.com/socialcopsdev/rLandsat,它提供了一种使用 ESPA API 和 sat-api 在 R 中请求 landsat 图像的简单方法。

包中的一个例子,

# 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())
于 2018-07-17T09:53:37.500 回答