python处理遥感图像_科学网—python处理Landsat系列影像的一些总结（1） - 温兆飞的博文...

1、Landsat数据(.tar.gz)的自动解压脚本# -*- coding: utf-8 -*-

"""

Created on Wed Mar 21 14:10:21 2018

@author: wenzhaofei

"""

import tarfile    #网上下的Landsat数据都是压缩包(.tar.gz)，因此需调用tarfile包(其他类似的解压缩包也可以)

folder = 'g:/test/'    #存放数据的文件夹

indata = 'LT51240392010355BKT00.tar.gz'    #需被解压的影像压缩包

tar = tarfile.open(folder+indata) # 打开影像压缩包

names = tar.getnames()    #获取影像压缩包里包含的文件名

for name in names:

tar.extract(name,path=folder) #开始解压(这里解压的是所有文件，也可以只解压其中的某一个文件，指定文件名就好)

tar.close()

2、Landsat波段影像等栅格数据的剪切(裁剪)

# -*- coding: utf-8 -*-

"""

Created on Wed Mar 21 14:10:21 2018

@author: wenzhaofei

"""

import subprocess

from osgeo import ogr    #需先按照gdal包(在这里:https://www.lfd.uci.edu/~gohlke/pythonlibs/,搜索"gdal")

folder = 'g:/test/'    #存放数据的文件夹

inraster = folder+'LT51240392010355BKT00_B2.TIF'    #需被剪切的栅格数据路径

outraster = inraster.replace('.TIF', '_cliped.TIF') #被剪切的栅格数据存放路径

inshape = folder+'clipping_area.shp'    #剪切的范围

subprocess.call(['gdalwarp', '-srcnodata', '-dstnodata', '-cutline', \

inshape, inraster, outraster, '-crop_to_cutline', '-overwrite'])

#这里，'-cutline'是关键参数，'-crop_to_cutline'可以让输出的范围为剪切范围的外包矩形，其他的参数参考官方解释

(http://www.gdal.org/gdalwarp.html)。

转载本文请联系原作者获取授权，同时请注明本文来自温兆飞科学网博客。

链接地址：http://blog.sciencenet.cn/blog-365459-1105168.html

上一篇：猜想——构建民航客货机遥感体系

下一篇：暴力解决"应用程序无法启动，因为应用程序的并行配置不正确。"