Question

有什么方法可以按日期自动向美国地质调查局提出要求获取Landsat 8图像的请求吗？实际上，我是通过Earthexplorer Web界面向USGD提交批量请求来手动执行此操作。另外，我可以使用R代码从Landsat 8图像中提取信息吗？（我知道已经有Landsat包了）

Answer 1

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

这是一个例子，我在那里读取.xml数据中的方位角

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

希望这对你有用

Answer 2

你有没有看过美国地质勘探局的ESPA API？它允许您请求元数据，构建订单等。他们的回购中有一个jupyter notebook demo，可以让您对功能有所了解。

您可能还会对USGS/EROS Inventory Service API感兴趣，以获取有关下载采集的元数据。

Answer 3

您应该检出rLandsat软件包https://github.com/socialcopsdev/rLandsat，该软件包提供了一种使用ESPA API和sat-api来请求R中的Landat图像的简便方法。

包装中的示例

# 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())

Landsat图片请求

3 个答案: