完整word版单像空间后方交会程序报告.doc

单像空间后方交会程序报告

指导老师：刘老师
班级：测绘101
姓名：尚锋
学号：19号

应用程序的主入口部分的代码：

usingSystem;
usingSystem.Collections.Generic;
usingSystem.Linq;
usingSystem.Windows.Forms;
namespace单像空间后方交会
{
staticclassProgram
{
///<summary>
///应用程序的主入口点。
///</summary>
[STAThread]
staticvoidMain()
{
Application.EnableVisualStyles();
Application.SetCompatibleTextRenderingDefault(false);
Application.Run(newForm1());
}
}
}

方法解算类（通用）部分的代码：

usingSystem;
usingSystem.Collections.Generic;
usingSystem.Linq;
usingSystem.Text;
namespace单像空间后方交会
{
classTongyong
{
structimage_point//一个像点结构，包含像点坐标和地面点坐标
{
publicdoublex;
publicdoubley;
publicdoubleX;
publicdoubleY;
publicdoubleZ;
}
privatedoublef;//主距
privatedoubleu;//u为外方位元素，下面5个相同
privatedoublew;
privatedoublek;
privatedoubleXs;
privatedoubleYs;
privatedoubleZs;
privateimage_point[]p=newimage_point[4];//四个控制点
privatedouble[]R=newdouble[9];//旋转矩阵
privatedouble[]a=newdouble[8];//像点坐标近似值
privatedouble[,]A=newdouble[8,6];//误差方程式系数
privatedouble[]L=newdouble[8];//误差方程式常数项
privateintcount=0;//统计代次数
publicTongyong(doubleg,double[]q)//构造函数，初始化各变量,单位m
{
f=g;
for(inti=0;i<4;i++)
{
intj=i*5;
p[i].x=q[j];
p[i].y=q[j+1];
p[i].X=q[j+2];
p[i].Y=q[j+3];
p[i].Z=q[j+4];
}
doubleave=0,sum=0;//求比例尺分母
for(inti=0;i<3;i++)
{
for(intj=i+1;j<4;j++)
{
sum+=Math.Sqrt(Math.Pow(p[i].X-p[j].X,2)+Math.Pow(p[i].Y-p[j].Y,2))/Math.Sqrt(Math.Pow(p[i].x-p[j].x,2)+Math.Pow(p[i].y-p[j].y,2));
}
}
ave=sum/6;
u=0;//给定外方位元素的初始值,角度均设置为0
w=0;
k=0;
Xs=(p[0].X+p[1].X+p[2].X+p[3].X)/4;//Xs为四个控制点X的平均值，Ys类似
Ys=(p[0].Y+p[1].Y+p[2].Y+p[3].Y)/4;
Zs=(p[0].Z+p[1].Z+p[2].Z+p[3].Z)/4+ave*f;
}
privatedoublesin(doublem)//正弦，为简化而写,下同
{
returnMath.Sin(m);
}
privatedoublecos(doublem)
{
returnMath.Cos(m);
}
privatevoidcalcos()//计算旋转矩阵
{
R[0]=cos(u)*cos(k)-sin(u)*sin(w)*sin(k);
R[1]=-cos(u)*sin(k)-sin(u)*sin(w)*cos(k);
R[2]=-sin(u)*cos(w);
R[3]=cos(w)*sin(k);
R[4]=cos(w)*cos(k);
R[5]=-sin(w);
R[6]=sin(u)*cos(k)+cos(u)*sin(w)*sin(k);
R[7]=cos(u)*sin(w)*cos(k)-sin(u)*sin(k);
R[8]=cos(u)*cos(w);
}
privatevoid