GF数据预处理

最新推荐文章于 2024-05-08 17:54:20 发布

今天也要加油吖！

最新推荐文章于 2024-05-08 17:54:20 发布

阅读量1.4k

点赞数 27

文章标签：学习方法课程设计笔记经验分享

本文链接：https://blog.csdn.net/YYQX212/article/details/134928586

版权

高分系列卫星介绍与数据下载 - 知乎 (zhihu.com)

如何在中国资源卫星应用中心下载GF1、GF4、GF6、CB04、CB04A、HJ1A、HJ1B、HJ2A、HJ2B遥感影像数据_哔哩哔哩_bilibili

一、数据查询

自然资源卫星遥感云服务平台 (sasclouds.com)

中国资源卫星应用中心 (cresda.com)

-1陆地观测卫星数据服务 (cresda.cn)

二、GF-1

2.1GF-1 PMS数据预处理

下载的是2m/8m分辨率的PMS数据

PMS预处理参考视频地址：101-102-处理专题：高分一号PMS WFV数据处理_哔哩哔哩_bilibili

GF-1WFV数据16m分辨率

WFV预处理参考视频：高分一号二号数据预处理_哔哩哔哩_bilibili

三、GF-2

样例数据下载

网盘链接

高分二号(GF-2)卫星样例数据下载高分二号样例数据_免费高速下载|百度网盘-分享无限制

国家科技资源共享服务平台下载4景免费数据。

网址：Dataset of GF-2 satellite images (2017)-National Tibetan Plateau Data Center (tpdc.ac.cn)

GF-2预处理

参考地址：

ENVI5.3下高分二号（GF2）数据预处理 - ENVI-IDL技术殿堂 - 博客园 (cnblogs.com)

GF2遥感影像预处理全流程 - smsmdmx - 博客园 (cnblogs.com)

高分二号影像数据预处理及裁剪过程 - 知乎 (zhihu.com)

1.多光谱图像预处理

辐射定标

一般来讲，辐射定标就是将图像的数字量化值（DN）转化为辐射亮度值或者反射率或者表面温度等物理量的处理过程。

1.打开ENVI，使用国产卫星扩展工具打开MSS影像，启动File→Open As→China Satellites→GF2，选择MSS.xml文件打开；

辐射定标需要高分二号传感器的绝对辐射定标系数，右键单击View etadata，查看元数据信息，这里软件已经自动正确读取两个绝对定标系数；

2.在Toolbox中，Radiometric Correction→Radiometric Calibration，在File Selection中选择待处理影像，点击OK；弹出Radiometric Calibration对话框，Calibration type确认为Radiance，单击Apply FLAASH Setting，设置输出路径与文件名，点击OK开始执行；

大气校正

FLAASH大气校正需要影像的中心波长信息，ENVI5.3.1能自动识别GF2数据的头文件信息，可通过View Metadata→Spectral查看；

在Toolbox中，双击Radiometric Correction→Atmospheric Correction Module→FLAASH Atmospheric Correction工具启动FLAASH模块；

(1). 设置输入与输出文件信息：

Input Radiance Image：输入辐射定标之后的数据；Radiance Scale Factors选择第二项Use single scale factor for all bands。因为在辐射定标时已经设置了单位，所以这里不用改，如果没有设置，这里应该填10

Output Reflectance File为输出的大气校正文件

Output Directory for FLAASH Files为输出的临时文件夹

Rootname for FLAASH Files：设置大气校正其他输出结果的根文件名

(2). 设置传感器及图像信息：

Sensor Type：传感器类型，这里选择Multispectral→UNKNOWN→MSI;

Ground Elevation：成像区域平均高度，

Ground Elevation的一种估测方法：1.打开需要估计区域对应的图像, 2.选择Flie->Open World Data->Elevation(GMTED2010)，打开ENVI自带的全球900m分辨率的DEM数据, 3.在ToolBox中，选择/Statistics/Compute Statistics, 4.选择GMTED2010.jp2数据，打开Stats Subset选项, 5.选择File，选择估计区域对应图像, 6.选择Image，查看选择后结果。红色区域为估计区域, 7.mean即为估计值，单位为m

(3). 大气模型和气溶胶模型：需要根据经纬度、影像区域以及成像时间选择

Atmospheric model：Mid-Latitude Summer；

Aerosol Model：Urban；

Aerosol Retrieval：None；

光谱响应函数设置Multispectral Settings：

选择GUI设置方式，点击“Filter Function File”，Open>New File,加载GF2的光谱响应函数（文件后缀为sli，需自主要PMS1还是PMS2 ，GF2-PMS1的光谱响应函数gf2_pms1_mss.sli，GF2-PMS2的光谱响应函数gf2_pms2_mss.sli）；位于ENVI安装目录：*\Program Files\Exelis\ENVI53\resource\filterfuncs下：

方法一：Advanced Settings：这里大部分都可以保持默认设置，但由于是多光谱数据，故需将Modtran Resolution设置为15cm-1.

方法二：点击Advancced Settings，其中Use Tiled Processing修改为No，确认后回到主面板点击Apply运行；如运行过程中出现下图错误，则返回Advancced Settings，将Use Tiled Processing修改为Yes，Tile Size(Mb)设置为100Mb-1024Mb，对应本计算机内存8GB-64GB，确认后返回主面板点击Apply运行。

所有设置完成之后，点击Apply执行大气校正，完成后会得到反演的能见度和水汽柱含量；

选择Display>Profiles>Spectral查看大气校正前后同一地物波谱曲线变化

正射校正

先进行辐射定标和大气校正，然后进行几何校正。因为几何校正的时候会重采样，重采样会改变像元值，从而影响辐射定标和大气校正的结果。

点击File→Open打开大气校正后(或者原始的多光谱)的影像，View Metadata查看其元数据信息，可以看到ENVI很好地识别了数据的RPC信息；有了RPC信息之后，就可以基于这些RPC信息进行正射校正；

在Toolbox中，启动/Geometric Correction/Orthorectification/RPC Orthorectification Workflow工具，打开正射校正流程化工具；

在File Selection面板中，Input File选择经过大气校正的多光谱数据，DEM File会默认选择全球分辨率为900米的DEM数据，我们这里保持默认（如果有更高分辨率的DEM数据，可以替换此数据），点击Next；

在RPC Refinement面板中，有四个选项卡可以选择。

如果有实测的或从其他途径获取的控制点数据，可以在该面板中进行添加，添加后在Statistics选项卡中可以看到相应的误差统计信息；

1) 切换到Advanced选项卡，修改输出像元大小Output Pixel Size为4米，重采样方法Image Resampling选择三次卷积法，其他参数保持默认；

2) 切换到Export选项卡，选择输出文件格式，设置输出路径及文件名，点击Finish；完成正射校正

2.全色波段影像处理

辐射定标

打开ENVI，启动File→Open As→China Satellites→GF2，选择PAN.xml文件打开；

在Toolbox中，打开/Radiometric Correction/Radiometric Calibration，

由于多光谱FLAASH大气校正的结果为扩大了10000倍的反射率数据，为了让融合图像效果好，需要将全色数据与多光谱数据的像元值变成一致。这里使用辐射定标工具将全色数据定标为大气表观反射率，并扩大10000倍。

参数设置如下：Calibration Type: Reflectance，全色影像定标为大气表观反射率；Output Data Type：Uint；Scale Factor：10000；

正射校正

在Toolbox中，启动/Geometric Correction/Orthorectification/RPC Orthorectification Workflow工具,

Input File选择上一步辐射定标后的文件，DEM File默认不变，点击Next

Advanced面板中注意Output Pixel Size为1m,

Export面板中在Out Filename中设置输出名称及路径，点击Finish，完成正射校正

影像融合

NNDiffuse Pan Sharpening方法

NNDiffuse Pan Sharpening工具要求输入的多光谱和全色数据的空间分辨率是整数倍的（本例正射校正时分别将多光谱的全色的分辨率重采样为4米和1米，就是为了方便该工具的使用）

在Toolbox中，启动/Image Sharpening/NNDiffuse Pan Sharpening, 低精度栅格：多光谱高精度栅格：全色设置输出路径

Gram-Schmidt Pan Sharpening方法

在Toolbox中，启动/Image Sharpening/Gram-Schmidt Pan Sharpening,

选择低分辨率影像：MSS, 选择高分辨率影像：PAN,

设置参数如图，设置输出路径。点击ok执行