#include <fstream.h>
#include <math.h>
#define nr 300
#define nc 300
typedef unsigned char BYTE;
//-------------------------------------------------------
// f9.cpp
// image shaping using 3x3 Laplacian
// MIAT Lab, Kaohsiung, Taiwan
// CHEN Ching-Han, 2001.10.09
//-------------------------------------------------------
void Laplacian(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];

   Laplacian(ima1, ima2);

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

}

void Laplacian(BYTE **InIm, BYTE **OutIm)
{
   int   r,c, y,x, Sum;
   static int Mask [3][3] = { { 0,  1, 0},
                              { 1, -4, 1},
                              { 0,  1, 0} };

   for (r=0; r<nr; r++) for (c=0; c<nc; c++)  OutIm [r][c] = 0;

   for (r=1; r<nr-1; r++)
   {
      for (c=1; c<nc-1; c++)
      {
         Sum = 0;
         for (y=-1; y<=1; y++) for (x=-1; x<=1; x++)
           Sum += InIm [r+y] [c+x] * Mask [y+1] [x+1];
         OutIm [r][c] = Sum/9;
      }
   }
}

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