辐射校正:像元DN值—辐照度—表观反射率—地表反射率
狭义的辐射校正是指将原始图像的DN值转化为有物理量纲的观测目标大气上界的辐射亮度值或者表观反射率;广义的除了狭义的步骤外还包括去除大气影响,也就是说还包括大气校正的过程,最终需要计算出地表目标真实的反射亮度或者是反射率。
大气校正就是将辐射亮度或者表观反射率转换为地表反射率,目的是消除大气散射、吸收、反射引起的误差,是辐射校正的一个方面。
狭义的辐射校正:
输入:像元DN值、辐射校正参数
输出:像元星上辐照度、表观反射率
一、由像元DN值到辐照度
利用影像头文件中的记录的辐射校正参数可方便地计算出地物在大气顶部的星上辐照度,式中DN是原始影像象元灰度值;gain和bias是传感器对应的增益和偏差值,可从头文件中获取,计算公式如下:
二、由辐照度到表观反射率
式中,
ρ* 是地物大气层顶表观反射率;
d是日地天文单位距离,一般取值1;
E0太阳辐照度,每个波段对应不同的E0;
θ是太阳天顶角,每景影像获取时的太阳天顶角也可在头文件中找到。
广义的辐射校正:
除了以上两个步骤之外,还有以下操作:
输入值:像元的气溶胶光学厚度AOD,像元星上辐照度值L*、查找表LUT
输出值:像元地表反射率
原理
假设地面为朗伯面,则在传感器处接受到的单像元光谱辐射亮度可以表示为:
式中,
L* 为在传感器处接收到的单个像元的辐射亮度;
ρ为像元地表反射率;
S为大气球面反照率;
La* 为大气程辐射进入传感器的辐射亮度;
Fd为从太阳到地面的下行辐射;
T为从地面到传感器的上行辐射。
求地表反射率,需要对上式进行变形:(下面这两个公式要表达的意思一样)
进一步理解
对遥感影像进行大气纠正,由辐射亮度L* 反推出地表反射率ρ一般分两步:
第一步是获取相应的大气参数即La*、Fd、T、S。
(Fd和T一直是乘积的形式组合在一起,因此可以视为一个参数,记为F);
第二步是根据上述经验公式反推出地表反射率ρ。
注意:
第一步中,确定气溶胶光学厚度值后,就可以确定这三个大气参数了(La*、F、S)。如果直接通过求解辐射传输方程来获取这几个值非常困难,因为需要对一个没有解析解的方程进行微积分运算 。在实际计算中 ,一般采取单点迭代算法或高斯-塞德儿算法 ,但这种计算非常耗时,而且不适合在大范围的研究(如遥感影像的大气纠正)中应用 ;另外一个方法是简化辐射传输方程求近似解 ,如二流近似算法(two-stream)和四流近似算法(four-stream),求近似解的方法一般效率较高,但精度有限;目前比较有效的方法被称为查找表法、在进行大气纠正之前利用辐射传输方程如6S模型或 Modtran模型等建立好一套查找表,在对影像进行大气纠正过程中,调用这些表进行插值运算求得参数 La * 、F、S的值。 具体的算法流程可以参考:赵春江等, 基于6S模型的遥感影像逐像元大气纠正算法. 光学技术, 2007. 33(1): 第11-15页.
辐射校正:像元DN值—辐照度—表观反射率—地表反射率相关推荐
- 计算GF1和GF6的表观反射率和投影问题
记录一下怕自己忘了 主要是之前一种用classic做出来的一直没有坐标信息 很烦 于是想了此下策 前提是要装了中国卫星的插件比较方便 首先用ENVI将需要的波段转成envi标准格式 不用classic ...
- 夜间灯光数据dn值_中国区域DMSP/OLS夜间灯光影像的校正
DMSP/OLS夜间灯光影像由于传感器未经过星上定标,导致长时间序列的各期影像之间缺乏连续性和可比性,且由于光谱分辨率的限制,使得大量的城市中心像元出现DN值为63的集聚饱和现象等问题,因此学者们对以 ...
- 【技术类】【ArcGIS对国产卫星的支持2:高分一号卫星】篇5、实时表观反射率产品生产
表观反射率就是指大气层顶的反射率,辐射定标的结果之一,又称视反射率.英文表示为:apparent reflectance (=地表反射率+大气反射率) 分发商提供的数据基本都是灰度值(DN值),必须转 ...
- 辐射亮度与表观反射率的转换
文章目录 DN值 辐射亮度 表观反射率 转换公式 亮度温度 其他资料 补充相关资料 DN值 就是影像灰度值 辐射亮度 略 表观反射率 大气层顶的反射率 转换公式 DN值转辐射亮度: Lλ=Gain∗D ...
- IDL调用6S文件生成数据,及表观反射率有关二维图像
;;=============================================================================================== ;; ...
- 遥感中的表观反射率、反射率、反照率等之间的区别与联系
遥感反射率的定义:地物表面反射能量与到达地物表面的入射能量的比值. 遥感表观反射率的定义:地物表面反射能量与近地表太阳入射能量的比值. 大气校正就是将辐射亮度或者表观反射率转换为地表实际反射率,目的是 ...
- 反射率(表观反射率、地表反射率)
反射率是物体表面所能反射的辐射量和它所接受的辐射量的比值,一般在[0,1]范围,有的时候为了储存方便而扩大一定的倍数,如放大一万倍[0,10000].一些材料可以通过他们的反射光谱来识别,因此,为了更 ...
- 利用Python和GDAL把MODIS的HDF格式表观反射率产品MOD09A1转换为TIFF格式并重投影到与LANDSAT8-OLI一致
解释 1.MOD09A1是表观反射率产品,是HDF格式并且投影和常用的Landsat数据不一致.为了统一反演时使用,想要将该产品转换为TIFF格式并重投影.特别是要批量的进行转换和重投影. 2.采用P ...
- (三)Landat_5 TM 遥感影像辐射标定以及大气表观辐射率计算
前言 一.辐射定标 辐射定标:将记录的原始DN值转换为辐射值,也就是将图像的亮度灰度值转换为绝对辐射亮度. 如下便是通过imtool函数查看Landsat_5 TM band5的原始亮度灰度值,其整体 ...
- 一体式表观高光谱仪Triplet AOP
一体式表观高光谱仪Triplet AOP是一款超便携.智能化.高性价比的一体式表观高光谱仪,安装了两个高光谱辐亮度传感器和一个高光谱辐照度传感器.配有小型自稳平台,可以保持"水面之上法&qu ...
最新文章
- Python 基础 - Day 5 Learning Note - 模块 之 标准库:xml (9)
- python画直方图代码-Python使用Plotly绘图工具,绘制直方图
- 线段树3——一些例题的题解
- 51单片机 自动重装载值计算
- CSS深入理解流体特性和BFC特性下多栏自适应布局
- 机器学习与分布式机器学习_机器学习的歧义
- python教程循环语句,Python基础教程之循环语句(for、while和嵌套循环)
- leetcode212. 单词搜索 II
- 声压级 matlab,语音信号处理教程(二)声音的声压级和响度
- vue中rules校验是验证首字符_小白也能秒懂Vue源码中那些精细设计(选项处理)...
- html十css十javascript,html代码,css代码,JavaScript代码, | 重蔚自留地
- webservice python开发接口_基于Python的Webservice开发(四)-泛微OA的SOAP接口
- Java 中的接口有什么作用?为什么不直接用类?
- ubuntu-Linux系统读取USB摄像头数据(uvc)
- 服务器监控报警系统软件设计,Monitor监控报警系统
- HTML 微信 获取微信ID,在微信浏览器里获取Cookie,获取微信用户的openId
- 各种颜色对应的RGB值
- 卷积 对图像进行卷积操作 卷积神经网络
- VueJs中的reactive函数
- iphoneX的safari浏览器操作栏隐藏时兼容底部Home键
热门文章
- NanoPi M1 Plus:Allwinner H3 搭建编译环境、编译主线Uboot、编译主线Linux Kernel 4.14、加Matrix2‘8屏、加Qt
- java excel转pdf有乱码_word,ppt,txt,excel转换pdf 解决中文乱码
- 华为BFD的配置规范
- C语言 二进制文件读写实例讲解
- 【April Fools Day Contest 2014F】【愚人节脑洞 OEIS大法】000001 输出序列某一项
- 旅通软件:旅行社管理系统怎么选?
- linux 联机游戏下载,星露谷物语多人联机版
- html语言代码 输入文字,html语言教程文字代码:
- 《用户故事与敏捷方法》 笔记
- MindManager思维导图 PRDBRD写作方法与技巧 产品规划模板下载