資源簡介
中南大學(xué)本科的攝影測量與遙感課程設(shè)計(jì),空間前交-后交程序設(shè)計(jì),附有詳細(xì)的代碼。

代碼片段和文件信息
#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è)文件信息
- 上一篇:力控組態(tài)信號(hào)燈
- 下一篇:二維碼生成源代碼
評(píng)論
共有 條評(píng)論