第四章遥感数据预处理-影像校正

遥感图像预处理遥感成像过程真实影像-----畸变影像遥感影像畸变辐射畸变：指遥感传感器在接收来自地物的电磁波辐射能时，电磁波在大气层中传输和传感器测量中受到遥感传感器本身特性、地物光照条件（地形影响和太阳高度角影响）以及大气作用等影响，而导致的遥感传感器测量值与地物实际的光谱辐射率的不一致。

几何畸变：当原始图像上各地物的几何位置、形状、尺寸、方位等特征与在参照系统中的表达要求不一致时，就产生几何畸变。

遥感影像辐射校正遥感辐射畸变辐射畸变的来源辐射畸变的主要来源传感器影响依赖于太阳辐射波长大气影响地形影响常见的遥感传感器系统误差随机坏像元行或列缺失行或列条纹随机坏像元的判定和消除由传感器的噪声或磁带等部件的误码率造成的，其特点是孤立和分散的，往往和周围的亮度值有明显的差别，并且彼此不相关。斑点可以通过将图像像元亮度值同它的邻近像无亮度值进行比较来判定。行列缺失的判定和消除在扫描图像中出现的与辐射信息无关的线条噪声，其表现为图像上的部分扫描行或线段的亮度值不反映地物的辐射，并与上下的亮度截然不同。部分数据缺失LandsatMSS传感器的衰减，所引起的条带噪声参考文献ZhuXL,LiuDS,ChenJ.AnewgeostatisticalapproachforfillinggapsinLandsatETMplusSLC-offimages[J].RemoteSensingofEnvironment.2012,124:49-60.ChenJ,ZhuX,VogelmannJE,GaoF,JinS.AsimpleandeffectivemethodforfillinggapsinLandsatETM+SLC-offimages[J].RemoteSensingofEnvironment.2011,115(4):1053-64.原始影像影像条带消除之后的效果遥感图像的大气校正地物发射的电磁能量在经过大气层时，大气对电磁波的吸收和散射对遥感图像造成影响，其中主要是大气散射，大气校正即消除大气散射对辐射失真的影响。包括大气的消光（吸收和辐射）、天空光照射、路径辐射。绝对校正示意图经验线法与卫星扫描同步进行野外波谱测试，将地面测量结果与卫星影像对应像元亮度值进行回归分析，回归方程为：常用大气辐射校正模型辐射校正算例辐射校正算例-续辐射校正算例-续相对辐射校正适用范围：归一化单时相遥感影像不同波段间的强度将多时相影像遥感数据各个波段的强度归一化到某一选定的标准影像上直方图校正方法

从图像像元亮度值中减去一个辐射偏置量（LP），辐射偏置量等于图像直方图中最小的辐射亮度值。

前提（假设）：水体（或阴影）等物体的灰度值为0，大气散射导致图像上这些物体的灰度值不为0（辐射偏置量）

暗物体法（Dark-objectmethod）回归分析法原理：假设图像中某像元反射率已知或固定，利用这些像元反射率和各波段反射率之间的线性关系来对整幅图像进行归一化方法：在不受大气影响的波段（如TM5或7）和待校正的某一波段图像中，选择一系列目标，将每个目标的两个待比较的波段灰度值提取出来进行回归分析，建立线性回归方程。Y：待校正波段的图像亮度值X：不受大气影响波段的图像亮度值怎么求a,b？辐射校正回顾辐射定标：传感器探测值的标定过程，用以确定传感器入口处的准确辐射值。辐射校正：消除或改正遥感图像成像过程中附加在传感器输出辐射能量中的各种噪声。大气校正：指大气散射校正，即消除大气散射对辐射失真的影响。包括大气的消光（吸收和辐射）、天空光照射、路径辐射。遥感几何校正中心投影定义：凡空间任意点A（物点）与一固定点S（投影中心）连成的直线或延长线（即中心光线）被一个平面（像平面）所截，则此直线与平面的交点a（像点）称为A点的中心投影。

从投影上而言，航空像片（正片）的位置，等于以投影中心为圆心，以焦距f为半径，将P旋转至P’（下图），P’即为正像的位置。投影距离的影响中心投影：投影距离不同或焦距不同则像片的比例尺也不同。垂直投影：投影距离不同与像片比例尺无关。（不存在焦距）投影倾斜面的影响中心投影：投影面的倾斜造成同一个像片不同部位比例尺的差异。垂直投影：不存在投影面的倾斜。地形起伏的影响中心投影：地形起伏造成像点位移。垂直投影：不存在像点位移。卫星遥感数据几何畸变几何校正方法

①图像空间坐标变换；首先建立图像像点坐标（行、列号）和物方（或参考图）对应点坐标间的映射关系，解求映射关系中的未知参数，然后根据映射关系对图像各个像素坐标进行校正；②确定各像素的灰度值（灰度内插）。

图像几何校正的基本方法是先建立几何校正的数学模型；其次利用已知条件确定模型参数；最后根据模型对图像进行几何校正。通常分两步：数字图像几何纠正的主要处理过程准备工作输入原始数字影像建立纠正变换函数确定输出影像范围像元坐标变换像元亮度值重采样输出纠正后的图像准备工作：图像、地图、大地测量资料、平台轨道参数、传感器参数、控制点的选择；（具体内容可选）纠正变换函数建立：输入和输出图像间的坐标变换关系；如多项式法、共线方程法等。遥感几何校正-校正模型校正模型通用校正模型

几何多项式DLT模型(DirectLinearTransformation)RPC模型(RationalFuntion)物理校正模型 共线方程多项式模型

根据影像变形情况选取校正模型，不同模型需要控制点数目不同，一阶多项式几何校正(理论最小值)：3个控制点；二阶多项式几何校正(理论最小值)：6个控制点；三阶多项式几何校正(理论最小值)：10个控制点；四阶多项式几何校正(理论最小值)：15个控制点；五阶多项式几何校正(理论最小值)：21个控制点；DLT模型

可以认为是共线方程的变化式x=(a1+a2X+a3Y+a4Z)/(c1+c2X+c3Y+c4Z)y=(b1+b2X+b3Y+b4Z)/(c1+c2X+c3Y+c4Z)RationalPolynomialCamera(RPC)model r’=fr(φ’,λ’,h’)/gr(φ’,λ’,h’)c’=fc(φ’,λ’,h’)/gc(φ’,λ’,h’)x’

…normalisedco-ordinatesx’=(x-x_offset)/x_scaleRationalpolynomialfunctionsusually,3rdorderpolynomialsareused

fr=a1+a2φ’+a3λ’+a4h’+a5φ’λ’+a6λ’h’+a7φ’h’+a8λ’2+a9φ’2+

+a10h’2+a11φ’λ’h’+a12λ’3+a13φ’2λ’+a14λ’h’2+a15φ’λ’2+

+a16φ’3+a17φ’h’2+a18λ’2h’+a19φ’2λ’+a20h’3gr=b1+b2φ’+b3λ’+b4h’+b5φ’λ’+b6λ’h’+b7φ’h’+b8λ’2+b9φ’2+

+b10h’2+b11φ’λ’h’+b12λ’3+b13φ’2λ’+b14λ’h’2+b15φ’λ’2+

+b16φ’3+b17φ’h’2+b18λ’2h’+b19φ’2λ’+b20h’3similarly,RPCcoefficientsc1,…,c20,d1,…,d20infunctionsfcandgc传感器校正物理模型航空影像--共线方程SPOT卫星影像IKONOS/QuickBird卫星影像共线方程XZYXsYsDSxzyXYXA-XsYA-YsNZAZA-ZsZsa地面点(XA-Xs,YA-Ys,ZA-Zs)像点（X,Y,Z）共线方程

影像坐标、地面坐标以及外方位参数之间的关系：

x,y:影像坐标；X,Y,Z:地面坐标；XsYsZs：摄影中心的地面坐标。f:像片主距。a,b,c由三个角元素定义的3×3旋转矩阵的系数，目的将影像坐标换成地面坐标系统。遥感几何校正-控制点选取控制点的选取在影像上可以分辨并能在地图上精确定位的地表位置控制点的选取来源地形图几何校正后的数字正射影像GPS野外测量获取GCP坐标控制点的选择地面控制点（GCP，GroundControlPoint）:一些特定的像元，其地图坐标或其它输出坐标为已知人工地物线性地物交叉点不易随时间变化的目标大比例尺的图像：道路交叉点、机场跑道、建筑物小比例尺的图像：城区、一些线性地物交叉点（河流、道路）分布：较均匀分布与图像范围内，保证足够数量1.多项式纠正法的精度与地面控制点（GCP）的精度、分布、数量及纠正范围有关；GCP的位置精度越高，则几何纠正的精度越高；GCP的个数不少于多项式的系数个数；适当增加GCP的个数，可以提高几何纠正的精度。20-30个GCP，一般可以满足需求2GCP分布应尽可能在整幅图像内均匀分布，否则在GCP密集区精度较高，在GCP分布稀疏区出现较大误差控制点的选择配准点的选取控制点选取的标准校正模型约定控制点数目控制点残差评价控制点选取好坏遥感几何校正-数据重采样影像重采样影像重采样影像重采样原始影像灰度表面最近邻内插法双线性内插法三次内插法像素灰度内插法效果比较遥感配准及镶嵌

影像配准

(Matching)是将同一地区的不同特性的相关影像(如不同传感器，不同日期，不同波段或传感器在不同位置获取的同一地区地物)在几何上互相匹配，即实现影像与影像间地理坐标及像元空间分辨率上的统一。

相对配准

选择某一卫星数据作为参考图象，将其他卫星数据与之配准，简称图像对图像的配准。

绝对配准

将所有的图像分别校正到地图坐标系下。精纠正后SPOT配准后TMDRGSPOT与TM影像的配准已有高精度纠正结果时：在对已有纠正图像进一步检查的基础上，以纠正好的、几何精度高的数据为底图，进行配准纠正。监测范围较小时：高几何分辨率全色SPOT先几何校正到地

图坐标系，再将TM配准到全色SPOT上。

监测范围较大时：先将高质量的某一时相TM数据进行几何纠

正，象元采样成10米的地面分辨率，作为标准影象地图；再

将其余卫星数据以景为单元与其进行相对配准。快速融合结果局部放大对配准后数据进行快速融合处理，目视检查融合图是否有重影现象。配准精度检查配准不合格，有重影的影像配准精度检查透明叠加检查利用图像处理软件的连接和分层交互显示功能目视检查数据配准精度。配准精度检查影像镶嵌对若干幅互为邻接的影像通过几何镶嵌、色调调整、去重叠等处理，拼接为统一的数字影像=>++++++++++++++++++++++++影像镶嵌的目标空间位置镶嵌－空间位置的一致色调调整－辐射信息的一致几何位置发生错位的影像几何错位改正后的影像空间位置拼接没有经过色彩调节的拼接影像色彩拼接→经过色彩调节的拼接影像色彩拼接68多源信息融合-影像复合

将不同类型传感器获得的同一地区的数据进行空间配准后，将各数据中的优势或互补性有机结合起来产生新数据的技术过程。目的是保留光谱信息和提高几何特征。TMSPOT融合结果69多源信息复合1、遥感信息复合1）不同传感器的遥感信息复合（以TM与SPOT图像为例）步骤（1）配准：先完成配准,使两幅图像所对应的地物吻合，分辨率一致。（2）复合：进行彩色合成，因该方法的效果比较明显。方法一：TM432每个波