快速排序中值滤波器.docVIP

快速排序中值滤波器实验报告 学生姓名： 黄永情 学 号： 6103314018 专业班级： 数字媒体技术141 实验项目名称 快速排序（非线性）中值滤波器 实验目的 作为期中实验测试，更好的巩固前面所学知识，自己动手动脑。掌握中值滤波器的构建方法，对图像进行噪声处理，并会中值滤波器对图像进行噪声去除。 实验背景知识 中值滤波是对一个滑动窗口内的诸像素灰度值排序，用其中值代替窗口中心象素的原来灰度值，它是一种非线性的图像平滑法，它对脉冲干扰级椒盐噪声的抑制效果好，在抑制随机噪声的同时能有效保护边缘少受模糊。 实验原理 对于一幅图像的像素矩阵，取以目标象素为中心的一个子矩阵窗口，这个窗口可以是3*3 ，5*5 等根据需要选取，本实验采用3*3（可改med变量值为5*5 7*7），对窗口内的象素灰度排序，取中间一个值作为目标象素的新灰度值。窗口示例如ooooxoooo上面x为目标象素，和周围o组成3*3矩阵pixelM，然后对这9个元素的灰度进行快速排序，以排序后的中间元素pixelM[4]为x的新灰度值，如此就完成对象素x的中值滤波，再迭代对其他需要的象素进行滤波即可。 实验步骤 先对图像进行椒盐噪声处理（同类还有高斯噪声） 中值滤波器处理噪声污染的图像，分别改变信噪比，对比图像处理效果。 通过从图像中的某个采样窗口取出奇数个数据进行排序 用排序后的中值取代要处理的数据即可 实验数据及处理结果 椒盐噪声： void CImageProcessingDoc::OnPepperNoisePollution() { // TODO: 在此添加命令处理程序代码 if (m_sImage.IsNull()) { AfxMessageBox(_T(请打开要处理的图像！)); return; } CDialogNoiseParameters mdlg;//对话框输入m_pa m_pb参数，调节信噪比 if (mdlg.DoModal() != IDOK) return; int x; int y; int sWidth = m_sImage.GetWidth(); int sHeight = m_sImage.GetHeight(); int sRowBytes = m_sImage.GetPitch(); BYTE *sImageData = (BYTE *)m_sImage.GetBits(); //生成伪随机种子 srand((unsigned)time(NULL)); //在图像中加椒盐噪声 for (y = 0; y sHeight; y++) for (x = 0; x sWidth; x++) { int srd = rand() % 100; if (srd ((int)(100 * (1 - mdlg.m_pb)))) { //图像中当前点置为黑(噪声） sImageData[y*sRowBytes + x * 3] = 0; sImageData[y*sRowBytes + x * 3 + 1] = 0; sImageData[y*sRowBytes + x * 3 + 2] = 0; } else if (srd ((int)(100 * (mdlg.m_pa)))) { sImageData[y*sRowBytes + x * 3] = 255; sImageData[y*sRowBytes + x * 3 + 1] = 255; sImageData[y*sRowBytes + x * 3 + 2] = 255; } } UpdateAllViews(NULL); } 处理图像效果： 对话框参数输入： 中值滤波器： void CImageProcessingDoc::OnMidFiltering() { // TODO: 在此添加命令处理程序代码 if (m_sImage.IsNull()) { AfxMessageBox(_T(请打开要处理的图像！)); return; } int x, y, My, Mx, x1, y1;//Mx,My为矩阵最大值，x1,y1为矩阵第一个访问点 int med=3;//矩阵大小，这里设为3*3，可改5,7,9 int Num = 0;//统计的矩阵中个数 int Numx[100], Numy[100];//保存每个像素X