资源描述：
《摄影测量编程实验报告.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、TONGJIUNIVERSITY摄影测量编程实验报告摄影测量编程实验项目名称孙海丽姓名061168学号土木工程学院测绘系学院（系）测绘工程专业叶勤指导教师2009年4月30日日期项目设计报告纸目录┊┊┊┊实验报告内容-------------------------------2┊┊数字内定向-------------------------2┊┊空间后方交会-----------------------4┊┊空间前方交会-----------------------8┊┊同名点量测成图--------------------10┊装心得体会---------
2、-------------------------12┊┊┊┊┊订┊┊┊┊┊线┊┊┊┊┊┊┊┊┊┊┊┊┊共12页第1页项目设计报告纸摄影测量编程实验报告及小结实验报告内容┊┊一实验目的：┊掌握数字内定向、空间后方交会、空间前方交会的原理，利用计算机编程语言实┊┊现双向解析计算的空间后交—前交法，并将计算结果利用AutoCAD进行简单成图。┊┊┊二、仪器用具及已知数据文件：┊┊计算机XP系统，编程软件（MATLAB），通用图像处理软件Photoshop，AutoCAD2004┊┊立体像对：左片1505.tif,右片：1504.tif┊装内定向数据文件：1505.
3、kb,1504.kb。内定向答案:1505photo.io,1504photo.io┊后方交会数据文件：1504.txt，1505.txt。后方交会答案（外定向参数）文件：┊┊1504_ao.txt，1505_ao.txt┊┊九对同名点坐标：tongmingdian.txt订┊┊三、实验内容：┊┊1.数字内定向┊线内定向的意思就是将含系统误差的框标坐标改正到正确的框标坐标，实验目的就┊┊是求内定向参数。试验中给出四个框标量测像素坐标和理论坐标，并要求用仿射变换┊公式：┊┊x=h0+h1*i+h2*j┊┊y=k0+k1*i+k2*j┊通过最小二乘平差求取内定向参
4、数：┊┊x方向：h0,h1,h2┊┊y方向：k0,k1,k2┊具体编程思路实现步骤：思路：因要求的是h0,h1,h2，k0,k1,k2，故可以将这六个待求参数当做误差方程中的X，则将得到一个有六个参数的误差方程，方程的系数A对于每个框标来说是：共12页第2页项目设计报告纸(1ij0000001ij)常数阵由x,y构成。即可平差求定内定向参数。┊┊具体步骤：┊┊1.数据文件的读取，我在读入的时候存入一个16*1维的数组，并分别赋给一个结构┊数组point，其具有四个属性：i,j,x,y。┊┊2.给系数矩阵Ｂ和常数项矩阵L赋值。（权阵为单位阵）┊┊3.利用平差知识
5、：X=inv(B'*B)*(B'*L)，可解求内定向参数。┊4.求改正数Ｖ，并做误差评定。┊┊5.输出结果。┊装实验数据：┊┊数据格式为：┊i,j┊┊x,y订┊具体数据见1505.kb,1504.kb┊平差得┊┊1505：┊线内方位元素：┊┊-121.642020.027980.00061115.321830.00060-0.02798┊平差后的x,y,坐标为：┊┊-109.99443110.00185┊┊-124.65253112.17801┊109.99343110.00315┊┊-118.49793112.30978┊┊精度：┊sx=0.004979sy
6、=0.0025521504：内方位元素：共12页第3页项目设计报告纸-121.766060.027980.00052115.125950.00052-0.02798平差后的x,y,坐标为：-109.99365109.99766┊┊-124.78592111.99113┊┊109.99265110.00734┊-118.63116112.10449┊┊精度：┊┊sx=0.005679sy=0.001197┊┊2.空间后方交会┊┊单幅影像的空间后方交会：利用已知地面控制点数据以及相应像点坐标，根据共装┊线方程，反求影像的外方位元素。┊数学模型：共线条件方程式┊┊a
7、(X−X)+b(Y−Y)+c(Z−Z)1s1s1s┊x=−fa(X−X)+b(Y−Y)+c(Z−Z)订3s3s3s┊┊a2(X−Xs)+b2(Y−Ys)+c2(Z−Zs)y=−f┊a(X−X)+b(Y−Y)+c(Z−Z)3s3s3s┊┊求解过程：线┊1）获取已知数据：┊┊a.像片比例尺1/m：┊┊由像片上两点间距离与相应地面点间距之比求得。┊b.平均航高H：┊┊H=m*f┊┊其中，n为已知控制点数。┊c.内方位元素x0,y0,f：┊┊d.地面摄影测量坐标Xtp,Ytp,Ztp2）量测控制点的像点坐标：像点坐标及相应地面摄测坐标如下表所示：数据格式：xyXtpY
8、tpZtp共12页第4页项目设计报告纸
