#include?
#include?
#include?
#include?
#include?
using?namespace?std;

#define?w?640
#define?h?480
#define?pi?3.14159
#define?t?32
#define?l?600
#define?d?500

#define?beta?2*pi*9*15.5/532

//參考平面
double?a0[h][w]?=?{0};
double?a1[h][w]?=?{0};
double?a2[h][w]?=?{0};
double?a3[h][w]?=?{0};

//被測物
double?b0[h][w]?=?{0};
double?b1[h][w]?=?{0};
double?b2[h][w]?=?{0};
double?b3[h][w]?=?{0};

double?det1[h][w]?=?{0};
double?det2[h][w]?=?{0};

int?k0[h][w]?=?{0};	//參考平面階次
int?k1[h][w]?=?{0};	//被測物階次

double?ph0[h][w]?=?{0};	//相對相位
double?ph1[h][w]?=?{0};	//參考平面絕對相位
double?ph2[h][w]?=?{0};	//被測物絕對相位
double?d_ph[h][w]?=??{0};?//相位差

double?z[?h?*w?]?=?{0};	//z軸坐標
int?main（）
{

	//CvMat*?mat?=?cvCreateMat（whCV_8UC1）;
	//cvZero（mat）;//將矩陣置0
	//參考平面投影光柵
	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			a0[i][j]?=?128+127*cos（2*pi/t*j）;
			a1[i][j]?=?128+127*cos（2*pi/t*j+0.5*pi）;
			a2[i][j]?=?128+127*cos（2*pi/t*j+pi）;
			a3[i][j]?=?128+127*cos（2*pi/t*j+1.5*pi）;
		}		
	}
	//被測物體投影光柵
	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			b0[i][j]?=?a0[i][j];
			b1[i][j]?=?a1[i][j];
			b2[i][j]?=?a2[i][j];
			b3[i][j]?=?a3[i][j];
		}		
	}

	for?（int?i?=?150;?i?	{
		for?（int?j?=?250;?j?		{
			b0[i][j]?=?128+127*cos（2*pi/t*j?+?pi/10）;
			b1[i][j]?=?128+127*cos（2*pi/t*j?+?pi/10?+?0.5*pi）;
			b2[i][j]?=?128+127*cos（2*pi/t*j?+?pi/10?+?pi）;
			b3[i][j]?=?128+127*cos（2*pi/t*j?+?pi/10?+?1.5*pi）;
		}		
	}

	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			ph0[i][j]=atan（（a3[i][j]-a1[i][j]）/（a0[i][j]-a2[i][j]））;
			
		}		
	}
	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			det1[i][j]=atan（（a2[i][j]-a1[i][j]）/（a0[i][j]-a1[i][j]））+pi/4;
		}		
	}

	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			det2[i][j]?=?atan（（（1-cos（beta））*（a0[i][j]-a2[i][j]））/（（sin（beta））*（2*a1[i][j]-a0[i][j]-a2[i][j]）））;
		}		
	}
	for?（int?i?=?0;?i?	{
		for?（int?j?=?1;?j?		{
			if（（det2[i][j]?-?det2[i][j-1]）?>=??pi/2）
				k0[i][j]?=?k0[i][j-1]?-1;
			if（abs（det2[i][j]?-?det2[i][j-1]）?				k0[i][j]?=?k0[i][j-1];
			if（（det2[i][j]?-?det2[i][j-1]）?<=??-pi/2）
				k0[i][j]?=?k0[i][j-1]?+1;
		}		
	}
	//參考平面絕對相位
	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			ph1[i][j]?=?det2[i][j]?+?k0[i][j]?*?pi;
		}		
	}

	//被測物解相位
	//被測物相對相位
		for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			ph0[i][j]=atan（（b3[i][j]-b1[i][j]）/（b0[i][j]-b2[i][j]））;
			
		}		
	}
	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			det1[i][j]=atan（（b2[i][j]-b1[i][j]）/（b0[i][j]-b1[i][j]））+pi/4;
		}		
	}

	for?（int?i?=?0;?i?	{
		for?（int?j?=?0;?j?		{
			det2[i][j]?=?atan（（（1-cos（beta））*（b0[i][j]-b2[i][j]））/（（sin（beta））*（2*b1[i][j]-b0[i][j

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2013-10-14?21:07??mat\
?????目錄???????????0??2013-10-14?20:56??mat\Debug\
?????文件??????114688??2013-10-14?21:05??mat\Debug\mat.exe
?????文件??????752168??2013-10-14?21:05??mat\Debug\mat.ilk
?????文件?????1788928??2013-10-14?21:05??mat\Debug\mat.pdb
?????目錄???????????0??2013-10-14?20:58??mat\mat\
?????目錄???????????0??2013-10-14?21:05??mat\mat\Debug\
?????文件???????17282??2013-10-14?21:05??mat\mat\Debug\CL.read.1.tlog
?????文件?????????578??2013-10-14?21:05??mat\mat\Debug\CL.write.1.tlog
?????文件????????1166??2013-10-14?21:05??mat\mat\Debug\cl.command.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link-cvtres.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link-cvtres.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link-rc.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link-rc.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.1120-cvtres.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.1120-cvtres.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.1120-rc.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.1120-rc.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.1120.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.1120.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2232-cvtres.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2232-cvtres.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2232-rc.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2232-rc.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2232.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2232.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2504-cvtres.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2504-cvtres.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2504-rc.read.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2504-rc.write.1.tlog
?????文件???????????2??2013-10-14?21:05??mat\mat\Debug\link.2504.read.1.tlog
............此處省略101個文件信息