遥感图像处理2PPT学习教案

《遥感图像处理2PPT学习教案》由会员分享，可在线阅读，更多相关《遥感图像处理2PPT学习教案（57页珍藏版）》请在装配图网上搜索。

1、会计学12George F. SmootJohn C. Mather第1页/共57页3for their discovery of the blackbody form and anisotropy of the cosmic microwave background radiation新浪：以表彰他们发现了黑体形态和宇宙微波背景辐射的扰动现象；163：他们的研究成果是发现了黑体结构以及宇宙背景辐射的微波各向异性。自译：发现了宇宙微波背景辐射的黑体形式和各向异性。 回望宇宙的婴儿时代第2页/共57页4第四章第四章 遥感图像处理遥感图像处理1）目标物的电磁波特性）目标物的电磁波特性2）信息的获取

2、（传感器和遥感平台）信息的获取（传感器和遥感平台）3）信息的接收）信息的接收4）信息的处理）信息的处理5）信息的应用）信息的应用第3页/共57页5遥感图像处理遥感图像处理光学影像光学影像光学原理和光学处理方法光学原理和光学处理方法数字影像数字影像数字图像校正数字图像校正数字图像增强数字图像增强多源信息复合多源信息复合第4页/共57页61 数字图像数字图像2 辐射校正辐射校正3 几何校正几何校正数字图像校正数字图像校正第5页/共57页7汤国安汤国安科学出版社科学出版社第6页/共57页8模拟图像与数字图像模拟图像与数字图像光学图像（模拟图像）连续图像照片数字图像离散图像计算机存储模拟图像与数字图像

3、之间的转换，称为模/数（A/D）转换。第7页/共57页9数字图像数字图像指能够被计算机存储、处理和使用的图像指能够被计算机存储、处理和使用的图像“离散化离散化”、二维矩阵：每个元素的取值是图像连续变化的灰度的离散整数值、二维矩阵：每个元素的取值是图像连续变化的灰度的离散整数值第8页/共57页10第9页/共57页11 模拟图像模拟图像 （灰度、颜（灰度、颜色的连续分色的连续分布布） 采采样样像元号像元号 像元像元 第10页/共57页121、离散化的单元像元值的涵义？代表什么？2、在计算机内是如何存储的（以什么数据类型）？3、图像所占存储空间如何计算？第11页/共57页13 用平面直角坐标系表示一

4、幅灰度范围为用平面直角坐标系表示一幅灰度范围为0-n的数字图像像元灰度分布状态，横轴表示灰度级，纵轴表示某一灰度级（或范围）的像元个数占像元总数的百分比。的数字图像像元灰度分布状态，横轴表示灰度级，纵轴表示某一灰度级（或范围）的像元个数占像元总数的百分比。 通过灰度直方图可以直观地了解图像特征，以确定图像增强方案并了解图像增强后的效果通过灰度直方图可以直观地了解图像特征，以确定图像增强方案并了解图像增强后的效果图像直方图图像直方图第12页/共57页14不同反差特征的图像不同反差特征的图像第13页/共57页15反差增强反差增强第14页/共57页16遥感图像处理系统遥感图像处理系统计算机（主机）计

5、算机（主机）图像输入输出设备图像输入输出设备磁带机、数字化器等；打印机、绘图仪、激光图像记录仪磁带机、数字化器等；打印机、绘图仪、激光图像记录仪专用处理设备专用处理设备图像计算机、阵列处理机图像计算机、阵列处理机外存设备外存设备磁盘、磁带、光盘磁盘、磁带、光盘显示器显示器软件部分软件部分系统软件、应用软件（图像处理软件）系统软件、应用软件（图像处理软件）第15页/共57页171 数字图像数字图像2 辐射校正辐射校正3 几何校正几何校正数字图像校正数字图像校正第16页/共57页18辐射校正辐射校正 由于由于传感器传感器响应特性和响应特性和大气大气的吸收、散的吸收、散射及其它随机因素影响，导致图像

6、模糊射及其它随机因素影响，导致图像模糊失真，造成图像分辨率和对比度相对下失真，造成图像分辨率和对比度相对下降。这些都需要通过辐射校正复原。降。这些都需要通过辐射校正复原。 包括：包括：系统辐射校正、大气校正系统辐射校正、大气校正第17页/共57页19大气影响的定量分析大气影响的定量分析0E0L地面上单位面积的辐照度为地面上单位面积的辐照度为 cos0EE假定地表面是朗伯体，其表面为漫反射，则某方向物体的亮度为假定地表面是朗伯体，其表面为漫反射，则某方向物体的亮度为 cos00ERERL系统增益系数因子系统增益系数因子 ，进入传感器的亮度值，进入传感器的亮度值 Scos00SERL无大气影响状况

7、无大气影响状况第18页/共57页20大气影响的定量分析大气影响的定量分析cos01STETRLDESTRL2PDSLSETRSTETRLcos0受入射方向透过率受入射方向透过率 和反射方向透过率和反射方向透过率 的影响，进入传感器的亮度值为：的影响，进入传感器的亮度值为：TT漫入射辐照度为ED，反射后进入传感器的亮度值为：散射光直接进入传感器的辐射（程辐射亮度LP）第19页/共57页21大气影响的定量分析大气影响的定量分析PDSLSETRSTETRLcos0大气的影响主要是减少了图像的对比度，使原始信号和背景信号都增加了因子。大气的影响主要是减少了图像的对比度，使原始信号和背景信号都增加了因子

8、。第20页/共57页22大气影响的粗略校正大气影响的粗略校正 通过比较简便的方法去掉上式中的通过比较简便的方法去掉上式中的LP，从而改，从而改善图像质量。善图像质量。 PDSLSETRSTETRLcos0问题：精确校正？问题：精确校正？第21页/共57页23大气影响的粗略校正大气影响的粗略校正（1 1）直方图最小值去除法）直方图最小值去除法（2 2）回归分析法）回归分析法 第22页/共57页24（1）直方图校正方法）直方图校正方法 从图像像元亮度值中减去一个辐射偏置量（从图像像元亮度值中减去一个辐射偏置量（LP），辐射偏置量等于图像直方图中最小的辐射亮度，辐射偏置量等于图像直方图中最小的辐射亮

9、度值。值。 前提（假设）：水体（或阴影）等物体的灰度值前提（假设）：水体（或阴影）等物体的灰度值为为0，大气散射导致图像上这些物体的灰度值不，大气散射导致图像上这些物体的灰度值不为为0（辐射偏置量）（辐射偏置量） 暗物体法（暗物体法（Dark-object method）第23页/共57页25第24页/共57页26（2 2）回归分析法）回归分析法原理原理：在遥感图像上大山的阴影区或深大水体区域，各个波段的反射为零。同时，大气散射主要影响短波部分，波长较长的波段几乎不受影响，因此可用其校正其它波段数据。：在遥感图像上大山的阴影区或深大水体区域，各个波段的反射为零。同时，大气散射主要影响短波部分，

10、波长较长的波段几乎不受影响，因此可用其校正其它波段数据。方法方法：在不受大气影响的波段（如：在不受大气影响的波段（如TM5或或7）和待校正的某一波段图像中，选择一系列目标，将每个目标的两个待比较的波段灰度值提取出来进行回归分析，建立线性回归方程。）和待校正的某一波段图像中，选择一系列目标，将每个目标的两个待比较的波段灰度值提取出来进行回归分析，建立线性回归方程。XbaYY：待校正波段的图像亮度值：待校正波段的图像亮度值X：不受大气影响波段的图像亮度值：不受大气影响波段的图像亮度值第25页/共57页27（2 2）回归分析法）回归分析法7TMbiaiTMi怎么求怎么求a,b？第26页/共57页28

11、1 数字图像数字图像2 辐射校正辐射校正3 几何校正几何校正数字图像校正数字图像校正第27页/共57页29遥感影像变形的原因遥感影像变形的原因（1 1）遥感平台位置和运动状态变化）遥感平台位置和运动状态变化第28页/共57页30遥感影像变形的原因遥感影像变形的原因（2 2）地形起伏影响）地形起伏影响第29页/共57页31遥感影像变形的原因遥感影像变形的原因（3 3）地球表面曲率的影响）地球表面曲率的影响 像点位置的移动像元对应于地面宽度的不等第30页/共57页32遥感影像变形的原因遥感影像变形的原因（3 3）地球表面曲率的影响）地球表面曲率的影响 第31页/共57页33遥感影像变形的原因遥感影

12、像变形的原因（4 4）大气折射的影响）大气折射的影响 第32页/共57页34遥感影像变形的原因遥感影像变形的原因（5 5）地球自转的影响）地球自转的影响 第33页/共57页35几何畸变校正几何畸变校正第34页/共57页36几何畸变校正几何畸变校正 数字纠正数字纠正：通过计算机对图象每个象元逐个：通过计算机对图象每个象元逐个地解析纠正处理完成，可以较精确改正线性地解析纠正处理完成，可以较精确改正线性和非线性变形误差。和非线性变形误差。 包括两个方面（基本环节）：包括两个方面（基本环节）： （1 1）像）像元坐标变换；元坐标变换； （2 2）像元灰度值重新计算（重采样）像元灰度值重新计算（重采样）

13、。第35页/共57页37准备工作输入原始数字影像建立纠正变换函数建立纠正变换函数确定输出影像范围确定输出影像范围像元坐标变换像元坐标变换像元亮度值重采样像元亮度值重采样输出纠正后的图像准备工作准备工作：图像、地图、大地测量资料、平台轨道参数、传感器参数、控制点的选择；（具体内容可选）纠正变换函数建立纠正变换函数建立：输入和输出图像间的坐标变换关系；如多项式法、共线方程法等第36页/共57页38数字图象几何校正方法数字图象几何校正方法(函数的建立函数的建立)多项式纠正法多项式纠正法：常用的方法：常用的方法基本原理基本原理：不考虑成像的空间几何过程，而直接对图象变形的本身进行数学模拟。把遥感图象的

14、总体变形看作是平移、缩放、旋转、仿射、偏扭、弯曲及更高次的基本变形的综合作用结果，因此，将纠正前后图象相应点间的坐标关系用一适当的多项式表达不考虑成像的空间几何过程，而直接对图象变形的本身进行数学模拟。把遥感图象的总体变形看作是平移、缩放、旋转、仿射、偏扭、弯曲及更高次的基本变形的综合作用结果，因此，将纠正前后图象相应点间的坐标关系用一适当的多项式表达一般取一般取2次，可以满足精度要求次，可以满足精度要求需要地面控制点需要地面控制点： 分布、数量分布、数量第37页/共57页39第38页/共57页40n 先求出原始图像四个角点（a, b, c, d）在纠正后图像中的对应点（a, b, c, d）

15、的坐标（Xa,Ya）(Xb,Yb) (Xc,Yc) (Xd,Yd)；然后求出最大值和最小值。n 为了将该边界范围转换为计算机纠正后的存储数组空间，须在其中划分出网格，每个网格代表一个输出像元；须根据精度要求确定输出像元的地面尺寸；确定输出像元的地面尺寸；n 地面坐标与输出图象坐标关系表达第39页/共57页41lX1 = min (Xa, Xb, Xc, Xd)lX2 = max (Xa, Xb, Xc, Xd)lY1 = min (Ya, Yb, Yc, YXd)lY2 = max (Ya, Yb,Yc, Yd)第40页/共57页42坐标变换坐标变换确定原始图象和纠正后图象间的坐标变换关系确定

16、原始图象和纠正后图象间的坐标变换关系（x, y）- (u, v) (u行数行数, v列数列数, 均为整数均为整数)有两种方案：有两种方案：直接纠正方案和间接纠正方案直接纠正方案和间接纠正方案直接纠正方案直接纠正方案：从原始图象阵列出发，依次对：从原始图象阵列出发，依次对 其中每一个象元其中每一个象元P(x，y)分别计算其在输出（纠正后）图象的坐标分别计算其在输出（纠正后）图象的坐标P(X，Y)，并计算，并计算P(X，Y)的灰度值的灰度值;间接纠正方案间接纠正方案：从空白图象阵列出发，依次计算每个象元：从空白图象阵列出发，依次计算每个象元P(X，Y)在原始图象中的位置在原始图象中的位置P(x，y

17、)，然后把该点的灰度值依次计算后返送给，然后把该点的灰度值依次计算后返送给P(X，Y). 二者间并无本质差别，二者间并无本质差别，互为逆变换互为逆变换第41页/共57页43),(ppXpyxFX ),(ppYpyxFY ),(PPxpYXfx ),(PPypYXfy （xp ,yp）（）（XP,YP）分别是任意一个像元在原始图像和纠正后图像中的坐标。）分别是任意一个像元在原始图像和纠正后图像中的坐标。直接间接间接第42页/共57页44直接纠正方法直接纠正方法:从原始图像，依次对每个像元根据变换函数 F（），求得它在新图像中的位置。并将灰度值付给新图像的对应位置上。间接纠正法间接纠正法: :从新

18、图像中依次每个像元，根据变换函数 f () 找到它在原始图像中的位置，并将图像的灰度值赋予新图像的像元。二者并无本质区别，互为逆变换第43页/共57页45第44页/共57页46第45页/共57页47第46页/共57页48第47页/共57页49原始图像原始图像纠正后图像（最邻近插纠正后图像（最邻近插值）值）第48页/共57页50双线性内插法双线性内插法利用利用X方向和方向和Y方向进行三次插值方向进行三次插值(4个邻点个邻点)第49页/共57页51双线性插值双线性插值原始图像纠正（双线性插值）第50页/共57页52第51页/共57页53第52页/共57页54原始图像原始图像几何纠正（三次卷积）几何

19、纠正（三次卷积）第53页/共57页55地面控制点地面控制点（GCP，Ground Control Point）: 一些特定的象元，其地图坐标或其它输出坐标为已知一些特定的象元，其地图坐标或其它输出坐标为已知人工地物人工地物 线性地物交叉点线性地物交叉点不易随时间变化的目标不易随时间变化的目标大比例尺的图象大比例尺的图象：道路交叉点、机场跑道、建筑物：道路交叉点、机场跑道、建筑物小比例尺的图象小比例尺的图象：城区、一些线性地物交叉点（河流、道路）：城区、一些线性地物交叉点（河流、道路）分布分布：较均匀分布与图象范围内，保证足够数量：较均匀分布与图象范围内，保证足够数量第54页/共57页56控制点

20、的选择控制点的选择注意问题注意问题：1. 多项式纠正法的精度与地面控制点（多项式纠正法的精度与地面控制点（GCP）的精度、分布、数量及纠正范围有关；）的精度、分布、数量及纠正范围有关；GCP的位置精度越高，则几何纠正的精度越高；的位置精度越高，则几何纠正的精度越高；GCP的个数不少于多项式的系数个数；适当增加的个数不少于多项式的系数个数；适当增加GCP的个数，可以提高几何纠正的精度。的个数，可以提高几何纠正的精度。20-30个个GCP，一般可以满足需求，一般可以满足需求2 GCP分布应尽可能在整幅图像内均匀分布，否则在分布应尽可能在整幅图像内均匀分布，否则在GCP密集区精度较高，在密集区精度较高，在GCP分布稀疏区出现较大误差分布稀疏区出现较大误差第55页/共57页57最少的控制点数量最少的控制点数量第56页/共57页