#include?“stdafx.h“
#include?“ImageProcessing.h“
#include?“stdio.h“

/************************************************************************/
/*?圖像傅里葉變換???????????????????????????????????????????????????????*/
/************************************************************************/
void?FFT_2D（IplImage?*?src?CvMat?*?dst）
{
	CvMat?*?real?=?cvCreateMat（dst->rows?dst->cols?CV_64FC1）;
	CvMat?*?imag?=?cvCreateMat（dst->rows?dst->cols?CV_64FC1）;
	CvMat?*?cplx?=?cvCreateMat（dst->rows?dst->cols?CV_64FC2）;
	CvMat?temp;

	cvGetSubRect（real?&temp?cvRect（00?src->width?src->height））;
	cvScale（src?&temp?1.0?0.0）;

	if?（dst->cols?>?src->width）
	{
		cvGetSubRect（real?&temp?cvRect（src->width?0?dst->cols?-?src->width?src->height））;
		cvZero（&temp）;
	}
	
	cvZero（imag）;
	cvMerge（real?imag?NULL?NULL?cplx）;

	cvReleaseMat（&real）;
	cvReleaseMat（&imag）;

	cvDFT（cplx?dst?CV_DXT_FORWARD?src->height）;

	cvReleaseMat（&cplx）;
}

/************************************************************************/
/*?相位相關的圖像配準???????????????????????????????????????????????????*/
/************************************************************************/
BOOL?MethodPhaseCorrelation（IplImage?*?image1?IplImage?*?image2?int?&?iOffX?int?&?iOffY）
{
	int?iDftX?=?image1->width;
	int?iDftY?=?image1->height;

	CvMat?*?mat1?=?cvCreateMat（iDftY?iDftX?CV_64FC2）;
	CvMat?*?mat2?=?cvCreateMat（iDftY?iDftX?CV_64FC2）;

	FFT_2D（image1?mat1）;
	FFT_2D（image2?mat2）;

	cvMulSpectrums（mat1?mat2?mat1?CV_DXT_MUL_CONJ）;
	
	for?（int?i?=?0;?i?		for?（int?j?=?0;?j?		{
			double?&?re?=?（（double*）（mat1->data.ptr?+?mat1->step?*?i））[j?*?2];
			double?&?im?=?（（double*）（mat1->data.ptr?+?mat1->step?*?i））[j?*?2?+?1];
			double?norm?=?sqrt（re?*?re?+?im?*?im）;

			re?/=?norm;
			im?/=?norm;
		}
	cvDFT（mat1?mat1?CV_DXT_INV_SCALE?mat1->rows）;

	double?fMax?=?-DBL_MAX;

	for?（int?i?=?0;?i?		for?（int?j?=?0;?j?		{
			double?re?=?（（double*）（mat1->data.ptr?+?mat1->step?*?i））[j?*?2];
			if?（re?>?fMax）
			{
				fMax?=?re;
				iOffX?=?j;
				iOffY?=?i;
			}
		}

//?	FILE?*?fp1;
//?	
//?	fopen_s（&fp1?“data.txt“?“w“）;
//?
//?	for?（int?i?=?0;?i?//?	{
//?		for?（int?j?=?0;?j?//?		{
//?			double?re?=?（（double*）（mat1->data.ptr?+?mat1->step?*?i））[j?*?2];
//?			double?im?=?（（double*）（mat1->data.ptr?+?mat1->step?*?i））[j?*?2?+?1];
//?			
//?			fprintf_s（fp1?“%.4lf\t“?re）;
//?		}
//?		fprintf_s（fp1?“\n“）;
//?	}
//?	fclose（fp1）;

	cvReleaseMat（&mat1）;
	cvReleaseMat（&mat2）;

	return?TRUE;
}

/************************************************************************/
/*?圖像配準?????????????????????????????????????????????????????????????*/
/************************************************************************/
BOOL?ImageRegistration（int?iMethod?IplImage?*?image1?IplImage?*?image2?int?&?iOffX?int?&?iOffY）
{
	if?（iMeth

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????6560??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\Debug\BuildLog.htm

?????文件??????18546??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\Debug\ImageProcessing.obj

?????文件?????786432??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\Debug\ImageProcessing.pch

?????文件???????9687??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\Debug\stdafx.obj

?????文件?????224256??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\Debug\vc90.idb

?????文件?????159744??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\Debug\vc90.pdb

?????文件???????3118??2009-11-05?18:54??Phase?Correlation\ImageProcessing\ImageProcessing\ImageProcessing.cpp

?????文件???????4269??2009-11-05?15:25??Phase?Correlation\ImageProcessing\ImageProcessing\ImageProcessing.vcproj

?????文件???????1377??2009-11-05?18:56??Phase?Correlation\ImageProcessing\ImageProcessing\ImageProcessing.vcproj.Ker.Administrator.user

?????文件???????1258??2009-11-05?14:45??Phase?Correlation\ImageProcessing\ImageProcessing\ReadMe.txt

?????文件????????302??2009-11-05?14:45??Phase?Correlation\ImageProcessing\ImageProcessing\stdafx.cpp

?????文件????????371??2009-11-05?14:45??Phase?Correlation\ImageProcessing\ImageProcessing\stdafx.h

?????文件????????765??2009-11-05?14:45??Phase?Correlation\ImageProcessing\ImageProcessing\targetver.h

?????文件????????911??2009-11-05?14:45??Phase?Correlation\ImageProcessing\ImageProcessing.sln

????..A..H.?????18432??2009-11-05?18:56??Phase?Correlation\ImageProcessing\ImageProcessing.suo

?????文件????????617??2009-11-05?16:35??Phase?Correlation\Include\ImageProcessing.h

?????文件??????29008??2009-11-05?18:54??Phase?Correlation\Lib\ImageProcessing.lib

?????目錄??????????0??2009-11-05?18:57??Phase?Correlation\ImageProcessing\ImageProcessing\Debug

?????目錄??????????0??2009-11-05?18:57??Phase?Correlation\ImageProcessing\ImageProcessing

?????目錄??????????0??2009-11-05?18:57??Phase?Correlation\ImageProcessing

?????目錄??????????0??2009-11-05?18:57??Phase?Correlation\Include

?????目錄??????????0??2009-11-05?18:57??Phase?Correlation\Lib

?????目錄??????????0??2009-11-05?18:57??Phase?Correlation

-----------?---------??----------?-----??----

??????????????1265653????????????????????23