1 前言
众所周知,Landsat8数据可在USGS以及地理空间数据云中下载得到。但是笔者经试验也发现现有下载资源的一些不足。
(1)USGS速度较慢。笔者网络环境电信100M,但在官网下载landsat8数据,速度一直维持在15KB左右。
(2)地理空间数据云。这个平台大家基本都知道,但是下载landsat8时还遇到一个大问题,可能我们需要的数据可以检索到,但是不提供下载。
2 如何使用
遇到以上说的情况,笔者经一番搜索,找到了另外的方法下载。
Landsat8有镜像数据在Google Storage上,所以国外老哥用python写了个脚本,可以直接从Google服务器上拿到数据。
2.1 必备环境
理论上来说Windows、Linux、Mac都可以。但是笔者试验过后,发现这个脚本所需的依赖比较难装,Win和Mac环境都失败了,只有Linux成功(Ubuntu 16.04)。
另外因为是从google storage下载,所以还需要一个FQ环境。笔者的科学上网是在路由器里设置好的,全局FQ。读者如果要用,请自行准备个科学上网的环境。并且,如果不是我这种路由器直接fq的,可能还会有坑,就是普通的FQ是pc端开代理,但是在Linux的shell中也开代理这个还得费点功夫。所以写到这里,发现我这个方法还是有点坑的。后续看看有没有办法改进吧。
推荐环境:
Ubuntu 16.04 系统
Python环境
科学上网环境
2.2 下载工具的部署
部署很简单,几行代码搞定。
第一行代码
‘sudo apt-get install python-pip python-numpy python-scipy libgdal-dev libatlas-base-dev gfortran libfreetype6-dev
’
第二行代码
‘pip install landsat-util’
如果遇到安装失败的问题,请仔细看下报错代码,然后搜索下,安装相关依赖就可以解决,具体因为Linux使用不同的源、不同的版本进行维护,很难说会有一种特定解决办法。
2.3 工具的使用
2.3.1 检索
本工具可以使用多种检索方式。
(1)用path &row搜索
示例代码如下:
‘landsat search --cloud 4 --start "january 1 2014" --end "january 10 2014" -p 009,045’
(2)经纬度查询
‘landsat search --lat 38.9004204 --lon -77.0237117’
(3)经纬度带geojson输出的查询
‘landsat search --lat 38.9004204 --lon -77.0237117 --geojson | geojosnio’
第三种我没用过,大家可以自己尝试。
2.3.2 下载
检索完毕就可以下载。下载方式是直接给出他的某一景ID,笔者推荐大家先去数据云上在地图查询出自己想要的某一景ID然后再用。
现在以我需要的ID:LC80090452014008LGN00为例。
那么下载就是这行代码:
‘landsat download LC80090452014008LGN00’
还可以下载指定波段,如:
‘landsat download LC80090452014008LGN00 --bands 432’
同时下载多景:
‘landsat download LC80090452014008LGN00 LC80090452015008LGN00 LC80090452013008LGN00’
下载好的影像存放在Home\Landsat文件夹下。
经本人试用,下载速度从90K会一直上升到800K/s,当然这也取决你的科学上网环境,大概半小时左右下好一景的数据。如果说你不想用数据云或者usgs那么这个工具是很好的选择了。
2.3.3 影像处理
本工具还带了个image processing的模块,看官网介绍主要是做pansharpen用。我是在envi里直接做了。有兴趣读者可以自行尝试下。
最后给出本工具的文档地址,更多功能在此探索:
‘https://pythonhosted.org/landsat-util’
最后
姚尧同学的投稿是没有图的,如果说非要加上个图的话,那只能是我送支付宝红包的图了,这个图可以免费领红包,最高99元,每天领一次
,当然你领一个,我也会得到一个哦,欢迎扫一扫哦!