注册 登录
  • 注册时,本站名称为:RGB空格3S博客,注意中间的空格。
  • 2018.10.14本人利用业余时间接遥感数据处理与编程小活,欢迎骚扰,QQ:853060844
  • 2018.8.14本站全面接入google广告
  • 本人闲暇时间(周末及其他闲暇时间)接遥感图像处理与IDL编程小活,欢迎骚扰,qq:853060844
  • 2017.2.14今天收到45条恶意评论(全是外文),故评论时请填写必要信息,匿名评论全部拉黑,迫不得已而为之
  • 2017.1.27,2017年春节及至,我谨代表本人祝大家新春快乐,本人年终总结文章请访问:2016年终总结
  • 为防止恶意转载,本站全面禁止复制,并添加图片水印:RGB 3S博客www.ixxin.cn。
  • 本站正式更名为RGB 3S博客,本站将撤消所有非3S内容,其将转移到新博客江湖时代

Python脚本速下载Landsat8

Python admin 3753次浏览 已收录 0个评论
[隐藏]

快速下载Landsat 8 数据

姚尧 2017.12.11 QQ:2841460972

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元,每天领一次,当然你领一个,我也会得到一个哦,欢迎扫一扫哦!


xxin blog , 版权所有丨如未注明 , 均为原创丨本网站采用BY-NC-SA协议进行授权 , 转载请注明Python脚本速下载Landsat8
喜欢 (3)
支付宝[853060844@qq.com]
分享 (0)
admin
关于作者:
坐标山科大遥感系小鲜肉一枚。
发表我的评论
取消评论
表情 贴图 加粗 删除线 居中 斜体 签到

Hi,您需要填写昵称和邮箱!

  • 昵称 (必填)
  • 邮箱 (必填)
  • 网址