前言

由于科研工作需要使用到国产的高分影像数据，并且是时序的数据，因此需要大规模的处理GF类型的数据。目前在网上很少看到有可利用的批量处理高分影像数据的工具，大多数是使用ENVI工具直接处理，不适合科研工作。因此写一份高分影像批处理的学习流程，以供自己回顾和大家参考。

高分数据

高分数据简明

高分专项目前已经发射了7颗遥感卫星。
高分一号卫星，高分辨率相机全色分辨率2米、多光谱分辨率8米，并携带16米分辨率幅宽800千米的宽幅相机，是典型的高时间分辨率（宽覆盖）遥感卫星；
高分二号卫星，携带了全色分辨率优于1米、多光谱分辨率优于4米的高分辨率相机，实现了较高的空间分辨率；
高分三号卫星，携带C波段多极化合成孔径雷达的，是中国第一颗民用合成孔径雷达卫星，空间分辨率最高可达1米，主要用于对海洋目标的监测；
高分四号卫星，是静止轨道光学遥感卫星，携带全色分辨率50米的凝视相机，在静止轨道上对7000千米X7000千米的区域进行实时观测；
高分五号卫星，是中国第一颗民用高光谱观测卫星，高光谱相机分辨率能力可达10米，它还将携带气溶胶探测仪、温室气体监测仪等探测设备；
高分六号卫星，属于高分一号卫星的后继，它仍是一颗全色2米分辨率、多光谱8米分辨率的普查卫星；
高分七号卫星，属于立体测绘卫星，将为国家测绘地理信息局、国家基础地理信息中心等部门提供高分辨率的空间立体测绘数据。

目前整体的思路以我拿到的高分数据为例子作学习，其他类型大家可以共同补充。

高分手册

首先作为用户，我们可以先查阅高分中心的用户手册，其用户手册的网站如下：高分应用综合信息服务共享平台 (gaofenplatform.com)http://gaofenplatform.com/contents/558/1303.html

网址指向的是GF6型的用户手册，如需其他类型切换即可。根据手册内容，我发现我拿到的数据为L1型数据，也就是只是经过相对辐射校正（即消除的传感器之间的相互差异），并没有经过大气校正、辐射定标以及几何校正。因此我们影像处理的任务就包括上述的内容，最终得到L2级用户数据。

但是，很遗憾，手册没有告诉我们如何处理L1级数据，需要自己去探索。

高分数据目录

打开一张GF影像的压缩包，发现其中具有一下内容：

不同类型的GF可能略有不同，基本包括：
辐射校正数据文件（*.tiff），其中宽幅相机分三块存储
RPC参数文件（*.rpb）
景元信息文件（*.xml）
浏览图（*.jpg）
拇指图（*_thumb.jpg）

标红的信息是使用的重点，其他信息无关紧要。

整体思路

初步考虑的思路如下：

第一步：将不同格式的高分数据组织好，便于主程序调用

第二步：对影像进行解压

第三步：对影像进行辐射定标

第四步：对影像进行大气校正

第五步：对影像进行几何校正

第六步：保存新文件，删除解压文件

高分数据的整理与解压

高分数据的整理

由于高分数据的格式可能没有对齐，可能有压缩包，也可能没有压缩包，因此该步骤旨在对其格式，得到一个压缩包路径的列表

    primaryFilelist = iglob(os.path.join(Input_dir,"GF*L1*[1-9]"))gzFilelist = glob(os.path.join(Input_dir,"GF*L1*.gz"))for item in primaryFilelist:if os.path.isdir(item):gzpath = glob(os.path.join(Input_dir,item,"GF*L1*.gz"))if len(gzpath)==1:gzFilelist.append(gzpath[0])else:continueelse:continue

这样就得到 Input_dir 路径下所有的压缩文件，存储在gzFilelist中.

解压缩

# 解压缩原始文件
def untar(fname, dirs):# print("文件路径",fname)try:t = tarfile.open(fname)except Exception as e:print("文件%s打开失败" % fname)t.extractall(path=dirs)

利用上述两段代码稍作组织即可实现批量解压缩的功能。

高分影像批处理第一回——数据格式分析与整理相关推荐

高分影像批处理第二回——辐射定标与大气校正
辐射定标遥感影像的辐射定标就是将传感器接收的无效的DN(digital number)值转化为我们所需要的辐射亮度值或者大气表观反射率.要达到这个目的,我们需要知道传感器的增益和偏移值.高分影像的g ...
基于ArcGIS与高分影像进行绿地变化分析
1. 需求现在有某区域的高分影像和对应范围的土地利用现状数据,现在的耕地面积比二调大了很多,需要把没有备案的土地清查出来.简单的说就是要找出哪些绿地是新增的或者由其他用地类型转化而来的. 2.总体技 ...
第一回开篇 D3D渲染流程简介
第一回开篇 D3D渲染流程简介 http://developer.178.com/201004/65293187592.html 开发这个3D engine已经两年半了,从06年8月刚开始统计的4万 ...
如何借助现有股票量化交易平台编写策略和回测分析
每个交易日的股票都会上涨或者下跌,在这个过程中笔者们偶尔会想针对部分股票进行股价的涨跌幅进行监控,或者自动进行交易,在这个需求前提下,现有券商.股票分析软件都会带有机器人自动交易策略功能,大部分都需要 ...
批量处理|基于ENVI的国产高分影像批量融合工具
从本文开始,介绍ENVI软件常用的遥感影像处理功能,并采用C#+IDL混合编程方式,实现ENVI常用功能的批量处理,对应的ENVI采用5.3版本,其他版本不保证能正常使用,尤其是5.3以下版本. 1. ...
C#《原Csharp》第一回买书未成炁自生惶惶回屋遇老翁
这是我在博客园首创的一个随笔性质的小说,旨在用还算风趣通俗的话来传达C#的魅力,在写的过程中发现有不少人有"转载"我的文章且未标明出处(是不是因为咕咕太久了?),所以本着打不过就加 ...
用Python编程借助现有量化平台编写股票交易策略和回测分析
一.简介接下来笔者会给大家介绍股票交易中的常见几种交易策略实现思路和源码编写过程,如果大家听说过量化交易这个词语的话,对其中的交易策略或许了解过,大概意思就是在股票.加密货币或者金融衍生物在价格的波 ...
第一课逆向分析基础知识
看雪软件安全论坛 > 初学者园地 > 『伴你成长』 > 『资料导航』加密与解密解密分析入门基础知识欢迎您,youcou 您上次访问的时间是: 2009-09-23,23:53 ...
[你必须知道的css系列]第一回：丰富的利器2：CSS选择符之子选择符、相邻选择符...
对了,接下来要讲的属性选择符,相邻选择符,子对象选择符可能大家稍微有点陌生了,这当然也是有原因的,因为IE6及以下的浏览器并不支持这几个选择符,而大多数从事这方面工作的技术人员,多数时候还是主要考虑 ...
Python 优化第一步: 性能分析实践使用cporfile+gprof2dot可视化
拿来主义: python -m cProfile -o profile.pstats to_profile.py gprof2dot -f pstats profile.pstats |dot -Tp ...

高分影像批处理第一回——数据格式分析与整理

前言