《图像几何校正与辐射校正.ppt》由会员分享,可在线阅读,更多相关《图像几何校正与辐射校正.ppt(77页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、 遥感数字图像校正与增强1本章主要内容1.数字图像基础2.图像辐射校正3.几何校正原理与方法4.图像增强处理28.3 几何校正遥感图象几何畸变遥感图象几何畸变遥感图象几何纠正方法遥感图象几何纠正方法3为什么要进行几何纠正为什么要进行几何纠正遥感调查分析结果遥感调查分析结果:一般是要求能满足量测和定位 要求的各类专题地图。利用多源数据多源数据进行计算机自动分类、地物特征的变化 监测等应用处理时,必须保证不同图像间的几何一致 性。几何误差的存在几何误差的存在4v定义:遥感图像上各地物的几何位置、形状、尺寸、方位等特征与在参照系统中的表达要求不一致时,即说明遥感图像发生了几何畸变。v遥感图像的总体变
2、形(相对于地面真实形态而言)是平移、缩放、旋转、偏扭、弯曲及其他变形综合作用的结果。8.1 遥感图像的几何变形v按照畸变的性质划分几何畸变可分为系统性畸变和随机性畸变。一、基本概念5系统性畸变系统性畸变是指遥感系统造成的畸变,这种畸变一般有一定的规律性,并且大小事先能够预测,例如扫描镜的结构方式和扫描速度等造成的畸变。随机性畸变随机性畸变是指大小不能事先预测、其出现带有随机性质的畸变,例如地形起伏造成的随地而异的几何偏差。随机畸变遥感图像几何畸变系统畸变6几何校正几何校正就是要校正成像过程所造成的各种几何畸变。几何校正分为两种:几何粗校正几何粗校正和和几何精校正几何精校正。几何粗校正几何粗校正
3、是针对引起畸变原因而进行的校正,这种畸变按照比较简单和相对固定的几何关系分布在图像中的,校正时只需将传感器原校准数据、遥感平台的位置以及卫星运行姿态等一系列测量数据代入理论校正公式即可。几何粗校正主要校正系统畸变。几何精校正几何精校正是利用控制点进行的几何校正,它是用一种数学模型来近似描述遥感图像的几何畸变过程,并利用畸变的遥感图像与标准地图之间的一些对应点(即控制点(GCP)求得这个几何畸变模型,然后利用此模型进行几何畸变校正,这种校正不考虑畸变的具体形成原因,而只考虑如何利用畸变模型来校正遥感图像。几何粗纠正几何精纠正随机畸变遥感图像几何畸变系统畸变7几何畸变几何畸变遥感器本身引起的畸变
4、外部因素引起的畸变 处理过程中引起的畸变二、几何变形的类型二、几何变形的类型 8.2 遥感图像的几何变形根据畸变产生的原因:根据畸变产生的原因:81、传感器本身引起的畸变、传感器本身引起的畸变传传感感器器本本身身引引起起的的几几何何畸畸变变与与遥遥感感器器的的结结构构、特特性性和和工作方式不同而异。这些因素主要包括:工作方式不同而异。这些因素主要包括:1)透镜的辐射方向畸变像差;透镜的辐射方向畸变像差;2)透镜的切线方向畸变像差;透镜的切线方向畸变像差;3)透镜的焦距误差;透镜的焦距误差;4)透镜的光轴与投影面不正交;透镜的光轴与投影面不正交;5)图像的投影面非平面;图像的投影面非平面;6)探
5、测元件排列不整齐;探测元件排列不整齐;7)采样速率的变化;采样速率的变化;8)采样时刻的偏差采样时刻的偏差;9)扫描镜的扫描速度变化扫描镜的扫描速度变化。二、几何变形的类型二、几何变形的类型9MSS举例举例:例如扫描形式成像的MSS,产生的几何畸变主要是由于扫描镜的非线性振动和其它一些偶然因素引起的。在地面上影响可达395米。全景畸变:10全景畸变全景畸变的图形变化情况的图形变化情况112、外部因素引起的畸变外部因素引起的畸变遥感平台位置和运动状态变化的影响地形起伏的影响地球表面曲率的影响大气折射的影响地球自转的影响二、几何变形的类型二、几何变形的类型121)遥感平台位置和运动状态变化的影响n
6、航高:当平台运动过程中受到力学因素影响标,或者说卫星运行的轨道本身就是椭圆的。航高始终发生变化,而传感器的扫描视场角不变,从而导致图像扫描行对应的地面长度地面长度发生变化。航高越向高处偏离,图像对应的地面越宽。2、外部因素引起的畸变13航速:卫星的椭圆轨道本身就导致了卫星飞行速度的不均匀,其他因素也可导致遥感平台航速的变化。航速快时,扫描带超前,航速慢时,扫描带滞后,由此可导致图像在卫星前进方向上(图像上下方向)的位置错动。俯仰:遥感平台的俯仰变化能引起图像上下方向的变化,即星下点俯时后移,仰时前移,发生行间位置错动.14翻滚:遥感平台姿态翻滚是指以前进方向为轴旋转了一个角度。可导致星下点在扫
7、描线方向偏移,使整个图像的行向翻滚角引起偏离的方向错动。偏航:指遥感平台在前进过程中,相对于原前进航向偏转了一个小角度,从而引起扫描行方向的变化,导致图像的倾斜畸变152 2)地形起伏的影响)地形起伏的影响n当地形存在起伏时,会产生局部像点的位移,使原本应是地面点的信号被同一位置上某高点的信号代替。由于高差的原因,实际像点P距像幅中心的距离相对于理想像点P0。距像幅中心的距离移动了r高差引起的像点位移高差引起的像点位移2、外部因素引起的畸变163)地球表面曲率的影响)地球表面曲率的影响n 地球是球体,严格地球是球体,严格说是椭球体,因此说是椭球体,因此地球表面是曲面。地球表面是曲面。这一曲面的
8、影响主这一曲面的影响主要表现在两个方面,要表现在两个方面,一是像点位置的移一是像点位置的移动,二是像元对应动,二是像元对应于地面宽度的不等。于地面宽度的不等。2、外部因素引起的畸变174)大气折射的影响)大气折射的影响n 大气对辐射的传播产生折射。由于大气的密度分布从下向上越来越小,折射率不断变化因此折射后的辐射传播不再是直线而是条曲线从而导致传感器接收的像点发生位移r.n大气折射影响大气折射影响2、外部因素引起的畸变185)地球自转的影响)地球自转的影响n 卫星前进过程中,传感器对地面扫描获得图像时,地球自转影卫星前进过程中,传感器对地面扫描获得图像时,地球自转影响较大,会产生影像偏离。因为
9、卫星自北向南运动,这时地球响较大,会产生影像偏离。因为卫星自北向南运动,这时地球自西向东自转。相对运动的结果,使卫星的星下位置逐渐产生自西向东自转。相对运动的结果,使卫星的星下位置逐渐产生偏离。偏离方向如下图所示,所以卫星图像经过校正后成为图偏离。偏离方向如下图所示,所以卫星图像经过校正后成为图c c的形态。的形态。地球自转引起偏离2、外部因素引起的畸变19遥感图像再处理过程中产生的误差,主遥感图像再处理过程中产生的误差,主要是由于处理设备产生的噪声引起的。要是由于处理设备产生的噪声引起的。传输、复制传输、复制 光学光学 数字数字3、处理过程中引起的畸变20n遥感图像的几何校正按照处理方式分为
10、光学纠正和数字纠正。n遥感图像的几何纠正就是将含有畸变的图像纳入到某种地图投影。对地面覆盖范围不大的单幅图像,一般以正射投影方式使其改正到地球切平面上。n光学纠正主要用于早期的遥感图像的处理中,现在的应用已经不多。除了对框幅式的航空照片(中心投影)可以进行比较严密的纠正以外,对于大多数动态获得的遥感影像只能进行近似的纠正。n主要介绍数字图像的几何精纠正。三三遥感图像的几何校正方法遥感图像的几何校正方法1、基本概念21两个基本环节:像元坐标变换像元坐标变换和像元灰度值重采样像元灰度值重采样校正思路校正思路(技术流程技术流程):22n校正前的图像,由于某种几何畸校正前的图像,由于某种几何畸变,图像
11、中变,图像中像元点间所对应的地面像元点间所对应的地面距离并不相等距离并不相等。n校正后的图像是由校正后的图像是由等间距的网格等间距的网格点点组成的,且以地面为标准,符合组成的,且以地面为标准,符合某种投影的均匀分布某种投影的均匀分布,图像中格网图像中格网的交点可以看作是像元的中心;的交点可以看作是像元的中心;校正的最终目的是确定校正后图像的行列数值,然后找到新图像中每一像元的亮度值。两个基本环节:像元坐标变换和像元灰度值重采样23(1)(1)确定输入图像和输出图像的坐标变换关系确定输入图像和输出图像的坐标变换关系 (2)(2)确定新的图像的边界确定新的图像的边界(3)(3)确定新确定新图图像的
12、分辨率像的分辨率(4)(4)灰度的重采灰度的重采样样2、校正过程24 数字图象几何纠正:通过计算机对离散结构的数字图像中的每一个像元逐个进行纠正处理的方法。这种方法能够精确地改正动态扫描图像所具备的各种误差。基本原理:利用图像坐标和地面坐标(另一图像坐标、地图坐标等)之间的数学关系,即输入图像和输出图像间的坐标转换关系实现几何校正。(1)确定输入图像和输出图像的坐标变换关系坐标变换关系2、校正过程25几何校正包括两个方面的内容几何校正包括两个方面的内容图像空间像元坐标的变换图像空间像元坐标的变换变换后的标准图像空间的各像元灰度值的计算。变换后的标准图像空间的各像元灰度值的计算。26直直接接纠纠
13、正正方方法法:从原始图像阵列出发,按行列的顺序依 次 对 每 个 原 始 图 像 像 元 点 位 用 变 换 函 数 F(x,y)(正解变换公式)求得它在新图像中的位置,并将该像元灰度值移置到新图像的对应位置上。间间接接纠纠正正法法:从空白的新图像阵列出发,按行列的顺序依次对新图像中每个像元点位用变换函数f(X,Y)(反解变换公式)求其在原始图像中的位置,然后把算得的原始图像点位上的灰度值赋予空白新图像相应的像元。纠正方法纠正方法:27图中(xp,yp)(XP,YP)分别是任意一个像元在原始图像和纠正后图像中的坐标。直接法(正解):间接法(反解):282.3.5遥感数字图像的多项式纠正遥感数字
14、图像的多项式纠正n多项式纠正的基本思想:图像的变性规律可以看作是多项式纠正的基本思想:图像的变性规律可以看作是平移、缩放、旋转、仿射、偏扭、弯曲等形变的合成。平移、缩放、旋转、仿射、偏扭、弯曲等形变的合成。一般的公式为:一般的公式为:2)直接法直接法1)间接法间接法利利用用有有限限的的控控制制点点的的已已知知坐坐标标,解解求求多多项项式式的的系系数数,确确定定变变换换函函数。然后将各个像元带入多项式进行计算,得到纠正后的坐标。数。然后将各个像元带入多项式进行计算,得到纠正后的坐标。29实际计算时常采用二元二次多项式:实际计算时常采用二元二次多项式:在这个方程组中有12个系数,需列12个方程才能
15、解出,因此需要6个已知的对应点,即这6个点的(u,v)与(x,y)均已知,这些已知坐标的对应点称为控制点(GCP)6个点只是解算方程组的理论最低数,实际工作中为提高校正精度需大量增加控制点数,这时就有了多余条件,可采用最小二乘法求解。3、遥感数字图像的多项式纠正30控制点的选取1)数目确定控制点数目的最低限是按未知系数的多少来确定的。求二次多项式有12个系数,需要12个方程(6个控制点)。依次类推,三次多项式至少需要10个控制点,n次多项式,控制点的最少数目为(n+1)(n+2)2。实际工作中,在条件允许的情况下,控制点数的选取都要大于最低数很多。3、遥感数字图像的多项式纠正311)1)表征空
16、间位置的可靠性,道路交叉点,标志物,表征空间位置的可靠性,道路交叉点,标志物,水域水域的边界的边界,山顶,小岛中心,机场等。,山顶,小岛中心,机场等。同名控制点要在图像上均匀分布;同名控制点要在图像上均匀分布;清楚辨认;清楚辨认;数量应当超过多项式系数的个数(数量应当超过多项式系数的个数((n+1)*(n+2)/2)。)。当控制点的个数超过多项式的系数个数时,采用最小当控制点的个数超过多项式的系数个数时,采用最小2乘法进行系数的确定,使得到的系数最佳。乘法进行系数的确定,使得到的系数最佳。2.3.5遥感数字图像的多项式纠正遥感数字图像的多项式纠正(续(续1)控制点的选择原则控制点的选择原则:3
17、2(2)确定新的图像的边界确定新的图像的边界n纠正后图像和原始图像的形状、大小、方向都不一样。纠正后图像和原始图像的形状、大小、方向都不一样。所以在纠正过程的实施之前,必须首先确定新图像的所以在纠正过程的实施之前,必须首先确定新图像的大小范围。大小范围。33(2)确定新的图像的边界(续确定新的图像的边界(续1 1)lX1=min(Xa,Xb,Xc,Xd)lX2=max(Xa,Xb,Xc,Xd)lY1=min(Ya,Yb,Yc,YXd)lY2=max(Ya,Yb,Yc,Yd)34(3)确定新图像的分辨率确定新图像的分辨率n目的是确定新图像目的是确定新图像宽度宽度和和高度高度;n根据精度要求,在新
18、图像的范围内,划分网格,根据精度要求,在新图像的范围内,划分网格,每个网格点就是一个像元。每个网格点就是一个像元。n新图像的行数新图像的行数M(Y2-Y1)/Y+1;n新图像的列数新图像的列数N(X2-X1)/XX+1;n新图像的任意一个像元的坐标由它的行列号唯新图像的任意一个像元的坐标由它的行列号唯一确定。一确定。35(4 4)灰度的重采样灰度的重采样n纠纠正正后后的的新新图图像像的的每每一一个个像像元元,根根据据变变换换函函数数,可可以以得得到到它它在在原原始始图图像像上上的的位位置置。如如果果求求得得的的位位置置为为整数,则该位置处的像元灰度就是新图像的灰度值。整数,则该位置处的像元灰度
19、就是新图像的灰度值。n如果位置不为整数如果位置不为整数,则有几种方法:则有几种方法:n1)最近邻法最近邻法n2)双线性内插法双线性内插法n3 3)三次卷积法)三次卷积法362.3.4 2.3.4 灰度的重采样(续灰度的重采样(续1 1)1)最近邻法:最近邻法:距离实际位置最近的像元的灰度值作为距离实际位置最近的像元的灰度值作为输出图像像元的灰度值;输出图像像元的灰度值;37原始图像纠正后图像(最邻近插值)最邻近法纠正效果最邻近法纠正效果最邻近法纠正效果最邻近法纠正效果382.3.4 2.3.4 灰度的重采样(续灰度的重采样(续2 2)2)双线性法)双线性法:以实际位置临近的以实际位置临近的4个
20、像元值,确定输出个像元值,确定输出像元的灰度值。公式为:像元的灰度值。公式为:式中,g(m,n)为输出像元灰度值gi为邻近点i的灰度值pi为邻近点对投影点的权重(pi=1/di,di表示邻近点到投影点的距离,最近者权重最大39双线性插值效果原始图像纠正(双线性插值)402.3.4 2.3.4 灰度的重采样(续灰度的重采样(续3 3)3 3)三次卷积法)三次卷积法以实际位置临近的以实际位置临近的16个像元值,确定输个像元值,确定输出像元的灰度值。公式为:出像元的灰度值。公式为:三次样条三次样条函数式中,g(m,n)为输出像元灰度值gi为邻近点i的灰度值pi为邻近点对投影投影点的权重(pi=1/d
21、i,di表示邻近点到投影点的距离,最近者权重最大)41原始图像几何纠正(三次卷积)三次卷积三次卷积三次卷积三次卷积法处理效果法处理效果法处理效果法处理效果422.3.4 2.3.4 灰度的重采样(续灰度的重采样(续4 4)几种采样方法的优缺点:几种采样方法的优缺点:1)最近邻法:算法简单且保持原光谱信息不变;缺)最近邻法:算法简单且保持原光谱信息不变;缺点是几何精度较差,图像灰度具有不连续性,边界出点是几何精度较差,图像灰度具有不连续性,边界出现锯齿状。现锯齿状。2)双线性插值:计算较简单,图像灰度具有连续性且双线性插值:计算较简单,图像灰度具有连续性且采样精度比较精确;缺点是细节丧失采样精度
22、比较精确;缺点是细节丧失3)三次卷积法:计算量大,图像灰度具有连续性且采三次卷积法:计算量大,图像灰度具有连续性且采样精度比较精确样精度比较精确43数字图象的纠正过程数字图象的纠正过程总上所述:总上所述:纠正的函数可有多种选择:多项式方法、共线方程纠正的函数可有多种选择:多项式方法、共线方程方法、随机场内插方法等等。其中多项式方法的应用最方法、随机场内插方法等等。其中多项式方法的应用最为普遍。为普遍。44四四几何校正类型几何校正类型n从影像到地图的校正从影像到地图的校正n地图控制点获取方法:地图控制点获取方法:简单的直尺、数字化地图、校正后的数字简单的直尺、数字化地图、校正后的数字正射影像、正
23、射影像、GPSn从影像到影像的校正从影像到影像的校正n同名地物进行匹配同名地物进行匹配n混合混合45五、几何校正操作过程五、几何校正操作过程n选择合适的参考平面图选择合适的参考平面图n采集地面控制点采集地面控制点n评价全部控制点的误差评价全部控制点的误差RMSE,确定最优几何,确定最优几何校正系数校正系数n利用坐标变换和灰度重采样方法输出校正后图利用坐标变换和灰度重采样方法输出校正后图像像46开始开始显显示示图图形形文件文件启启动动几几何何校校正正模型模型采采集集地地面面控控制制点点计算转换模型计算转换模型图图像像重重采采样样检检验验校校正正结结果果结束结束 几何校正流程图几何校正流程图47补
24、充补充:最小二乘法最小二乘法n最小二乘法最早称为回归分析法。由著名的英国生物最小二乘法最早称为回归分析法。由著名的英国生物学家、统计学家道尔顿(学家、统计学家道尔顿(F.GalltonF.Gallton)所创。所创。(道尔顿研究英国男子中父亲们的身高与儿子们的身道尔顿研究英国男子中父亲们的身高与儿子们的身高之间的关系时,创立了回归分析法。高之间的关系时,创立了回归分析法。)n探索变量之间关系最重要的方法,用以找出变量之间探索变量之间关系最重要的方法,用以找出变量之间关系的具体表现形式。关系的具体表现形式。n后来,回归分析法从其方法的数学原理后来,回归分析法从其方法的数学原理误差平方误差平方和最
25、小(和最小(二乘是平方的意思二乘是平方的意思)出发,改称为最小二乘)出发,改称为最小二乘法。法。482.3.6最小二乘法(续最小二乘法(续1)1.为了精确地描述为了精确地描述Y Y与与X X之间的关系,必须使用这两个变量的所有观察值,之间的关系,必须使用这两个变量的所有观察值,才不至于以才不至于以“点点”概面(作到同步与全面)。概面(作到同步与全面)。2 2Y Y与与X X之间是否是直线关系(用协方差或相关系数衡量)?若是,将用之间是否是直线关系(用协方差或相关系数衡量)?若是,将用一条直线描述它们之间的关系。一条直线描述它们之间的关系。3 3在在Y Y与与X X的散点图上画出直线的方法很多。
26、的散点图上画出直线的方法很多。任务?任务?找出一条能够最好地描述找出一条能够最好地描述Y Y与与X X(代表所有点)之间关系的直(代表所有点)之间关系的直线。线。4 4什么是最好?什么是最好?找出判断找出判断“最好最好”的原则。的原则。最好指的是找这么一条直线,使得所有点到该直线的纵向距离的和(平最好指的是找这么一条直线,使得所有点到该直线的纵向距离的和(平方和)最小。方和)最小。5 5三种距离三种距离49点到直线的距离点到直线的距离点到直线的垂直线的长度。点到直线的垂直线的长度。横向距离横向距离点沿(平行)点沿(平行)X X轴方向到直线的距离。轴方向到直线的距离。纵向距离纵向距离点沿(平行)
27、点沿(平行)Y Y轴方向到直线的距离。也就是轴方向到直线的距离。也就是实际观点的实际观点的Y Y坐标减去根据直线方程计算出来的坐标减去根据直线方程计算出来的YY的拟合值。的拟合值。2.3.6最小二乘法(续最小二乘法(续2)50n纵向距离是纵向距离是Y Y的实际值与拟合值之差,差异大拟合不好,的实际值与拟合值之差,差异大拟合不好,差异小拟合好。差异小拟合好。n将所有纵向距离平方后相加,即得误差平方和。所以,将所有纵向距离平方后相加,即得误差平方和。所以,“最好最好”直线就是使误差平方和最小的直线。直线就是使误差平方和最小的直线。n运用求极值的原理,将求最好拟合直线问题转换为求运用求极值的原理,将
28、求最好拟合直线问题转换为求误差平方和最小的问题。误差平方和最小的问题。2.3.6最小二乘法(续最小二乘法(续3)512.3.6最小二乘法(续最小二乘法(续3)522.3.6最小二乘法(续最小二乘法(续4)最小二乘法的性质最小二乘法的性质:之一:即拟合直线过之一:即拟合直线过y和和x的平均数点。的平均数点。之二:之二:残差与自变量残差与自变量x的乘积和等于的乘积和等于0,即残差与自变量,即残差与自变量x两两者不相关。者不相关。53遥感图像的辐射纠正遥感图像的辐射纠正1为什么要进行辐射纠正?为什么要进行辐射纠正?1)传感器本身的特性,传感器本身的特性,2)大气对于电磁辐射的衰减;(散射、)大气对于
29、电磁辐射的衰减;(散射、反射和吸收)反射和吸收)3)地形因子的影响地形因子的影响阴影阴影4)其它生态环境因子其它生态环境因子形成形成“同物异谱,异物同谱同物异谱,异物同谱”现象。现象。图图像像不不能能全全部部真真实实地地反反映映不不同同地地物物地地特特征征,影影响响了了数数字字图图像的质量。像的质量。54552 辐射纠正类型n1)遥感器纠正:遥感器的设计n2)大气辐射纠正:n3)地形辐射纠正:需要DEM56大气校正后,图像的亮度显著降低;对比度增强572.1遥感探测器系统误差校正遥感探测器系统误差校正n比较常见的遥感系统引入的辐射误差有:比较常见的遥感系统引入的辐射误差有:n随机坏像元随机坏像
30、元n行起始行起始 /终止问题终止问题n行或列缺失行或列缺失n行或列部分缺失行或列部分缺失n行或列条纹行或列条纹582.1遥感探测器系统误差校正遥感探测器系统误差校正n随机坏像元随机坏像元n使用简单的阈值法判断坏像元(比如等于使用简单的阈值法判断坏像元(比如等于0 0或小于或小于0 0)n使用周围使用周围8 8像元的均值代替坏像元像元的均值代替坏像元592.1遥感探测器系统误差校正遥感探测器系统误差校正n行起始行起始/终止问题终止问题n造成扫描线错位,不是系统性的,需要手工进行移造成扫描线错位,不是系统性的,需要手工进行移位。位。602.1遥感探测器系统误差校正遥感探测器系统误差校正n行或列缺失
31、行或列缺失n行或列部分缺失行或列部分缺失n利用上下相邻两行的平均值代替缺失行的像元利用上下相邻两行的平均值代替缺失行的像元值值3445763538405880773856787860546076612.1遥感探测器系统误差校正遥感探测器系统误差校正n行或列条纹行或列条纹n确定误定标扫描行,即明确定误定标扫描行,即明显过亮或过暗的行显过亮或过暗的行n利用相同地物有相同的辐利用相同地物有相同的辐射等级的原理,采用线性射等级的原理,采用线性定标法调整所有坏扫描行定标法调整所有坏扫描行的辐射值。的辐射值。nDNDN好好 DNDN坏坏 A+B A+B622.2 大气辐射纠正方法n一、绝对大气辐射校正n1
32、辐射传输过程n2辐射传输模型n3经验线定标法n二、相对大气辐射校正n1直方图最小值去除法(DOS)n2回归法n三、坡度坡向影响校正n1、余弦法n2、C校正n3、统计经验校正631基于辐射传输过程的绝对辐射校正方法基于辐射传输过程的绝对辐射校正方法传感器记录的总辐射传感器记录的总辐射n根据传输过程根据传输过程n根据图像像元值根据图像像元值64遥感系统接收的多种程辐射遥感系统接收的多种程辐射研究目标反射研究目标反射邻域反射邻域反射1351,3,524太阳辐照度太阳辐照度天空漫射天空漫射辐照度辐照度652基于辐射传输模型的大气校正基于辐射传输模型的大气校正nATCOR(The ATmospheric
33、 CORrection)nFLAASH(Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes)nATREM(The ATmospheric REMovel)nACORN(The Atmospheric CORrection Now)662基于辐射传输模型的大气校正基于辐射传输模型的大气校正辐射传输模型所需参数辐射传输模型所需参数n遥感影像的经纬度遥感影像的经纬度n遥感数据采集的日期和时间遥感数据采集的日期和时间n影像获取高度影像获取高度n影像景的平均高度(如:海拔影像景的平均高度(如:海拔200m200m)n大气模式(如:
34、中纬度夏季,中纬度冬季,热带)大气模式(如:中纬度夏季,中纬度冬季,热带)n辐射定标后的影像数据(即数据单位必须是辐射定标后的影像数据(即数据单位必须是 )n各个波段的信息(波段均值和半幅全宽)各个波段的信息(波段均值和半幅全宽)n获取遥感数据时的当地大气能见度(如:获取遥感数据时的当地大气能见度(如:10km10km)672基于辐射传输模型的大气校正基于辐射传输模型的大气校正nATCOR(TheAtmosphericCORrection)n包括包括ACTOR2用于平坦地区用于平坦地区ACTOR3用于山区用于山区ACTOR4用于亚轨道遥感系统获取的数据用于亚轨道遥感系统获取的数据n用用MODT
35、RAN4+辐射传输代码计算大气校正函数辐射传输代码计算大气校正函数(程辐射、大气透射、太阳的直射、漫射辐射通量)(程辐射、大气透射、太阳的直射、漫射辐射通量)nERDAS软件扩展模块软件扩展模块682基于辐射传输模型的大气校正基于辐射传输模型的大气校正nFLAASH(FastLine-of-sightAtmosphericAnalysisofSpectralHypercubes)n用用MODTRAN4+辐射传输代码对影像逐像元地校正辐射传输代码对影像逐像元地校正大气中的水汽、氧气、二氧化碳、甲烷、臭氧、气大气中的水汽、氧气、二氧化碳、甲烷、臭氧、气溶胶等散射的影响。溶胶等散射的影响。nENVI
36、软件里自带软件里自带692基于辐射传输模型的大气校正基于辐射传输模型的大气校正ATREM(The ATmospheric REMovel)n采用6S辐射传输模型和用户指定的气溶胶模型计算分子散射n利用大气模式(温度、压力、水汽垂直分布)计算大气吸收n水汽总量通过水汽波段和三通道比值方法从高光谱数据中逐像元获得n得用以上值获得经验辐射校正后的表面反射率值703经验线性定标法经验线性定标法n经验线性定标公式经验线性定标公式nDNDN,像元灰度值,像元灰度值,像元表面反射率,像元表面反射率,K,波段,波段nA A通常称为增益、通常称为增益、B B称为偏移量称为偏移量n利用与图像同时获取的地面定标物体
37、的反射率,带入利用与图像同时获取的地面定标物体的反射率,带入公式得到公式得到A A和和B B,n定标物体通常是地面布置的白板和黑体,或者是沙堆定标物体通常是地面布置的白板和黑体,或者是沙堆和较深的平静湖水和较深的平静湖水71二、相对辐射校正1)直方图最小值去除法(DOS)基本思想:影像中总有一些地物,其辐射亮度或反射率接近0,则这些地物的像元值应该为零;但是实际影像中,这些地物的像元不为零;认为该值是大气散射、反射等进入的程辐射值。洁净并且具有一定深度和面积的水体深暗地形阴影72n1)直方图最小值去除法(DOS)作出影像直方图找出直方图最小值(影像最小值)每个像元值都减去本波段的最小值732)回归法基本思想:n以红外波段最低值校正可见光波段。随波长增长,大气散射作用减弱,红外波段受大气影响最小。可以把红外图象当作无散射的标准图象,将其他波段图像与之比较,其差值就是需要矫正的散射辐射值。n选择可见光和红外波段进行2维散点图,建立线性回归方程。742)回归法(续1)75三、地形辐射纠正n需要DEMn余弦校正n统计经验校正法nC校正n太阳天顶角,i太阳入射角,LH,水平面辐射,LT坡面辐射,c,调整因子坡坡遥感系统遥感系统i9076三、地形辐射纠正n简单的处理方法n比值法:有效消除阴影的影响。比值法:有效消除阴影的影响。77
限制150内