#include?“stdio.h“
#include?“math.h“
#include?“iostream.h“
#define?PI?3.1415926
#define?N?4

//求矩陣a的轉(zhuǎn)置矩陣b，a為m行、n列
void?transpose（double?*a?double?*b?int?m?int?n）;

//矩陣a乘以矩陣b，結(jié)果存儲(chǔ)在c中，a為m×n大小，b為n×l大小
void?multiply（double?*a?double?*b?double?*c?int?m?int?n?int?l）;

//求矩陣a的逆
int?inverse1（double?*a?int?n）;

//輸出m行、n列的矩陣a
void?shuchu（double?*a?int?m?int?n）;

//求矩陣a與b的差值矩陣，a、b、c大小均為m×n
void?subtrat（double?*a?double?*b?double?*c?int?m?int?n）;

double?R1[3][3]={0.0}R2[3][3]={0.0};?

double?L[8][1]={0.0};
double?A[8][6]={0.0}AT[6][8]={0.0}ATA[6][6]={0.0}AX[8][1]={0.0}
???????ATL[6][1]={0.0}M[6][6]={0.0}B[6][8]={0.0};
double?AM[8][1]={0.0}X[6][1]={0.0}V[8][1]={0.0}VT[1][8]={0.0}VTV[1][1]={0.0};
double?xy1[4][2]={0.0}xy2[4][2]={0.0};
double?_X[4][1]={0.0}??_Y[4][1]={0.0}??_Z[4][1]={0.0};
double?M01M02;
double?W[6][1]={123456};
int?n1（0）n2（0）;

void?main（）
{?
	printf（“???????????????????????????空間后交—前交程序???????????????????????????\n“）;
????printf（“???????????????????????????????????????????——測繪0902班0405090307??張超??\n“）;
	FILE?*fp?=?NULL;
	FILE?*fp1?=?NULL;
	if（（fp=fopen（“image.txt““r“））?==?NULL）
	{
		printf（“Open?file?error!“）;
		return;
	}
	if（（fp1=fopen（“ground.txt““r“））?==?NULL）
	{
		printf（“Open?file?error!“）;
		return;
	}
	
	//像點(diǎn)坐標(biāo)和地面點(diǎn)坐標(biāo)
	double?imagecontrol[4][4]={0.0};
	double?groundcontrol[4][3]={0.0};
??
	double?f=0.150000;
	
	long?ijr;


	for（i=0;?i<4;?i++）
	{
		for（j=0;?j<4;?j++）
		{
			fscanf（fp?“%lf“?&imagecontrol[i][j]）;
			imagecontrol[i][j]?/=?1000.0;
		}
		for（r=0;?r<3;?r++）
		{
			fscanf（fp1?“%lf“?&groundcontrol[i][r]）;
		}
	}
	fclose（fp）;
	fclose（fp1）;
	double?Ix1Iy1Gx1Gy1;
	double?m1=0.0;
	for（i=0;i<3;i++）
	{
		for（j=i+1;j<4;j++）
		{
			Ix1=imagecontrol[i][0]-imagecontrol[j][0];
			
			Iy1=imagecontrol[i][1]-imagecontrol[j][1];
			
			Gx1=groundcontrol[i][0]-groundcontrol[j][0];
			
			Gy1=groundcontrol[i][1]-groundcontrol[j][1];
			
			m1+=sqrt（（pow（Gx12）+pow（Gy12））/（pow（Ix12）+pow（Iy12）））;
			
		}
	}
	m1/=6;

?//???shuchu（?&imagecontrol[0][0]?4?4）;

/****************************************************************************************************************/
/****************************************************************************************************************/
/****************************************************************************************************************/
	printf（“\n-----------------------------左片外方位元素的計(jì)算---------------------------\n“）;

	double?Zs1（0.0）Xs1（0.0）Ys1（0.0）p1（0.0）w1（0.0）k1（0.0）;
	//1.計(jì)算攝影中心XsYs的初始近似值
	printf（“\n??????????????????????????攝影比例尺m1=%lf\n“m1）;
	for（i=0;i<4;i++）???????????????????
	{
		Xs1+=groundcontrol[i][0];
		Ys1+=groundcontrol[i][1];
		Zs1+=groundcontrol[i][2];
	}
	Xs1/=4.0;
	Ys1/=4.0;
	Zs1/=4.0;
	Zs1+=m1*f;
	//計(jì)算旋轉(zhuǎn)矩陣R
	do???????
	{
		
		R1[0][0]=cos（p1）*cos（k1）-sin（p1）*sin（w1）*sin（k1）;
		R1[0][1]=（-1）

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2015-07-01?14:12??程序\
?????文件???????20225??2011-11-15?22:03??程序\001.cpp
?????目錄???????????0??2015-07-01?09:47??程序\Debug\
?????文件?????????772??2015-07-01?09:47??程序\Debug\cl.command.1.tlog
?????文件????????1318??2015-07-01?09:47??程序\Debug\CL.read.1.tlog
?????文件?????????308??2015-07-01?09:47??程序\Debug\CL.write.1.tlog
?????文件???????46006??2011-11-14?10:34??程序\Debug\HFJH.obj
?????文件??????278593??2012-05-26?10:59??程序\Debug\SY0.exe
?????文件??????321100??2012-05-26?10:59??程序\Debug\SY0.ilk
?????文件?????????146??2015-07-01?09:47??程序\Debug\SY0.lastbuildstate
?????文件????????1182??2015-07-01?09:47??程序\Debug\SY0.log
?????文件???????????0??2015-07-01?09:47??程序\Debug\SY0.unsuccessfulbuild
?????文件???????19456??2015-07-01?09:47??程序\Debug\vc100.idb
?????文件???????45056??2015-07-01?09:47??程序\Debug\vc100.pdb
?????文件???????41984??2012-05-26?10:59??程序\Debug\vc60.idb
?????文件???????61440??2012-05-26?10:59??程序\Debug\vc60.pdb
?????文件?????????116??2011-10-26?16:52??程序\ground.txt
?????文件?????????128??2011-10-26?16:56??程序\image.txt
?????文件?????????153??2011-11-14?21:23??程序\image1.txt
?????文件?????????143??2011-11-14?21:23??程序\image2.txt
?????目錄???????????0??2015-07-01?08:30??程序\ipch\
?????目錄???????????0??2015-07-01?08:30??程序\ipch\sy0-a4a86840\
?????文件????????4248??2011-10-30?13:54??程序\SY0.dsp
?????文件?????????529??2011-10-30?13:51??程序\SY0.dsw
?????文件???????58368??2012-05-26?11:01??程序\SY0.ncb
?????文件???????48640??2012-05-26?11:01??程序\SY0.opt
?????文件?????????880??2012-05-26?10:59??程序\SY0.plg
?????文件?????1789952??2015-07-01?14:12??程序\SY0.sdf
?????文件?????????872??2015-07-01?09:47??程序\SY0.sln
?????文件???????11776??2015-07-01?14:12??程序\SY0.suo
?????文件????????6084??2015-07-01?08:29??程序\SY0.vcxproj
............此處省略2個(gè)文件信息