前方交会实验报告(含VB程序代码)

1立体像对前方交会实验报告一、实验目的在掌握前方交会原理的基础上，自己编写前方交会程序，在计算机上调试，输出计算结果并对计算结果进行检验。通过上机调试程序加强动手能力的培养，通过对实验过程的掌握以及对实验结果的分析，增强学生综合运用所学知识解决实际问题的能力。二、实验仪器计算机，VB6.0三、实验数据1.模拟像片一对：左片号23，右片号24；2.航摄机主距：f=150mm；3.左片23号片外方位元素：𝜑=−0°25′00″𝜔=−1°00′00″𝑘=−0°10′00″Xs=103007.006117Ys=139998.994849Zs=4801.9989994（m）右片24号片外方为元素：𝜑=1°39′59″𝜔=−0°10′00″𝑘=0°40′00″Xs=106002.023762Ys=140005.002780Zs=4797.009648（m）待求像点坐标如下表：四、实验内容利用所给立体像对两张像片的内、外方位元素，编写空间前方交会程序，根据所给像对中若干同名像点在左右像片上的坐标，解求其对应的地面点的物方坐标，实现空间前方交会的过程。五、实验成果程序流程图：23号片24号片点号X（mm）Y(mm)X(mm)Y(mm)50.8642.652-100.2511.443897.0842.933-4.4250.32041.082-76.227-102.695-79.61860.64082.003-99.90481.754795.207-75.521-5.334-78.443996.79783.77-3.50979.5872程序设计界面：开始分别写入两张相片的外方位元素写入两张像片同名像点坐标求解转矩阵计算像点像空间辅助坐标计算投影系数计算地面点地面摄影测量坐标求地面点坐标（Ｘ，Ｙ，Ｚ）结束3程序运行界面：4运行结果：（注：表上显示地面点坐标依次是：7,9,4,6,5,8）附：excel进行角度转换：六、程序如下：Dimf#,x1#,y1#,x2#,y2#,i%,j%,u1#,u2#,v1#,v2#,w1#,w2#,fai1#,kab1#,omg1#,fai2#,kab2#,omg2#Dima12#,a13#,b11#,b12#,b13#,c11#,c12#,c13#,a21#,a22#,a23#,b21#,b22#,b23#,c21#,c22#,c23#Dimn1#,n2#,bu#,bv#,bw#Dimxs1#,xs2#,ys1#,ys2#,zs1#,zs2#5PrivateSubCommand1_Click()fai1=Val(Text1.Text):kab1=Val(Text3.Text):omg1=Val(Text2.Text)fai2=Val(Text11.Text):kab2=Val(Text10.Text):omg2=Val(Text7.Text)f=Val(Text13.Text)xs1=Val(Text4.Text):xs2=Val(Text9.Text):ys1=Val(Text5.Text):ys2=Val(Text8.Text):zs1=Val(Text6.Text):zs2=Val(Text12.Text)EndSubPrivateSubCommand3_Click()x1=Val(InputBox(输入23片坐标x的值))y1=Val(InputBox(输入23片坐标y的值))x2=Val(InputBox(输入24片坐标x的值))y2=Val(InputBox(输入24片坐标y的值))EndSubPrivateSubCommand2_Click()Text1.Visible=FalseText2.Visible=FalseText3.Visible=FalseText4.Visible=FalseText5.Visible=FalseText6.Visible=FalseText7.Visible=FalseText8.Visible=FalseText9.Visible=FalseText10.Visible=FalseText11.Visible=FalseText12.Visible=FalseText13.Visible=FalseLabel1.Visible=FalseLabel2.Visible=FalseLabel3.Visible=FalseLabel4.Visible=FalseLabel5.Visible=FalseLabel6.Visible=FalseLabel7.Visible=FalseLabel8.Visible=FalseLabel9.Visible=FalseLabel10.Visible=FalseLabel11.Visible=FalseLabel12.Visible=FalseLabel13.Visible=FalseLabel14.Visible=Falsea11=Cos(fai1)*Cos(kab1)-Sin(fai1)*Sin(omg1)*Sin(kab1)a12=-Cos(fai1)*Sin(kab1)-Sin(fai1)*Sin(omg1)*Cos(kab1)a13=-Sin(fai1)*Cos(omg1)b11=Cos(omg1)*Sin(kab1)6b12=Cos(omg1)*Cos(kab1)b13=-Sin(omg1)c11=Sin(fai1)*Cos(kab1)+Cos(fai1)*Sin(omg1)*Sin(kab1)c12=-Sin(fai1)*Sin(kab1)+Cos(fai1)*Sin(omg1)*Cos(kab1)c13=Cos(fai1)*Cos(omg1)a21=Cos(fai2)*Cos(kab2)-Sin(fai2)*Sin(omg2)*Sin(kab2)a22=-Cos(fai2)*Sin(kab2)-Sin(fai2)*Sin(omg2)*Cos(kab2)a23=-Sin(fai2)*Cos(omg2)b21=Cos(omg2)*Sin(kab2)b22=Cos(omg2)*Cos(kab2)b23=-Sin(omg2)c21=Sin(fai2)*Cos(kab2)+Cos(fai2)*Sin(omg2)*Sin(kab2)c22=-Sin(fai2)*Sin(kab2)+Cos(fai2)*Sin(omg2)*Cos(kab2)c23=Cos(fai2)*Cos(omg2)u1=a11*x1+a12*y1-a13*fv1=b11*x1+b12*y1-b13*fw1=c11*x1+c12*y1-c13*fu2=a21*x2+a22*y2-a23*fv2=b21*x2+b22*y2-b23*fw2=c21*x2+c22*y2-c23*fbu=xs2-xs1bv=ys2-ys1bw=zs2-zs1n1=(bu*w2-bw*u2)/(u1*w2-u2*w1)n2=(bu*w1-bw*u1)/(u1*w2-u2*w1)Dimxx#,yy#,yy1#,yy2#,zz#xx=Fix((xs1+u1*n1)*1000+0.5)/1000yy1=Val(Text5.Text)+v1*n1yy2=Val(Text8.Text)+v2*n2yy=Fix(((yy1+yy2)/2)*1000+0.5)/1000zz=Fix((Val(Text6.Text)+w1*n1)*1000+0.5)/1000Form1.Print地面坐标为:Form1.Printxx=;xx,yy=;yy;zz=;zzEndSub七、实验心得通过本次前方交会实验，熟悉掌握了前方交会原理，并利用计算机程序对其进行了解决，收益颇多！