第三章遥感数字图像增强处理

遥感数字图像增强处理•图像增强的目的改变图像的灰度等级,提高图像的对比度;消除边缘和噪声,平滑图像;突出边缘或线状地物,锐化图像;合成彩色图像;压缩图像数据量,突出主要信息等•图像增强的方法空间域增强(点运算,邻域运算)频率域增强彩色增强多图像代数运算多光谱图像增强直方图变换累积直方图1789101114526714121534786911214788984591112108101115161013136916131210灰度级012345678910111213141516累积像元数02458101317232833374043454749累积比例值0.000.040.080.10.160.200.270.350.470.570.670.760.820.880.920.961.00为最大灰度级）；，，，，＝为灰度级（1L1L...210k)()(0kxhxckjjk直方图均衡化直方图均衡化是将原图像的直方图通过变换函数变为均匀的直方图，然后按均匀的直方图修改原图像，从而获得一幅灰度分布均匀的新图像。可以证明，这个变换函数就是累积直方图。变换式：为最大灰度值。的概率为原图像任一灰度出现为像元数，变换后新图像的值1L;hN)(1a0kbkjjaakbxxhNLx对一幅图像进行直方图均衡化的具体步骤如下：（1）统计原图像每一灰度级的像元数和累积像元数（2）根据变换函数式计算每一灰度级xa均衡化后对应的新值，并对其四舍五入取整，得到新灰度级xb（3）以新值代替原灰度值，形成均衡化后的新图像。（4）根据原图像像元统计值对应找到新图像像元统计值，作出新直方图。原灰度级xa像元统计值累积像元统计值变换后值新灰度级新像元统计值0000001220.7142241.313151.6214382.63552103.3363134.24374175.66486237.58695289.1951053310.81151143712.11241234013.11331334314.01431424514.71541524715.3151624916162直方图规定化直方图规定化是指使一幅图像的直方图变成规定形状的直方图而对图像进行变换的增强方法.规定的直方图可以是一幅参考图像的直方图,通过变换,使两幅图像的亮度变化规律尽可能地接近;规定的直方图也可以是特定函数形式的直方图,从而使变换后图像的亮度变化尽可能地服从这种函数分布.直方图规定化设T(xa)为原图像直方图均衡化的变换函数,G(yc)为参考图像直方图均衡化的变换函数,变换后的灰度值均为Zb,由上述可知)]([)()()()()(1100abckjjccbkjjaabxTGzGyyhcyGZxhaxTZ直方图规定化直方图规定化的具体步骤如下:(1)做出原图像的直方图(2)做出原图像的累积直方图Zb=T(xa),对原图像进行均衡化变换(3)做出参考图像的直方图或确定参考直方图(4)做出参考图像累积直方图zb=G(yc),进行均衡化变换(5)对于原图像中的每一灰度级xa的累积值zb,在参考累积直方图中找到对应的累积值G(yc);(6)以新值yc替代原灰度值xa,形成均衡化后的新图像.(7)根据原图像像元统计值对应找到新图像像元统计值,做出新直方图。10911121011109101112111211911101213141312121012141513111314161614151412121516151413121216131415参考图像原灰度级xa像元统计值ha(xa)累积像元统计值T(xa)对应参考累积像元值G(yc)新灰度级yc新像元统计值hd(xd)0000001／1620.040.069/1652／1620.080.063／1610.100.064／1630.160.1810/1655／1620.200.186／1630.270.3311/1677／1640.350.338／1660.470.5112/16119／1650.570.5110／1650.670.6713/16511／1640.760.8214/16712／1630.820.8213／1630.880.9215/16514／1620.920.9215／1620.961.00141211.00灰度级yc9/1610/1611/1612/1613/1614/1615/161像元统计值hc(xc)36798754累积像元统计值G(yc)0.060.180.330.510.670.820.921.00空域增强-点运算•线性变换y=a•x+b单一线性、分段线性、特定区域或覆盖类型的局部线性空域增强-点运算•非线性变换指数扩展xb=a•ebxa+c扩展高亮度区对数扩展xb=a•log(xa+1)+c扩展低亮度区空域增强-点运算•非线性变换•三角函数扩展xb=a•sin(b•xa)+c扩展中部亮度区•xb=a•tan(b•xa)+c扩展两端亮度区空域增强-邻域增强•邻域对于图像中的某个像元f（x,y），把以像元为中心一定距离内的像元集合Axy=｛x±p，y±q｝（p,q取任意整数）叫做该像元的邻域。(x,y)(x,y)4-邻域8-邻域建立在离散卷积基础上的卷积运算是在空间域上对图像进行邻域变换的运算。选定一个卷积函数，又称为模板，它实际上是一个M*N大小的小图像。图像的卷积运算是运用模板来实现的，如图，从图像的左上角开始，将一个给定大小的模板，逐行逐列依次放在图像的每一个像元位置上，计算两者之间对应各点的乘积并求和，以和数作为中心像元的输出值，从而产生新的图像。卷积运算如果假设卷积函数为H(m,n)，其大小为M,N。对应图像窗口内灰度值为f(m,n)，则卷积运算可通过下式表示。MmNnnmHnmfyxg11),(),(),(MmNnMmNnnmHnmHnmfyxg1111),(),(),(),(在实际应用中，经常使图像窗口与模板像元的灰度值对应相乘再相加，相加的总和再除以模板内所有值的和作为中心像元新的灰度值。模板运算的公式为平滑•均值平滑设待处理图像f（x,y）有N行N列，平滑后的图像为g（x,y）。均值平滑指对原图像每个像元在以它为中心的邻域A内取平均值，作为该像元新的灰度值。即AjijifMyxg),(),(1),(均值滤波器模板010101010111101111均值滤波器Mean5x5原始图象Mean11x11中值滤波器在邻域平均法中，是将n×n局部区域中的灰度的平均值作为区域中央象元的灰度值。而在中值滤波中，是把局部区域中灰度的中央值作为区域中央象元的值。如，在3×3区域内进行中值滤波，是将区域内9个灰度值按由小到大排列，从小的一方开始的第5个值即为中央象元的值。用这种非线性的滤波，比邻域平均法可以在很大的程度上防止边缘的模糊。)),((),(yxfmedianyxg351012162481055683746745833081019试用1*3和3*3的窗口对此进行中值滤波空间锐化为了突出边缘和轮廓、线状目标信息，可以采用锐化的方法。锐化可使图像上边缘与线状目标的反差提高，因此也称为边缘增强。平滑通过积分过程使得图像边缘模糊，图像锐化则通过微分使图像边缘突出、清晰。空间域图像锐化1.梯度法最常用的微分方法是梯度法。1梯度法•矢量微分----梯度二元函数f(x,y)在坐标点(x,y)处的梯度向量的定义：yfxfffyxfGyx'')],([梯度的幅度（模）为各分量的平方和再求平方根：22)()()],([yfxfyxfG空间域图像锐化222121y)gradf(x,),(,),(ttyyxftxyxft设用绝对值可得到以下近似的结果21),(ttyxgradf•连续域的微分----离散域的差分对于数字图像，连续导数形式可以用求差来近似表示)1,(),(),1(),(21jifjiftjifjift21),(ttjigradf梯度算法对应的模板为：10-101-100t1＝t2=常见的梯度算子模板100-101-10RobertsRoberts梯度Roberts梯度采用交叉差分的方法。设：)1,(),1()1,1(),(21yxfyxftyxfyxftt1=t2＝常见的梯度算子模板-101-101-101-1-1-1000111Prewitt-101-202-101-1-2-1000121SobelPrewitt和Sobel梯度t1＝t2＝t1＝t2＝66612121212666121212126661212121266612121212666666666666666666666从上面可知，Sobel算子、Prewitt算子不像普通梯度算子那样用两个像素之差值，而用两列或两行加权和之差值，其优点为：由于引入了平均因素，因而对图像中的随机噪声有一定的平滑作用由于它是相隔两行或两列的差分，故边缘两侧的象元得到了增强，边缘显得粗而亮。适用情况梯度算子法取值式说明保留背景，强调边缘或：：为规定的亮度级研究边缘灰度级的变化，但不受背景影响：为规定的亮度级只对边缘位置感兴趣、：为规定的亮度级其它),,(),(),(),(yxfTyxfyxfyxg其它),,(),(),(yxfTyxfLyxgggL其它,),(),(),(bLTyxfyxfyxgbL其它,),(),(bgLTyxfLyxggLbL水平、垂直梯度Roberts梯度Sobel梯度或检测垂直边界：检测水平边界：检测对角线边界：图像锐化2.定向检测：当有目的地检测某一方向的边、线或纹理特征时，可选择特定的模板卷积运算作定向检测。常用的模板为：1011011011h1211211212h1110001113h1112221114h或0111011105h1101010116h2111211127h1121212118h图像锐化3.拉普拉斯算子•拉普拉斯算子处理是常用的边缘增强处理算子，它是各向同性的二阶导数22222yfxff对数字图像来讲，的二阶偏导数可表示为：),(2)1,()1,(),(),(2),1(),1(),1(),(),(),1(),(),1(),(2222jifjifjifyyxfjifjifjifjifjifjifjifjifjifxyxfxx),(yxf•为此，拉普拉斯算子为：),(4)1,()1,(),1(),1(22222jifjifjifjifjifyfxff以模板形式表示为：010141010可见数字图像在某点的拉普拉斯算子，可以由中心像素点灰度级值和邻域像素灰度级值通过加减运算来求得。拉普拉斯锐化：用原图像的值减去模板运算结果的整倍数，即：),(),(),(2yxfkyxfyxg拉普拉斯算子114111h115112h彩色增强①伪彩色密度分割：把一幅黑白图像的不同灰度按一定的函数关系变换成彩色，得到一幅彩色图像的方法。密度分割法是对单波段黑白遥感图像按灰度分层，对每层赋予不同的色彩，使之变为一幅彩色图像。例如：灰度范围为0～15的赋为一种颜色；灰度范围15～20赋为一种颜色；…，依此类推，生成一幅彩色图像。密度分割中的彩色是人为赋予的，与地物的真实色彩毫无关系，因此也称为伪彩色。②彩色合成与假彩色合成彩色合成时，波段如何选择，可通过求最佳指数OIF进行选择：最佳指数大