8位位图的中值滤波
for(int i=0;i<nHeight;i++)
{
for(int j=(select-1)/2;j<nWidth-(select-1)/2;j++)
{
//取得临近的像素
int m=-(select-1)/2;
for(int n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j+m);
m++;
}
//对取得的像素进行排序
for(int k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置中值为其点的像素值
*(image+lLineBytes*i+j)=data[(select-1)/2];
}
}
24位位图的中值滤波
for(int i=0;i<nHeight;i++)
{
for(int j=((select-1)/2)*3;j<(nWidth-(select-1)/2)*3;j++)
{
//对B分量进行中值滤波
int m=-(select-1)/2;
//取得j两旁的B分量的值
for(int n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j-2+m);
m++;
m++;
m++;
}
//对取得的B分量进行由小到大的排序
for(int k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置B分量的值为排序后的序列的中间值
*(image+lLineBytes*i+j)=data[(select-1)/2];
//对G分量进行中值滤波
j++;
m=-(select-1)/2;
//取得j两旁的G分量的值
for(n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j-2+m);
m++;
m++;
m++;
}
//对取得的G分量进行由小到大的排序
for(k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置G分量的值为排序后的序列的中间值
*(image+lLineBytes*i+j)=data[(select-1)/2];
//对R分量进行中值滤波
j++;
m=-(select-1)/2;
//取得j两旁的R分量的值
for(n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j-2+m);
m++;
m++;
m++;
}
//对取得的R分量进行由小到大的排序
for(k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置R分量的值为排序后的序列的中间值
*(image+lLineBytes*i+j)=data[(select-1)/2];
}
}
滤波效果说明
由于对24位位图进行中值滤波的话,会改变RGB各分量的值,所以图形的颜色是会发生变化的。但是对于8位的位图来说,由于都是灰度的颜色所以不会有特别明显的变化,而且滤波的窗口选的越大,对应的滤波效果的模糊会比较明显。
8位位图的滤波效果图片
24位位图的滤波效果图片
完整的中值滤波程序
void CICETIMDlg::OnBtnMedianfilter()
{
// TODO: Add your control notification handler code here
if(m_dib.GetBitCount()==0)
{
AfxMessageBox("请先打开位图文件");
return;
}
int select;
CMedianFilterDlg dlg=new CMedianFilterDlg();
if(dlg.DoModal()==IDOK)
{
UpdateData();
select=dlg.m_select;
}
else
{
return;
}
int nHeight=m_dib.GetHeight();
int nWidth=m_dib.GetWidth();
BYTE* image=m_dib.GetDibData();
long lLineBytes=m_dib.GetLineBytes();
BYTE* data=new BYTE[select];
BeginWaitCursor();
int nBitCount=m_dib.GetBitCount();
if(nBitCount==8)
{
for(int i=0;i<nHeight;i++)
{
for(int j=(select-1)/2;j<nWidth-(select-1)/2;j++)
{
//取得临近的像素
int m=-(select-1)/2;
for(int n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j+m);
m++;
}
//对取得的像素进行排序
for(int k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置中值为其点的像素值
*(image+lLineBytes*i+j)=data[(select-1)/2];
}
}
}
else if(nBitCount==24)
{
for(int i=0;i<nHeight;i++)
{
for(int j=((select-1)/2)*3;j<(nWidth-(select-1)/2)*3;j++)
{
//对B分量进行中值滤波
int m=-(select-1)/2;
//取得j两旁的B分量的值
for(int n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j-2+m);
m++;
m++;
m++;
}
//对取得的B分量进行由小到大的排序
for(int k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置B分量的值为排序后的序列的中间值
*(image+lLineBytes*i+j)=data[(select-1)/2];
//对G分量进行中值滤波
j++;
m=-(select-1)/2;
//取得j两旁的G分量的值
for(n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j-2+m);
m++;
m++;
m++;
}
//对取得的G分量进行由小到大的排序
for(k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置G分量的值为排序后的序列的中间值
*(image+lLineBytes*i+j)=data[(select-1)/2];
//对R分量进行中值滤波
j++;
m=-(select-1)/2;
//取得j两旁的R分量的值
for(n=0;n<select;n++)
{
data[n]=*(image+lLineBytes*i+j-2+m);
m++;
m++;
m++;
}
//对取得的R分量进行由小到大的排序
for(k=0;k<select;k++)
{
for(int kk=k+1;kk<select;kk++)
{
if(data[kk]<data[k])
{
int a;
a=data[k];
data[k]=data[kk];
data[kk]=a;
}
}
}
//设置R分量的值为排序后的序列的中间值
*(image+lLineBytes*i+j)=data[(select-1)/2];
}
}
}
else
{
AfxMessageBox("暂时只能处理8或24位位图");
return;
}
ShowImage(m_dib,"中值滤波");
EndWaitCursor();
}