ENVI拓展插件—焦糖计划：预处理高分系列卫星光学数据（更新至20.03.07）

近期做了理想状态的更新并改成了非插件的白鸽，需要在IDL环境下调用sav或源码来处理数据，目前已发布v1.0正式版本，如有需要敬请移步至https://blog.csdn.net/desertsTsung/article/details/104456051

 

说明：焦糖计划（Project JiaoTang）支持一步预处理多景高分一/二/四/六号数据，包括GF1-PMS/WFV，GF2-PMS，GF4-PMI，GF6-PMS/WFV。

特点：只用输入原始压缩文件即可得到预处理的结果，中间不用人为操作，也不用人为删除中间文件。

预处理流程：解压，正射校正，（配准），（裁剪），辐射定标，融合（GF1/2/6-PMS），（快速大气校正）。

版本要求：ENVI5.3/IDL8.5及更高。

最新下载地址：蓝奏云免费下载（0522版本），CSDN积分下载（0522版本）

问题汇总：1）拓展最好在本地使用，在服务器运行概率性出现未知错误；2）处理过程会生成临时文件，如果处理中断可以自行删除临时文件夹中envi开头的文件，默认为路径是C:\Users\用户名\AppData\Local\Temp；3）如果默认的临时文件夹剩余空间过小，可以在ENVI中设置其他路径，方法：左上角File->Preferences->Directories->Temporary Directories；

使用简介：安装拓展后（下载文件含安装说明），重启ENVI，该拓展位于Toolbox / Extensions / Optical GaoFen Auto Process。参数TGZ Files（必填）：原始未经解压文件；参数DEM：用于正射校正的DEM；参数Shapefile：用于裁剪的矢量；参数Registration：是否配准；参数Fusion Method：融合方法；参数QUAC：是否快速大气校正；参数Divide 10k：是否在QUAC基础上除以10000；参数Display Result：是否在屏幕显示结果；参数Output：输出路径。推荐中高配置电脑。

单击界面左下角问号图标可以查看具体的填写说明：

处理结果弹窗（GF6-PMS，裁剪区域为局部邓州市）：

局部处理结果（5%拉伸，标准假彩色）：

 

更新日志：

---

由于准备考研，预计2020年上半年才会第十二次更新（如果上岸的话），以下是预更新内容：

1.大量代码重构：1）可能会有一个非插件版的程序，如果可以的话，两个版本尽量集成到一个sav中；2）调整输入为文件夹输入，批处理更方便；3）直接从文件名读取数据信息，不再依赖xml文档对象；4）小概率会支持TOA的生产；5）优化保存结果的逻辑，之前版本的思路是先输出再复制，如果结果文件过大，会多消耗时间和计算机性能；6）优化删除临时文件的逻辑，之前版本的思路是按照参数删除，可能后续会改为按照临时文件的时间属性去删除；7）生成time.log格式的处理日志[感觉这个很有必要]；8）非插件版可能会有非常简单的进度显示。

2.数据相关：1）移除高分四号；2）新增高分一号BCD星；3）更准确的高分六号光谱信息；4）如果赶得上的话，会有2019年的定标系数；5）准备将中心波长和定标系数移出程序外，改为从内置txt(或JSON或xml)文件中读取，用户可以自行更改满足需求。

3.这一路上得到了很多师兄师姐的无私帮助，也有各位使用者的支持和信赖，与你们握手，感谢。

 

2019/05/22第十一次更新：（点击下载）

优化全色与多光谱分辨率为非整数倍的NNDiffuse融合效果。

 

2019/05/08第十次更新：（点击下载）

1.自主重构矢量裁剪部分的代码，移除掩膜操作；

2.精简解压部分代码。

 

2019/04/27第九次更新：（点击下载）

1.新增对高分六号的支持。对于三段式的WFV数据，会根据裁剪区域（如果有的话）进行选择性处理及拼接；

2.支持多文件批处理，输入与输出为n:n，输入与裁剪区域为n:1或n:n；

3.更改输入为原始tar.gz压缩文件；

4.修复shapefile进行裁剪出现位置不正确的问题，考虑到实际情况，删除对xml格式ROI文件的支持；

5.更改辐射定标逻辑，对于未公布系数的影像，采用最近的定标系数；

6.更改配准的逻辑顺序；

7.更改添加波长及其单位的时间段；

8.更改QUAC后除以10000的方式；

9.更改ENVI的COMMON声明为DEFSYSV；

10.删除处理完成后冗余的步骤提示；

11.彩蛋：做了一个小功能，可以通过错误的逻辑（即不选择quac并选择除以10000），来检查此插件的新版本。不过今年可能不会更新了；

12.源码移至GitHub：https://github.com/desertstsung/project_jiaotang。

 

2019/04/08第八次更新：（点击下载（内含源码））

1.新增自主选择项，优化体验；

2.帝都好冷。

 

2019/03/23第七次更新--紧急修复版本2：（点击下载（内含源码））

1.完全修复了不输入矢量会报错的情况；

2.其他代码结构优化调整。

3.加入的矢量文件必须全部在单景影像中。最近研究研究怎么在相交的情况下进行裁剪。

 

2019/03/21第七次更新--紧急修复版本：（点击下载（内含源码））

1.修复了不输入矢量会报错的情况；（感谢使用者“云中剑”的反馈）

2.给各位使用者带来了不便，敬请谅解。

 

2019/03/16第七次更新：（点击下载（内含源码））/*此版本存在BUG，不推荐下载*/

1.支持ESRI Shapefile格式矢量数据进行裁剪操作；

2.更新了较为完善的数据填写说明，点击界面左下角的问号图标可以查看；

3.修复矢量数据未关闭的情况；

4.界面添加版本号；

5.单独PRO/FUNC配置单独文件，便于后期更新；

6.将之前的两个sav文件合并为一个；

7.感谢每一位使用者的信赖和支持。

 

2019/02/01第六次更新：（点击下载（内含源码））

1.更新2018年定标系数；

2.完成后提示预处理步骤；

3.删除中间文件不再用ENVI打开后删除；

4.增加FUNCTION处理，增加主程序可读性；

5.更改XML文件结点的读取方式，不再每个结点单独一个FUNCTION；

 

2018/12/24第五次更新：（点击下载（内含源码））

1.融合方法选择栏回归，支持自定义两种融合方法：NNDiffusePanSharpening和GramSchmidtPanSharpening，分辨率不再降低。

2.愿你度过一个幸福祥和的平安夜。

 

2018/12/23第四次更新：（点击下载（内含源码））

1.上个版本中，如果选择NNDiffuse Pan Sharpening融合方法处理一整副影像（GF1-PMS和GF2-PMS），结果会因为文件过大而降低分辨率（一小幅范围则不会出现这种情况），所以此版本不再提供融合方法选择，默认为Gram-Schmidt Pan Sharpening融合方法，虽然融合速度会比前者慢很多很多，但是可以保证结果分辨率与全色一致。

2.更新辐射定标机制，对定标系数未予以公布的数据不进行定标操作，上个版本的思路是将定标系数设置为1和0进行定标。

3.支持自定义输入xml格式ROI文件进行裁剪，减少工作量。需要预先在Data Manager中打开ROI数据。

4.在最后快速大气校正的基础上除以10000，使结果为地表反射率。

 

2018/12/06第三次更新：（点击下载（内含源码））

1.更改名字为焦糖计划（Project JiaoTang），源于为食品着色的焦糖。

2.更改传感器和年份读取为从XML文件读取，更合理。

3.支持自定义融合方法（GS和NNDiffuse），个人感觉两种方法各有所长，前者在纹理方面与全色最接近，后者在色调方面与多光谱最接近。

4.修复DEM数据未关闭的情况。

 

2018/12/05第二次更新：（点击下载（包含源码））

1.支持高分一号，高分二号，高分四号的所有光学遥感数据，即高分一号二号四号的PMS传感器和高分一号的WFV传感器。如果是高分一号WFV数据或高分四号PMS数据，则无需输入全色影像（PAN Image）。高分四号IRS数据预处理步骤较少，未作考虑。

2.最低版本要求为5.3，不再是5.3.1。

3.上个版本主攻高分二号，忽略了高分一号的中心波长，本次更新进行了添加（包括高分四号）。中心波长采用Palmer提出的矩方法求得（下式），下载压缩包里也包含每个传感器的中心波长和每年的定标系数。

4.自动获取年份和传感器类型，不再手动输入。

5.输入影像必须为原始未更改影像，不能在File Select面板中进行Spatial Subset裁剪，否则会报错找不到影像的年份和传感器，从而终止进程。但是可以先对原始影像进行裁剪，再用裁剪结果替换原始影像就可以达到预处理某一区域的效果了。

6.读取XML文件的方法来源于麻辣GIS，表示感谢。

 

 

 

 

 

不知道好久的第一次更新：

之前的一篇博客写的太匆忙了，有很多需要修改或新增的，于是有了新版本（下载地址），要求最低版本为ENVI5.3.1。

建议电脑内存为8G以上，4G内存运行可能会现卡在配准界面的情况

为什么最低版本是ENVI5.3.1：相较于之前赤裸的源码，这次加入了对Metadata里一些属性的判断（HasTag函数判断，ENVI5.3.1引入），添加了背景忽略值，既能使NNDiffuse融合结果得到提升，也可以使最后的结果更加美观。

预处理过程：正射校正，辐射率定标，全色和多光谱的配准，影像融合，快速大气校正

更改地方：（1）由于2018年的定标系数资源卫星应用中心还未公布，所以只支持2017年之前的高分一号和高分二号PMS传感器数据，后续随官网更新；（2）融合前将多光谱的储存方式改为了BIL，使融合效率提高；（3）更改融合方法为NNDiffusePanSharpening，并和快速大气校正交换顺序；（4）对临时文件进行删除；（5）更新为ENVI拓展，交互界面更人性化。