第四章遥感数字图像处理教材

1、14.1.1 4.1.1 颜色视觉颜色视觉光和颜色光和颜色 电磁波谱中电磁波谱中0.38-0.760.38-0.76 m m的波段能引起人的视觉，的波段能引起人的视觉，称为称为可见光可见光, , 简称光简称光 可见光加上紫外和红外部分来自原子和分子的发可见光加上紫外和红外部分来自原子和分子的发光辐射，称为光学辐射光辐射，称为光学辐射 -光学遥感光学遥感 光与色有什么不同？光与色有什么不同？ 光是色的源泉，色是光的表现。光是色的源泉，色是光的表现。“五光十色五光十色”24.1.1 颜色视觉颜色视觉(一一)亮度对比亮度对比 亮度对比：视场中对象与背景的亮度差与背景亮亮度对比：视场中对象与背景的亮度

2、差与背景亮度之比度之比 C=（L对象对象-L背景）背景）/L背景背景 选择适宜的对象及背景的亮度，可以提高对比，选择适宜的对象及背景的亮度，可以提高对比，从而提高视觉效果。从而提高视觉效果。 两个或多个对象间的对比，即两个或多个对象间的对比，即C= L对象对象/L对象对象3亮度对比例子亮度对比例子高度对比度合适亮度对比度过量44.1.1 颜色视觉（一）颜色视觉（一）颜色对比颜色对比 颜色对比：在视场中，相邻区域的不同颜色相互影响。 受视觉影响很大 (所有颜色都是对某段波长有选择地反射而对其它波长吸收的结果) 两种颜色互相影响的结果，使每种颜色向其影响色的补色变化。两种颜色的边界，对比现象更加明

3、显 在色度学中，当两种颜色混合产生白色或灰色时，这两种颜色称为互补色。 绿是品红的补色，蓝是黄的补色 人眼能分辨一百多种不同颜色，彩色图象能表现更丰富的信息量564.1.1 颜色视觉颜色视觉(二二)颜色的性质颜色的性质 颜色描述对遥感图象很重要，颜色变换是遥感图颜色描述对遥感图象很重要，颜色变换是遥感图象处理的重要方法。象处理的重要方法。 所有颜色都是对某段波长有选择地反射而对其它所有颜色都是对某段波长有选择地反射而对其它波长吸收的结果波长吸收的结果 颜色的性质由明度、色调、饱和度来描述颜色的性质由明度、色调、饱和度来描述 明度（Lightness, 强度, Intensity） 色调（Hue

4、） 饱和度(Saturation)74.1.1颜色视觉颜色视觉(二二)颜色的性质：颜色的性质：明明度度(Lightness) 明度：是人眼对光源或物体明亮程度的感觉。与明度：是人眼对光源或物体明亮程度的感觉。与电磁波辐射亮度的概念不同，明度受人的视觉感电磁波辐射亮度的概念不同，明度受人的视觉感受性和经验影响。受性和经验影响。 物体反射率越高，明度就越高。白比灰，黄比红物体反射率越高，明度就越高。白比灰，黄比红 光源：亮度越大，明度越高光源：亮度越大，明度越高 黑白图象：灰度、灰阶黑白图象：灰度、灰阶8灰度（明度）、灰阶灰度（明度）、灰阶94.1.1 颜色视觉颜色视觉(二二)颜色的性质：色颜色的

5、性质：色调调(Hue)色彩彼此间相互区分的特性。可见光波段的不同波色彩彼此间相互区分的特性。可见光波段的不同波长刺激人眼产生了不同色彩（红长刺激人眼产生了不同色彩（红紫）的感觉。紫）的感觉。 从红到紫是可见光谱上存在的颜色，每种颜色对从红到紫是可见光谱上存在的颜色，每种颜色对应一个波长值，是应一个波长值，是光谱色光谱色。 有时刺激人眼的光波不是单一波长，而是有时刺激人眼的光波不是单一波长，而是一些波一些波长的组合长的组合，也可构成一些颜色，但它们找不到对，也可构成一些颜色，但它们找不到对应的波长值，不叫光谱色。应的波长值，不叫光谱色。104.1.1 颜色视觉颜色视觉(三三)颜色的性质：饱颜色的

6、性质：饱和度和度(Saturation) 彩色纯洁的程度，即光谱中波长段是否窄，频率是否单一彩色纯洁的程度，即光谱中波长段是否窄，频率是否单一的表示。对于光源，发出的若是单色光就是最饱和的彩的表示。对于光源，发出的若是单色光就是最饱和的彩色，如激光。对于物体颜色，如果物体对光谱反射性有很色，如激光。对于物体颜色，如果物体对光谱反射性有很高的选择性，只反射很窄的波段则饱和度高。高的选择性，只反射很窄的波段则饱和度高。 如果光源或物体反射光在某种波长中混有许多其它波长的如果光源或物体反射光在某种波长中混有许多其它波长的光或混有白光则饱和度变低。光或混有白光则饱和度变低。 黑白色只用明度描述，不用色

7、调、饱和度描述。黑白色只用明度描述，不用色调、饱和度描述。11彩色彩色 IH S124.1.1 颜色视觉（三）颜色视觉（三）颜色立体（球颜色立体（球形）形） 球形明度球形明度-色度色度-饱和度色彩坐标系饱和度色彩坐标系 垂直坐标轴代表亮度垂直坐标轴代表亮度 圆周代表色度，色度值从红色中点处的圆周代表色度，色度值从红色中点处的0开始，开始，沿着球体圆周的逆时针方向增加沿着球体圆周的逆时针方向增加 饱和度描述了色彩的纯度饱和度描述了色彩的纯度134.1.1 颜色视觉（三）颜色视觉（三）颜色立体（锥颜色立体（锥形）形）中间垂直轴代表明度，中间水平面的圆周代表色调、圆周上的半径大中间垂直轴代表明度，中

8、间水平面的圆周代表色调、圆周上的半径大小代表饱和度小代表饱和度144.1.1 颜色视觉（三）颜色视觉（三）颜色立体颜色立体(孟孟赛尔颜色立体赛尔颜色立体)154.1.2 加色法与减色法加色法与减色法(一一)：颜色相加原：颜色相加原理理 三原色：若三种颜色，其中的任一种都不能由其三原色：若三种颜色，其中的任一种都不能由其余二种颜色混合相加产生，三种颜色按一定比例余二种颜色混合相加产生，三种颜色按一定比例混合，可以形成各种色调的颜色，则称之为三原混合，可以形成各种色调的颜色，则称之为三原色。实验证明，红、绿、蓝三种颜色是最优的三色。实验证明，红、绿、蓝三种颜色是最优的三原色。原色。 由三原色混合，

9、可以产生其他颜色，称为加色由三原色混合，可以产生其他颜色，称为加色法。互补色：若两种颜色混合产生白色或灰色法。互补色：若两种颜色混合产生白色或灰色 , 这两种颜色就称为互补色。这两种颜色就称为互补色。16174.1.2 加色法与减色法加色法与减色法(二二)：色度图：色度图颜色相加的原理可以用色度图来表现，比加色法示意图更接近实际情颜色相加的原理可以用色度图来表现，比加色法示意图更接近实际情况。因为每种波长的光都可以用红绿蓝三原色相加产生。况。因为每种波长的光都可以用红绿蓝三原色相加产生。对任何一种颜色的光，当匹配的各波长光谱能量相同（等能光谱）对任何一种颜色的光，当匹配的各波长光谱能量相同（等

10、能光谱）时，都可以推算出其所需要的红绿蓝三原色的数量值。时，都可以推算出其所需要的红绿蓝三原色的数量值。所有光谱色混合时，即形成等能光谱中的白光，且白光时由相同数量所有光谱色混合时，即形成等能光谱中的白光，且白光时由相同数量的红绿蓝三原色组成。的红绿蓝三原色组成。设光的总能量为设光的总能量为1，则白光由三原色各，则白光由三原色各1/3构成，即红构成，即红=绿绿=蓝蓝=1/3白，白，红红+绿绿+蓝蓝=白白根据上述原理，设计了色度图根据上述原理，设计了色度图18 左图代表色度图上的色左图代表色度图上的色彩分布彩分布 其中，其中，x: 红色的比例；红色的比例； y: 绿色的比例；蓝色比绿色的比例；蓝

11、色比例可由例可由x+y+z=1导出导出 弧形曲线代表光谱，线弧形曲线代表光谱，线上每一点代表一种波长上每一点代表一种波长（nm）和光谱颜色，）和光谱颜色，曲线曲线包围的部分及直包围的部分及直线部分代表非光谱线部分代表非光谱色。色。19色度图可以粗略推算两种色度图可以粗略推算两种颜色相混合得到的中间颜色相混合得到的中间色，色，M和和N两种颜色按照两种颜色按照一定比例合成，一定能得一定比例合成，一定能得到到MN连线上的中间色连线上的中间色K（只要比例合适，（只要比例合适，MN上上其他的点也可得到）其他的点也可得到）连接连接C点与色度图曲线内点与色度图曲线内的一点，可得该点的光的一点，可得该点的光谱

12、，例如连接谱，例如连接C与与K，可，可得得K点的光谱色点的光谱色(0.573um)该点距离该点距离C点的远近反映点的远近反映了了C点的饱和度点的饱和度过过C点的直线与边缘交于点的直线与边缘交于两个点，则两个点对应的两个点，则两个点对应的颜色一定是互补色，两者颜色一定是互补色，两者混合可产生白光。混合可产生白光。204.1.2加色法与减色法加色法与减色法(三三)：减色法：减色法 减色法是从自然光（白光）中，减去一种或二种基色光而减色法是从自然光（白光）中，减去一种或二种基色光而生成色彩的方法。生成色彩的方法。 颜料本身的色彩是由于本身选择性地吸收了入射自然光中颜料本身的色彩是由于本身选择性地吸收

13、了入射自然光中一定波长的光，反射出未被吸收的色光而呈现出本身的色一定波长的光，反射出未被吸收的色光而呈现出本身的色彩。彩。 例如，黄色颜料是由于本身吸收了自然光中的蓝色光，反例如，黄色颜料是由于本身吸收了自然光中的蓝色光，反射出未被吸收的红光和绿光叠加混合的结果；品红颜料是射出未被吸收的红光和绿光叠加混合的结果；品红颜料是由于吸收了自然光中的绿光，反射出红光和蓝光相加的结由于吸收了自然光中的绿光，反射出红光和蓝光相加的结果；同样，青颜料是由于吸收了自然光中的红光成分，反果；同样，青颜料是由于吸收了自然光中的红光成分，反射蓝光与绿光的结果。射蓝光与绿光的结果。214.1.2加色法与减色法加色法与

14、减色法(三三)：减法三原色 减法三原色：指加法三原色的补色，即黄、品红和青色。减法三原色：指加法三原色的补色，即黄、品红和青色。用白光由红、绿、蓝三色组成这种理想模型来理解，可以用白光由红、绿、蓝三色组成这种理想模型来理解，可以认为黄色，是减去蓝色的的红绿组合；同样地，品红色是认为黄色，是减去蓝色的的红绿组合；同样地，品红色是减去绿色的红蓝组合，青色是减去红色的蓝绿组合。这减去绿色的红蓝组合，青色是减去红色的蓝绿组合。这样，黄、品红、青便是减色法的三原色。样，黄、品红、青便是减色法的三原色。 将彩色涂料的三色叠加时，由于光线依次通过减红、减将彩色涂料的三色叠加时，由于光线依次通过减红、减绿、减

15、蓝层就成黑色。只有当涂料浓度不够，减色不彻底绿、减蓝层就成黑色。只有当涂料浓度不够，减色不彻底时才会出现灰白色，但这仍是减色法而不是加色法。时才会出现灰白色，但这仍是减色法而不是加色法。22补充：遥感数字图像的表现形式补充：遥感数字图像的表现形式黑黑白、彩色白、彩色 遥感影像常以照片（或数字影像）来表现遥感影像常以照片（或数字影像）来表现: 单波段或全色波段表现为单波段或全色波段表现为黑白图像黑白图像， 三波段组合表现为三波段组合表现为彩色影像彩色影像 彩色：彩色： 真彩色真彩色(true color)：（三波段组合）：（三波段组合） 假彩色假彩色(false color)：（三波段组合）：（

16、三波段组合） 伪彩色伪彩色(pseudo color) 伪彩色 (pseudocolor): 灰度图像的彩色表示或显示23单波段单波段三波段三波段241. 数字图像基础2. 图像辐射校正3. 几何校正原理与方法4. 图像增强处理25数字图像概念：指能够被计算机存储、处理和使用的图像 。“离散化”、二维矩阵是数字图像的特点，每个像元的取值是图像连续变化的灰度的离散整数值。26数字图像数字图像 数字图像的性质与特点数字图像的性质与特点 什么是数字图像？什么是数字图像？ 模拟图像：普通像片那样的灰度级及颜色连续变化的图像 数字图像：而是以数字形式表示的遥感影像。包括把模拟图像分割成同样形状的小单元，

17、以各个小单元的平均亮度值或中心部分的亮度值作为该单元的亮度值进行数字化的图像。 把前一部分的空间离散化处理叫把前一部分的空间离散化处理叫采样采样（sampling），）， 而 后 一 部 分 的 亮 度 值 的 离 散 化 处 理 叫而 后 一 部 分 的 亮 度 值 的 离 散 化 处 理 叫 量 化量 化（quantization），）， 以 上 两 种 过 程 结 合 起 来 叫 图 像 的以 上 两 种 过 程 结 合 起 来 叫 图 像 的 数 字 化数 字 化（digitization）。）。 27 数字图像又称数字化图像，是一种以二维数组数字图像又称数字化图像，是一种以二维数组（

18、矩阵）形式表示的图像。或者称为相应区域（矩阵）形式表示的图像。或者称为相应区域内地物电磁辐射强度的二维分布。内地物电磁辐射强度的二维分布。 该数组由对连续变化的空间图像作等间距抽样该数组由对连续变化的空间图像作等间距抽样所产生的抽样点所产生的抽样点像元（像素）组成。像元（像素）组成。 28采样采样采样的原理29量化量化量化的概念30 像元的量值，通常为抽样区间内连续变化的景像元的量值，通常为抽样区间内连续变化的景物的均值化量值，一般为亮度值或灰度值，它物的均值化量值，一般为亮度值或灰度值，它们的最大、最小值区间代表该数字图像的动态们的最大、最小值区间代表该数字图像的动态范围。范围。 31 0=

19、f(x, y)=G x belong to 0, xmax, y belong to 0, ymax 式中式中, G为灰度值的上界为灰度值的上界. 因此因此, 一幅图像可用一幅图像可用M(行行)N(列列)的矩阵函数表示的矩阵函数表示: 1 2 7 12 34 56 3 4 56 12 3 4 23 23 5 6 8 9 45 8 6 5 12 6 45 34 5 34 34 23 9 9 23 34 34 5 6 45 67 5632 在遥感图像处理中在遥感图像处理中, 既需要将光学图像转换为数既需要将光学图像转换为数字图像进行计算机处理字图像进行计算机处理,也需要将处理后的数字图也需要将处理

20、后的数字图像变成光学图像输出。光学图像为模拟量，数字像变成光学图像输出。光学图像为模拟量，数字图像又称数字量，它们之间的转换称模图像又称数字量，它们之间的转换称模/数转换，数转换，记作记作A/D, 反之称数反之称数/模转换，记模转换，记D/A 33 以一景以一景Landsat 4 或或5的的TM数字图像为例，共数字图像为例，共7个波段，其中个波段，其中6个波段有个波段有6166行，行，6166列，列，代表地面约代表地面约185 km x 185 km 的范围。每一个的范围。每一个波段约有波段约有6166x6166=38M个像元，则个像元，则7个波段个波段共共(6x38M)+24M = 252M

21、个像元，即需要个像元，即需要252M字节存储空间才能存下一景全部的字节存储空间才能存下一景全部的TM数数据。据。34 遥感卫星地面站（气象卫星接受站）提供计算兼容遥感卫星地面站（气象卫星接受站）提供计算兼容的数字磁带，输入计算机图像处理系统，形成数字的数字磁带，输入计算机图像处理系统，形成数字图像。图像。 记录在胶片上的影像可在专用设备上进行数字记录在胶片上的影像可在专用设备上进行数字化。也可以使用胶片，即透明正片，用一束强度固化。也可以使用胶片，即透明正片，用一束强度固定的光束扫描，透射光用光电增强管进行量度，通定的光束扫描，透射光用光电增强管进行量度，通过抽样和量化而成为一连串的数字，可以

22、存储在磁过抽样和量化而成为一连串的数字，可以存储在磁盘上或记录磁带上，成为数字图像文件。盘上或记录磁带上，成为数字图像文件。 35数字图像的性质与特点数字图像的性质与特点 由于传感器上探测元件的灵敏度直接影响有效由于传感器上探测元件的灵敏度直接影响有效量化的级数，因此，不同传感器提供的有效量量化的级数，因此，不同传感器提供的有效量化的级数是不同的。化的级数是不同的。传感器类型传感器类型卫星名称卫星名称有效量化级数有效量化级数信息量信息量/bitMSSLandsat646TMLandsat2568HRV(S)Spot2568HRV(PA)Spot646AVHRRNOAA102410常用的遥感数字

23、图像有效量化级数。36数字图像数字图像373839特点：(与光学图像相比) 便于计算机处理 数字存储方式 获取、传输、分发等过程中没有信息损失 抽象性强（表示在计算机里是很抽象 的“数据流”） 数字形式表示，便于建立分析模型，进行计算机解译和采用遥感图像专家系统。40二维数组表示：位置隐含（象素坐标），属性用亮度表示；遥感数字图像根据波段分： 二值数字图像：0-1 一个bit存储 单波段数字图像：0-255 一个字节（8位） 彩色数字图像： 0-255 0-2550-255 三个数字层 多波段数字图像：0-255 0-255 多个数字层41非数字式遥感图像（通常为航空摄影像片），进行数字图像处

24、理、解译前，需转换为数字图像方法：黑白扫描/彩色扫描扫描时需注意：（1）扫描的空间分辨率 一般300dpi（像片）/ 600dpi（负片）可满足要求（2）灰度级：0-255（黑白）/ RGB（彩色）42分辨率设置过程: （1）空间采样 （2）属性量化如何设置分辨率43例：将一张1:50000的航空图像扫描成分辨率是2米的数字化图(1pix=2m) 。l 1:50000 lcm=500m 1cm内要有250个pix l pix边长=1cm/250pix=0.004cm=0.001575inch (1cm=0.3937inch) 635pix/inch44图像直方图图像直方图 用平面直角坐标系表示

25、一幅灰度范围为0-n（0-255）的数字图像像元灰度分布状态，横轴横轴表示灰度级，纵轴纵轴表示某一灰度级（或范围）的像元个数占像元总数的百分比。 通过灰度直方图可以直观地了解图像特征，以确定图像增强方案并了解图像增强后的效果45数字图像直方图：以每个像元为单位，表示图像中各亮度值或以每个像元为单位，表示图像中各亮度值或亮度值区间像元出现的频率的分布图。亮度值区间像元出现的频率的分布图。直方图的作用：直观地了解图像的亮度值分布范围、峰值的直观地了解图像的亮度值分布范围、峰值的位置、均值以及亮度值分布的离散程度。直方图的曲线位置、均值以及亮度值分布的离散程度。直方图的曲线可以反映图像的质量差异。可

26、以反映图像的质量差异。正态分布：反差适中，亮度分布均匀，层次丰富，图像正态分布：反差适中，亮度分布均匀，层次丰富，图像质量高。质量高。偏态分布：图像偏亮或偏暗，层次少，质量较差。偏态分布：图像偏亮或偏暗，层次少，质量较差。小结 46直方图直方图47l图像校正：包括辐射校正、几何校正。l增强处理：增强图像中的有用信息，利于识别分析。 包括彩色增强、直方图增强、图像运算、邻域增强、频率彩色增强、直方图增强、图像运算、邻域增强、频率域增强、信息融合域增强、信息融合等。l图像变换：消除干扰和滤掉噪声，提高图像质量。l信息提取：图像分类（监督分类、非监督分类、神经网络分类、模糊分类）、空间信息提取、光谱

27、信息提取。数字图像处理的主要内容48辐射校正辐射校正l 由于由于传感器传感器响应特性和响应特性和大气大气的吸收、散射及其它的吸收、散射及其它随机因素影响，导致图像模糊失真，造成图像分随机因素影响，导致图像模糊失真，造成图像分辨率和对比度相对下降。这些都需要通过辐射校辨率和对比度相对下降。这些都需要通过辐射校正复原。正复原。l 包括：包括：系统辐射校正、大气校正系统辐射校正、大气校正4950系统辐射校正系统辐射校正l光学摄影机内部辐射误差校正镜头中心和边缘透射光的强度不一致，造成图像上不同位置的同一类地物有不同的灰度值。l光电扫描仪内部辐射误差校正A）光电转换误差；B）探测器增益变化引起的误差5

28、1大气校正大气校正概念：消除主要由大气散射、吸收引起的辐射误差的处理过程。l何时需要进行大气校正？定量信息提取；不同时相间的定量比较；不同波段间的运算；52大气校正方法大气校正方法1 公式法（较准确）；简单的相对大气校正方法：2 回归分析法；3 直方图校正法；53大气影响的粗略校正大气影响的粗略校正 （1）直方图最小值去除法）直方图最小值去除法 从图像像元亮度值中减去一个辐射偏置量（从图像像元亮度值中减去一个辐射偏置量（LP），辐射偏），辐射偏置量等于图像直方图中最小的辐射亮度值。置量等于图像直方图中最小的辐射亮度值。前提（假设）：水体（或阴影）等物体的灰度值为前提（假设）：水体（或阴影）等物

29、体的灰度值为0，大气，大气散射导致图像上这些物体的灰度值不为散射导致图像上这些物体的灰度值不为0（辐射偏置量）（辐射偏置量） 暗物体法（暗物体法（Dark-object method）54（1）直方图最小值去除）直方图最小值去除法法55（2）回归分析法）回归分析法 原理：在遥感图像上大山的阴影区或深大水体原理：在遥感图像上大山的阴影区或深大水体区域，各个波段的反射为零。同时，大气散射区域，各个波段的反射为零。同时，大气散射主要影响短波部分，波长较长的波段几乎不受主要影响短波部分，波长较长的波段几乎不受影响，因此可用其校正其它波段数据。影响，因此可用其校正其它波段数据。 方法：在不受大气影响的波

30、段（如方法：在不受大气影响的波段（如TM5或或7）和待校正的某一波段图像中，选择一系列目和待校正的某一波段图像中，选择一系列目标，将每个目标的两个待比较的波段灰度值提标，将每个目标的两个待比较的波段灰度值提取出来进行回归分析，建立线性回归方程。取出来进行回归分析，建立线性回归方程。XbaYY Y：待校正波段的图像亮度值：待校正波段的图像亮度值X X：不受大气影响波段的图像亮度值：不受大气影响波段的图像亮度值56(2)回归分析法回归分析法57公式法与卫星扫描同步进行野外光谱测量，将地面测量结果与卫星影象对应像元亮度值进行回归分析.回归方程为：Ri = a+bCiLAi =CiK+Lmin ， L

31、Ai = a1Ri +b1K=（Lmax-Lmin）/CmaxLAi 为卫星观测值（辐射亮度，辐射度）Ri 为地面反射率Ci灰度值58回归分析法回归分析法l 原理：大气散射主要影响短波部分，波长较长的波段几乎不受影响，因此可用其校正其它波段数据。l 方法：在不受大气影响的波段（如TM5或7）和待校正的某一波段图像中，选择一系列目标，将每个目标的两个待比较的波段灰度值提取出来进行回归分析，建立线性回归方程。XbaYY：待校正波段的图像亮度值X：不受大气影响波段的图像亮度值59回归法分析法回归法分析法7TMbiaiTMi60直方图校正方法直方图校正方法l 直方图校正方法： 从图像象元亮度值中减去一

32、个辐射偏置量，辐射偏置量等于图像直方从图像象元亮度值中减去一个辐射偏置量，辐射偏置量等于图像直方图中最小的亮度值；图中最小的亮度值； 前提（假设）：水体（或阴影）等物体的灰度值为前提（假设）：水体（或阴影）等物体的灰度值为0，大气散射导致，大气散射导致图像上这些物体的灰度值不为图像上这些物体的灰度值不为0（辐射偏置量）；（辐射偏置量）； 暗物体法暗物体法（Dark-object method）；）；615.2.0 基本术语基本术语 图像配准图像配准Registration 图像纠正（校正）图像纠正（校正）-Rectification 或者叫做或者叫做 Geo-referencing 地理参照地

33、理参照. 图像地理编码图像地理编码 Geo-coding 图像正射投影校正图像正射投影校正 Ortho-rectification 地面控制点地面控制点GCPs (Ground Control Points) 最邻近重采样最邻近重采样 Nearest neighbor resampling 双线性重采样双线性重采样 Bilinear resampling 立方卷积重采样立方卷积重采样Cubic convolution resampling 均方根误差均方根误差 Root Mean Square Error，RMSE625.2.1 为何要进行几何校正？ 几何畸变：遥感图像的几何位置上发生变化，产

34、生几何畸变：遥感图像的几何位置上发生变化，产生诸如行列不均匀，像元大小与地面大小对应不准诸如行列不均匀，像元大小与地面大小对应不准确，地物形状不规则变化等变形。确，地物形状不规则变化等变形。63遥感影像变形（几何畸变）的原因：遥感影像变形（几何畸变）的原因：遥感平台位置和运动状态变化的影响：遥感平台位置和运动状态变化的影响： 航高、航速、航高、航速、俯仰、翻滚、偏航。俯仰、翻滚、偏航。地形起伏的影响：产生像点位移。地形起伏的影响：产生像点位移。地球表面曲率的影响：一是像点位置的移动；二是像地球表面曲率的影响：一是像点位置的移动；二是像元对应于地面宽度不等，距星下点愈远畸变愈大，对元对应于地面宽

35、度不等，距星下点愈远畸变愈大，对应地面长度越长。应地面长度越长。大气折射的影响：产生像点位移。大气折射的影响：产生像点位移。地球自转的影响：产生影像偏离。地球自转的影响：产生影像偏离。p总之，几何畸变是：平移、缩放、旋转、偏扭、弯曲总之，几何畸变是：平移、缩放、旋转、偏扭、弯曲等共同作用的结果。等共同作用的结果。645.2.2 几何畸变校正p基本思路：把存在几何畸变的图像，校正成符合某种地图把存在几何畸变的图像，校正成符合某种地图投影的图像，且要找到新图像中每一像元的亮度值。投影的图像，且要找到新图像中每一像元的亮度值。p具体步骤 1 1）定位。计算校正后每一点所对应原图中的位置）定位。计算校

36、正后每一点所对应原图中的位置 2 2）定）定量。计算每一点的亮度值量。计算每一点的亮度值p计算方法 1 1）建立两图像像元点之间的对应关系；）建立两图像像元点之间的对应关系； 2 2）求出原图所对应点的亮度：最近邻法、双）求出原图所对应点的亮度：最近邻法、双 线性内插法、三次卷积内插法。线性内插法、三次卷积内插法。6566图像灰度重采样方法1) 最近邻法：最近邻法：距离实际位置最近的像元的灰度值作为输出图像像元的灰度值；672）双线性内插法双线性内插法：以实际位置临近的4个像元值，确定输出像元的灰度值。公式为： 4141432144332211,iiiiinmpgpppppgpgpgpgpg）

37、（ G2 G3 G4 G1 图像灰度重采样方法683）三次卷积内插法：三次卷积内插法：以实际位置临近的16个像元值，确定输出像元的灰度值。公式为： 161161,iiiiinmpgpg）（ PSF 三次样条函数 sinc函数图像灰度重采样方法69p控制点的选取1. 数目的确定:数量应当超过多项式系数的个数（(n+1)*(n+2)/2）2. 选择的原则：易分辨、易定位的特征点：道路的交叉口、水库坝址、河流弯曲点等特征变化大的地区应多选些尽可能满幅均匀选取5.2.2 几何畸变校正70数字图像几何校正过程纠正函数可有多种选择：多项式方法、共线方程方法、随机场内插方法等。其中多项式方法的应用最为普遍。

38、 准准备备 工工作作 输输入入原原始始图图像像 建建立立校校正正函函数数 确确定定输输出出图图像像的的范范围围 逐逐个个像像元元进进行行几几何何变变化化 灰灰度度的的重重采采样样 输输出出校校正正后后的的图图像像 效效果果 评评价价 71 问题：几何校正过程主要由哪两个步骤组成？分问题：几何校正过程主要由哪两个步骤组成？分别可采取哪些方法予以实现？别可采取哪些方法予以实现？Step1：变换后像元坐标位置的计算：变换后像元坐标位置的计算Step2：变换后像元灰度值（亮度值、：变换后像元灰度值（亮度值、DN）的计算）的计算725.3 数字图像增强技术简介 必要性：图像目视效果不好、信息不突出、要进

39、必要性：图像目视效果不好、信息不突出、要进行有用信息的自动或半自动识别提取时。行有用信息的自动或半自动识别提取时。 目的：提高图像质量、突出有用信息、为图像判目的：提高图像质量、突出有用信息、为图像判读和进一步分析处理服务。读和进一步分析处理服务。73常见数字图像增强处理技术 直方图变换直方图变换 基于频率域和空间域的滤波基于频率域和空间域的滤波 彩色变换（单波段假彩色密度分割，多波段。）彩色变换（单波段假彩色密度分割，多波段。） 图像运算图像运算 多光谱变换多光谱变换 数据融合数据融合74几何校正 几何畸变：遥感图像的几何位置上发生变化，产几何畸变：遥感图像的几何位置上发生变化，产生诸如行列

40、不均匀，像元大小与地面大小对应不生诸如行列不均匀，像元大小与地面大小对应不准确，地物形状不规则变化等变形。准确，地物形状不规则变化等变形。 几何畸变是平移、缩放、旋转、偏扭、弯曲等作几何畸变是平移、缩放、旋转、偏扭、弯曲等作用的结果。用的结果。75几何校正几何校正从具有几何变形的图像中消除变形的过程。一般步骤：1 确定校正方法2 确定校正公式3 验证校正方法4 对原始输入图像进行重采样76遥感图像几何畸变来源遥感图像几何畸变来源误差来源：l 内部误差：由于传感器自身的性能、结构等因素造成；l 外部误差：传感器以外的个因素造成，如地球曲率、地形起伏、地球旋转等 传感器成像几何形态影响 传感器外方

41、位元素变化畸变 地球自转的影响 地球曲率的影响 其它77传感器成像几何形态影响传感器成像几何形态影响l 成像几何特点：成像几何特点：中心投影、全景投影、斜距投影、平行投影；中心投影、全景投影、斜距投影、平行投影；l 全景投影：每条扫描线相当于中心投影，扫描视场角越全景投影：每条扫描线相当于中心投影，扫描视场角越大，边缘变形越大大，边缘变形越大全景畸变；全景畸变；l 斜距投影变形：侧视雷达；斜距投影变形：侧视雷达；78传感器外方位元素变化畸变传感器外方位元素变化畸变l 传感器外方位元素变化：指决定遥感平台姿态的6个自由度：三轴方向（X，Y，Z）和姿态角（，），其中任何一个发生变化，都会给遥感图像

42、带来不同变形，这种畸变是成像瞬间的综合影响。对不同的传感器，畸变可能不尽相同。l 航高、航速、俯仰、翻滚、偏航79中心投影中心投影传感器外方位元素变化畸变传感器外方位元素变化畸变80多中心投影多中心投影例如MSS,TM等传感器外方位元素变化畸变传感器外方位元素变化畸变81地球自转引起的误差地球自转引起的误差82地球曲率引起的误差地球曲率引起的误差83遥感图像几何校正的原理遥感图像几何校正的原理l 包括光学校正和数字纠正包括光学校正和数字纠正l 数字纠正数字纠正：通过计算机对图像每个象元逐个地解析纠正处：通过计算机对图像每个象元逐个地解析纠正处理，可以较精确改正线性和非线性变形误差。理，可以较精

43、确改正线性和非线性变形误差。l 包括两个方面（基本环节）：包括两个方面（基本环节）： （1 1）象元坐标变换；象元坐标变换； （2 2）象元灰度值重新计算（重采样）。）象元灰度值重新计算（重采样）。84遥感图像几何校正的原理遥感图像几何校正的原理85数字图像几何纠正的主要处理过程数字图像几何纠正的主要处理过程准备工作输入原始数字影象建立建立纠正纠正变换变换函数函数确定确定输出输出影象影象范围范围像元像元坐标坐标变换变换像元像元亮度亮度值重值重采样采样输出纠正后的图像准备工作准备工作：图像、地图、大地测量资料、平台轨道参数、传感器参数、控制点的选择；（具体内容可选）纠正变换函数建立纠正变换函数建

44、立：输入和输出图像间的坐标变换关系；如多项式法、共线方程法等86坐标变换坐标变换确定原始图像和纠正后图像间的坐标变换关系（x， y）（u， v）（u行数， v列数， 均为整数）有两种方案：直接纠正方案直接纠正方案：从原始图像阵列出发，依次对 其中每一个象元P（x，y）分别计算其在输出（纠正后）图像的坐标P（X，Y），并计算P（X，Y）的灰度值;间接纠正方案间接纠正方案：从空白图像阵列出发，依次计算每个象元P（X，Y）在原始图像中的位置P（x，y），然后把该点的灰度只依次计算后返送给P（X，Y）. 二者间并无本质差别，互为逆变换87坐标关系坐标关系),(ppXpyxFX ),(ppYpyxFY

45、),(PPxpYXfx ),(PPypYXfy （xp ，yp）（XP，YP）分别是任意一个像元在原始图像和纠正后图像中的坐标。直接间接88坐标变换坐标变换89图像灰度值的重采样图像灰度值的重采样l 几何校正过程中，由于校正前后图像的像元大小可能变几何校正过程中，由于校正前后图像的像元大小可能变化、象元点位置的相对变化等，不能简单用原图像象元灰化、象元点位置的相对变化等，不能简单用原图像象元灰度值代替输出象元灰度值。度值代替输出象元灰度值。l 对应的坐标值可能不是整数对应的坐标值可能不是整数l 因此，需要插值（重采样）因此，需要插值（重采样）l 三种插值方法：三种插值方法：最邻近法、双线性法、

46、三次卷积法最邻近法、双线性法、三次卷积法90最邻近法（最邻近法（nearest neighbor）l 用距离投影点（采样点）最近象元灰度值代替用距离投影点（采样点）最近象元灰度值代替输出象元灰度值输出象元灰度值l 简单、省时，保留原始图像的值简单、省时，保留原始图像的值l 边缘出现锯齿状边缘出现锯齿状91最邻近法原理最邻近法原理92最邻近法最邻近法l 优点：优点： 1 保留大量原始灰度值，保留大量原始灰度值，没有经过平滑处理，对没有经过平滑处理，对于区分植被类型、识别于区分植被类型、识别线性特征等有重要意义线性特征等有重要意义 2 简易、省时简易、省时 3 分类前使用分类前使用 4 适合于专题

47、文件适合于专题文件93原始图像纠正后图像（最邻近插值）最邻近法纠正效果最邻近法纠正效果94双线性内插法双线性内插法 取取(z，y)点周围的点周围的4邻点，在邻点，在y方向方向(或或x方向方向)内插内插二次，再在二次，再在x方向方向(或或y方向方向)内插一次，得到内插一次，得到(x，y)点的亮度值点的亮度值f(x，y)，该方法称双线性内插法，该方法称双线性内插法.95双线性（双线性（bilinear）l 优点优点 1 较平滑，没有锯较平滑，没有锯齿状齿状 2 与最邻近法相与最邻近法相比，空间信息更准比，空间信息更准确些确些 3 常用于改变象元常用于改变象元大小时，如数据融大小时，如数据融合合l缺

48、点：缺点：象元值被平均化，某象元值被平均化，某些地物边缘更平滑，某些地物边缘更平滑，某些极值可能丢失些极值可能丢失96双线性插值效果原始图像纠正（双线性插值）97三次卷积法（三次卷积法（cubic convolution）l 取与投影点邻近的取与投影点邻近的16个象元灰度值（个象元灰度值（44） ，计算输，计算输出象元的灰度值出象元的灰度值 取与计算点(x，y)则围相邻的16个点，与双向线性内插类似，可先在在某一入向上内插，如先在x方向上，每4个值依次内插4次，再根据这四个计算结果在y方向上内插，得到f(x，y)。每一组4个样点组成一个连续内插函数。实际上是一种卷积运算，也叫三次卷积内插。98

49、三次卷积法三次卷积法l 优点：优点： 1 与其它重采样方法相比，均与其它重采样方法相比，均值和标准偏差与原始象元的相值和标准偏差与原始象元的相一致一致 2 改变象元大小时使用（改变改变象元大小时使用（改变幅度更大时）幅度更大时） 3 可以锐化图像、平滑噪声，可以锐化图像、平滑噪声，实际的效应与数据有关实际的效应与数据有关99原始图像几何纠正（三次卷积）三次卷积法处理效果三次卷积法处理效果100l增强处理：增强图像中的有用信息，利于识别分析。 包括彩色增强、直方图增强、图像运算、邻域增强、频率彩色增强、直方图增强、图像运算、邻域增强、频率域增强、信息融合域增强、信息融合等。遥感图像处理图象增强 图象增强是指对图像的某些特征，如边缘、轮廓、对比度进行强调或尖锐化，以便于显示、观察或进一步地