图像增强与复原

图像增强与复原第一页，共三十七页，编辑于2023年，星期二图像在生成和传输过程中常受到各种噪声的干扰和影响，使图像质量下降。为了抑制噪声改善图像质量，必须对图像进行平滑处理，可以在空域或频域中进行。数字图像的平滑技术分为两大类：一、全局处理，即对噪声图像的整体或大的块进行校正以得到平滑的图像。缺点：计算量大。二、局部算子，仅对某一像素的局部小邻域的一些像素加以运算。优点：计算效率高。第二页，共三十七页，编辑于2023年，星期二2.1模板操作模板（Template）又称卷积（Ｃonvolution）,可以理解为矩阵1111●

1111192124●

2121116

1112●

1111110第三页，共三十七页，编辑于2023年，星期二模板操作设图像中某像素的灰度值为f(x,y),它的邻域为3×3，点集总数为9，则平滑后这点的灰度值为

f’(x,y)=1/9∑

f(i,j)03这种类似点象矩阵的表示方法，称为模板（template）。中间的黑点表示中心元素，即用哪个点作为处理后的元素。1111●

111119第四页，共三十七页，编辑于2023年，星期二举例原图为：11111222223333344444模板为：1·

001操作后图像为：1.51.51.51.5×2.52.52.52.5×3.53.53.53.5×

×

×

×

×××表示边界上无法进行模板操作的点，通常复制原图的灰度1/2第五页，共三十七页，编辑于2023年，星期二模板操作实现的实际上就是邻域运算（NeighborhoodOperation）即，某个像素点的结果不仅和本像素灰度有关，而且和其邻域的值有关。一般常用的模板并不大，如3×3，3×5。很多专用的图像处理系统，用硬件来完成模板运算，以提高处理速度。平滑模板的思想是通过一点和周围8个点的平均来去除突然变化的点，从而滤掉一定的噪声，其代价是图像有一定的模糊。常用的模板有BOX模板和高斯模板：第六页，共三十七页，编辑于2023年，星期二1111●

1111192124●

2121116BOX模板高斯模板BOX模板对所有9个点都一视同仁，没有考虑各点位置的影响。高斯模板进行改进，距离越近的点，加权系数越大第七页，共三十七页，编辑于2023年，星期二构造卷积函数１、当在图象上逐个像素的移动模板时，图象边界数据无法进行卷积操作方法：在进行卷积运算前，先复制原图象，保留边界数据。需要解决的问题：２、在进行卷积操作时，求出的新值可能远远超出像素值的有效范围，因此需要最卷积结果按比例放缩方法：对卷积操作的结果除以放缩比例（非０）３、当模板中有负数时，卷积结果有可能是负数方法：一、负数取０二、取绝对值第八页，共三十七页，编辑于2023年，星期二以３Ｘ３模板为例构造模板类publicclassMatrix3x3:ImageInfo{inttopLeft=0,topMid=0,topRight=0;intmidLeft=0,center=1,midRight=0;intbottomLeft=0,bottomMid=0,bottomRight=0;intscale=1;intkernelOffset=0;／／用于增加或减少图象的亮度topLefttopMidtopRightmidLeftcenter.

midRightbottomLeftbottomMidbottomRightscale1第九页，共三十七页，编辑于2023年，星期二以３Ｘ３模板为例构造模板类

///<summary>///初始化窗口所有点为同一权值///</summary>///<paramname="degree">权值</param>

publicvoidInit(intdegree){topLeft=topMid=topRight=midLeft=center=midRight=bottomLeft=bottomMid=bottomRight=degree;}//endofInit进行卷积转换的代码第十页，共三十七页，编辑于2023年，星期二2.2平滑处理的几种方法2.2.1邻域平均法(near-areaaveragemethod)2.2.2中值滤波2.2.3多图象平均法第十一页，共三十七页，编辑于2023年，星期二2.2.1邻域平均法(near-areaaveragemethod)原理：给定一幅N×N的图像f(x,y)，取围绕(x,y)点的预设邻域内的几个像素点（不含(x,y)点）的灰度级平均值作为新图像中该点(x,y)的灰度级，并对N×N个像素点都这样做，由此构成新图像g(x,y)。简单邻域平均法：

x,y=0，1，…N-1,s是(x,y)点邻域内点的坐标（不包含点(x,y)）的集合，M是集合s内坐标点的总数。第十二页，共三十七页，编辑于2023年，星期二OOOO*OOOOOOOOOOOOOOOO*OOOOOOOOOOOOOO*OOOOOO*OOOO3×3邻域5×5邻域4邻域8邻域4邻域={(x+1,y),(x,y+1),(x-1,y),(x,y-1)}几种邻域举例第十三页，共三十七页，编辑于2023年，星期二对于3×3邻域OOOO*OOOO其模板为：1111●

111119OOOO*OOOO对于8邻域其模板为：1110●

111118第十四页，共三十七页，编辑于2023年，星期二加权邻域平均法：举例:高斯模板2124●

2121116=1

∑

a(m,n)=第十五页，共三十七页，编辑于2023年，星期二阈值邻域平均法：T是预先设定的阈值，当某些点的灰度值与其邻域点的灰度平均值之差不超过阈值T时，仍保留这些点的灰度值。当某些点的灰度值与其邻域点的灰度平均值差别较大时，这些点必然是噪声，这时再取其邻域平均值作为这些点的灰度值。第十六页，共三十七页，编辑于2023年，星期二2.2.2中值滤波在一定条件下，中值滤波可以克服线性滤波器所带来的图像细节模糊而且对滤除脉冲干扰及图像扫描噪声非常有效,但是对一些细节多，特别是点、线、尖顶细节较多的图像则不宜采用中值滤波的方法。中值滤波也是一种局部平均平滑方法。

中值滤波的目的是在保护图像边缘的同时去除噪声。第十七页，共三十七页，编辑于2023年，星期二

中值滤波是一种非线性的图像平滑方法，它对一个滑动窗口内的诸像素灰度排序，用其中值代替窗口中心像素(x,y)原来的灰度（若窗口中有偶数个像素，则取两个中间值的平均）原理例如，若窗口长度为5，窗口中像素的灰度值分别为80、90、200、120、110。首先按从小到大排序，结果为80、90、110、120、200，其中间位置上的值为110。于是原来窗口正中的灰度值200就由110代替。如果200是一个噪声的尖峰，则将被滤除。然而如果它是一个信号，那么此法处理的结果将会造成信号损失。第十八页，共三十七页，编辑于2023年，星期二例题由一个序列为{0，3，4，0，7}，当窗口m=5时，试分别求出采用中值滤波和简单邻域平均法的结果。中值滤波：重新排序后为{0，0，3，4，7}，则中值滤波的结果{0，3，3，0，7}简单邻域平均法：（5邻域）（0+3+4+0+7）/5=2.8，则平滑滤波的结果{0，3，2.8，4，7}第十九页，共三十七页，编辑于2023年，星期二几种信号进行中值滤波示例1（m=5）原信号中值滤波a）阶跃b）斜坡可以看到中值滤波不影响阶跃函数和斜坡函数，因而对图象的边缘由保护作用。第二十页，共三十七页，编辑于2023年，星期二原信号中值滤波c）单脉冲几种信号进行中值滤波示例2（m=5）d）双脉冲对于持续周期小于窗口尺寸的1/2的脉冲将进行滤除，因而可能损坏图像中的某些细节。持续周期大于窗口１／２会怎样第二十一页，共三十七页，编辑于2023年，星期二

中值滤波窗口也可以是二维的，而且可以由各种不同形状、如线、方形、圆形、十字形等。在实际使用窗口时，窗口的尺寸一般先取3，再取5，依次增大，直到滤波效果满意位置。对于由缓变的较长轮廓线物体的图像，采用方形或圆形窗口较合适；对于包含尖顶角物体的图像，采用十字形窗口较合适。第二十二页，共三十七页，编辑于2023年，星期二示例000000000000000011100001610000111000000000000000000001110011100111000000原图处理后的图原图中的6和周围的灰度相差很大，是一个噪声点，采用长度为3的窗口进行中值滤波，得到处理后的图，噪声点被去除。第二十三页，共三十七页，编辑于2023年，星期二课堂习题000111001/32/311001/43/411000111原图分别利用BOX模板、高斯模板和中值滤波法对其进行处理。BOX模板：1/3*（11.1）高斯模板：1/4*（12.1）从原图中可以看出左边的区域灰度值低，右边灰度值高，中间有一条明显的边界，称这类图像为台阶（step）。第二十四页，共三十七页，编辑于2023年，星期二习题结果BOX模板处理后：

01/32/31

1/91/32/38/91/121/32/311/1201/32/31高斯模板处理后：01/43/411/121/32/311/121/165/1611/1615/1601/43/41中值滤波处理后：001101/32/31

01/43/41

0011从处理结果可看出，应用平滑模板，图像平滑了，但使边界模糊了。二中值滤波很好的保持了原来的边界。第二十五页，共三十七页，编辑于2023年，星期二2.2.3多图像平均法

多图像平均法是利用对同一景物的多幅图像取平均来消除噪声产生的高频部分，在图象采集中常应用这种方法来去除噪声。第二十六页，共三十七页，编辑于2023年，星期二课堂习题有已幅图像由于受到干扰，图中有若干各亮点（灰度值为255）如下图。试问此类图像如何处理，并将处理后的图像画出来。11187422552333332554333332554633452558234678第二十七页，共三十七页，编辑于2023年，星期二2.3平滑的实现利用模板操作可以实现邻域平均法，可以用在光滑边缘过于清晰或对比度过于强烈的区域，产生模糊效果来柔化边缘．以ＢＯＸ为例．第二十八页，共三十七页，编辑于2023年，星期二

///<summary>///对图像进行平滑处理///</summary>///<paramname="b">位图流</param>///<returns></returns>publicBitmapSmooth(Bitmapb){//111//111//111/9Matrix3x3m=newMatrix3x3();m.Init(1);m.Scale=9;returnm.Convolute(b);}//endofSmooth第二十九页，共三十七页，编辑于2023年，星期二中值滤波经过中值滤波后，图象中的随机噪声将被有效的消除．原因：亮度值发生随机突变的像素，经过排序后，要么排在队首要么排在队尾，而中心像素的新值是取自队列最中间的那个像素的值．第三十页，共三十七页，编辑于2023年，星期二第三十一页，共三十七页，编辑于2023年，星期二中值滤波编程思路１、获取用户指定的滤波窗口２、取出当前像素和周围ＮＸＮ个像素，对他们进行排序３、取得中值４、将中值作为亮度值输出设:滤波窗口长度为len;用数组sequence存放窗口中的值;取得的中值为count.写出中值滤波的算法程序第三十二页，共三十七页，编辑于2023年，星期二boolisMovable=true;intt=0;while(isMovable){isMovable=false;for(inti=1;i<len;i++){if(sequence[i-1]>sequence[i]){