#include <fstream.h>
#include <math.h>
#include "array.h"
//-------------------------------------------------------
// f7.cpp
// image enhencement using accumulate histogram density function
// MIAT Lab, Kaohsiung, Taiwan
// CHEN Ching-Han, 2001.08.26
//-------------------------------------------------------
void enhencement(uc2D &ima0, uc2D ima1);

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];
   int ii,jj,count;

   enhencement(ima0, ima1);

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

}


void enhencement(uc2D &ima0, uc2D ima1)
{
	long ImaSize=ima0.nr*ima0.nc;
   int histo[256],table[256];
	float accpbhisto[256];
	for(int i=0;i<256;i++)
   {
     histo[i]=0;
     table[i]=0;
     accpbhisto[i]=0.0;
   }
	for(int i=0;i<ima0.nr;i++)for(int j=0;j<ima0.nc;j++)
     histo[ima0.m[i][j]]++;
   accpbhisto[0]=float(histo[0])/float(ImaSize);
	for(int i=1;i<256;i++)
	{
		accpbhisto[i]=accpbhisto[i-1]+float(histo[i])/float(ImaSize);
	}
	for(int i=0;i<256;i++)table[i]=char(accpbhisto[i]*256.);
	for(int i=0;i<ima0.nr;i++)for(int j=0;j<ima0.nc;j++)
     ima1.m[i][j]=table[ima0.m[i][j]];
}