图像的中值滤波

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();
}