快速中值滤波及c语言实现

《快速中值滤波及c语言实现》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、快速中值滤波及c语言实现学生姓名：刘勇学号：6100410218专业班级：数媒101【摘要】本文讨论了用c语言在微机上实现中值滤波及快速算法，在程序设计的过程中充分考虑到程序运行的时间复杂度和空间复杂度的问题．解决了由于图像太大而内存不够的问题，运用对程序运行时的方法，得出在PENTIUM-S100MHz上中值滤渡的一般算法运行4．23秒．而快速算法运行258秒。【关键词】c语言；中值滤波；快速算法1引言中值滤波是涂基发明的一种非线性信号处理技术，对抑制图像的噪声非常有效，在二维形式下，中值滤渡

2、器是一个古有奇数个像素的滑动窗口，窗口正中的象素的灰度值用窗口内各个象素的中值代替窗口的中值为窗口中象素按大小顺序排列后处于中间位置的象素；本文讨论中值滤的一般算法并比较其运算速度。2用C语言实现算法的若干问题在设计算法编制程序的时候，我们充分考虑到程序运行的时间复杂度和空间复杂度问题，在解决问题的前提下，使算法尽量简单，使程序运行占有的空间尽量的小，这样来减少不必要的时问浪费和空间浪费，从而太大的提高程序执行的效率。首先考虑到的内存问题。由于在本文算法中用的图像是512+5128bit，这就存

3、在一个内存不够大一整幅图像不能一次性调入的问题。为了解受此问题，可以只开辟一个3"512的缓冲区n，将原图像采用分批调入缓冲区，使内存不够的问题得到了圆满的解决。另外为了对中值滤波的快速算法和普通算法进行精确的比较，采用对程序运行计时的方法，并精确计算每个算法运行的时间，使得出的结论更可靠。3中值滤波算法的C语言程序实现本算法采用对开辟的3*512的缓冲区从左到右依次形成一个3*3的窗口．然后将此3*3的窗口放人一个一维数组中，调用求中值子函数．通过排序得出中值，当此中值不等于窗口中间位置的象素

4、时．用此中值来代替窗VI中间位置的象素灰度值．若此缓冲区处理完毕后，将缓冲区的第一行存入新建的文件中，将第二、第三行分别向上移动一行，若存人新建的文件中的行数小于或等于511(即这样处理的行数小于或等于511)，则从原文件中调入一行作为缓冲区第三行，按上述方法进行直到处理的总行数等于511为止，最后，将缓冲区的第二、三行存人新建的文件，程序流程框图如图14中值滤波快速算法的C语言程序实现本算法充分利用了上一次处理的结果．采用迭代，逐次逼近的方法得到本次的中值，在一行处理完毕后转人下一行也采用走S

5、型的方法．这样除第一个窗口采用了一伏排序得到中值外，其它的窗口都利用上伏的窗口的象素删除无用的3个象素后再加人新的3个象素，利用迭代的方法得到本次窗口的中值．这样太大地提高了程序执行的效率。4．1算法的解释首先是开辟一个3*512的缓冲区a，在初始化缓冲区时考虑到时间复杂度的问题，所以只初始化了第二、三行，而对第一行只初始化了前三个象素，这样便在缓冲区中可以得到一个3*3的窗口，对此窗口进行排序求中值后得出第一个窗口的中间象素的值．将文件指针定位在2*512处。然后开始循环，当处理的行数小于或等

6、于511时，将缓冲区a中的第二、三行分别向上移动一行变为第一、二行，从文件中读人512个字节作为缓冲区的第三行，并用行数模2的方法设置方向标志k．当k为0时，从左向右移动窗口．当k为1时．从右向左移动窗口．而每一窗口都利用上次的窗口的像素删除无用的3个象素后再加入新的3个象素．利用迭代的方法从上次的中值得到本次的中值。当处理完一行后将缓冲区的第一行存入新建的文件中，最后将缓冲区的第二、三行存入文件中。4．2算法代码//ImageProcessingDoc.cpp:implementationof

7、theCImageProcessingDocclass//#include"stdafx.h"#include"ImageProcessing.h"#include"ImageProcessingDoc.h"#include"GreyRatio.h"#include#definePI(acos(0.0)*2)#ifdef_DEBUG#definenewDEBUG_NEW#undefTHIS_FILEstaticcharTHIS_FILE[]=__FILE__;#endif////

8、///////////////////////////////////////////////////////////////////////////CImageProcessingDocIMPLEMENT_DYNCREATE(CImageProcessingDoc,CDocument)BEGIN_MESSAGE_MAP(CImageProcessingDoc,CDocument)//{{AFX_MSG_MAP(CImageProcessingDoc)ON_COMMAND(ID_HISTOGRA