OPENCV中值滤波

#include "stdio.h"
#include "cv.h"
#include "highgui.h"
uchar FilterMid(IplImage* image,int k);
int main()
{
 IplImage *img=0,*img2;
 int height,width,step,channels;
 uchar *data;
 img=cvLoadImage("11.jpg",0);
 img2=cvCloneImage(img);
 if(!img)
 {
  printf("can not load image file:%s\n","tln.jpg");
 }
 height=img->height;
 width=img->width;
 step=img->widthStep;
 channels=img->nChannels;
 data=(uchar*)img->imageData;
 printf("processing a %dx%d image with %d channels\n",height,width,channels);
 cvNamedWindow("before",CV_WINDOW_AUTOSIZE);
 cvNamedWindow("after",CV_WINDOW_AUTOSIZE);
    FilterMid( img,3);
 cvShowImage("after",img);
 cvShowImage("before",img2);
 cvWaitKey(0);
 cvReleaseImage(&img);
 return 0;
}


uchar FilterMid(IplImage* image,int k)
{
int x,y,i,j;
 IplImage *image2=cvCreateImage(cvSize(image->width,image->height),IPL_DEPTH_8U, 1);
 uchar *ImagePix=(uchar *)image2->imageData;
 int m=(k-1)/2;
 for(x=m;x<image->height-m;++x)
  for(y=m;y<image->width-m;++y)
  {
   uchar PixArray[100];
   int t=0;
   for(i=-m;i<m+1;++i)
    for(j=-m;j<m+1;j++)
    {
      PixArray[t++]=((uchar *)image->imageData)[(x+i)*image->widthStep+y+j];
    }
    for(i=0;i<k*k-1;++i)
     for(j=0;j<k*k-i-1;++j)
     {
      if(PixArray[j]>PixArray[j+1])
      {
       uchar k=PixArray[j];
      PixArray[j]=PixArray[j+1];
      PixArray[j+1]=k;
      }
     }
    
     ImagePix[x*image->widthStep+y]=PixArray[(k*k-1)/2];
  }
  image=cvCloneImage(image2);
  cvReleaseImage(&image2);
  return 0;
}