文档介绍：《摄影测量学》实验报告 学期12-13学年第1学期专业地理信息系统班级1011学号1020209125姓名马明炎任课教师江标初苏州科技学院环境科学与工程学院空间信息与测绘工程系实验一摄影坐标系变换实验地点:C1—201日期:2012年9月19日目的通过影像外方位角元素q、w、k成功的完成像空间坐标与像辅助坐标之间的变换,熟悉并掌握两者之间的变换公式。完成这两种坐标系的相互转换。要求界面设计可以参考老师上课所给的界面。转换过程中,一个坐标系可以按某种顺序依次地旋转三个角度即可变换为另一个同原点的坐标系。通过VB程序来实现两种坐标系之间的转换。方法与详细步骤1、先在电脑里打开VB,,在界面上用lable、mand设置如下界面:2、如上界面所示,9像点坐标为(x,y,z),影像外方位元素中的角元素为(q,w,k),像空间辅助坐标为(X,Y,Z),根据书上第28页的公式:a1=Cos(a)*Cos(c)-Sin(a)*Sin(b)*Sin(c)a2=-Cos(a)*Sin(c)-Sin(a)*Sin(b)*Cos(c)a3=-Sin(a)*Cos(b)b1=Cos(b)*Sin(c)b2=Cos(b)*Cos(c)b3=-Sin(b)c1=Sin(a)*Cos(c)+Cos(a)*Sin(b)*Sin(c)c2=-Sin(a)*Sin(c)+Cos(a)*Sin(b)*Cos(c)c3=Cos(a)*Cos(b)(其中a=q/180*,b=w/180*,c=k/180*,因为我们输入text中的是实数,而不是度数,所以要进行转换)然后就可以求出我们所要求的像空间辅助坐标中所必要的元素。再根据书本第27页的公式:X=a1*x+a2*y+a3*zY=b1*x+b2*y+b3*zZ=c1*x+c2*y+c3*z这样就可以求出所要求的像空间辅助坐标了。3、按照以上所给的公式,用VB编写代码,代码内容如下:mand1_Click()DimaAsSingle,bAsSingle,cAsSingleDima1AsSingle,a2AsSingle,a3AsSingleDimb1AsSingle,b2AsSingle,b3AsSingleDimc1AsSingle,c2AsSingle,c3AsSinglea=*=*=*=Cos(a)*Cos(c)-Sin(a)*Sin(b)*Sin(c)a2=-Cos(a)*Sin(c)-Sin(a)*Sin(b)*Cos(c)a3=-Sin(a)*Cos(b)b1=Cos(b)*Sin(c)b2=Cos(b)*Cos(c)b3=-Sin(b)c1=Sin(a)*Cos(c)+Cos(a)*Sin(b)*Sin(c)c2=-Sin(a)*Sin(c)+Cos(a)*Sin(b)*Cos(c)c3=Cos(a)*Cos(b)=a1*+a2*+a3*=b1*+b2*+b3*=c1*+c2*+c3*()=""=""=""=""=""=""=""=""=""()EndEndSub4、运行程序代码,得出正确结果。实验成果实验二单幅影像空间后方交会实验地点C1—201日期9月26日一、目的掌握单幅影像的空间后方交会方法二、要求用VB进行编程,利用像片上三个以上控制点的像点坐标及相应的物方坐标,反算像片外方位元素,实现单幅影像的空间后方交会。三、方法与详细步骤1、将界面设计好后,并在相应的text内输入原始数据。如下图所示:2、计算和确定初值3、组成旋转矩阵R4、计算(x)、(y)和Lx,Ly5、逐点组成误差方程式并法化6、解法方程,求未知数改正数7、计算改正后的外方位元素8、检查未知数改正数是否小于限差9、输出最终结果四、实验成果1、mand1_Click()Dimx(4)AsDouble,y(4)AsDoubleDima(8,6)AsDoubleDimXd(4)AsDouble,Yd(4)AsDouble,Zd(4)AsDoubleDimnAsInteger,iAsIntegerDimXy(4)AsDoubl