图像处理(增强)

1、第4章 遥感图像处理w4.1 光学原理与光学处理w4.2 数字图像概念w4.3 数字图像的校正w4.4 数字图像增强w4.5 多元信息复合w4.6 遥感图像计算机分类4.4 数字图像增强（Image Enhancement ）w一一. 对比度变换（对比度变换（contrast conversion ） 也称为辐射增强，或直方图拉伸w 对比度变换：通过改变图像像元的亮度值来改变对比度，从而改善图像质量的图像处理方法。w 直方图：以横坐标代表像元的亮度值，纵坐标代表某一亮度或亮度间隔的像元数的统计图。从直方图的形态可以判断图像的质量。一. 对比度变换一. 对比度变换w1. 线性变换w线性变换：按比

2、例扩大原始亮度值的范围。将过于集中的像元分布区域扩展开来，增加图像表现的层次性。 a1 ，a2 b1，b2变换函数为直线方程：111212)(baxaabbxab一. 对比度变换w1. 线性变换xbb2线性变换分段线性变换xbxaxab1a1a2b0b1b2b3a0a1a2a3一. 对比度变换w2.非线性变换 变换函数为非线性函数，即对暗区与亮区不同比例地扩展，常用的有：n指数变换指数函数n对数变换对数函数n查表法1. 对比度变换w2.非线性变换指数变换（突出亮区，压制暗区）对数变换（突出暗区，压制亮区）xbb2xab1a1a2xbb2xab1a1a2一. 对比度变换1. 对比度变换w 3.

3、直方图变换w A.直方图均衡（Histogram equalization） 1. 对比度变换1. 对比度变换w 3. 直方图变换w A.直方图标准化（Histogram normalization ）4.4 数字图像增强w二. 空间滤波（邻区法处理）w 空间滤波：通过处理像元与其周围相邻像元的关系，突出图像上某些特征的图像处理方法。w 图像卷积运算：选定一“模板”（即滤波器）作为卷积函数，将模板对应像元内的灰度值与模板像元亮度值对应相乘后相加，然后将计算结果作为窗口中心的像元的灰度值。 卷积运算v 空间滤波处理：对数字图像中各行或各列所构成的复式波形曲线，通过滤波算子将其分解成不同频率的简单

4、波形曲线，根据需要选择适宜的频率波形曲线，重新构成新的图像。亮度值像元高通滤波低通滤波 二. 空间滤波（Spatial filtering ）二. 空间滤波w1.平滑（低通滤波）：消除变化过大的区域，使亮度平缓或去掉不必要的“噪声”点。n均值滤波l卷积运算公式：),(1),(11nmMNjirMmNn 二. 空间滤波l模板（33）：11111111191),(nmt11110111181),(nmt或均值滤波均值滤波均值滤波前均值滤波后中值滤波n卷积运算法则：取窗口内所有像元亮度值的中值为中间像元的亮度值。13模板均值滤波与中值滤波w 图象亮度阶梯状变化时均值平滑w 图象中有突出亮点噪声中值平

5、滑二. 空间滤波w2.锐化（高通滤波）：突出图像中亮度变化大的细节。如突出边缘和线性地物等。经锐化后，原遥感图像只保留边缘的信息。二. 空间滤波w罗伯特梯度：=t1+ t2 10011t01102t二. 空间滤波w索伯尔梯度：=t1+ t2 1210001211t1012021012t二. 空间滤波w拉普拉斯算法010141010),(nmt二. 空间滤波w3.定向检测：n检测垂直边界101101101),(nmt121121121或二. 空间滤波w3.定向检测：n检测水平边界111000111),(nmt111222111或二. 空间滤波w3.定向检测：n检测对角线边界 t（m，n）=，01

6、1101110110101011112121211211121112 彩色合成TM432图 地物边界形状信息 水平线性地物信息 垂直线性地物信息4.4 数字图像增强w三. 彩色变换w人眼对色彩有较强的识别力，因此往往将黑白影像处理后赋予各种颜色。三. 彩色变换三. 彩色变换w1.单波段彩色变换：彩色密度分割n将一幅图像的亮度值分成若干个等级，每个等级以不同的颜色来表示。 DEM数字高程图像 伪彩色坡度分割图 伪彩色高度分割图三. 彩色变换2.多波段彩色变换：彩色合成 RGB Color Composite n真彩色合成（natural color composite. ）： TM 321n假彩

7、色合成（false color composite ）：l标准假彩色合成：TM 432l其他：TM 453，TM742 彩色合成TM321图钦州市 彩色合成TM432图钦州市 彩色合成TM742图钦州市TM1TM4TM7754TM2TM5321543645TM3TM6432742643423南宁NDVI 2000与2006年 两时相假彩色合成图红色-植被指数减少 青色-植被指数上升 灰色、黑色-植被指数不变R:1998NDVIG:2006NDVIB:2006NDVI三. 彩色变换w3.HLS变换：hue，lightness，saturationnHLS变换：将RGB模式转换为HLS模式。nH：

8、0360红、黄、绿、青、蓝、品nL： 01 黑 白nS： 01 极不纯 纯l S=0时，H无定义，为黑或灰或白l S=1时，L=0.5 三. 彩色变换w对于RGB：三种颜色混合相加。其亮度分别为LR，LG，LB. 设： Lmax=max（ LR，LG，LB） Lmin=min（ LR，LG，LB）三. 彩色变换w对于HLS：w明度L=（ Lmax + Lmin ）/2 当Lmax= Lmin 时，即LR=LG=LB， 色调为灰色，H无定义，S=0 三. 彩色变换w饱和度S=（Smax-Smin）/（Smax+Smin） （L0.5） 或 =（Smax-Smin）/(（1-Smax)+(1-Sm

9、in）) （L0.5） 三. 彩色变换w色调H =60(HG-HB) /H （当HR=HRrmax） 黄品或 =60(2+(HB-HR) /H ) （当HG=HGrmax ） 青黄或 =60(4+(HR-HG) /H （当HB=HBmax ） 品青 式中： H =（ H max - H min ）4.4 数字图像增强w四. 图像运算（Arithmetic Operations）利用两幅或多幅已经完成空间配准的不同波段的影像，进行加、减、乘、除四则混合运算，达到提取或去掉某些信息的目的。四. 图像运算n1.差值运算: fD(x,y)=f1(x,y) f2(x,y)l红外波段红光波段:突出植被信息

10、l事件发生前后相减:动态监测l错行(列)后相减:边缘增强，四. 图像运算n2.比值运算： fR(x,y)=f1(x,y) / f2(x,y)l植被指数：近红外/红 或 (近红外-红) /(近红外+红) 用于突出植被特征,提取植物类别或估算生物量。四. 图像运算n2.比值运算：lTM1/TM2,可以消除地形的影响。亮度值TM1TM2TM1/TM2阳坡砂岩28430.65阴坡砂岩22340.651998年南宁NDVI 2006年南宁NDVI 0.54.4 数字图像增强w五. 多光谱变换n对遥感图像实行线性变换。n多光谱空间就是一个n维坐标系，每个坐标轴代表一个波段，每个像元点的位置可表示为一个n维

11、向量 X，每个分量xi 为该点在第i个坐标轴上的投影，即亮度值。五. 多光谱变换lK-L变换主成分变换（Principal Component Analysis ）lK-T变换缨帽变换（Tasseled Cap Analysis ） K-L变换（ 主成分分析）v K-L变换，也称主成分变换。它是对多光谱图像X进行K-L变换，得到矩阵A，而后进行线性组合，产生新的光谱图像Y，公式为： Y = A X 式中:X原图像像元的光谱矢量； Y变换后主分量空间的光谱矢量； A变换矩阵v经过K-L变换后，新图像Y中的各yi之间具有最小相关性。 K-L变换（ 主成分分析）v多光谱图像中，各波段数据间存在相关的

12、关系，利用主成分分析法将原图像中所含的大部分信息用虚拟的少数波段来表示。x1x2y1在多波段数据彩色合成中，能够同时被感知仅限于三个波段RGB，因此通常利用主成分分析法将多波段数据压缩到三个波段。 K-L变换（ 主成分分析）x1x2y1y2 K-L变换（ 主成分分析）w 对n维的原始图像向量X，作如下的线性变换:YAX An n正交矩阵Yn维新的图像特征向量其中: Y中n（n）个特征为“有效特征”，剩余的n-n 个为“噪声特征”。A的确定：通过计算X的协方差矩阵的特征值和特征向量求得。 K-L变换（ 主成分分析） 经变换后，新向量的各分量之间具有最小相关性。第一主成分集中了最大的信息量，前三个

13、主成分集中了原向量的大部分的信息量，可以用前几个主成分向量代替原向量，从而实现了降维处理。当特征值ii+1，便可舍去yi+1特征。 K-T变换（缨帽变换）v K-T变换，也称缨帽变换。也是一种线性变换。主要针对MSS、TM图像数据进行，其变换公式为： Y = B X 式中:X原图像像元的光谱矢量； Y变换后主分量空间的光谱矢量； B变换矩阵。对于MSS变换矩阵为44， TM变换矩阵为66。 K-T变换（缨帽变换）B的确定：经长期研究，由专家提出。的确定：经长期研究，由专家提出。对于对于TM数据：数据： B= 0.30 0.28 0.47 0.56 0.51 0.19 -0.28 -0.24 -

14、0.54 0.72 0.08 -0.18 0.15 0.20 0.33 0.34 -0.71 -0.46 -0.82 -0.08 0.44 -0.06 0.20 -0.28 -0.33 -0.05 0.11 0.19 -0.44 0.81 0.11 -0.90 0.41 0.06 -0.03 0.02 TM1 TM2 TM3 TM4 TM5 TM7 K-T变换（缨帽变换） TM数据变换后矢量分量为： Y=(y1,y2 ,y3,y4,y5,y6 )T 各分量的物理意义：其中y1为亮度分量，反映总体的亮度值； y2为绿度分量，TM4-TM3-TM2-TM1； y3为湿度分量, TM1+TM2+TM

15、3+TM4-TM5-TM6； K-T变换（缨帽变换）（y1,y2 ）：植被视面；（y1 ,y3）：土壤视面；（y3 ,y2）：过渡区视面； 在由y1 , y2 , y3构成的三维坐标体系中，农作物的生长过程的形态，很像一顶带穗的帽子，所以也把K-T变换形象地称为缨帽（ Tasseled Cap ）变换。K-T变换（缨帽变换）B1B2B3B4IgvIsbK-T变换（缨帽变换）w经变换后，：n植被的特征点集中在帽空间中，不同植被周期，特征点沿垂直于帽的轴线发生变化。n土壤地物特征点主要集中在帽的底部，沿过特征空间原点的一线轴射线散布。TM543合成图主成分变换123合成图缨帽变换-123合成图基于

16、主成分变换的分类图融合分类图第4章 遥感图像处理w4.1 光学原理与光学处理w4.2 数字图像概念w4.3 数字图像的校正w4.4 数字图像增强w4.5 多元信息复合w4.6 遥感图像计算机分类4.5 多元信息复合w1、 遥感信息的复合v 遥感信息复合： 将不同来源的遥感数据或其它数据统一到同一地理坐标系中，采用一定的算法生成新的信息或合成图像的过程。 1、 遥感信息的复合HLS变换信息融合变换信息融合DEM信息融合信息融合基于HLS变换的多源影像融合w 基本原理：n将低分辩率的多光谱图像进行重采样，使其图像大小和高分辨率全色波段的图像大小一致。n对采样得到的多光谱图像进行HLS正变换处理，得到色调H、饱和度S和亮度L三个分量。n用高分辨率的全色波段替换亮度I，再进行HLS逆变换