遥感图像几何纠正 下载本文

内容发布更新时间 : 2024/4/28 8:04:45星期一 下面是文章的全部内容请认真阅读。

第三章 遥感图像的几何纠正

教学目标

1、使学生了解引起遥感图像几何畸变的原因及其进行遥感图像几何纠正的必要性;

2、使学生掌握理解进行遥感几何纠正的原理与方法;

3、要求学生通过本章的学习可以熟练使用ENVI进行遥感数字图像的几何纠正;

教学内容:

1、几何纠正的概念 2、引起几何畸变的原因

3、遥感图像几何纠正的原理与方法

4、在ENVI中进行影像到影像的配准实践

一、遥感图像几何纠正的概念:

由于搭载传感器的平台(如飞机,卫星)的姿态,速度等的不稳定,以及地球曲率,空气折射等的影响,形成的图像常有畸变,几何纠正即消除遥感图像中所包含的几何畸变的过程。通常有两个叫法:

1、图像配准(Registration):同一区域里一幅图像(基准图像)对另一幅图像的校准,以使两幅图像中的同名像元配准。

2、图像校正(Rectification):借助一组地面控制点,对一幅图像进行地理坐标的校正。

在应用之前,进行遥感图像几何纠正是很必要的。第一,对遥感原始图像进行几何变形纠正后,才能对图像信息进行各种分析,制作满足两侧和定位要求的各类地球资源及环境的遥感专题图;第二,当应用不同遥感方式、不同光谱范围以及不同成像时间的各种同一区域复合图像数据来进行计算机自动分类、地物特征的变化监测或其他应用处理时,必须进行图像间的几何配准,保证各不同图像间的几何一致性;第三,利用遥感图像进行地形图测图或更新对遥感图像的几何纠正提出了更严格的要求。

几何纠正的类型:

1)把畸变图像往地形图上配 2)把畸变图像往底图上配 3)把畸变图图像相互之间配

二、引起遥感图像几何畸变的原因

1、遥感平台位置和运动状态变化引起的畸变 无论是飞机还是卫星,运动过程中都会由于种种原因产生飞行姿势的变

化(如:航高、航速、仰俯、翻滚、偏航等)从而引起图像变形; 2、地形起伏引起的几何畸变

当地形存在起伏时,会产生局部像点的位移,使原本应是地面点的信号

被同一位置上某一高点的信号所代替。由于高差的原因,实际像点距像幅中

心的距离相对于理想像点距像幅中心的距离移动了一点。 3、地球表面曲率引起的几何畸变

地球是椭球体,因此其表面是曲面,这一曲面的影响主要体现在两个方

面,一是像点位置的移动,二是像点对应于地面宽度不等。当扫描角较大时,影响尤为突出,造成边缘景物在图像显示时被压缩。 4、大气折射引起的几何畸变

大气对电磁辐射的传播产生折射。由于大气的密度分布从下向上越来越

小,折射率不断变化,因此折射后的辐射传播不再是直线而是一条曲线,从而导致传感器接收的像点发生位移。 5、地球自转引起的几何畸变

卫星前进过程中,传感器对地面扫描获取影像时,地球自转影响较大,

会产生影像偏离。因为多数卫星在轨道运行的降段接收图像,即卫星自北向南运动,这时地球自西向东自转。相对运动的结果,使卫星的星下位置逐渐产生偏离。

三、遥感图像几何纠正的原理与方法

通常,遥感图像接收部门在收到图像时,会根据卫星轨道方程将卫星的位置、姿态、轨道、及扫描特征作为时间的函数加以计算,以获得传感器的位置和姿态,并在此基础上进行几何粗纠正。当用户拿到图像时,进行的几何纠正是精纠正,即根据已知地面控制点坐标和图像中的对应像元坐标建立含有畸变影像和要纠正的地理参考之间的关系,并利用该关系式对整幅影像进行几何纠正。

遥感图像几何纠正的过程:

准备工作输入原始数字图像建立纠正变换函数确定输出影像范围像元几何位置变换像元的灰度重采样输出纠正数字影像

如上图示,在精纠正的时候最重要的两个过程是:像元几何位置变换和灰度重采样。

一)像元几何位置变换

找到一种数学关系,建立变换前图像坐标(x,y)与变换后图像坐标(u,v)的关系,通过每一个变换后图像像元的中心位置计算出变换前对应的像元坐标点,通常,整数uv的像元点在原图像坐标系中一般不在整数点上,即不在原图像像元的中心。

得到变换关系后,对图像逐行逐点进行计算在原始图像中的位置。

建立几何变换关系的方法:多项式法。通常,图像的几何变形可以看作是挤

压、扭曲、缩放、偏移以及更高次基本变形的综合作用的结果。这些变形均可用多项式模型来概括。

多项式纠正模型: 1)一次多项式 Xw=a0+a1*Xb+a2*Yb Yw=b0+b1*Xb+b2*Yb 2)二次多项式

Xw=a0+a1*Xb+a2*Yb+a3*Xb^2+ a3*Yb^2+a5*Xb*Yb Yw=b0+b1*Xb+b2*Yb+b3*Xb^2+ b3*Yb^2+b5*Xb*Yb 3)三次多项式

Xw=a0+a1*Xb+a2*Yb+a3*Xb^2+

a3*Yb^2+a5*Xb*Yb+a6*Xb^3+a7*Yb^3+a8*Xb^2*Yb+ a9*Xb*Yb^2

Yw=b0+b1*Xb+b2*Yb+b3*Xb^2+

b3*Yb^2+b5*Xb*Yb+b6*Xb^3+b7*Yb^3+b8*Xb^2*Yb+ b9*Xb*Yb^2

求多项式系数:

通过寻找控制点(GCP)去求系数,有几个未知数就至少要几对控制点!一次至少要3对,二次至少要6对,三次至少要10对。有个经验公式: GCPmin=(t+1)(t+2)/2

二)灰度重采样:

由纠正后像元坐标求取的原图像坐标系中的坐标一般不在整数点上,即不在原图像像元的中心,因此需要根据其周围像元的亮度值求取纠正后像元的亮度值。方法有以下几种:

1)最近邻法:将最近像元的亮度值直接赋给计算后的像元。

(k,l(x,y)

2)双线性内插法: