遥感 数字影像 处理

1. 遥感数字影像处理 薛云 湖南城市学院城市建设系 2008.9.1

2. 遥感技术的目的-- 是为了获得地物的几何属性和物理属性 原始的遥感图像并不能地提供实现这个目的所需的准确而完备的条件。为了实现这个目的，原始遥感影像需要经过图像处理，来消除成像过程中的误差，改善图像质量。从这讲开始介绍图像处理。 遥感图像处理包括以下几个阶段：图像的校正（预处理），图像的变换，图像的增强，图像的分类。所采用的手段有：光学图像处理和数字图像处理两种方法。 -1-

3. 第二讲 影像校正 1 数字影像的性质和特点 2 影像校正 -1-

4. 1 数字影像的性质与特点 1.1模拟影像与数字影像 1.2 数字影像的特点 1.3 多波段数字影像的数据格式 -1-

5. 1 数字影像的性质与特点 1.1模拟影像与数字影像 • 模拟影像：普通像片那样的灰度级及颜色连续变化的影像 • 数字影像：把模拟影像分割成同样形状的小单元，以各个小单元的平均亮度值或中心部分的亮度值作为该单元的亮度值进行数字化的影像。 把前一部分的空间离散化处理叫采样（sampling）， 而后一部分的亮度值的离散化处理叫量化（quantization）， 以上两种过程结合起来叫影像的数字化（digitization）。 -1-

6. 采样 采样的原理 -1-

7. 量化 量化的概念 -1-

8. 1.2 数字影像的特点 • 便于计算机处理与分析：计算机是以二进制方式处理各种数据的。采用数字形式表示遥感影像，便于计算机处理。因此，与光学影像处理方式相比，遥感数字影像是一种适于计算机处理的影像表示方法。 • 影像信息损失低：由于遥感数字影像是用二进制表示的，因此在获取、传输和分发过程中，不会因长期存储而损失信息，也不会因多次传输和复制而产生影像失真。而模拟方法表现的遥感影像会因多次复制而使影像质量下降。 • 抽象性强：尽管不同类别的遥感数字影像，有不同的视觉效果，对应不同的物理背景，但由于它们都采用数字形式表示，便于建立分析模型，进行计算机解译和运用遥感影像专家系统。 -1-

9. 1.3 多波段数字影像的数据格式 (1) BSQ格式（band sequential） 各波段的二维影像数据按波段顺序排列。 -1-

10. （2）BIL方式（band interleaved by line） 对每一行中代表一个波段的光谱值进行排列，然后按波段顺序排列该行，最后对各行进行重复。 -1-

11. （3）BIP方式（band interleaved by pixel） 在一行中，每个像元按光谱波段次序进行排列，然后对该行的全部像元进行这种波段次序排列，最后对各行进行重复。 -1-

12. 2 影像校正 从具有畸变的影像中消除畸变的处理过程叫影像校正，包括：由于遥感成像过程中多种因素影响，致使遥感图像质量的衰减。遥感图像数据的校正处理就是消除遥感图像因辐射度失真、大气消光和几何畸变等造成的图像质量的衰减。遥感图像质量衰减产生的原因和作用结果都不相同，因此一般采用不同的校正处理方法。 包括： 辐射校正 大气校正 几何校正 -1-

13. 2 影像校正 2.1 辐射校正 2.2 大气校正 2.3 几何校正 -1-

14. 2.1 辐射校正（radiometric correction） • 辐射校正：针对遥感影像辐射失真或辐射畸变进行的影像校正。由于这种校正是通过纠正辐射亮度的办法来实现的，因此称作辐射校正。 2.1.1 引起辐射畸变的因素： (1)遥感器的灵敏度特性 (2)太阳高度及地形 (3)大气吸收、散射等。 -1-

15. 2.1.2 辐射校正方法 （1）由遥感器的灵敏度特性引起的畸变校正 • 由光学系统的特性引起的畸变校正：在使用透镜的光学系统中，例如在摄像面中，存在着边缘部分比中心部分发暗的现象（边缘减光）。如果以光轴到摄象面边缘的视场角为θ，则理想的光学系统中某点的光量与cosnθ几乎成正比，利用这一性质可以进行校正（cosnθ校正）。 • 由光电变换系统的特性引起的畸变校正：由于光电变换系统的灵敏度特性通常有很高的重复性，所以可以定期地在地面测定其特性，根据测定值进行校正。 -1-

16. （2）太阳高度及地形等引起的畸变校正 • 视场角和太阳角的关系引起的亮度变化的校正：太阳光在地表反射、扩散时，其边缘更亮的现象叫太阳光点（sun spot），太阳高度高时容易产生。太阳光点与边缘减光等都可以用推算阴影曲面的方法进行校正。阴影曲面是指在影像的明暗范围内，由太阳光点及边缘减光引起的畸变部分。 • 地形倾斜的影响校正：当地形倾斜时，经过地表扩散、反射再入射到遥感器的太阳光的辐射亮度就会依倾斜度而变化。可以采取用地表的法线矢量和太阳光入射矢量的夹角进行校正的方法，以及对消除了光路辐射成分的影像数据采用波段间的比值进行校正的方法等。 -1-

17. 注：由遥感器引起的误差或由太阳高度引起的误差，一般在数据生产过程中由生产单位根据遥感器参数进行校正，而不需要用户进行自行处理。用户应该考虑大气影响引起的辐射畸变。注：由遥感器引起的误差或由太阳高度引起的误差，一般在数据生产过程中由生产单位根据遥感器参数进行校正，而不需要用户进行自行处理。用户应该考虑大气影响引起的辐射畸变。 -1-

18. 2.2 大气校正 进入大气的太阳辐射会发生反射、折射、吸收、散射和透射。其中对传感器接收影响较大的是吸收和散射。为消除由大气的吸收、散射等引起失真的辐射校正，称作大气校正。 -1-

19. 2.2.1 影响遥感影像辐射失真的大气因素 （1）大气的消光（吸收和散射） （2）天空光（大气散射）照射 （3）路径辐射 -1-

20. 2.2.2 大气校正方法 常用的大气校正方法有两类: (1)一类为基于理论模型的方法，该方法必须建立大气辐射传递方程，在此基础上近似地求解。 (2)另一类方法为基于经验或统计的方法，如回归分析方法。 利用大气辐射传输方程来建立大气校正模型在理论上是可行的。实现精确的大气校正，必须找到每个波段像元亮度值和地物反射率的关系。这需要知道模型中成像时刻气溶胶的密度、水汽的浓度等大气参数。在现实中，一般很难得到这些数据，需要专门的观测来准确地测量这些数据，因此其方法应用受到一定限制 -1-

21. 以可见光和近红外为例：在辐射传输过程中，到达地表面的总辐射能量主要是太阳直射辐照和天空漫射辐照之和，由于地表目标的反射是各向异性的，因此在遥感器观测方向的目标反射辐射经大气散射和吸收之后进入遥感器视场，这一部分能量中含有目标信息。在太阳入射辐射中，有一部分能量在未到达地面之前就被大气散射和吸收了，其中有一部分散射能量进入了遥感器视场，这一部分能量（通常称之为程辐射）中不含有任何目标信息。另外，由于周围环境的存在，入射到环境表面的辐射被其反射后有一部分经过大气散射后而进入遥感器视场，另一部分又被大气反射到目标表面，再被目标表面反射和大气透过进入遥感器视场。电磁波辐射与地球大气作用理论如此复杂，因此，从遥感的应用角度出发，有必要对其进行简化，以便于理解掌握。若忽略了大气折射、湍流和偏振，并假设天空是均朗伯散射，而且假定地表面假设为均质平坦的朗伯体（而该假设条件即各向同性天空辐照度和地面朗伯面反射在遥感研究中也广泛应用）。以可见光和近红外为例：在辐射传输过程中，到达地表面的总辐射能量主要是太阳直射辐照和天空漫射辐照之和，由于地表目标的反射是各向异性的，因此在遥感器观测方向的目标反射辐射经大气散射和吸收之后进入遥感器视场，这一部分能量中含有目标信息。在太阳入射辐射中，有一部分能量在未到达地面之前就被大气散射和吸收了，其中有一部分散射能量进入了遥感器视场，这一部分能量（通常称之为程辐射）中不含有任何目标信息。另外，由于周围环境的存在，入射到环境表面的辐射被其反射后有一部分经过大气散射后而进入遥感器视场，另一部分又被大气反射到目标表面，再被目标表面反射和大气透过进入遥感器视场。电磁波辐射与地球大气作用理论如此复杂，因此，从遥感的应用角度出发，有必要对其进行简化，以便于理解掌握。若忽略了大气折射、湍流和偏振，并假设天空是均朗伯散射，而且假定地表面假设为均质平坦的朗伯体（而该假设条件即各向同性天空辐照度和地面朗伯面反射在遥感研究中也广泛应用）。 -1-

22. 严格地说，去除大气影响是将公式 中的附加项和附加因子求出，最终求出地物反射率R，从而恢复遥感影像中地面目标的真实面目。当大气透过率变化不大时，有时只要去掉含ED和Lp的数据项就可修正影像的亮度，使影像中像元之间的亮度变化真正反映不同像元地物反射率之间的变化关系。这种对大气影响的纠正是通过纠正辐射亮度的办法实现的，因此也称作辐射校正。 -1-

23. 大气影响辐射纠正 精确的校正公式需要找出每个波段像元亮度值与地物反射率的关系。为此需得到卫星飞行时的大气参数，以求出透过率Tθ、Tφ等因子。如果不通过特别的观测，一般很难得到这些数据，所以，常常采用一些简化的处理方法，只去掉主要的大气影响，使影像质量满足基本要求。 -1-

24. 粗略校正指通过比较简便的方法去掉式 中的Lp，即程辐射度，从而改善影像质量。式中还有漫入射因子ED及其他如透过率等影响，这些因子都作为地物反射率的因子出现，直接相减不易去除，常用比值法或其他校正方法去除。严格地说，程辐射度的大小与像元位置有关，随大气条件、太阳照射方向和时间变化而变化，但因其变化量微小而忽略。可以认为，程辐射度在同一幅影像的有限面积内是一个常数，其值的大小只与波段有关。 -1-

25. 大气影响的回归分析法纠正 假定某红外波段，存在程辐射为主的大气影响，且亮度增值最小，接近于零，设为波段a。现需要找到其他波段相应的最小值，这个值一定比a波段的最小值大一些，设为波段b，分别以a，b波段的像元亮度值为坐标，作二维光谱空间，两个波段中对应像元在坐标系内用一个点表示。由于波段之间的相关性，通过回归分析在众多点中一定能找到一条直线与波段b的亮度Lb轴相交，且 -1-

26. 回归分析法 -1-

27. 回归分析法 是斜率： 和 分别为a、b波段亮度的平均值。 是波段a中的亮度为0处波段b中所具有的亮度。可以认为就是波段b的程辐射度。校正的方法是将波段b中每个像元的亮度值减去 ，来改善影像，去掉程辐射。 同理依次完成其他波段的校正。 -1-

28. 2.3 几何校正 2.2.1 影像几何畸变的因素 遥感影像成像过程中所造成的各种几何畸变称为几何校正。影响影像几何畸变的因素主要包括以下三类： （1）遥感器的内部畸变：由遥感器结构引起的畸变，如遥感器扫描运动中的非直线性等。 （2）遥感平台的运行状态：包括由于平台的高度变化、速度变化、轨道偏移及姿态变化引起的图像畸变。 （3）地球本身对遥感影像的影响：包括地球的自转、高程的变化、地球曲率、大气折射等引起的图像畸变。 （4）影像投影面及地图投影法的选取 -1-

29. （1）遥感器的内部畸变： 透镜的辐射方向畸变像差 透镜的切线方向畸变像差 透镜的焦距误差 透镜的光轴与投影面的非正交性 图象投影面的非平面性 探测元件排列的不整 采样速率的变化 采样时刻的偏差 扫描镜的扫描速度的变化 -1-

30. （2）遥感平台运动状态变化 航高：当平台运动过程中受到力学因素影响，产生相对于原标准航高的偏离，或者说卫星运行的轨道本身就是椭圆的。航高始终发生变化，而传感器的扫描视场角不变，从而导致影像扫描行对应的地面长度发生变化。航高越向高处偏离，影像对应的地面越宽 -1-

31. 航速：卫星的椭圆轨道本身就导致了卫星飞行速度的不均匀，其他因素也可导致遥感平台航速的变化。航速快时，扫描带超前，航速慢时，扫描带滞后，由此可导致影像在卫星前进方向上（影像上下方向）的位置错动。航速：卫星的椭圆轨道本身就导致了卫星飞行速度的不均匀，其他因素也可导致遥感平台航速的变化。航速快时，扫描带超前，航速慢时，扫描带滞后，由此可导致影像在卫星前进方向上（影像上下方向）的位置错动。 -1-

32. 俯仰：遥感平台的俯仰变化能引起影像上下方向的变化，即星下点俯时后移，仰时前移，发生行间位置错动。俯仰：遥感平台的俯仰变化能引起影像上下方向的变化，即星下点俯时后移，仰时前移，发生行间位置错动。 -1-

33. 翻滚：遥感平台姿态翻滚是指以前进方向为轴旋转了一个角度。可导致星下点在扫描线方向偏移，使整个影像的行向翻滚角引起偏离的方向错动。翻滚：遥感平台姿态翻滚是指以前进方向为轴旋转了一个角度。可导致星下点在扫描线方向偏移，使整个影像的行向翻滚角引起偏离的方向错动。 -1-

34. 偏航：指遥感平台在前进过程中，相对于原前进航向偏转了一个小角度，从而引起扫描行方向的变化，导致影像的倾斜畸变。偏航：指遥感平台在前进过程中，相对于原前进航向偏转了一个小角度，从而引起扫描行方向的变化，导致影像的倾斜畸变。 -1-

35. （3）地球本身对遥感影像的影响 地形起伏的影响 当地形存在起伏时，会产生局部像点的位移，使原来本应是地面点的信号被同一位置上某高点的信号代替。由于高差的原因，实际像点P距像幅中心的距离相对于理想像点P0距像幅中心的距离移动了△r。 高差引起的像点位移 -1-

36. 地表曲率的影响 地球是球体，严格说是椭球体，因此地球表面是曲面。这一曲面的影响主要表现在两个方面，一是像点位置的移动，当选择的地图投影平面是地球的切平面时，使地面点P0相对于投影平面点P有一高差△h。 像点位移 -1-

37. 二是像元对应于地面宽度的不等。由于传感器通过扫描取得数据，在扫描过程中每一次取样间隔是星下视场角的等分间隔。如果地面无弯曲，在地面瞬时视场宽度不大的清况下，L1，L2，L3，…的差别不大。但由于地球表面曲率的存在，对应于地面的P1，P2，P3，…，显然P3-P1 ＞L3-L1，距星下点越远畸变越大，对应地面长度越长。 像元对应于地面宽度的不等 -1-

38. 大气折射的影响 大气对辐射的传播产生折射。由于大气的密度分布从下向上越来越小，折射率不断变化，因此折射后的辐射传播不再是直线而是一条曲线，从而导致传感器接收的像点发生位移 N P 大气折射的影响 -1-

39. (a)获得影像 （b）实际对应的 地面位置 （c）影像变形 地球自转的影响 卫星前进过程中，传感器对地面扫描获得影像时，地球自转影响较大，会产生影像偏离。因为多数卫星在轨道运行的降段接收影像，即卫星自北向南运动，这时地球自西向东自转。相对运动的结果，使卫星的星下位置逐渐产生偏离。偏离方向如图所示，所以卫星影像经过校正后成为图C的形态。 地球自转引起偏离

40. 2.2.2 几何校正的方法 从具有几何畸变的影像中消除畸变的过程。也可以说是定量地确定影像上的像元坐标（影像坐标）与目标物的地理坐标（地图坐标等）的对应关系（坐标变换式）。 -1-

41. （1）基本思路 校正前的影像看起来是由行列整齐的等间距像元点组成的，但实际上，由于某种几何畸变，影像中像元点间所对应的地面距离并不相等（图a）。校正后的影像亦是由等间距的网格点组成的，且以地面为标准，符合某种投影的均匀分布（图b），影像中格网的交点可以看作是像元的中心。校正的最终目的是确定校正后影像的行列数值，然后找到新影像中每一像元的亮度值。 -1-

42. （2）几何校正的方法 系统性校正： 当知道了消除影像几何畸变的理论校正公式时，可把该式中所含的与遥感器构造有关的校准数据（焦距等）及遥感器的位置、姿态等的测量值代入到理论校正式中进行几何校正。该方法对遥感器的内部畸变大多是有效的。可是在很多情况下，遥感器的位置及姿态的测量值精度不高，所以外部畸变的校正精度也不高。 -1-

43. 非系统性校正： 利用实地测量的地物的真实坐标值，寻找实测值与畸变之后的图像之间的函数关系，从而得到几何校正的方法。利用控制点的影像坐标和地图坐标的对应关系，近似地确定所给的影像坐标系和应输出的地图坐标系之间的坐标变换式。坐标变换式经常采用1次、2次等角变换式，2次、3次投影变换式或高次多项式。坐标变换式的系数可从控制点的影像坐标值和地图坐标值中根据最小2乘法求出。 -1-

44. 复合校正： 实际工作中常常将两种方法结合起来。 把理论校正式与利用控制点确定的校正式组合起来进行校正。① 分阶段校正的方法，即首先根据理论校正式消除几何畸变（如内部畸变等），然后利用少数控制点，根据所确定的低次校正式消除残余的畸变（外部畸变等）；② 提高几何校正精度的方法，即利用控制点以较高的精度推算理论校正式中所含的遥感器参数、遥感器的位置及姿态参数。 -1-

45. （3）几何精校正的步骤 几何畸变有多种校正方法，但常用的是一种通用的精校正方法，适合于在地面平坦，不需考虑高程信息，或地面起伏较大而无高程信息，以及传感器的位置和姿态参数无法获取的情况时应用。有时根据遥感平台的各种参数已做过一次校正，但仍不能满足要求，就可以用该方法作遥感影像相对于地面坐标的配准校正，遥感影像相对于地图投影坐标系统的配准校正，以及不同类型或不同时相的遥感影像之间的几何配准和复合分析，以得到比较精确的结果。 -1-

46. 几何精校正概括为两个步骤： 第一步是构建一个模拟几何畸变的数学模型，以建立原始畸变图像空间与标准图像空间的某种对应关系，实现不同图像空间中像元位置的变换； 第二步是利用这种对应关系把原始畸变图像空间中全部像素变换到标准图像空间中的对应位置上，完成标准图像空间中每一像元亮度值的计算。 具体步骤： -1-

47. 第一步 选择几何校正计算模型 -1-

48. 第二步 选择几何校正采点模式 -1-

49. 第三步 采集地面控制点（GCP） 第四步 采集地面检查点

50. 第五步 影像重采样 -1-