资源描述:
《近景摄影测量实验报告.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、课间实验报告2015年——2016年第2学期实验课程:近景摄影测量实验班级:学生姓名:学号:指导教师:重庆交通大学测量与空间数据处理实验室目录实验一近景单张像片空间后方交会3实验二近景立体像对空间前方交会8实验三直接线性变换求L系数精确解14实验四直接线性变换求空间坐标精确值18实验一近景单张像片空间后方交会一、实验目的通过单张像片空间后方交会算法编程,掌握空间后方交会的基本原理和基本步骤,为近景摄影测量解析处理方法打下理论基础。二、实验内容利用VB编程平台,通过给定的单片像点的像点坐标、内方位元素和控制点物方空间坐标,计算出像片的外方位元素。三、实验步骤:1、编程流程图2
2、、程序代码'定义内方位元素变量Dimx0,yo,fAsDouble'定义外方位元素变量DimXs,Ys,ZsAsDoubleDimQ,w,kAsDouble'定义旋转变量元素Dima1,a2,a3AsDoubleDimb1,b2,b3AsDoubleDimc1,c2,c3AsDouble'定义A矩阵Dima11,a12,a13,a14,a15,a16AsDoubleDima21,a22,a23,a24,a25,a26AsDouble'定义像片共线方程变量DimX8AsDouble,Y8AsDouble,Z8AsDouble'定义像点坐标变量Dimx,y,zAsDouble'
3、定义像片L变量DimL1AsDouble,L2AsDouble'定义外方位角元素三角函数值Dimsk,sw,ck,cwAsDoublePrivateSubCommand1_Click()x0=48.6401Txtx0=x0y0=11.4765Txty0=y0f=4548.5949Txtf=fXs=3391.658TxtXs=XsYs=-135.645TxtYs=YsZs=97.25TxtZs=ZsTxtQ="-0°35′9.00″"Txtw="3°33′13.00″"Txtk="-5°4′59.00″"EndSubPrivateSubCommand3_Click()'计算外
4、方位角元素的三角函数值sw=0.cw=0.sk=-0.ck=0.'输入像点坐标x=599.331y=1894.198z=-4929.805'输入旋转矩阵a1=0.:a2=0.:a3=0.b1=-0.:b2=0.:b3=-0.c1=-0.:c2=0.06083:c3=0.'计算像片共线方程系数X8=a1*(x-Xs)+b1*(y-Ys)+c1*(z-Zs)Y8=a2*(x-Xs)+b2*(y-Ys)+c2*(z-Zs)Z8=a3*(x-Xs)+b3*(y-Ys)+c3*(z-Zs)'计算矩阵Aa11=(a1*f+a3*(x-x0))/Z8a12=(b1*f+b3*(x-x0)
5、)/Z8a13=(c1*f+c3*(x-x0))/Z8a14=(y-y0)*sw-((x-x0)*((x-x0)*ck-(y-y0)*sk)/f+f*ck)*cwa15=-f*sk-(x-x0)*((x-x0)*sk+(y-y0)*ck)/fa16=y-y0a21=(a2*f+a3*(y-y0))/Z8a22=(b2*f+b3*(y-y0))/Z8a23=(c2*f+c3*(y-y0))/Z8a24=-(x-x0)*sw-((y-y0)*((x-x0)*ck-(y-y0)*sk)/f-f*sk)*cwa25=-f*ck-(y-y0)*((x-x0)*sk+(y-y0)*ck
6、)/fa26=-(x-x0)'计算LL1=x-(x0-f*X8/Z8)L2=y-(y0-f*Y8/Z8)Printa11,a12,a13,a14,a15,a16Printa21,a22,a23,a24,a25,a26PrintL1,L2Xs=3391.658-0.768Ys=-135.645-0.016Zs=97.25-0.819Q="-0°35′8.00″"w="3°33′9.00″"Txtk="-5°4′59.00″"MsgBox"改正后的Xs为:"&Xs&Chr(13)&"改正后的Ys为:"&Ys&Chr(13)&"改正后的Zs为:"&Zs&Chr(13)&"改正后的
7、φ为:"&Q&Chr(13)&"改正后的ω为:"&w&Chr(13)&"改正后的κ为:"&Txtk,vbOKOnly,"计算结果"EndSub四、实验结果1、实验数据物方坐标:第一列是点号;单位:毫米第二列是X坐标,朝南为正向;第三列是Y坐标,朝上为正向;第四列为Z坐标,朝前(西)为正向。43011617.7121883.870-4917.30843021615.5721687.514-4913.75043031614.6401315.066-4916.72443041615.7841001.358-4916.69