数字图像处理- 图像平滑与锐化资料讲解.ppt

《数字图像处理- 图像平滑与锐化资料讲解.ppt》由会员分享，可在线阅读，更多相关《数字图像处理- 图像平滑与锐化资料讲解.ppt（46页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、数字图像处理-图像平滑与锐化1数字图像处理图像平滑与尺度空像平滑与尺度空间2数字图像处理图像平滑像平滑低通滤波法均值滤波法中值滤波法高斯滤波法3数字图像处理低通低通滤波法波法低通滤波器，就是通过低频信号且抑制高频信号的滤波器。一个理想的低通滤波器，就是能够完全剔除高于截止频率的所有信号，并且保留所有低于截止频率的信号。在实时应用中，实际滤波器可以通过将信号延长一小段时间以得到未来的一小部分信号信息，这样就可以近似的实现理想滤波器。巴特沃斯滤波器切比雪夫滤波器椭圆滤波器4数字图像处理巴特沃斯巴特沃斯滤波器波器5数字图像处理切比雪夫切比雪夫滤波器波器切比雪夫滤波器是在通带或阻带上频率响应幅度等波纹

2、波动的滤波器。切比雪夫滤波器在过渡带比巴特沃斯滤波器衰减要快，但是频率响应的幅频特性不如巴特沃斯滤波器平坦。在通带波动的切比雪夫滤波器的数学表达为在阻带波动的切比雪夫滤波器的数学表达为6数字图像处理椭圆滤波器波器7数字图像处理低通低通滤波器波器巴特沃斯滤波器通带波动下的切比雪夫滤波器阻带波动下的切比雪夫滤波器椭圆函数滤波器8数字图像处理傅里叶傅里叶变换9数字图像处理傅里叶傅里叶变换示意示意图10数字图像处理傅里叶傅里叶变换示意示意图11数字图像处理快速傅里叶快速傅里叶变换12数字图像处理基于傅里叶基于傅里叶变换域的低通域的低通滤波波涉及4种图像初始、中间或最终结果，和三个主要处理步骤。图像结果

3、包括原始图像、原始变换域、滤波后的变换域和滤波后的图像；处理步骤包括傅里叶正变换、低通滤波和傅里叶逆变换。将原始图像转化为复数数组进行傅里叶正变换，输出得到原始变换域；对原始变换域进行频域上的处理，如低通滤波，得到滤波后的变换域；将滤波后的变换域进行傅里叶逆变换，便得到最终的复数数组，取复数的实部信息便可以得到低通滤波后的图像结果。13数字图像处理基于傅里叶基于傅里叶变换域的低通域的低通滤波波低通滤波，就是将高频信息剔除，只保留低频信息，对于空域图像结果而言就是图像变得模糊（或者说平滑）了，因为抑制了图像的边缘信息；而对于频域中间结果图像而言，就是图像四周的信息被抑制了，看起来就是中间低频信息

4、所占区域变小了。14数字图像处理傅里叶正傅里叶正变换和逆和逆变换的的总源源码框架框架傅里叶变换是基于复数空间的，所以需要一个基本的数据结构来表示复数，为typedef struct double m_re;double m_im;complex;，其中m_re为实部，m_im为虚部。int CImageProcessing:forward_fft2d(complex*array,int rows,int cols)return(fft2d(array,rows,cols,FFT_FORWARD);int CImageProcessing:inverse_fft2d(complex*array,

5、int rows,int cols)return(fft2d(array,rows,cols,FFT_INVERSE);15数字图像处理从原始从原始图像向复数数像向复数数组转化化最直观的做法，就是将每一个像素的灰度值作为每一个复数的实部，实际上这个做法是完全正确的。但是，如何设置复数的虚部呢？至此，读者应该能想到，设置为0。这样，原始图像的总像素个数就是复数数组的维数，每个复数的实部就对应像素的灰度，虚部设为0。16数字图像处理CTArray CImageProcessing:Image_to_complex(const CTMatrix&gray_image,long&new_height,

6、long&new_width)long image_height=gray_image.Get_height();long image_width =gray_image.Get_width();new_height=image_height;new_width=image_width;if(!(power_of_2(image_height)new_height=long(pow(2.0f,fastlog2(image_height);if(!(power_of_2(image_width)new_width=long(pow(2.0f,fastlog2(image_width);CTArr

7、ay array_of_complex(new_height*new_width);for(int row=0;row new_height;row+)for(int column=0;column new_width;column+)array_of_complex row*new_width+column.m_re=gray_image row column;array_of_complex row*new_width+column.m_im=0;return array_of_complex;17数字图像处理从复数数从复数数组得到最得到最终的的结果空域果空域图像像对于最后的复数数组向空域

8、图像转换的问题，只要取复数的实部作为图像像素的灰度值就可以了。CTMatrix CImageProcessing:Complex_to_image(const CTArray&complex_array,long image_height,long image_width)CTMatrix gray_image(image_height,image_width);ASSERT(complex_array.GetDimension()=image_height*image_width);for(int row=0;row image_height;row+)for(int column=0;co

9、lumn image_width;column+)gray_image row column =BYTE(complex_array row*image_width+column.m_re);return gray_image;18数字图像处理将将频域中域中间结果果转化化为显示示图像像对于一般的频率信息的显示，则需要取复数的幅度作为显示图像的像素值，并且考虑到灰度像素的值域（0-255）问题，需要除以一个常量（如100）来使得人眼能够看到频率的变化信息。CTMatrix CImageProcessing:Display_frequency(const CTArray&complex_array

10、,long image_height,long image_width)CTMatrix gray_image(image_height,image_width);ASSERT(complex_array.GetDimension()=image_height*image_width);for(int row=0;row image_height;row+)for(int column=0;column image_width;column+)gray_image row column =BYTE(sqrt(complex_array row*image_width+column.m_re*c

11、omplex_array row*image_width+column.m_re+complex_array row*image_width+column.m_im*complex_array row*image_width+column.m_im)/100.0);return gray_image;19数字图像处理CTArray CImageProcessing:Low_pass_filter(CTArray original_signal)long dimension=original_signal.GetDimension();double threshold=0;for(int ind

12、ex=0;index threshold)threshold=magnitude;threshold/=100;for(int index=0;index dimension;index+)double magnitude=sqrt(original_signal index.m_re*original_signal index.m_re +original_signal index.m_im*original_signal index.m_im);double eplon=1.0/sqrt(1+(threshold/magnitude)*(threshold/magnitude);origi

13、nal_signal index.m_re*=eplon;original_signal index.m_im*=eplon;return original_signal;20数字图像处理基于傅里叶基于傅里叶变换域的低通域的低通滤波框波框架架原始图像到复数数组的转化傅里叶正变换低通滤波傅里叶逆变换复数数组到结果图像的转化21数字图像处理基于傅里叶基于傅里叶变换域的低通域的低通滤波框波框架架/1）原始图像到复数数组的转化CTArray array_of_complex=CImageProcessing:Image_to_complex(m_image.Get_gray_image(),new_h

14、eight,new_width);/2）傅里叶正变换complex*pointer_of_complex=new complex array_of_complex.GetDimension();array_of_complex.ExportTo(pointer_of_complex);CImageProcessing:forward_fft2d(pointer_of_complex,new_height,new_width);array_of_complex.ImportFrom(pointer_of_complex);/3）低通滤波array_of_complex=CImageProcess

15、ing:Low_pass_filter(array_of_complex);/4）傅里叶逆变换array_of_complex.ExportTo(pointer_of_complex);CImageProcessing:inverse_fft2d(pointer_of_complex,new_height,new_width);array_of_complex.ImportFrom(pointer_of_complex);delete pointer_of_complex;/5）复数数组到结果图像的转化CTMatrix filter_image=CImageProcessing:Complex

16、_to_image(array_of_complex,new_height,new_width);22数字图像处理高斯高斯滤波法波法23数字图像处理高斯高斯滤波法波法24数字图像处理高斯高斯滤波法波法25数字图像处理高斯高斯滤波波结果示例果示例26数字图像处理均均值滤波法波法均值滤波法，就是取一个图像区域中的平均值作为滤波后的像素值，其非常类似于高斯滤波，也是属于基于空域的模板卷积的方法。均值滤波所对应的卷积模板中的各个权重都是一样的，在实际计算中，也根本不需要进行模板卷积操作，只需要进行各个像素累加取平均就行了。均值滤波，是属于线性的窗口滤波，是典型的低通滤波方法，其基本的思想就是对任何一个

17、元素取其邻域中的平均值作为滤波响应。27数字图像处理均均值滤波法波法英语单词发音均值滤波模板示意图28数字图像处理均均值滤波法波法一维均值滤波实现示意图均值滤波结果示例29数字图像处理中中值滤波法波法中值滤波，是一种非线性的窗口滤波方法，可以在很容易去除噪声的基础上保留边缘信息。中值滤波的基本思想就是，对于图像中的任意一个像素，从其邻域中的所有像素中挑选与其它像素最为相似的像素。30数字图像处理混合中混合中值滤波波混合中值滤波，也是一种非线性的滤波算法，相对简单的中值滤波其具有很好的边缘保留特性。混合中值滤波的基本思想是，对图像中的任意一个像素运用多次简单中值滤波，每次中值滤波所使用的邻域形状

18、会发生变化，然后取多次中值结果中的中值作为最后的滤波响应。具体来说，可以先运用十字模板进行简单中值滤波，然后再运用叉形模板进行滤波，最后取上述两个中值与原始像素本身进行比较，并提取它们之间的中值作为最后的结果。31数字图像处理混合中混合中值滤波波32数字图像处理33数字图像处理34数字图像处理中中值滤波波结果示例果示例35数字图像处理图像像锐化化图像锐化，顾名思义，就是让图像看起来更尖锐一点，更清晰一点。图像锐化的目的，则是让图像的内部对比度增大，使得用户能够看起来更加清晰，其存在一个副作用，就是可能会添加用户不感兴趣的信息，即噪声。从噪声的角度来说，图像平滑是起抑制作用的，而图像锐化是起促进

19、作用的；从对比度的角度来说，图像平滑是起缩小作用的，而图像锐化是起放大作用的。36数字图像处理基于傅里叶基于傅里叶变换域的高通域的高通滤波框波框架架基本思路就是，首先对原始图像进行傅里叶正变换，得到原始傅里叶变换域；然后对原始变换域进行高通滤波操作，滤波后的直接结果还是变换域；最后对滤波后的变换域再进行傅里叶逆变换，得到滤波处理后的图像，即为最终结果。37数字图像处理巴特沃斯低通与高通巴特沃斯低通与高通滤波器示意波器示意图38数字图像处理基于巴特沃斯基于巴特沃斯滤波器的低通波器的低通滤波波与高通与高通滤波波低通滤波高通滤波39数字图像处理基于基于圆形形滤波器的高通波器的高通滤波波CTArray

20、 CImageProcessing:High_pass_filter(CTArray original_signal,long image_height,long image_width)long dimension=original_signal.GetDimension();long radius=min(image_height,image_width)/7;for(int row=0;row image_height;row+)for(int column=0;column image_width;column+)if(row-image_height/2)*(row-image_he

21、ight/2)+(column-image_width/2)*(column-image_width/2)radius*radius)original_signal row*image_width+column.m_re=0;original_signal row*image_width+column.m_im=0;return original_signal;40数字图像处理基于基于圆形形滤波器的低通波器的低通滤波与高波与高通通滤波波41数字图像处理基于基于圆形形滤波器的高通波器的高通滤波波结果果示例示例42数字图像处理差差值滤波法波法与像素邻域的平均值相反的概念，就是像素邻域内的差值；也就

22、是说，可以取像素邻域内的差值来作为像素的最后响应值，以实现基于空域的图像锐化，此方法就叫作差值滤波法。43数字图像处理各种各种图像像锐化化结果的果的对比比原始图像巴特沃斯滤波器圆形滤波器差值滤波器44数字图像处理本章小本章小结图像平滑与锐化，是两个截然相对的概念。图像平滑，就是让图像变得模糊；图像锐化，就是让图像变得尖锐。从频域角度来说，图像平滑就对应着低通滤波法（保留低频信息），而图像锐化则对应着高通滤波法（保留高频信息）。对于低通滤波法，本章阐述了基于傅里叶变换域的低通滤波，并介绍了快速傅里叶变换原理。从空域角度来说，图像平滑使用了三种滤波方法，分别为高斯滤波法、均值滤波法和中值滤波法，而图像锐化则使用了差值滤波法。45数字图像处理此此课件下件下载可自行可自行编辑修改，修改，仅供参考！供参考！感感谢您的支持，我您的支持，我们努力做得更好！努力做得更好！谢谢46数字图像处理