#include <fstream.h>
#include <math.h>
#define nr 300
#define nc 300
typedef unsigned char BYTE;
//-------------------------------------------------------
// f12.cpp
// Convolution filtering using Median fiter
// MIAT Lab, Kaohsiung, Taiwan
// CHEN Ching-Han, 2001.10.09
//-------------------------------------------------------
void Median(int WinSize, BYTE **InIm, BYTE **OutIm);
BYTE ** array2D(int nnr,int nnc);
void main()
{

   ifstream in1("01.raw",ios::binary);
   ofstream out1("test.raw",ios::binary);
   BYTE **ima1,**ima2;
   ima1=array2D(nr,nc);
   ima2=array2D(nr,nc);

   for(int i=0;i<nr;i++)for(int j=0;j<nc;j++)in1>>ima1[i][j];

   Median(5, ima1, ima2);

   for(int i=0;i<nr;i++)for(int j=0;j<nc;j++)
    	out1<<ima2[i][j];

}

void Median(int WinSize, BYTE **InIm, BYTE **OutIm)
{
   int   r,c, y,x, i,j, n, Area;
   BYTE  Buf;
   BYTE  *Lst;

   n = (WinSize-1) >> 1;
   Area = (2*n+1) * (2*n+1);
   Lst = new BYTE[Area];

   for (r=0; r<nr; r++)for (c=0; c<nc; c++)OutIm[r][c] = 0;
   for (r=n; r<nr-n; r++)
   {
     for (c=n; c<nc-n; c++)
     {
       i=0;
       for (y=-n; y<=n; y++)for (x=-n; x<=n; x++){Lst [i] = InIm [r+y] [c+x];i++;}
/*---------- bubble sort -------------------*/
       for (i=0; i<Area-1; i++)for (j=Area-1; i<j; j--)if (Lst[j-1] > Lst[j])
       {
           Buf      = Lst[j-1];
           Lst[j-1] = Lst[j];
           Lst[j]   = Buf;
       }
       OutIm [r][c] = Lst [Area/2];  // get madian value
     }
   }
}

BYTE ** array2D(int nnr,int nnc)
{
  BYTE **m;
  m=new BYTE*[nnr];
  for(int i=0;i<nnr;i++)m[i]=new BYTE[nnc];
  for(int i=0;i<nnr;i++)for(int j=0;j<nc;j++)m[i][j]=0;
  return m;
}