#include <fstream.h>
#include <math.h>
#define nr 300
#define nc 300
typedef unsigned char BYTE;
//-------------------------------------------------------
// f8.cpp
// Convolution filtering using moving average
// MIAT Lab, Kaohsiung, Taiwan
// CHEN Ching-Han, 2001.10.09
//-------------------------------------------------------
void Average (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];

   Average(5, ima1, ima2);

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

}

void Average(int WinSize, BYTE **InIm, BYTE **OutIm)
{
   int   r,c, y,x, n, Area;
   long  Sum;

   n = (WinSize-1) >> 1;
   Area = (2*n+1) * (2*n+1);

   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++)
      {
         Sum = 0;
         for (y=-n; y<=n; y++)for (x=-n; x<=n; x++)Sum += InIm [r+y] [c+x];
         OutIm [r][c] = (BYTE) (Sum/Area);
      }
   }
}

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