引言:高分一号(GF-1)卫星搭载了两台2m分辨率全色/8m分辨率多光谱相机(PMS),四台16m分辨率多光谱相机(WFV)。
处理流程:

- 全色图像很难进行大气校正,所以一般在定量遥感中不使用全色图像。
- 本流程中包括了大气校正,如果不涉及到遥感定量方面的应用,可不进行大气校正,也就是只包括正射校正和图像融合。
(一)先安装两个需要的插件
将素材包中的…\常用ENVI功能扩展工具\植被增强工具envi_vegetation_enhance.sav;以及数据101-处理专题:高分一号PMS数据处理里面ENVI_China_Satellites_Support_V5.3_9文件拷贝到本地ENVI安装目录extensions中,重启ENVI
注:扩展工具也可以在中国资源卫星ENVI扩展工具商店下载
(二)打开数据并查看详细信息
1)打开数据(查看解压的GF1号卫星PMS数据:分为MSS和PAN)

2)ENVI界面里面使用安装的扩展工具打开


3)查看打开的全色波段数据(2m分辨率)

查看打开的多光谱波段数据(8m分辨率)

4)以多光谱数据为例查看浏览源信息
双击图层,或者鼠标右击,打开View Metadata(浏览源信息)

可以查看坐标系,时间等各种源信息

(三)多光谱大气校正——辐射定标
1)在Toolbox中,打开/Radiometric Correction/Radiometric Calibration,对多光谱图像辐射定标为辐射量度图像
辐射定标:DN值转化为辐射亮度值
原理:Rad=gain(增益)*DN(图像中的信息)+offset(偏移)

2)在工具框中选择多光谱数据

空间裁剪:可以自行选择裁剪所需区域

波段裁剪:可以自行选择裁剪区域的所需波段数目

3)打开辐射定标面板,进行设置
整个面板设置是为了后面的大气校正

点击Apply FLAASH ,选择存储路径,点击

得到多光谱辐射定标的结果

4)查看辐射定标的结果与原始多光谱信息的变化

即进行了一个线性变换

(四)多光谱大气校正——FLAASH大气校正
1)在Toolbox中,打开/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction

2)在大气校正面板中进行设置参数
a)输入辐射影像:进行辐射定标的多光谱影像
- 输出大气校正位置:GF1_PMS2_MSS_Dat_Ref.dat
- 输出临时文件位置:建议放在存储空间较大的非系统盘

因为在辐射定标时已经设置了单位,所以这里不用改,如果没有设置,这里应该填10

存储面板参数已经设置好

b)进行设置经纬度(在这里已经默认被系统识别)

如果系统无法识别,可用识别工具打开,让红十字落在大致影像中心位置即可,读取消息框中的经纬度即可

c)设置传感器类型、高度、海拔以及像元尺寸
这里传感器类型、高度、像元尺寸已经自动识别

海拔查询操作如下
打开Open World Data——Elevation

打开Statistics——Compute Statistics

点击Stats Subset——打开统计面板

点击添加框中的File,选择辐射定标后的影像

自动识别辐射定标的尺寸大小

点击确定后,在Stats Subset中已经选择好统计的信息

进行统计

统计好信息,海拔大约为2km

d)进行设置时间
系统已经默认读取

如果时间读取错误,可根据GF1号的详细信息进行查询

e)进行大气模型的选择
这里可以点击Help进行查询

如下这张表,根据时间和经纬度可查询特定的大气模型

f)气溶胶模型设置

气溶胶繁衍参数:是否可进行KT变换——打开多光谱设置
Defaults中选择陆地,发现KT变换没有改变(陆地;海陆;海洋)

所以气溶胶繁衍参数——选择None

不同传感器具有不同的波谱文件(这里已经自行构建,但其它大多数传感器需要手动构建波谱文件函数)

g)进行高级设置

设置完成后点击Apply

大气校正完成

3)大气校正结果没有显示,将结果进行标准假彩色显示


4)进行查询对比光谱曲线


打开大气校正与辐射定标查看变化
蓝波段吸收率下降,绿波段有个反射波峰,而红光有个波谷后直线上升

(五)多光谱正射校正
1)在Toolbox中,启动/Geometric Correction/Orthorectification/RPC Orthorectification Workflow工具

2)输入文件里面选择:大气校正后的影像

DEM File里面选择:文件数据自带的MSS的dem数据

设置完成

在下一步后的面板中Advanced里面现在输出分辨率8m,选择方法双线内插

设置输出路径:GF1_PMS2_MSS_Rad_Ref_rpcortho.dat

纠正完成

纠正完成后查看详细信息发现影像地图的坐标系统发生了变化

(六)全色大气校正——辐射定标
1)在Toolbox中,打开/Radiometric Correction/Radiometric Calibration

打开工具框选择全色波段影像

由于多光谱FLAASH大气校正的结果为扩大了10000倍的反射率数据,为了让融合图像效果好,需要将全色数据与多光谱数据的像元值变成一致。这里使用辐射定标工具将全色数据定标为大气表观反射率,并扩大10000倍
2)在辐射定标工具框中对全色波段进行设置
辐射定标工具将全色数据定标为大气表观反射率,并扩大10000倍

(七)全色波段——正射校正
1)在Toolbox中,启动/Geometric Correction/Orthorectification/RPC Orthorectification Workflow工具

输入文件:辐射定标后的全色影像

DEM File:解压数据里面的PAN的dem

2)在下一步后的工具框中进行设置
Advanced:输出像元大小为2m,方法为双线内插

设置输出路径:FG1_PMS2_PAN_Rad_repcortho

得到正射纠正后的影像

(八)图像融合
1)在Toolbox中,启动/Raster Management/Convert Interleave,将多光谱图像的储存顺序由BSQ转成BIP。提高融合速度。

2)查看影像发现有黑框,这会影响融合效果,当消除

双击图层打开消息信息。点击Edit Metadata(编辑源信息)

在工具框中点击Add,添加Data Ignore Value(数据忽略值)

在Data Ignore Value里面设置0

去除黑框完成

对全色影像采取同样的方法进行去除

3)转化数据格式
打开Convert Interleave,选择正射纠正后的多光谱影像

设置对话框参数


重新打开

4)在Toolbox中,启动/Image Sharpening/NNDiffuse Pan Sharpening

在工具框中输入相应参数
低精度栅格:多光谱
高精度栅格:全色

设置输出路径,保存

融合完成

5)对正射纠正后的影像和融合的影像进行对比
正射纠正的影像

融合的影像

发现融合后的图像保留了多光谱与全色数据的优点,图像质量高
也可以通过查看工具进行比较查看
点击 多光谱正射纠正影像,打开Display in Portal




也可以通过波谱曲线查看
先打看融合图像的光谱曲线,然后点击Additional Profiles

选择多光谱正射纠正后的影像,可导入此影像光谱曲线

两条曲线显示在一个矩形面板中
点击城市的某个地方,波谱曲线差别很大

点击植被区域,发现差别很小,融合后影像较低了吸收的光谱曲线


这里再介绍一个植被增强工具(安装的扩展工具)
点击Vegetation Enhance

选择正射纠正的多光谱数据为例(融合图像数据太大)

在植被增强面板中进行设置,这里选择Low增强就可,可以根据需要进行设置

进行植被增强的影像

发现已经很明显提取出来了影像的植被

原理:NIR*0.2+green*0.8
进行植被增强和没有增强的多光谱正射纠正的影像对比

