相对辐射校正和绝对辐射校正

基于物理模型的绝对辐射校是利用一系列参数(例如,卫星过境时的地物反射率,大气的能见度,太阳天顶角和卫星传感器的标定参数等)将遥感图像进行校正的方法。仪器引起的误差畸变一般在数据生产过程中由生产单位根据传感器参数进行了校正。对于用户来所,绝对辐射校正的方法主要是辐射传输模型法,该方法校正精度较高,它是利用电磁波在大气中的辐射传输原理建立起来的模型对遥感图像进行大气校正的方法。由于有不同的不同的假设条件和适用的范围,因此产生很多可选择的大气较正模型,例如6S模型、LOWTRAN模型、MODTRAN模型、ATCOR模型等。

基于统计模型的相对辐射校正,主要包括不变目标法、黑暗像元法与直方图匹配法等等。不变目标法假定图像上存在具有较稳定反射辐射特性的像元,并且可确定这些像元的地理意义,那么就称这些像元为不变目标,这些不变目标在不同时相的遥感图像上的反射率将存在一种线性关系。当确定了不变目标以及它们在不同时相遥感图像中反射率的这种线性关系,就可以对遥感图像进行大气校正。黑暗像元法的基本原理就是在假定待校正的遥感图像上存在黑暗像元区域、地表朗伯面反射、大气性质均一,忽略大气多次散射辐照作用和邻近像元漫反射作用的前提下,反射率很小的黑暗像元由于大气的影响,而使得这些像元的反射率相对增加,可以认为这部分增加的反射率是由于大气程辐射的影响产生的。利用黑暗像元值计算出程辐射,并代入适当的大气校正模型,获得相应的参数后,通过计算就得到了地物真实的反射率。直方图匹配法是指如果确定某个没有受到大气影响的区域和受到大气影响的区域的反射率是相同的,并且可以确定出不受影响的区域,就可以利用它的直方图对受影响地区的直方图进行匹配处理。此外,还有很多基于统计模型的方法,如有人提出利用小波变换的遥感图像相对辐射校正方法。该方法对源图像小波变换域的低频成分实施辐射变换,并保持高频成分不变,重构的图像具有保持高频信息的特性,因而能够较好地保留原图像中由于地物变化引起的辐射差异;也有人利用主成分分析法把遥感图像中有用的信息和大气影响噪音区分开来。

大气辐射传输模型6S

1986年,法国Université des

Sciences et Technologies de

Lille(里尔科技大学)大气光学实验室Tanré等人为了简化大气辐射传输方程,开发了太阳光谱波段卫星信号模拟程序5S(SIMULATION

OF THE SATELLITE SIGNAL IN THE SOLAR

SPECTRUM),用来模拟地气系统中太阳辐射的传输过程并计算卫星入瞳处辐射亮度。1997年,Eric

Vemote对5S进行了改进,发展到6S(SECOND SIMULATION OF THE SATELLITE SIGNAL IN

THE SOLAR

SPECTRUM),6S吸收了最新的散射计算方法,使太阳光谱波段的散射计算精度比5S有所提高。

这种模式是在假定无云大气的情况下,考虑了水汽、CO2、O3和O2的吸收、分子和气溶胶的散射以及非均一地面和双向反射率的问题。6S是对5S的改进,光谱积分的步长从5nm改进到2.5nm,同5S相比,它可以模拟机载观测、设置目标高程、解释BRDF作用和临近效应,增加了两种吸收气体的计算(CO、N2O)。采用SOS

(successive order of scattering) 方法计算散射作用以提高精度。缺点是不能处理球形大气和limb

(临边)观测。

它其中主要包括以下几个部分:

(1)太阳、地物与传感器之间的几何关系:用太阳天顶角、太阳方位角、观测天顶角、观测方位角四个变量来描述;

(2)大气模式:定义了大气的基本成分以及温湿度廓线,包括7种模式,还可以通过自定义的方式来输入由实测的探空数据,生成局地更为精确、实时的大气模式,此外,还可以改变水汽和臭氧含量的模式;

(3)气溶胶模式:定义了全球主要的气溶胶参数,如气溶胶相函数、非对称因子和单次散射反照率等,6S中定义了7种缺省的标准气溶胶模式和一些自定义模式;

(4)传感器的光谱特性:定义了传感器的通道的光谱响应函数,6S中自带了大部分主要传感器的可见光近红外波段的通道相应光谱响应函数,如TM,MSS,POLDER和MODIS等;

(5)地表反射率:定义了地表的反射率模型,包括均一地表与非均一地表两种情况,在均一地表中又考虑了有无方向性反射问题,在考虑方向性时用了9种不同模型)。

这5个部分便构成了辐射传输模型,考虑了大气顶的太阳辐射能量通过大气传递到地表,以及地表的反射辐射通过大气到达传感器的整个辐射传输过程。

6S的输入参数主要有9个部分组成:

(1)几何参数

6S两种输入方法⑴

太阳和卫星的天顶角和方位角以及观测时间(月,日)。⑵

卫星的接收时间(月,日,年)、像素点数、升交点时间,由程序计算太阳和卫星的天顶角和方位角。特别注意的是这里的时间采用世界时且要精确到1/6秒。

(2)大气模式

6S给出几种可供选择的大气模式,热带、中纬度夏季、中纬度冬季、近极地夏季、近极地冬季、美国62标准大气也可自定义大气模式。

(3)气溶胶模式

三种选择:⑴

无气溶胶。⑵

自定义气溶胶模式。如,四种基本气溶胶的体积的加权平均;气溶胶的谱分布加光度计测量结果(光学厚度)和复折射指数;直接给出消光系数。⑶

提供的三种气溶胶模式大陆型,海洋型和乡村型。

(4)气溶胶浓度

两种选择:⑴

在550nm处的光学厚度 ⑵ 气象能见度(km)。故它也提供了两者的相互关系。

(5)地面高度

以千米为单位的地面海拔高度(设为负值)。

(6)探测器高度

-1000代表卫星测量,0为地基观测,飞机航测输入以千米为单位的负值。

(7)探测器的光谱条件

给出了常见卫星Meteosat,Goes,NOAA/AVHRR和HRV,Landsat

TM 和MSS,Modis Polder的每个通道的光谱响应函数,也可选择自定义。

(8)地表特性

可以选择地表均一或不均一,也可选择地表为郎伯体或双向反射。6S给出了九种比较成熟的BRDF模式供用户选择,也可自定义BRDF函数(输入个角度的反射率及入射强度)

(9)表观反射率

输入反射率或辐射亮度,同时也决定模式是正向还是反向工作。当RAPP0(辐射亮度)或-1

“湖泊水色遥感大气校正研究”是中科院南京地理与湖泊研究所委托南京大学GIS与遥感实验室的项目,受到中科院领域前沿项目(CXNIGLAS-A02-014)的资助。起止年限:2005.6—2007.6。

2.

研究内容

本研究通过卫星遥感影像获取湖泊水体的离水辐亮度或遥感反射率,即研究如何对湖泊水色遥感中的卫星影象进行大气校正。具体内容包括:

对现有的一般的大气校正方法/模式(如6S、LOWTRAN、MODTRAN等)与海洋水色遥感的大气校正方法/模式进行比较,并分析现有现有软件(如ENVI、PCI等)中的大气校正方法;

提出最适合湖泊水体的大气校正方法/流程/模式,用IDL或MATLAB编程实现。

研究的重点是内陆湖泊水体的遥感大气校正,其目的是找到一种适合内陆湖泊水体遥感大气校正的方法和流程。水体的遥感的目标是通过遥感手段来了解水体的物化参数,就目前而言,能从遥感波谱反射特征所获得的参数主要有叶绿素、悬浮物和黄色物质等。因此,研究这三种物质的光谱特征是水体遥感大气校正的前提。

理想的遥感模型是没有大气存在,地面为朗伯体,遥感接收到的光谱直接反映地面目标的状况。但实际情况中,太阳光在下行到达地面和上行到达传感器的过程中,受到大气的影响而衰减,大气校正就是要研究并消除大气条件对遥感的影响,在传感器位置恢复地面目标的光谱特征。

研究区选在太湖,太湖地处长江三角洲,面积36900平方公里,为流域第一大湖,又是长江中下游五大淡水湖之一。水面面积为2338.1km2。

3.

技术路线

水色遥感是快速获取大面积水体组分参数的重要手段。太阳辐射在太阳-水体目标-传感器的传输过程中,气体的吸收、气溶胶与分子的散射这两种大气过程影响了获得的水体反射率信息。由于水体在可见光和近红外波段的反射率很低,传感器接收到的信号只有5-15%来自水体。因此,利用遥感影像精确反演水体信息,大气校正是至关重要的。

大气校正就是在遥感图像的基础上,消除大气影响,获得地面反射率的过程。许多学者研究大气对辐射传输的影响,提出了不同的大气校正模型,主要包括6S,LOWTRAN/MODT-

RAN ,TURNER等。本文主要研究前两种模型。

6S

(Second Simulation of the Satellite Signal in the Solar

Spectrum)大气校正模型是Eric F. Vermote et

al.(1997)在5S模型的基础上发展起来的。6S模型可以很好地模拟太阳光在太阳-地面目标-传感器的传输过程中所受到的大气影响。相对于5S模型,6S模型考虑了地面目标的海拔高度、非朗伯平面的情况和新的吸收气体种类(CH4,N2O,CO)。通过采用the

art

approximation近似算法和SOS运算法则,提高了瑞利和气溶胶散射作用的计算精度。光谱步长提高到了2.5nm。6S模型建立在辐射传输理论基础之上,模型应用范围广,不受研究区特点及目标类型等的影响。

LOWTRAN是一种低分辨率(分辨率≥20cm-1)大气辐射传输模式。在20多年的发展过程中不断扩充和修订基础资料,改进算法,增加可计算的辐射传输结果,从原意义上的“低分辨率大气透过率计算模式”扩展到目前能导出复杂天气条件下多种辐射传输量的“低分辨率大气辐射传输计算模式”,提供了许多新的应用可能性,已被国际上许多应用专家广泛应用于各自的实际问题。MODTRAN(中光谱分辨率大气辐射传输模式)较之LOWTRAN不但提高了光谱分辨率,而且还包括了多次散射辐射传输精确算法——离散纵标法,对有散射大气的辐射传输如太阳短波辐射,比LOWTRAN中的二流近似算法有更高的精度和更大的灵活性。

由于基于辐射传输理论的大气校正需要相应的大气参数,而在一般研究中,这些参数一般很难获得,针对不同研究区的的大气校正方法有时是很实用的。本研究中也尝试使用黑暗象元法、不变目标法、直方图匹配法、大气阻抗植被指数法、主成分分析法等获取水体的遥感反射率或直接获取水体参数信息,并与基本基于辐射传输理论的大气校正结果进行比较

modtran matlab,大气辐射传输模型(6S,MODTRAN)相关推荐

  1. 浅谈经典大气辐射传输模型

    5S模型 1986年由法国里尔科技大学大气光学实验室Tanre D.等人研发了5S模型,5S模型不仅可以用来模拟太阳辐射传输过程,还可以计算辐亮度.该模型假设大气吸收和散射作用可以耦合. 6S模型 6 ...

  2. modtran matlab,modtran4软件及配套中文资料

    [实例简介] modtran是辐射传输模型,可用于计算大气透过率和辐射亮度,本资料中内附有中文的使用说明.参数说明及其相关资料. [实例截图] [核心代码] modtran_4(1) └── modt ...

  3. MODTRAN辐射传输模型使用笔记

    目录 MODTRAN辐射传输模型使用笔记 CARD1主要控制参数设置: CARD2主要控制参数设置: CARD3主要控制参数设置: CARD4主要控制参数设置: MODTRAN辐射传输模型使用笔记 第 ...

  4. 6S大气传输模型修改源码添加、自定义CASI传感器光谱响应

    6S大气传输模型编译以及修改源码添加.自定义CASI传感器光谱响应 在利用6S模型进行CASI影像数据的大气纠正时,发现6S模型中没有对应的光谱响应函数,又不想自己整个输进去,就查资料对源码进行了修改 ...

  5. 6S大气传输模型修改源码添加、自定义卫星光谱响应(以HJ-1B CCD为例)

    6S大气传输模型修改源码添加.自定义卫星光谱响应(以HJ-1B CCD为例) 最近要做国产卫星的大气校正,打算用6s模型模拟气溶胶的查找表,但是发现6s模型中没有国产卫星的相应光谱响应函数,只能在输入 ...

  6. 无线电能传输LCL-S拓扑/WPT MATLAB/simulink仿真模型 (模型左边为两电平H桥逆变器,LCL-S串联谐振,右边不可控整流结构)

    无线电能传输LCL-S拓扑/WPT MATLAB/simulink仿真模型 (模型左边为两电平H桥逆变器,LCL-S串联谐振,右边不可控整流结构) [控制方法]:有用滑模控制移相控制和用PI控制的两种 ...

  7. AI Earth 深度学习模型替换数值天气预报模型中的参数化方案-大气辐射传输方案

    1.背景 太阳辐射和热辐射是大气和海洋运动的最根本的驱动力.大气辐射传输过程实际上已经可以通过一种叫做LBLRTM的辐射模型精确计算,但是LBLRTM模型同时也最为耗时.因此,有各种各样的辐射传输参数 ...

  8. py6s 光谱响应函数_Windows7 64位环境6sv2.1大气传输模型修改源码添加国产高分卫星GF-1 GF-2光谱响应支持...

    下面开始添加国产卫星光谱响应的支持: 以下主要参考文章"6S大气传输模型修改源码添加.自定义卫星光谱响应(以HJ-1B CCD为例)"网址:http://blog.csdn.net ...

  9. matlab求迟滞,MATLAB PI迟滞模型问题

    之前我用matlab进行迟滞模型建模,误差非常大,我之前怀疑是自己的参数辨识出问题了,所以我用论文中的阈值向量和权重向量进行建模,输入数据也是用论文数据,但是出来的拟合的线误差非常大,和文献中相差太大 ...

最新文章

  1. 生信人写程序2. Editplus添加Perl, Shell, R模板和语法高亮
  2. GPTEE中定义的RSA的Algorithm Identifier详解
  3. [luoguP1773] 符文之语_NOI导刊2010提高(02)(DP)
  4. 如何高效搜索汪子熙公众号发表的文章
  5. vuex mapGetters使用
  6. pat1100. Mars Numbers (20)
  7. 随想录(构建自己的代码库)
  8. Low Power概念介绍<Level Shifter>
  9. [转载] python classmethod存在的意义_@classmethod和@staticmethod对初学者的意义?
  10. 截获3389远程登陆的密码
  11. 《领域驱动设计:软件核心复杂性应对之道(修订版)》—第2章 2.3节一个团队,一种语言...
  12. 记录一次bypy使用提示无法创建目录错误
  13. js_单例模式制作无限弹窗(3s内销毁)
  14. 有了这两本书,学习领域驱动设计会很容易
  15. TSP问题-简介与部分解法
  16. 微信公众号支付开发流程
  17. 运维工程师绩效考核表_IT运维绩效考核表
  18. easyui表格 序号如何进行自适应宽度
  19. mysql无参的存储过程_创建无参的存储过程(四十七)
  20. 金融python入门书籍推荐_学习金融工程,有哪些推荐的入门书籍?

热门文章

  1. 中国动力电池PACK行业商业模式创新及发展方向建议报告2022-2028年版
  2. 大部分公司可在两天内被渗透攻击
  3. 解析FeatureCollection(Geotools对geojson操作出现的问题)
  4. 【Jetpack】学穿:LiveData →
  5. 走进元宇宙:是什么、有什么、为什么?
  6. 我的择业思考:在AI最火的时候来到工业界!
  7. 分布式技术-Zookeeper
  8. Java面试题及答案2019版(上),springboot缓存技术
  9. 2019年北京摇号窗口办理地址、电话
  10. 内存不能为READ的原因