ENVI下高分一号PMS相机多光谱数据大气校正 envi地形校正

高分一号 PMS相机的多光谱数据拥有很好的空间分辨率(8米)和幅宽(60km),包括蓝、绿、红、近红外4个波段。本文介绍利用ENVIFLAASH工具完成高分一号PMS多光谱图像的大气校正,主要分为以下三步:

l数据定标

l波谱响应函数制作

lFLAASH大气校正

注:本文是在ENVI5.0下和ENVIclassic完成,其他版本基本类似。

1 数据定标

利用以下公式可将GF-1卫星各载荷的通道观测值计数值DN转换为卫星载荷入瞳处等效表观辐亮度数据。

式中:Gain为定标斜率;DN为卫星载荷观测值;Bias为定标截距。

定标参数参照表1,来源于中国资源卫星应用中心:http://www.cresda.com/n16/n1115/n1522/n2103/index.html

表1 高分一号星PMS相机的定标系数

卫星载荷

波段号

Gain

Bias

PMS1

PAN

0.1886

-13.127

Band1

0.2082

4.6186

Band2

0.1672

4.8768

Band3

0.1748

4.8924

Band4

0.1883

-9.4771

PMS2

PAN

0.1878

-7.9731

Band1

0.2072

7.5348

Band2

0.1776

3.9395

Band3

0.177

-1.7445

Band4

0.1909

-7.2053

WFV1

Band1

0.1709

-0.0039

Band2

0.1398

-0.0047

Band3

0.1195

-0.0030

Band4

0.1338

-0.0274

WFV2

Band1

0.1588

5.5303

Band2

0.1515

-13.642

Band3

0.1251

-15.382

Band4

0.1209

-7.985

WFV3

Band1

0.1556

12.28

Band2

0.1700

-7.9336

Band3

0.1392

-7.031

Band4

0.1354

-4.3578

WFV4

Band1

0.1819

3.6469

Band2

0.1762

-13.54

Band3

0.1463

-10.998

Band4

0.1522

-12.142

进行数据定标计算时候,我们可以有两种方法选择使用:

一、Bandmath

以下是利用Bandmath工具进行定标,这种方法灵活,不过需要一个一个波段定标。

1)有了绝对定标参数和定标公式,选择Toolbox/Band Ratio/BandMath工具很容易进行传感器定标。

2)由于是单个波段文件的定标,选择Toolbox/Raster Management/LayerStacking将定标后的单波段文件组合成一个多波段文件。

二、Apply Gain and Offset工具

我们可以使用Apply Gain and Offset工具。。

1)启动Toolbox/Radiometric Correction/Apply Gain andOffset。

2)在Gain and Offset Values中填入增益Gain值和偏移Offset。

3)选择输出路径和文件名执行定标。

图1Gain and Offset Values面板

计算好之后,继续执行以下步骤(原数据的储存顺序是BIP):

(1)选择Toolbox/Raster Management/ Edit ENVI Header->Edit Attributes->Wavelength,将每个波段的中心波长输入——b1(503nm),b2(576nm),b3(680nm),b4(810nm)

讨论:由于没有找到官方公布的中心波长信息,这里取与波谱响应中反射率为1的波长。

目前为止,影像数据准备工作完成了。

2 波谱响应函数

波段响应函数,英文名为spectral responsefunction(SRF)或者叫RelativeSpectral Response (RSR) ,与宽波段传感器出现。我们知道每一个波段都有一个波段范围,比如HJ-b1(475nm)波段为520-430纳米,实际上传感器的感光元件在这个波段范围内的每一点所感应的强度都是不一样。在成像中,原则上讲应该根据波段响应函数来进行加权平均,但由于处理起来比较麻烦,而且一般的精度要求不太高,所以大多数图像都是直接取了波段范围内的中点值来运算。波段响应函数是描述一定波长范围内(超出波段范围)的量子效应,当需要精确计算像元响应时候,比如大气校正反演真实地表反射率,就需要使用波谱响应函数。

任何传感器在设计时都会给出严格的波段响应函数,高分一号也不例外。在中国资源卫星应用中心可下载:http://www.cresda.com/n16/n1115/n1522/n2118/186272.html。与资源一号02c和资源三号一样,高分一号也是以.xls格式提供,值得注意的是,第三列和第四例后面有部分无值区,手动填入0。

将数据拷贝到文本文件中。如下图中为波谱响应函数数字表达的一部分,第一列表示波长,后面五列分别表示一个全色和4个多光谱对应波长的波谱响应值。

图2高分一号波谱响应值

在ENVI中,使用波谱曲线来描述波谱响应函数,也就是以波长作为x轴,波谱响应值作为y轴,存储格式为ENVI波谱库文件(.sli)。下面介绍ENVIclassic中的波谱响应函数的制作。

(1)启动ENVI Classic

(2)选择Window->Start New Plot Window,ENVIPlot Window窗口中,选择File->InputData->ASCII,如图3所示,自动将第一列作为X轴,后面3-6列作为Y轴。单击OK。

(3)如图所示,生成了4条曲线。选择Edit->DataParameters,更改每一条曲线的名称:b1,b2,b3,b4,便于区分。

(4)选择File->Save Plot As->SpectralLibrary,将波谱曲线保存为波谱库文件.sli。

图3导入ASCII文件

图4高分一号波谱响应函数的曲线表达

注:如果出现以下情况,原因是没有将多光谱波段1和波段2空白部分填补0。

图5错误情况

3 FLAASH大气校正

(1)工具箱/RadiometricCorrection/Atmospheric Correction Module/FLAASH AtmosphericCorrection打开FLAASH大气校正模块;

(2)点击InputRadiance Image,前面处理好的数据,在RadianceScale Factors面板中选择Usesingle scale factor for all bands,由于定标的辐射量数据与FLAASH的辐射亮度的单位相差10倍,所以在此Singlescale factor选择默认:10,单击OK;

注:定标后的单位是W⋅m−2⋅sr−1⋅μm−1,与FLAASH要求的单位(μW)/(cm2*nm*sr)相差10倍关系,因此在RadianceScale Factors中输入10缩放系数。

(3)设置输出文件及路径设置;

(4)传感器基本信息设置:

l成像中心点经纬度FLAASH自动从影像中获取。

l传感器高度(Sensor Altitude):645km

注:写这篇文档时没有得到官方正式数据,下面截图是公布之前填写的505,最终公布的数据为645km。

l像元大小(pixelSize):8m

l成像区域平均高度可以通过统计DEM数据获取

l成像时间:从数据头文件中读取(.XML),位置为:2013-09-0813:14:07;减去8换算GMT时间

(5)大气模型和气溶胶模型,根据经纬度和影像区域选择。

(6)气溶胶反演方法选择None(缺少短波红外),能见度设置为40km。

图6FLAASH大气校正参数设置

(7)单击MultispectralSetting按钮,在FilterFunction File 导入光谱响应曲线“gf-1pms.sli”,单击OK;

(8)单击AdvancedSettings,在高级设置中

lTile Size 默认的是Cashsize 的大小,手动改为50-100Mb(根据内存大小设定),单击OK;

l从xml文件中得到:97.0867

79.1983

讨论:一般中低分辨率的影像近似天顶角:180和方位角:0,即垂直观测。ENVIFLAASH的天顶角:90-180度,方位角:-180and 180,没有查到高分一号的角度说明,这里选择默认。

ENVI下高分一号PMS相机多光谱数据大气校正 envi地形校正

图7大气校正高级参数设置

(9)设置好后,在大气校正模块面板中,单击Apply。

(10)大气校正完成后,检查大气校正的结果,分别加载校正前后的图像(选择CIR假彩色方式加载,可以更好的识别植被),查看典型地物的大气校正前后的光谱曲线。

4 总结

从上面可以看到,对于未知的多光谱传感器的FLAASH大气校正,关键是获取波谱响应函数,而波谱响应函数一般都是公开的。类似的方法可用于高分一号宽幅传感器的FLAASH大气校正。

  

爱华网本文地址 » http://www.413yy.cn/a/25101017/350183.html

更多阅读

“润元昌大班章一号”开汤 润元昌普洱茶

这次有幸得到润元昌茶业的关照,寄来新茶样大班章一号二小盒子。哈哈,开心呀。就喜欢上这“大班章地图”盒子,也标注了润元昌“位于大班章茶区的中心,是目前大班章茶区范围内第一家也是目前唯一一家精致茶厂”。感觉不错。打开小包装取出

电视剧《天字一号》故事原型 张国福为原型的电视剧

电视剧《天字一号》故事原型最近热播的反特电视剧《天字一号》,取材于“炮击天安门”阴谋案和老电影《国庆十点钟》。更多的改编线索,来自于国庆十点钟。编剧比较大胆,把“炮击天安门”阴谋案的素材也用上了。这是一个经典的侦破案件

2008年中央一号文件 1987年中央一号文件

关于2008年中央一号文件跨越藩篱超越梦想——写在中央关于“三农”工作的第10个一号文件公布之际《中共中央国务院关于切实加强农业基础建设进一步促进农业发展农民增收的若干意见》30日正式公布。这是近5年来第5个关于“三农”问题

声明:《ENVI下高分一号PMS相机多光谱数据大气校正 envi地形校正》为网友緗約看齣分享!如侵犯到您的合法权益请联系我们删除