中南大学摄影测量空间前交-后交课程设计源代码

中南大学本科的摄影测量与遥感课程设计，空间前交-后交程序设计，附有详细的代码。

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

//求矩阵a的转置矩阵b，a为m行、n列
void transpose（double *a double *b int m int n）;

//矩阵a乘以矩阵b，结果存储在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;
	}
	
	//像点坐标和地面点坐标
	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-----------------------------左片外方位元素的计算---------------------------\n“）;

	double Zs1（0.0）Xs1（0.0）Ys1（0.0）p1（0.0）w1（0.0）k1（0.0）;
	//1.计算摄影中心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;
	//计算旋转矩阵R
	do       
	{
		
		R1[0][0]=cos（p1）*cos（k1）-sin（p1）*sin（w1）*sin（k1）;
		R1[0][1]=（-1）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           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个文件信息