#include?“StdAfx.h“
#include?“Gabor.h“
#include?
#include?

CGabor::CGabor（void）
{
}

CGabor::~CGabor（void）
{
}

#define		SIZE_L	32
#define		SIZE_W	16
#define		SIZE_L2	16
#define		SIZE_W2	8

#ifndef?PI
#define?PI?3.1415926535897932384626433832795
#endif

/*
?Function?Name	:?Frequency
?Functionality	:?Prepare?frequency?For?Gabor?Filter
				??
*/
int?CGabor::Frequency（unsigned?char?*?ucImg?double?*?fDirection?double?*?fFrequency?long?lWidth?long?lHeight）
{
	//?the?peak?point?of?the?sin?wave
	int		peak_pos[SIZE_L];
	int		peak_cnt;
	double	peak_freq;
	double	Xsig[SIZE_L];
//	double	pmin?=?0?pmax?=?0;

	//?directions
????double	dir?=?0.0;
????double	cosdir?=?0.0;
????double	sindir?=?0.0;
	double	maxPeak?minPeak;

	double	*frequency1?=?new?double[lWidth*lHeight];

	memset（fFrequency?0?sizeof（double）*lWidth*lHeight）;
	memset（frequency1?0?sizeof（double）*lWidth*lHeight）;

	long	xy;
	long	dk;
	long	uv;

	for（y=SIZE_L2;?y	{
		for（x=SIZE_L2;?x		{
			dir		=?fDirection[（y+SIZE_W2）*lWidth?+?（x+SIZE_W2）];
			cosdir	=?-sin（dir）;
			sindir	=?cos（dir）;

			//	calc?the?x-signature?X[0]...X[L-1]
			for（k=0;?k			{
				Xsig[k]?=?0.0;
				for（d=0;?d				{
					u?=?（long）（x?+?（d-SIZE_W2）*cosdir?+?（k-SIZE_L2）*sindir）;
					v?=?（long）（y?+?（d-SIZE_W2）*sindir?-?（k-SIZE_L2）*cosdir）;
					//?never?over?the?edge
					if（u<0）
						u?=0;
					else?if（u>lWidth-1）
						u?=?lWidth-1;
					if（v<0）
						v?=?0;
					else?if（v>lHeight-1）
						v?=?lHeight-1;

					Xsig[k]?+=?ucImg[u+v*lWidth];
				}
				Xsig[k]?/=?SIZE_W;
			}

			//?here?we?have?got?the?SIN?wave

			//?see?if?the?wave?exists?max?and?min?value?should?have?a?range
			maxPeak?=?minPeak?=?Xsig[0];
			for（k=0;?k			{
				if（minPeak>Xsig[k]）
					minPeak	=?Xsig[k];
				if（maxPeak					maxPeak	=?Xsig[k];
			}

			peak_cnt	=?0;
			if（（maxPeak?-?minPeak）>64）
			{
				for（k=0;?k				{
					if（（Xsig[k-1]?=?Xsig[k+1]））
						peak_pos[peak_cnt++]?=?k;
				}
			}

			//?mean?peak
			peak_freq	=?0.0;
			if（peak_cnt>=2）
			{
				for（k=0;?k				{
					peak_freq	+=?（peak_pos[k+1]?-?peak_pos[k]）;
				}
				peak_freq	/=?peak_cnt-1;
			}
			//?peak?frequency?should?be?limited
			if（peak_freq<3.0?||?peak_freq>25.0）
				frequency1[x+y*lWidth]?=?0.0;
			else
				frequency1[x+y*lWidth]?=?1.0/peak_freq;
		}
	}

	//?mean?filter?for?the?frequency
	for（y=SIZE_L2;?y	{
		for（x=SIZE_L2;?x		{
			k?=?x+y*lWidth;
			peak_freq?=?0.0;
			for（v=-2;?v<=2;?v++）
			{
				for（u=-2;?u<=2;?u++）
				{
					peak_freq	+=?frequency1[（x+u）+（y+v）*lWidth];
				}
			}
			fFrequency[k]?=?peak_freq/25;
		}
	}

	delete?frequency1;

	return?0;
}

/*
?Function?Name	:?ImageEnhance
?Functionality	:?image?enc