R 陆地卫星图像请求
我有没有办法在截止日期前自动向美国地质勘探局提出请求,以获取陆地卫星8号的图像?实际上,我是通过EarthXPlorer web界面向USGD提交批量请求来手动完成这项工作的。此外,是否有R代码可用于从Landsat 8图像中提取信息?(我知道已经有陆地卫星软件包了)我不知道下载陆地卫星图像的任何代码或软件包。但是为了获得关于它们的信息,我在plotKML包中使用read.metadata函数 下面是一个示例,我在其中读取.xml数据中的方位角R 陆地卫星图像请求,r,download,landsat,R,Download,Landsat,我有没有办法在截止日期前自动向美国地质勘探局提出请求,以获取陆地卫星8号的图像?实际上,我是通过EarthXPlorer web界面向USGD提交批量请求来手动完成这项工作的。此外,是否有R代码可用于从Landsat 8图像中提取信息?(我知道已经有陆地卫星软件包了)我不知道下载陆地卫星图像的任何代码或软件包。但是为了获得关于它们的信息,我在plotKML包中使用read.metadata函数 下面是一个示例,我在其中读取.xml数据中的方位角 library(plotKML) meta<
library(plotKML)
meta<-read.metadata(paste(DataPath,"LC82330852015268LGN00.xml",sep=""))
azi=as.numeric(meta$value[9])
库(plotKML)
meta你从美国地质调查局查过了吗?它允许您请求元数据、构建订单等。他们的repo中有一个可供您了解的功能
您可能还对获取采集元数据的软件包感兴趣。您应该查看rLandsat
软件包,该软件包提供了使用ESPA api和sat api在R中请求陆地卫星图像的简便方法
包中的一个示例
# 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())