數位影像的中值濾波

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