【教程】GF-1和GF-2PMS数据预处理全流程

GF-1和GF-2的PMS数据全部的预处理流程主要包括5个部分：（1）辐射定标（2）大气校正（3）正射校正（4）图像融合（5）影像配准。由于全色图像无法进行大气校正，在定量遥感中一般仅使用多光谱数据。本期所用数据为GF1_PMS1_E82.3_N34.2_20151103_L1A0001149149

1.多光谱辐射定标

在ENVI5.3.1版本中，直接支持高分一号PMS数据的辐射定标和大气校正。

（1）选择Open As>China Satellites>GF-1，选择

GF1_PMS1_E82.3_N34.2_20151103_L1A0001149149-PAN1.xml

和GF1_PMS1_E82.3_N34.2_20151103_L1A0001149149-MSS1.xml文件打开。

(2)在Toolbox中，打开/Radiometric Correction/Radiometric Calibration，选择多光谱数据。

(3)在Radiometric Calibration面板中，单击Apply FLAASH Settings按钮，几个参数自动选择符合FLAASH大气校正要求，包括定标类型（Radiance）、存储顺序（Interleave）和辐射亮度单位（Scale Factor）。

(4)选择数据路径和文件名，单击OK执行。

2.多光谱数据大气校正

(1)在Toolbox中，打开Radiometric Correction>Atmospheric Correction Module>FLAASH Atmospheric Correction。

(2)点击Input Radiance Image，前面辐射定标好的数据，在Radiance Scale Factors面板中选择Use single scale factor for all bands，由于定标的辐射量数据与FLAASH的辐射亮度的单位一致，所以在此Single scale factor选择：1，单击OK；由于使用Radiometric Calibration自动将定标后的辐射亮度单位调整为（μW）/（cm2*nm*sr），与FLAASH要求的一致，因此在Radiance Scale Factors中输入1。

（3）设置输出文件及路径设置；

（4）传感器基本信息设置：成像中心点经纬度FLAASH自动从影像中获取；传感器高度（Sensor Altitude）：645km；像元大小（pixel Size）：8m；成像区域平均高度可以通过统计DEM数据获取；成像时间：在图层管理中右键选View metadata，在Time选项中可以获取。

（5）大气模型和气溶胶模型，根据经纬度和影像区域选择（单击Help，找到经纬度和成像时间的对照表）。

（6）气溶反演方法选择None（缺少短波红外），能见度设置为40km（查看Help中说明）。

（7）多光谱设置面板按照默认参数

（8）打开Advanced settings面板，设置Use Tied Peocessing：No。

（9）单击Apply执行处理。

3.全色数据辐射定标

(1)在Toolbox中，打开/Radiometric Correction/Radiometric Calibration，选择全色数据。

(2)在Radiometric Calibration面板中，设置：Calibration Type: Reflectance；Output Data Type：Uint；Scale Factor：10000。

(3)选择输出路径和文件名。

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

4.多光谱和全色数据正射校正

高分一号的L1A级包括了RPC文件，在经过了辐射定标、大气校正等处理，ENVI会自动将RPC嵌入处理结果中，可以在图层管理中辐射定标或者大气校正结果图层右键选View metadata，RPC选项就是嵌入的RPC文件。可以直接使用Geometric Correction>Orthorectification>RPC Orthorectification Workflow工具进行正射校正。

(1)在Toolbox中，启动/Geometric Correction/Orthorectification/RPC Orthorectification Workflow工具。在File Selection中选择全色数据辐射定标结果。DEM使用ENVI自带DEM（可以选择使用更高精度的DEM）。

(2)在RPC Refinement步骤中，打开Advanced面板，设置Output Piexl Size：2。

(3)在Exports面板中，选输出路径和文件名。

(4)单击Finish执行处理。同样的方法对多光谱的大气校正结果进行正射校正，多光谱数据输出像元设置为8。

5.图像融合

为了提高融合速度，将多光谱数据的储存顺序由BSQ转成BIP。

(1)在Toolbox中，启动/Raster Management/Convert Interleave，选择上一步中多光谱正射校正结果。如下设置参数，Convert In Place：Yes，不生成新的文件。

(2)在Toolbox中，启动/Image Sharpening/NNDiffuse Pan Sharpening，分别选择多光谱和全色数据。

(3)选择输出路径和文件名。

(4)单击ok执行处理。

6.数据配准

在进行多期影像对比分析时，需要对多景影像进行配准操作，从而保证多景影像位置的一致性。在进行数据配准时需要以一景影像为基础，这里选择具有更高分辨率的GF-2数据为基准影像，对裁剪后的GF-1影像进行配准。

(1)在Toolbox中，启动Geometric Correction>Registration>Image Registration，Base Image File选择基准影像，Warp Image File 选择带配准影像，单击Next。

(2)在Toolbox中，启动Geometric Correction>Registration>Image Registration，Base Image File选择基准影像，Warp Image File 选择带配准影像，单击Next。

(3)再单击Next，选择呢Tie Points面板中的Show Table，去除误差较大的点，保证RMS误差小1。在Warping面板中，设置Output Pixel Size From 参数为Warp Image。单击Next，设置输出路径和文件名，单击OK执行

这样数据处理的全部流程就结束了。不过有些操作比如影像融合、地理配准等操作并不是必须的，可以根据自己的需要进行选择是否处理。

觉着有用的话，欢迎分享、点赞和关注生态与遥感应用公众号。

【教程】GF-1和GF-2PMS数据预处理全流程相关推荐

数据分析问题（异常值识别）中数据预处理部分流程（含2022年全国服务外包大赛实例）
博主个人理解的数据预处理主要包括个方面:读取文件 => 数据概览 => 缺失值填补 => 数据分布预览 => 衍生特征设计.这套流程在完成异常值识别时作为数据预处理时没有 ...
软件杯龙源风电赛题培训！千万分钟数据和全流程基线等你来战
‍‍"中国软件杯"大学生软件设计大赛是一项面向中国在校学生的公益性赛事,大赛由国家工业和信息化部.教育部.江苏省人民政府共同主办,是全国软件行业规格最高.最具影响力的国家级一类赛事 ...
数据预处理常用方法流程
文章目录数据读取和选择缺失数据找出缺失数据处理缺失数据丢弃观察值丢弃特征填充缺失数据数值变量分类变量替换缺失值异常值如何找出异常值数值变量直方图箱形图分类变量条形图 ...
音频数据建模全流程代码示例：通过讲话人的声音进行年龄预测
来源:DeepHub IMBA 本文约6100字,建议阅读10+分钟本文展示了从EDA.音频预处理到特征工程和数据建模的完整源代码演示. 大多数人都熟悉如何在图像.文本或表格数据上运行数据科学项目. ...
从实战中了解数据开发全流程——DataWorks OpenAPI实战
简介:DataWorks作为飞天大数据平台操作系统,历经11年发展,形成了涵盖数据集成.数据开发.数据治理.数据服务的一站式大数据开发治理平台.很多企业用户在使用产品的过程中希望他们的本地服务能够和阿 ...
3000字长文为你解读数据仓库与复杂业务数据建模全流程
近日,2022年个推TechDay"治数训练营"系列直播课第一期圆满举办.个推资深大数据研发工程师为大家深入浅出地介绍了数据仓库的前世今生以及数据建模的常用方法. 本文对" ...
机器学习笔记2：数据建模全流程
数学建模标准流程包括商业理解.数据理解.数据清洗.建立模型.部署应用留个流程,具体如下: 我们将数据建模标准流程转化为可执行可实施部分,那就是以下几个步骤: 商业理解即对建模的项目需求和目标进行综合 ...
在机器学习中，如何用Python进行数据预处理？
数据分析之路持续学习中- - - 近期学习了机器学习中的数据预处理章节,在此进行分享,欢迎大家讨论指正. 顺便说一下,这里我使用的软件是Anacnoda 3中已经安装好的Spyder 3,这个软件用起 ...
AlphaFold2源码解析(3)--数据预处理
AlphaFold2源码解析(3)–数据预处理数据预处理整体流程数据处理入口: feature_dict = data_pipeline.process( input_fasta_path=fas ...

【教程】GF-1和GF-2PMS数据预处理全流程

【教程】GF-1和GF-2PMS数据预处理全流程相关推荐

最新文章

热门文章