中值滤波实现

最新推荐文章于 2024-07-25 15:39:09 发布
最新推荐文章于 2024-07-25 15:39:09 发布
阅读量1.3w 收藏 23
点赞数 3
分类专栏： 图像处理 文章标签： 中值滤波

转自http://blog.csdn.net/foreverhuylee/article/details/22732003

1、什么是中值滤波？

中值滤波是对一个滑动窗口内的诸像素灰度值排序，用其中值代替窗口中心象素的原来灰度值，它是一种非线性的图像平滑法，它对脉冲干扰级椒盐噪声的抑制效果好，在抑制随机噪声的同时能有效保护边缘少受模糊。

中值滤波可以过滤尖峰脉冲。目的在于我们对于滤波后的数据更感兴趣。滤波后的数据保留的原图像的变化趋势，同时去除了尖峰脉冲对分析造成的影响。

    以一维信号的中值滤波举例。对灰度序列80、120、90、200、100、110、70，如果按大小顺序排列，其结果为70、80、90、10O、110、120、200，其中间位置上的灰度值为10O，则该灰度序列的中值即为100。一维信号中值滤波实际上就是用中值代替规定位置（一般指原始信号序列中心位置）的信号值。对前面所举的序列而言，中值滤波的结果是用中值100替代序列80、120、90、200、100、110、70中的信号序列中心位置值200，得到的滤波序列就是80、120、90、100、100、110、70。如果在此序列中200是一个噪声信号，则用此方法即可去除这个噪声点。

    二维中值滤波算法是：对于一幅图像的象素矩阵，取以目标象素为中心的一个子矩阵窗口，这个窗口可以是3*3 ，5*5 等根据需要选取，对窗口内的象素灰度排序，取中间一个值作为目标象素的新灰度值。窗口示例如ooooxoooo上面x为目标象素，和周围o组成3*3矩阵Array，然后对这9个元素的灰度进行排序，以排序后的中间元素Array[4]为x的新灰度值，如此就完成对象素x的中值滤波，再迭代对其他需要的象素进行滤波即可。

图像处理中，中值滤波的实现方法

1：通过从 图像中的某个 采样窗口取出奇数个数据进行排序
2： 用排序后的中值取代要处理的数据即可
中值滤波的算法实现过程，重点是排序，最常用的冒泡排序~~
把滤波区间的数据从小到大进行排序，然后取中值，（如果是奇数个数据，那么中值就只有一个了，如果偶数个数据，中值有两个，可以对两个数据再求平均）

下面是一个C语言实现中值滤波的函数：

  1. unsigned char GetMedianNum(int * bArray, int iFilterLen) 
  2. { 
  3.     int i,j;// 循环变量 
  4.     unsigned char bTemp; 
  5.      
  6.     // 用冒泡法对数组进行排序 
  7.     for (j = 0; j < iFilterLen - 1; j ++) 
  8.     { 
  9.         for (i = 0; i < iFilterLen - j - 1; i ++) 
  10.         { 
  11.             if (bArray[i] > bArray[i + 1]) 
  12.             { 
  13.                 // 互换 
  14.                 bTemp = bArray[i]; 
  15.                 bArray[i] = bArray[i + 1]; 
  16.                 bArray[i + 1] = bTemp; 
  17.             } 
  18.         } 
  19.     } 
  20.      
  21.     // 计算中值 
  22.     if ((iFilterLen & 1) > 0) 
  23.     { 
  24.         // 数组有奇数个元素，返回中间一个元素 
  25.         bTemp = bArray[(iFilterLen + 1) / 2]; 
  26.     } 
  27.     else 
  28.     { 
  29.         // 数组有偶数个元素，返回中间两个元素平均值 
  30.         bTemp = (bArray[iFilterLen / 2] + bArray[iFilterLen / 2 + 1]) / 2; 
  31.     } 
  32.  
  33.     return bTemp; 
  34. }  
unsigned char GetMedianNum(int * bArray, int iFilterLen)
{
	int i,j;// 循环变量
	unsigned char bTemp;
	
	// 用冒泡法对数组进行排序
	for (j = 0; j < iFilterLen - 1; j ++)
	{
		for (i = 0; i < iFilterLen - j - 1; i ++)
		{
			if (bArray[i] > bArray[i + 1])
			{
				// 互换
				bTemp = bArray[i];
				bArray[i] = bArray[i + 1];
				bArray[i + 1] = bTemp;
			}
		}
	}
	
	// 计算中值
	if ((iFilterLen & 1) > 0)
	{
		// 数组有奇数个元素，返回中间一个元素
		bTemp = bArray[(iFilterLen + 1) / 2];
	}
	else
	{
		// 数组有偶数个元素，返回中间两个元素平均值
		bTemp = (bArray[iFilterLen / 2] + bArray[iFilterLen / 2 + 1]) / 2;
	}

	return bTemp;
}
[cpp] view plain copy print ? 在CODE上查看代码片 派生到我的代码片
  1. 注：bArray 是一个整形指针，我们传入的一般是一个数组，用来存储待排序的数据 
  2.     iFilterLen 是滤波器的长度 
  3.    用在图像处理中时，由于像素的取值范围是0~255，刚好是unsigned char 的范围，所以函数的返回值是unsigned char，如果我们要处理的数是float型，或其他类型，返回值也可以更改~~返回值是bTemp，也即是我们想得到的中值  
注：bArray 是一个整形指针，我们传入的一般是一个数组，用来存储待排序的数据
    iFilterLen 是滤波器的长度
   用在图像处理中时，由于像素的取值范围是0~255，刚好是unsigned char 的范围，所以函数的返回值是unsigned char，如果我们要处理的数是float型，或其他类型，返回值也可以更改~~返回值是bTemp，也即是我们想得到的中值
[cpp] view plain copy print ? 在CODE上查看代码片 派生到我的代码片
  1. <span style="color: rgb(51, 51, 51);">下面是一个完整的C语言程序，用在图像处理中</span>  
<span style="color:#333333;">下面是一个完整的C语言程序，用在图像处理中</span>
[cpp] view plain copy print ? 在CODE上查看代码片 派生到我的代码片
  1. /*************************************************************************
  2. * 函数名称：
  3. *   MedianFilter()
  4. * 参数:
  5. *   int   iFilterH         - 滤波器的高度
  6. *   int   iFilterW         - 滤波器的宽度
  7. *   int   iFilterMX        - 滤波器的中心元素X坐标
  8. *   int   iFilterMY        - 滤波器的中心元素Y坐标
  9. * 说明:
  10. *   该函数对DIB图像进行中值滤波。
  11. ************************************************************************/ 
  12. #define iFilterW 1 
  13. #define iFilterH 1 
  14. #define iFilterMX 1 
  15. #define iFilterMY 1 
  16. #define WIDTHBYTES(bits)    (((bits) + 31) / 32 * 4) 
  17.  
  18.  
  19. unsigned char GetMedianNum(int * bArray, int iFilterLen); 
  20. void MedianFilter(unsigned char *pImg1,unsigned char *pImg,int nWidth,int nHeight) 
  21. {        
  22.     unsigned char   *lpSrc;                         // 指向源图像的指针  
  23.     unsigned char   *lpDst;                         // 指向要复制区域的指针 
  24.     int         aValue[iFilterH*iFilterW];          // 指向滤波器数组的指针 
  25.     int         i,j,k,l;                            // 循环变量  
  26.     int         lLineBytes;                         // 图像每行的字节数  
  27.     lLineBytes = WIDTHBYTES(nWidth * 8); 
  28.     for ( i=0;i<nWidth;i++,pImg++ ) 
  29.         (*pImg)=0; 
  30.     // 开始中值滤波 
  31.     // 行(除去边缘几行) 
  32.     for(i = iFilterMY; i < nHeight - iFilterH + iFilterMY + 1; i++) 
  33.     { 
  34.         // 列(除去边缘几列) 
  35.         for(j = iFilterMX; j < nWidth - iFilterW + iFilterMX + 1; j++) 
  36.         { 
  37.             // 指向新DIB第i行，第j个象素的指针 
  38.             lpDst = pImg + lLineBytes * (nHeight - 1 - i) + j; 
  39.              
  40.             // 读取滤波器数组 
  41.             for (k = 0; k < iFilterH; k++) 
  42.             { 
  43.                 for (l = 0; l < iFilterW; l++) 
  44.                 { 
  45.                     // 指向DIB第i - iFilterMY + k行，第j - iFilterMX + l个象素的指针 
  46.                     lpSrc = pImg1 + lLineBytes * (nHeight - 1 - i + iFilterMY - k) + j - iFilterMX + l; 
  47.                  
  48.                     // 保存象素值 
  49.                     aValue[k * iFilterW + l] = *lpSrc; 
  50.                 } 
  51.             } 
  52.              
  53.             // 获取中值 
  54.             * lpDst = GetMedianNum(aValue, iFilterH * iFilterW); 
  55.         } 
  56.     } 
  57.  
  58. } 
  59.  
  60. unsigned char GetMedianNum(int * bArray, int iFilterLen) 
  61. { 
  62.     int     i,j;            // 循环变量 
  63.     unsigned char bTemp; 
  64.      
  65.     // 用冒泡法对数组进行排序 
  66.     for (j = 0; j < iFilterLen - 1; j ++) 
  67.     { 
  68.         for (i = 0; i < iFilterLen - j - 1; i ++) 
  69.         { 
  70.             if (bArray[i] > bArray[i + 1]) 
  71.             { 
  72.                 // 互换 
  73.                 bTemp = bArray[i]; 
  74.                 bArray[i] = bArray[i + 1]; 
  75.                 bArray[i + 1] = bTemp; 
  76.             } 
  77.         } 
  78.     } 
  79.      
  80.     // 计算中值 
  81.     if ((iFilterLen & 1) > 0) 
  82.     { 
  83.         // 数组有奇数个元素，返回中间一个元素 
  84.         bTemp = bArray[(iFilterLen + 1) / 2]; 
  85.     } 
  86.     else 
  87.     { 
  88.         // 数组有偶数个元素，返回中间两个元素平均值 
  89.         bTemp = (bArray[iFilterLen / 2] + bArray[iFilterLen / 2 + 1]) / 2; 
  90.     } 
  91.      
  92.     return bTemp; 
  93. }  
/*************************************************************************
 * 函数名称：
 *   MedianFilter()
 * 参数:
 *   int   iFilterH			- 滤波器的高度
 *   int   iFilterW			- 滤波器的宽度
 *   int   iFilterMX		- 滤波器的中心元素X坐标
 *   int   iFilterMY		- 滤波器的中心元素Y坐标
 * 说明:
 *   该函数对DIB图像进行中值滤波。
 ************************************************************************/
#define iFilterW 1
#define iFilterH 1
#define iFilterMX 1
#define iFilterMY 1
#define WIDTHBYTES(bits)    (((bits) + 31) / 32 * 4)


unsigned char GetMedianNum(int * bArray, int iFilterLen);
void MedianFilter(unsigned char *pImg1,unsigned char *pImg,int nWidth,int nHeight)
{		
    unsigned char   *lpSrc;			                // 指向源图像的指针	
	unsigned char   *lpDst;		                 	// 指向要复制区域的指针
	int         aValue[iFilterH*iFilterW];		    // 指向滤波器数组的指针
	int			i,j,k,l;		                    // 循环变量	
	int			lLineBytes;		                    // 图像每行的字节数	
	lLineBytes = WIDTHBYTES(nWidth * 8);
	for ( i=0;i<nWidth;i++,pImg++ )
		(*pImg)=0;
	// 开始中值滤波
	// 行(除去边缘几行)
	for(i = iFilterMY; i < nHeight - iFilterH + iFilterMY + 1; i++)
	{
		// 列(除去边缘几列)
		for(j = iFilterMX; j < nWidth - iFilterW + iFilterMX + 1; j++)
		{
			// 指向新DIB第i行，第j个象素的指针
			lpDst = pImg + lLineBytes * (nHeight - 1 - i) + j;
			
			// 读取滤波器数组
			for (k = 0; k < iFilterH; k++)
			{
				for (l = 0; l < iFilterW; l++)
				{
					// 指向DIB第i - iFilterMY + k行，第j - iFilterMX + l个象素的指针
					lpSrc = pImg1 + lLineBytes * (nHeight - 1 - i + iFilterMY - k) + j - iFilterMX + l;
				
					// 保存象素值
					aValue[k * iFilterW + l] = *lpSrc;
				}
			}
			
			// 获取中值
			* lpDst = GetMedianNum(aValue, iFilterH * iFilterW);
		}
	}

}

unsigned char GetMedianNum(int * bArray, int iFilterLen)
{
	int		i,j;			// 循环变量
	unsigned char bTemp;
	
	// 用冒泡法对数组进行排序
	for (j = 0; j < iFilterLen - 1; j ++)
	{
		for (i = 0; i < iFilterLen - j - 1; i ++)
		{
			if (bArray[i] > bArray[i + 1])
			{
				// 互换
				bTemp = bArray[i];
				bArray[i] = bArray[i + 1];
				bArray[i + 1] = bTemp;
			}
		}
	}
	
	// 计算中值
	if ((iFilterLen & 1) > 0)
	{
		// 数组有奇数个元素，返回中间一个元素
		bTemp = bArray[(iFilterLen + 1) / 2];
	}
	else
	{
		// 数组有偶数个元素，返回中间两个元素平均值
		bTemp = (bArray[iFilterLen / 2] + bArray[iFilterLen / 2 + 1]) / 2;
	}
	
	return bTemp;
}
转自 http://blog.csdn.net/computer_liuyun/article/details/16822341

数字图像在其形成、传输记录的过程中往往会受到很多噪声的的污染，比如：椒盐噪声、高斯噪声等，为了抑制和消除这些随即产生的噪声而改善图像的质量，就需要去、对图像进行去噪处理，去噪也就是滤波处理。
中值滤波和同态滤波都是图像增强的方法，但是这两种方法是从不同的方式进行改善一副图片的质量。
中值滤波是图像平滑的一种方法 它是一种非线性平滑滤波技术，在一定条件下可以克服线性滤波带来的图像细节的模糊问题，特别是针对被椒盐噪声污染的图像。
所谓中值滤波，其中滤波就是前面讲的去噪，关键在于中值两字，中值从字面意思上讲就是中间的那个值也就是中心值。

举一个简单的例子：一维序列{0,3,4,0,7}，进行中值滤波排序后为{0,0,3,4,7}，则其中值为3。

数字图像是以二维图像来描述的，故对图像的滤波也就是对二维数据序列的滤波，这个二维序列相当于一个二维矩阵，里面元素的值就是每个像素点的像素。

中值滤波通常采用一个含奇数个点的滑动窗口，用窗口的中的灰度值的中值来代替中心点的灰度值，其实就是对这个窗口中的灰度值进行排序，然后将其中值赋值给中心点即可。常用的中值滤波窗口形状有线状、方形、圆形以及十字形等。

注：对每一个像素的m*n邻域进行计算，中值滤波对图像的边界用0做扩张，所以对边界可能会出现扭曲。

算法描述：

[1] 获得源图像的首地址及图像的宽和高

[2] 开辟一块内存缓冲区，用以暂存结果图像，并初始化为0

[3] 逐个扫描图像中的像素点，将其邻域各元素的像素值从小到大进行排序，将求得到的中间值赋值给目标图像中与当前点对应的像素点

[4] 循环步骤[3]，直到处理完源图像的全部像素点

[5] 将结果从内存缓冲区复制到源图像的数据区

PS: 个人觉得中值滤波原理非常简单，在做的过程中原理没什么难度，一直让我纠结的是怎么去读取一副图片的数据区(也就是这副图片的像素值)，对于我这个没接触过图像的水货来说困扰了我很久，当然利用matlab软件读取是很容易的事情，中值滤波也就其内嵌写好的函数调用就好了，但是为了锻炼一下自己，尝试用c写。

Bmp图像的存储、格式以及读取

现在讲解BMP的4个组成部分：

1.文件头信息块

0000-0001：文件标识，为字母ASCII码“BM”。
0002-0005：文件大小。
0006-0009：保留，每字节以“00”填写。
000A-000D：记录图像数据区的起始位置。各字节的信息依次含义为：文件头信息块大小，图像描述信息块的大小，图像颜色表的大小，保留（为01）。


2.图像描述信息块

000E-0011：图像描述信息块的大小，常为28H。
0012-0015：图像宽度。
0016-0019：图像高度。
001A-001B：图像的plane总数（恒为1）。
001C-001D：记录像素的位数，很重要的数值，图像的颜色数由该值决定。
001E-0021：数据压缩方式（数值位0：不压缩；1：8位压缩；2：4位压缩）。
0022-0025：图像区数据的大小。
0026-0029：水平每米有多少像素，在设备无关位图（.DIB）中，每字节以00H填写。
002A-002D：垂直每米有多少像素，在设备无关位图（.DIB）中，每字节以00H填写。
002E-0031：此图像所用的颜色数，如值为0，表示所有颜色一样重要。

3.颜色表

　　颜色表的大小根据所使用的颜色模式而定：2色图像为8字节；16色图像位64字节；256色图像为1024字节。其中，每4字节表示一种颜色，并以B（蓝色）、G（绿色）、R（红色）、alpha（32位位图的透明度值，一般不需要）。即首先4字节表示颜色号1的颜色，接下来表示颜色号2的颜色，依此类推。

4.图像数据区

　　颜色表接下来位为位图文件的图像数据区，在此部分记录着每点像素对应的颜色号，其记录方式也随颜色模式而定，既2色图像每点占1位（8位为1字节）；16色图像每点占4位（半字节）；256色图像每点占8位（1字节）；真彩色图像每点占24位（3字节）。所以，整个数据区的大小也会随之变化。究其规律而言，可的出如下计算公式：图像数据信息大小=（图像宽度*图像高度*记录像素的位数）/8。
　　然而，未压缩的图像信息区的大小。除了真彩色模式外，其余的均大于或等于数据信息的大小。这是为什么呢？原因有两个：
　　1.BMP文件记录一行图像是以字节为单位的。因此，就不存在一个字节中的数据位信息表示的点在不同的两行中。也就是说，设显示模式位16色，在每个字节分配两个点信息时，如果图像的宽度位奇数，那么最后一个像素点的信息将独占一个字节，这个字节的后4位将没有意义。接下来的一个字节将开始记录下一行的信息。
　　2.为了显示的方便，除了真彩色外，其他的每中颜色模式的行字节数要用数据“00”补齐为4的整数倍。如果显示模式为16色，当图像宽为19时，存储时每行则要补充4-(19/2+1)%4=2个字节（加1是因为里面有一个像素点要独占了一字节）。如果显示模式为256色，当图像宽为19时，每行也要补充4-19%4=1个字节。

由于bmp图象是从下至上存储的，所以我们不能进行直接顺序读取。详细的说，bmp图象存储区数据是从1078偏移字节开始。文件内第一个图象点实际上是对应图象（320*200）第200行的最左边的第一个点，而从1078开始的320个点则是图象最下面一行对应的点，之后的321个点是图象倒数第二行最左边的第一个点。这样，bmp文件最后一个字节对应的点是第一行最后边的点了。

BMP图像的读取

1 文件格式

典型的BMP图像文件由四部分组成。

1>位图文件头数据结构，它包含BMP图像文件的类型、显示内容等信息；

2>位图信息数据结构，它包含有BMP图像的宽、高、压缩方法以及定义颜色等信息；

3>调色板，这个部分是可选的，有些位图需要调色板，有些位图，比如真彩色图（24位的BMP）就不需要调色板；

4>位图数据，这部分的内容根据BMP位图使用的位数不同而不同，在24位图中直接使用RGB，而其它小于24位的使用调色板中颜色所引值。也就是说 1,4,8位图(2色图像为8字节；16色图像位64字节；256色图像为1024字节。)。

2 位图的读取

1 //读取文件头  

2 fread(&bf,sizeof(bf),1,fi); 

3 //读取文件信息头  

4 fread(&bi,sizeof(bi),1,fi); 

      在默认字节对齐方式下，这样读取BMP图像会出现错误。由于系统默认情况下，结构体以4字节对齐，这样，对于结构体中不够4个字节的字段，会自动填充，使其也占据4个字节。对于像BMP图像这样采用1字节保存的文件，读取就会出错。

      对于这个问题，解决办法：

I> 在读取之前使用

     #pragma pack(1)

     这条语句，将字节对齐方式设置为1字节对齐即可；

这样即可

如果以上的关于BMP的知识理解了，那么就不难实现中值滤波了。

代码如下：

[cpp] view plain copy print ?
  1. /*
  2.     filter.c
  3. */ 
  4. #include<stdio.h> 
  5. #include<stdlib.h> 
  6. #include <string.h> 
  7. typedef unsigned char BYTE;  
  8. //颜色表定义 
  9. typedef struct tagRGBQUAD { 
  10. BYTE rgbBlue;// 蓝色的亮度（值范围为0-255) 
  11. BYTE rgbGreen; // 绿色的亮度（值范围为0-255) 
  12. BYTE rgbRed; // 红色的亮度（值范围为0-255) 
  13. BYTE rgbReserved;// 保留，必须为0 
  14. } RGBQUAD; 
  15. /*
  16.     定义头文件型
  17. */ 
  18. #pragma pack(1)//将字节对齐方式设置为1字节对齐 
  19. typedef struct 
  20. { 
  21. unsigned char id1;//位图文件的类型，必须为BM(占用0-1字节)  
  22. unsigned char id2; 
  23. unsigned int filesize;//位图文件的大小，以字节为单位（2-5字节） 
  24. unsigned int reserved;// 位图文件保留字，必须为0(6-9字节) 
  25. unsigned int bitmapdataoffset;//位图数据的起始位置，以相对于位图(10-13字节) 
  26. unsigned int bitmapheadersize;//BMP头的大小，固定为40（14-17字节） 
  27. unsigned int width;//图片宽度；以像素为单位（18-21字节） 
  28. unsigned int height;//图片高度；以像素为单位（22-25字节） 
  29. unsigned short planes;//图片位面数，必须为1(26-27字节） 
  30. unsigned short bitperpixel;//每个像素所需的位数，每个像素所需的位数（28-29字节） 
  31.    //只能是以下几个数：1（双色），4(16色），8(256色）或24（真彩色）  灰度级 
  32. unsigned int compression;//是否压缩（30-33字节） 
  33. //只能是以下几个数：0（不压缩），1(BI_RLE8压缩类型），2(BI_RLE4压缩类型） 
  34. unsigned int bitmapdatasize;//位图的大小，以字节为单位（34-37字节） 
  35. unsigned int hresolution;//位图水平分辨率，每米像素数（38-41字节） 
  36. unsigned int vresolution;//位图垂直分辨率，每米像素数（42-45字节) 
  37. unsigned int colors;//位图实际使用的颜色表中的颜色数（46-49字节） 
  38. unsigned int importantcolors;//位图显示过程中重要的颜色数（50-53字节） 
  39. //unsigned int  bmiColors[1];//调色板；(54 - 57字节) 
  40. unsigned char palette[256][4];//调色板 占256*4=1024字节 
  41. }BMPheader;//总大小40+14+1024=1078字节 
  42. typedef struct 
  43. { 
  44. BMPheader* bmpheader ; 
  45. unsigned char* bitmapdata;//图片数据； 
  46. }BMPheaderfile; 
  47. /*
  48.     求文件长度的函数
  49. */ 
  50. long getfilesize(FILE *f) 
  51. { 
  52. long pos,len; 
  53. pos=ftell(f);//ftell函数用于得到文件指针当前位置相对于文件首的偏移字节数 
  54. fseek(f,0,SEEK_END);//fseek函数用于移动文件指针相对于SEEK_END的偏移量为0 
  55. len=ftell(f);//len就是文件的长度 
  56. fseek(f,pos,SEEK_SET);//将文件指针移动到原来的地方 
  57. return len; 
  58. } 
  59. /*
  60.     主函数
  61. */ 
  62. int main() 
  63. { 
  64. BMPheaderfile *output=(BMPheaderfile*)malloc(sizeof(BMPheaderfile));//定义一个输出指针 
  65. unsigned char *data=NULL; 
  66.     FILE *fpr,*fpw; 
  67. /*
  68.   打开文件
  69. */ 
  70. if((fpr=fopen("start.bmp","rb"))==NULL) 
  71. { 
  72. printf("cannot open this file"); 
  73. exit(0); 
  74. } 
  75. if((fpw=fopen("end1.bmp","wb"))==NULL) 
  76. { 
  77. printf("cannot wirte this file"); 
  78. exit(0); 
  79. } 
  80. long length=getfilesize(fpr); 
  81. printf("文件的长度为%ld\n",length); 
  82. printf("文件的头长度为%d\n",sizeof(BMPheader)); 
  83. data=(unsigned char*)malloc(length*sizeof(char));//分配空间 
  84. if(0==fread(data,1,length,fpr))//读文件，从fpr指向的文件中读出length到data所指的内存空间去 
  85. { 
  86. printf("read failed\n"); 
  87. exit(0); 
  88. } 
  89. fclose(fpr);//释放指针 
  90. output->bmpheader=(BMPheader*)malloc(sizeof(BMPheader)); 
  91. memcpy(output->bmpheader,data,sizeof(BMPheader));//从data中拷贝sizeof(BMPheader)大小到output->bmpheader 
  92. /*
  93. 打印出图像中头文件的信息
  94. */ 
  95.     int height=output->bmpheader->height; 
  96. int width=output->bmpheader->width; 
  97. printf("filesize is %d\n",output->bmpheader->filesize); 
  98. printf("该图像每个像素所需要的位数：%d\n",output->bmpheader->bitperpixel); 
  99. printf("height is %d\n",output->bmpheader->height); 
  100. printf("width is %d\n",output->bmpheader->width); 
  101. data=data+sizeof(BMPheader); 
  102. output->bitmapdata=data; 
  103.   /*
  104.    中值滤波算法(选择3×3的滑动窗口)
  105.   */ 
  106. unsigned char pixel[9]={0};//滑动窗口的像素值，初始为0 
  107. unsigned char mid;//中值 
  108. unsigned char temp;//中间变量 
  109. int flag; 
  110. int m,i,j,x,h,w,y; 
  111. for(j=1;j<height-1;j++) 
  112. { 
  113. for(i=1;i<width-1;i++) 
  114. { 
  115. //将3×3滑动窗口中的所有像素值放入pixel[m] 
  116. m=0; 
  117. for(y=j-1;y<=j+1;y++) 
  118. for(x=i-1;x<=i+1;x++) 
  119. { 
  120. pixel[m]=data[y*width+x]; 
  121. m=m+1; 
  122. } 
  123. //让一位数组pixel[9]进行降序排列 
  124. do  
  125. { 
  126. flag=0;//循环结束的标志 
  127. for(m=0;m<9;m++) 
  128. { 
  129. if(pixel[m]<pixel[m+1]) 
  130. { 
  131. temp=pixel[m]; 
  132. pixel[m]=pixel[m+1]; 
  133. pixel[m+1]=temp; 
  134. flag=1; 
  135. }//if 
  136. }//for 
  137. } while (flag==1); 
  138. mid=pixel[4]; 
  139. output->bitmapdata[width*j+i]=mid; 
  140. } 
  141. } 
  142. /*
  143. N=5的中值滤波线段状
  144. */ 
  145. /*
  146.     for(i=0;i<height;i++)
  147. {
  148. for(j=2;j<width-1;j++)
  149. {
  150. m=0;
  151. for(x=j-2;x<=j+2;x++)
  152. pixel[m++]=data[i*width+x];
  153. for(h=0;h<5;h++)
  154. for(w=h+1;w<5;w++)
  155. {
  156. if(pixel[h]>pixel[w])
  157. {
  158. temp=pixel[w];
  159. pixel[w]=pixel[h];
  160. pixel[h]=temp;
  161. }
  162. }
  163. output->bitmapdata[i*width+j]=pixel[2];
  164. }
  165. }
  166. */ 
  167. /*
  168. 保存图像文件
  169. */ 
  170. fseek(fpw,0,0);  //fseek(fpw,0,SEEK_SET) 
  171. fwrite(output->bmpheader,1,sizeof(BMPheader),fpw);//写入图像的头文件 
  172. fwrite(output->bitmapdata,1,length-sizeof(BMPheader),fpw);//写入图像的数据信息 
  173. fclose(fpw);//释放指针 
  174. return 0; 
  175. } 

数字图像在其形成、传输记录的过程中往往会受到很多噪声的的污染，比如：椒盐噪声、高斯噪声等，为了抑制和消除这些随即产生的噪声而改善图像的质量，就需要去、对图像进行去噪处理，去噪也就是滤波处理。
中值滤波和同态滤波都是图像增强的方法，但是这两种方法是从不同的方式进行改善一副图片的质量。
中值滤波是图像平滑的一种方法 它是一种非线性平滑滤波技术，在一定条件下可以克服线性滤波带来的图像细节的模糊问题，特别是针对被椒盐噪声污染的图像。
所谓中值滤波，其中滤波就是前面讲的去噪，关键在于中值两字，中值从字面意思上讲就是中间的那个值也就是中心值。

举一个简单的例子：一维序列{0,3,4,0,7}，进行中值滤波排序后为{0,0,3,4,7}，则其中值为3。

数字图像是以二维图像来描述的，故对图像的滤波也就是对二维数据序列的滤波，这个二维序列相当于一个二维矩阵，里面元素的值就是每个像素点的像素。

中值滤波通常采用一个含奇数个点的滑动窗口，用窗口的中的灰度值的中值来代替中心点的灰度值，其实就是对这个窗口中的灰度值进行排序，然后将其中值赋值给中心点即可。常用的中值滤波窗口形状有线状、方形、圆形以及十字形等。

注：对每一个像素的m*n邻域进行计算，中值滤波对图像的边界用0做扩张，所以对边界可能会出现扭曲。

算法描述：

[1] 获得源图像的首地址及图像的宽和高

[2] 开辟一块内存缓冲区，用以暂存结果图像，并初始化为0

[3] 逐个扫描图像中的像素点，将其邻域各元素的像素值从小到大进行排序，将求得到的中间值赋值给目标图像中与当前点对应的像素点

[4] 循环步骤[3]，直到处理完源图像的全部像素点

[5] 将结果从内存缓冲区复制到源图像的数据区

PS: 个人觉得中值滤波原理非常简单，在做的过程中原理没什么难度，一直让我纠结的是怎么去读取一副图片的数据区(也就是这副图片的像素值)，对于我这个没接触过图像的水货来说困扰了我很久，当然利用matlab软件读取是很容易的事情，中值滤波也就其内嵌写好的函数调用就好了，但是为了锻炼一下自己，尝试用c写。

Bmp图像的存储、格式以及读取

现在讲解BMP的4个组成部分：

1.文件头信息块

0000-0001：文件标识，为字母ASCII码“BM”。
0002-0005：文件大小。
0006-0009：保留，每字节以“00”填写。
000A-000D：记录图像数据区的起始位置。各字节的信息依次含义为：文件头信息块大小，图像描述信息块的大小，图像颜色表的大小，保留（为01）。


2.图像描述信息块

000E-0011：图像描述信息块的大小，常为28H。
0012-0015：图像宽度。
0016-0019：图像高度。
001A-001B：图像的plane总数（恒为1）。
001C-001D：记录像素的位数，很重要的数值，图像的颜色数由该值决定。
001E-0021：数据压缩方式（数值位0：不压缩；1：8位压缩；2：4位压缩）。
0022-0025：图像区数据的大小。
0026-0029：水平每米有多少像素，在设备无关位图（.DIB）中，每字节以00H填写。
002A-002D：垂直每米有多少像素，在设备无关位图（.DIB）中，每字节以00H填写。
002E-0031：此图像所用的颜色数，如值为0，表示所有颜色一样重要。

3.颜色表

　　颜色表的大小根据所使用的颜色模式而定：2色图像为8字节；16色图像位64字节；256色图像为1024字节。其中，每4字节表示一种颜色，并以B（蓝色）、G（绿色）、R（红色）、alpha（32位位图的透明度值，一般不需要）。即首先4字节表示颜色号1的颜色，接下来表示颜色号2的颜色，依此类推。

4.图像数据区

　　颜色表接下来位为位图文件的图像数据区，在此部分记录着每点像素对应的颜色号，其记录方式也随颜色模式而定，既2色图像每点占1位（8位为1字节）；16色图像每点占4位（半字节）；256色图像每点占8位（1字节）；真彩色图像每点占24位（3字节）。所以，整个数据区的大小也会随之变化。究其规律而言，可的出如下计算公式：图像数据信息大小=（图像宽度*图像高度*记录像素的位数）/8。
　　然而，未压缩的图像信息区的大小。除了真彩色模式外，其余的均大于或等于数据信息的大小。这是为什么呢？原因有两个：
　　1.BMP文件记录一行图像是以字节为单位的。因此，就不存在一个字节中的数据位信息表示的点在不同的两行中。也就是说，设显示模式位16色，在每个字节分配两个点信息时，如果图像的宽度位奇数，那么最后一个像素点的信息将独占一个字节，这个字节的后4位将没有意义。接下来的一个字节将开始记录下一行的信息。
　　2.为了显示的方便，除了真彩色外，其他的每中颜色模式的行字节数要用数据“00”补齐为4的整数倍。如果显示模式为16色，当图像宽为19时，存储时每行则要补充4-(19/2+1)%4=2个字节（加1是因为里面有一个像素点要独占了一字节）。如果显示模式为256色，当图像宽为19时，每行也要补充4-19%4=1个字节。

由于bmp图象是从下至上存储的，所以我们不能进行直接顺序读取。详细的说，bmp图象存储区数据是从1078偏移字节开始。文件内第一个图象点实际上是对应图象（320*200）第200行的最左边的第一个点，而从1078开始的320个点则是图象最下面一行对应的点，之后的321个点是图象倒数第二行最左边的第一个点。这样，bmp文件最后一个字节对应的点是第一行最后边的点了。

BMP图像的读取

1 文件格式

典型的BMP图像文件由四部分组成。

1>位图文件头数据结构，它包含BMP图像文件的类型、显示内容等信息；

2>位图信息数据结构，它包含有BMP图像的宽、高、压缩方法以及定义颜色等信息；

3>调色板，这个部分是可选的，有些位图需要调色板，有些位图，比如真彩色图（24位的BMP）就不需要调色板；

4>位图数据，这部分的内容根据BMP位图使用的位数不同而不同，在24位图中直接使用RGB，而其它小于24位的使用调色板中颜色所引值。也就是说 1,4,8位图(2色图像为8字节；16色图像位64字节；256色图像为1024字节。)。

2 位图的读取

1 //读取文件头  

2 fread(&bf,sizeof(bf),1,fi); 

3 //读取文件信息头  

4 fread(&bi,sizeof(bi),1,fi); 

      在默认字节对齐方式下，这样读取BMP图像会出现错误。由于系统默认情况下，结构体以4字节对齐，这样，对于结构体中不够4个字节的字段，会自动填充，使其也占据4个字节。对于像BMP图像这样采用1字节保存的文件，读取就会出错。

      对于这个问题，解决办法：

I> 在读取之前使用

     #pragma pack(1)

     这条语句，将字节对齐方式设置为1字节对齐即可；

这样即可

如果以上的关于BMP的知识理解了，那么就不难实现中值滤波了。

代码如下：

[cpp] view plain copy print ?
  1. /*
  2.     filter.c
  3. */ 
  4. #include<stdio.h> 
  5. #include<stdlib.h> 
  6. #include <string.h> 
  7. typedef unsigned char BYTE;  
  8. //颜色表定义 
  9. typedef struct tagRGBQUAD { 
  10. BYTE rgbBlue;// 蓝色的亮度（值范围为0-255) 
  11. BYTE rgbGreen; // 绿色的亮度（值范围为0-255) 
  12. BYTE rgbRed; // 红色的亮度（值范围为0-255) 
  13. BYTE rgbReserved;// 保留，必须为0 
  14. } RGBQUAD; 
  15. /*
  16.     定义头文件型
  17. */ 
  18. #pragma pack(1)//将字节对齐方式设置为1字节对齐 
  19. typedef struct 
  20. { 
  21. unsigned char id1;//位图文件的类型，必须为BM(占用0-1字节)  
  22. unsigned char id2; 
  23. unsigned int filesize;//位图文件的大小，以字节为单位（2-5字节） 
  24. unsigned int reserved;// 位图文件保留字，必须为0(6-9字节) 
  25. unsigned int bitmapdataoffset;//位图数据的起始位置，以相对于位图(10-13字节) 
  26. unsigned int bitmapheadersize;//BMP头的大小，固定为40（14-17字节） 
  27. unsigned int width;//图片宽度；以像素为单位（18-21字节） 
  28. unsigned int height;//图片高度；以像素为单位（22-25字节） 
  29. unsigned short planes;//图片位面数，必须为1(26-27字节） 
  30. unsigned short bitperpixel;//每个像素所需的位数，每个像素所需的位数（28-29字节） 
  31.    //只能是以下几个数：1（双色），4(16色），8(256色）或24（真彩色）  灰度级 
  32. unsigned int compression;//是否压缩（30-33字节） 
  33. //只能是以下几个数：0（不压缩），1(BI_RLE8压缩类型），2(BI_RLE4压缩类型） 
  34. unsigned int bitmapdatasize;//位图的大小，以字节为单位（34-37字节） 
  35. unsigned int hresolution;//位图水平分辨率，每米像素数（38-41字节） 
  36. unsigned int vresolution;//位图垂直分辨率，每米像素数（42-45字节) 
  37. unsigned int colors;//位图实际使用的颜色表中的颜色数（46-49字节） 
  38. unsigned int importantcolors;//位图显示过程中重要的颜色数（50-53字节） 
  39. //unsigned int  bmiColors[1];//调色板；(54 - 57字节) 
  40. unsigned char palette[256][4];//调色板 占256*4=1024字节 
  41. }BMPheader;//总大小40+14+1024=1078字节 
  42. typedef struct 
  43. { 
  44. BMPheader* bmpheader ; 
  45. unsigned char* bitmapdata;//图片数据； 
  46. }BMPheaderfile; 
  47. /*
  48.     求文件长度的函数
  49. */ 
  50. long getfilesize(FILE *f) 
  51. { 
  52. long pos,len; 
  53. pos=ftell(f);//ftell函数用于得到文件指针当前位置相对于文件首的偏移字节数 
  54. fseek(f,0,SEEK_END);//fseek函数用于移动文件指针相对于SEEK_END的偏移量为0 
  55. len=ftell(f);//len就是文件的长度 
  56. fseek(f,pos,SEEK_SET);//将文件指针移动到原来的地方 
  57. return len; 
  58. } 
  59. /*
  60.     主函数
  61. */ 
  62. int main() 
  63. { 
  64. BMPheaderfile *output=(BMPheaderfile*)malloc(sizeof(BMPheaderfile));//定义一个输出指针 
  65. unsigned char *data=NULL; 
  66.     FILE *fpr,*fpw; 
  67. /*
  68.   打开文件
  69. */ 
  70. if((fpr=fopen("start.bmp","rb"))==NULL) 
  71. { 
  72. printf("cannot open this file"); 
  73. exit(0); 
  74. } 
  75. if((fpw=fopen("end1.bmp","wb"))==NULL) 
  76. { 
  77. printf("cannot wirte this file"); 
  78. exit(0); 
  79. } 
  80. long length=getfilesize(fpr); 
  81. printf("文件的长度为%ld\n",length); 
  82. printf("文件的头长度为%d\n",sizeof(BMPheader)); 
  83. data=(unsigned char*)malloc(length*sizeof(char));//分配空间 
  84. if(0==fread(data,1,length,fpr))//读文件，从fpr指向的文件中读出length到data所指的内存空间去 
  85. { 
  86. printf("read failed\n"); 
  87. exit(0); 
  88. } 
  89. fclose(fpr);//释放指针 
  90. output->bmpheader=(BMPheader*)malloc(sizeof(BMPheader)); 
  91. memcpy(output->bmpheader,data,sizeof(BMPheader));//从data中拷贝sizeof(BMPheader)大小到output->bmpheader 
  92. /*
  93. 打印出图像中头文件的信息
  94. */ 
  95.     int height=output->bmpheader->height; 
  96. int width=output->bmpheader->width; 
  97. printf("filesize is %d\n",output->bmpheader->filesize); 
  98. printf("该图像每个像素所需要的位数：%d\n",output->bmpheader->bitperpixel); 
  99. printf("height is %d\n",output->bmpheader->height); 
  100. printf("width is %d\n",output->bmpheader->width); 
  101. data=data+sizeof(BMPheader); 
  102. output->bitmapdata=data; 
  103.   /*
  104.    中值滤波算法(选择3×3的滑动窗口)
  105.   */ 
  106. unsigned char pixel[9]={0};//滑动窗口的像素值，初始为0 
  107. unsigned char mid;//中值 
  108. unsigned char temp;//中间变量 
  109. int flag; 
  110. int m,i,j,x,h,w,y; 
  111. for(j=1;j<height-1;j++) 
  112. { 
  113. for(i=1;i<width-1;i++) 
  114. { 
  115. //将3×3滑动窗口中的所有像素值放入pixel[m] 
  116. m=0; 
  117. for(y=j-1;y<=j+1;y++) 
  118. for(x=i-1;x<=i+1;x++) 
  119. { 
  120. pixel[m]=data[y*width+x]; 
  121. m=m+1; 
  122. } 
  123. //让一位数组pixel[9]进行降序排列 
  124. do  
  125. { 
  126. flag=0;//循环结束的标志 
  127. for(m=0;m<9;m++) 
  128. { 
  129. if(pixel[m]<pixel[m+1]) 
  130. { 
  131. temp=pixel[m]; 
  132. pixel[m]=pixel[m+1]; 
  133. pixel[m+1]=temp; 
  134. flag=1; 
  135. }//if 
  136. }//for 
  137. } while (flag==1); 
  138. mid=pixel[4]; 
  139. output->bitmapdata[width*j+i]=mid; 
  140. } 
  141. } 
  142. /*
  143. N=5的中值滤波线段状
  144. */ 
  145. /*
  146.     for(i=0;i<height;i++)
  147. {
  148. for(j=2;j<width-1;j++)
  149. {
  150. m=0;
  151. for(x=j-2;x<=j+2;x++)
  152. pixel[m++]=data[i*width+x];
  153. for(h=0;h<5;h++)
  154. for(w=h+1;w<5;w++)
  155. {
  156. if(pixel[h]>pixel[w])
  157. {
  158. temp=pixel[w];
  159. pixel[w]=pixel[h];
  160. pixel[h]=temp;
  161. }
  162. }
  163. output->bitmapdata[i*width+j]=pixel[2];
  164. }
  165. }
  166. */ 
  167. /*
  168. 保存图像文件
  169. */ 
  170. fseek(fpw,0,0);  //fseek(fpw,0,SEEK_SET) 
  171. fwrite(output->bmpheader,1,sizeof(BMPheader),fpw);//写入图像的头文件 
  172. fwrite(output->bitmapdata,1,length-sizeof(BMPheader),fpw);//写入图像的数据信息 
  173. fclose(fpw);//释放指针 
  174. return 0; 
  175. } 

 

Andy20081251
关注
专栏目录
快速二维中值滤波算法的实现（灰度+RGB图）
qq_42974561的博客
12-14 3724
快速二维中值滤波算法的实现（灰度+RGB图）摘要算法流程及伪代码算法流程伪代码程序代码 摘要 本文主要完成对 Thomas S Huang, George J Yang, and Gregory J Yang. A Fast Two-Dimensional Median Filtering Algorithm[J]. IEEE Transactions on Acoustics Speech an...
图形学－降噪技术－2维中值滤波
qiuqiu173的专栏
01-25 1955
导读： 　　在图像采集过程中，由于电子设备的不稳定性（比如毛刺电压，电磁干扰等）会对获取的图像产生一些影响，这种影响叫噪声，降噪技术就是图像处理中消除这种噪声影响的技术。二维中值滤波就是其中一种方法。 　　中值滤波是一种局部图像平滑技术，属于非线性滤波，他可以是1维的也可以是2维的，因为图像是2维象素矩阵，所以这里使用2维中值滤波。 　　2维中值滤波算法是： 　　对于一幅图像的象素矩阵，取
ITK-中值滤波
最新发布
翟天保的博客
07-25 1221
使用ITK进行中值滤波
图像处理小tip——中值滤波的多种实现（包括快速中值滤波算法）
qq_48176859的博客
05-09 4298
文章目录中值滤波介绍opencv自带的中值滤波函数简单算法实现中值滤波快速算法实现中值滤波 中值滤波介绍 中值滤波介绍 中值滤波定义：将????(????为奇数)个数据按其值????????进行从大到小或者从小到大排列后得到一个有序序列:????0,????1, … ,????????−1,则????(????/2)称为中值。例如：有序序列10,11,12,13,14,15,16,17,18的???? = 9，有 [9/2] = 4，则中值为????4，即14。   根据以上可知，中值滤波就是以当前像素为
中值滤波的快速算法
LinJM-机器视觉
09-09 1万+
我想学过图像处理的人没有人会不知道中值滤波的，最早的时候我是在冈萨雷斯的图像处理课本[1]中学到的，后来在看Sonka的书[2]的时候又看到了中值滤波的介绍，下面我试着结合课本所学和网上的资料自己整理一篇中值滤波的介绍。 Jeremy Lin 中值滤波器是一种统计排序滤波器，由Tukey于1971年在文献[3]中提出。所谓的统计排序滤波器是一种非线性的空间滤波器，它的响应基于图像滤波器包围的图像区域中像素的排序，然后用统计排序结果决定的值代替中心像素的值。除了中值滤波器外，最大值滤波器和最小值滤波器也是统计
matlab中图像双边滤波_MATLAB--数字图像处理 中值滤波
weixin_39890517的博客
11-28 831
中值滤波概念： 中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替，让周围的像素值接近的真实值，从而消除孤立的噪声点。方法是用某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升（或下降）的为二维数据序列。二维中值滤波输出为g（x,y）=med{f(x-k,y-l),(k,l...
MATLAB的图像中值滤波?_matlab_中值滤波_
09-29
在MATLAB环境中，实现中值滤波是非常直观且便捷的。本篇文章将详细探讨MATLAB中的中值滤波操作及其相关知识点。 首先，我们需要理解中值滤波的基本原理。中值滤波是一种非线性的滤波方法，它不依赖于图像的统计特性...
均值滤波与中值滤波的应用_MATLAB中值滤波_均值滤波中值滤波的matlab代码_中值滤波_plandop_
09-30
在MATLAB中，我们可以利用内置函数实现均值滤波和中值滤波。这两种滤波方法各有特点，广泛应用于不同场景。下面我们将详细介绍这两种滤波方法以及它们在MATLAB中的实现。 **1. 均值滤波** 均值滤波是一种线性滤波...
图像去噪基于中值滤波实现图像去噪含Matlab源码 上传.zip
10-23
此项目适合本科及硕士级别的学生进行教研学习，因为图像处理是计算机视觉和信号处理课程中的常见主题，而中值滤波则是初学者可以理解和实现的一个基础去噪方法。 **图像示例**： 文件列表中的`image7.JPG`和`image3...
快速中值滤波及c语言实现.docx
07-03
"快速中值滤波及C语言实现" 快速中值滤波是一种非线性信号处理技术，对抑制图像噪声非常有效。在二维形式下，中值滤波器是一个奇数个像素的滑动窗口，窗口正中的像素灰度值用窗口内各个像素的中值代替。在本文中，...
中值滤波_高斯噪声_中值滤波_matlab图像处理_
09-29
3. **中值滤波**：中值滤波是一种非线性滤波方法，特别适合去除椒盐噪声。其基本原理是对每个像素点，用其邻域内像素值的中位数来代替原值。中值滤波器可以有效地保护图像边缘，因为它不依赖于平均值，而是在像素值...
中值滤波算法C语言实现MidValueFilter.rar
07-02
中值滤波算法C语言实现代码，采用VC++6.0环境验证，编写成了子函数形式（无main），原创函数代码，欢迎下载
图像处理中值滤波C/C++算法
06-02
该算法用于处理8位灰度值图像，使用N=3的值滤波的算法程序，可以很方便的转换成N=5、7、9的中值滤波算法。 经该算法处理带有椒盐噪声的图像后，可以看出与原噪声图像的不同，但是效果不是那么的明显。
中值滤波算法VC实现（代码）
04-26
中值滤波算法VC实现（代码），各位不要嫌小弟心黑 ，这10分花的绝对值，小弟费了好大劲才写出中并且加了详细注释，而且保证程序能够运行通过，效果与Maltab实现相差无几
中值滤波算法Quartus实现
03-04
中值滤波算法Quartus实现，一个计算机专业的毕业设计，Oh,shit! How can a students do this?
中值滤波算法
热门推荐
liuyun的专栏
11-19 4万+
数字图像在其形成、传输记录的过程中往往会受到很多噪声的的污染，比如：椒盐噪声、高斯噪声等，为了抑制和消除这些随即产生的噪声而改善图像的质量，就需要去、对图像进行去噪处理，去噪也就是滤波处理。 中值滤波和同态滤波都是图像增强的方法，但是这两种方法是从不同的方式进行改善一副图片的质量。 中值滤波是图像平滑的一种方法 它是一种非线性平滑滤波技术，在一定条件下可以克服线性滤波带来的图像细节的模糊问题，
中值滤波的一种快速计算方法
攻城狮凌风
06-24 2万+
        中值滤波(median filter)为常用图像滤波算法，用于处理椒盐噪声且能够保持边界，由于是非线性操作，传统的做法即对局部区域进行排序取中值的方法，较为耗时。        一种快速的计算法，即使用直方图取中位数代替排序，且利用相邻元素的直方图变动较少更新直方图，这样迭代获得滤波结果, 如下：Matlab源码如下：clc; clear all; close all; imag...
中值滤波算法及例程
嵌入式之斋
10-09 2276
在函数内部，我们通过遍历图像的每个像素，获取周围邻域像素的值，并使用np.median函数计算中值，然后将中值赋值给当前像素。中值滤波是一种常用的非线性图像滤波算法，它能够有效去除图像中的椒盐噪声（即孤立的亮或暗像素点），同时保持图像边缘和细节的清晰度。中值滤波的主要思想是使用一个滑动窗口，在窗口内对像素值进行排序，并将排序后的中间值作为中心像素的新值。中值滤波算法的关键在于选择合适的窗口大小，较小的窗口大小可以更好地保留图像细节和边缘特征，但可能无法有效去除较大的噪声；将中值赋给当前像素作为滤波后的值。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值