#include <fstream.h>
#include <math.h>
#include "array.h"
//-------------------------------------------------------
// f5.cpp
// Compute Lacal minimum image by a 3x3 mask
// MIAT Lab, Kaohsiung, Taiwan
// CHEN Ching-Han, 2001.08.24
//-------------------------------------------------------
void main()
{

   ifstream in1("01.raw",ios::binary);
	ofstream out1("test.raw",ios::binary);
   uc2D ima0,ima1;
   ima0.Initialize(300,300);
   ima1.Initialize(300,300);

   for(int i=0;i<ima0.nr;i++)for(int j=0;j<ima0.nc;j++)
    	in1>>ima0.m[i][j];
   for(int i=0;i<ima1.nr;i++)for(int j=0;j<ima1.nc;j++)ima1.m[i][j]=255;
   int ii,jj,count;

   for(int i=1;i<ima1.nr-1;i++)for(int j=1;j<ima1.nc-1;j++)
   {
// Get local minimum
     count=0;
     for(ii=-1;ii<=1;ii++)for(jj=-1;jj<=1;jj++)
       if(ima0.m[i][j]>ima0.m[i+ii][j+jj])count++;
     if(count==0) ima1.m[i][j]=0;
     else ima1.m[i][j]=255;
   }

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

}