第一篇:摄影测量学 作业一
现有一80KM*50KM的山地测区,要求地面像元分辨率达到40cm。设定航向重叠度为60%,旁向重叠度为35%。
1、选用摄影机,并依据所选的摄影机制定摄影方案,包括:
(1)所选摄影机的基本参数;
(2)确定航高;
(3)确定每条基线的长度;
(4)应布设几条航线?每条航线应布设多少条基线?
2、如果采用双像摄影测量(单模型)的方式来解算测区内加密点的坐标,给出一种解算方案及其基本原理,并说明控制点的布设方案。
要求:查阅资料,独立完成。温馨提示:请同学们重视“独立完成”的要求,讨论有益,拒绝抄袭!表达时,务请思路完整、清晰,并力求简洁。完成即随时上交,上交时限: 10月9日。
第二篇:摄影测量学作业三
1、简述低空摄影测量的发展背景及技术难点。
2、一个车载的移动测量系统(MMS)利用多传感器集成技术实现了快速的位置
服务,请分析MMS中各传感器系统的作用。
上交时间:下周(第九周)
周五。
第三篇:摄影测量学
摄影测量学:是对研究的对象进行摄影,根据所获得的构想信息,从几何方面和物理面加分析研究,从而对所摄影的对象本质提供各种资的一门学科。
航向重叠:供测图用的航测相片沿飞行方向上相邻像片的重叠。
4D产品:是指 DEM、DLG、DRG、DOM。
空间前方交会:由立体像对中两张像片的内、外方位元素和像点坐标来确定相应地面点的地面坐标的方法,称为空间前方交会。
点像空间后方交会:知道像片的内方位元素,以及三个地面点坐标和量测出的相应像点的坐标,就可以根据共线方程求出六个外方位元素的方法。
相对航高:摄影瞬间航摄飞机相对于某一索取基准面的高度。
相片纠正:将中心投影转换成正射投影时,经过投影变换来消除相片倾斜所引起的像点位移,使它相当于水平相片的构象,并符合所规定的比例尺的变换过程。
解析空中三角测量:是将建立的投影光束,单元模型或航带模型以及区域模型的数字模型,根据少数地面控制点,按最小二乘法原理进行平差计算,并求加密点地面坐标的方法。透视平面旋转定律:当物面和合面分别绕透视轴合线旋转后,只要旋转地角度相同,则投影射线总是通过物面和像面的统一相对应点。
外方位元素:用以确定摄影瞬间摄影机或像片空间位置,即摄影光束空间位置的数据。核面:通过摄影基线与任意物方点所作的平面称作通过该点的核面。
绝对定向元素:确定相对定向所建立的几何模型比例尺和恢复模型空间方位的元素。像主点:像片主光轴与像平面的交点。
立体像对:相邻摄站获取的具有一定重叠度的两张影像。
数字影像重采样:当欲知不位于采样点上的像素值时,需进行灰度重采样。
中心投影:所有投影光线均经过同一个投影中心。
摄影基线:相邻两摄站点之间的连线。
相对定向:恢复两张像片的相对位置和方位称为相对定向。
双像解析摄影测量:按照立体像对与被摄物体的几何关系,以数学计算方式,通过计算机解求被摄物体的三维空间坐标的方法,称为双像解析摄影测量。
xf
yfa1(XAXS)b1(YAYS)c1(ZAZS)a3(XAXS)b3(YAYS)c3(ZAZS)a2(XAXS)b2(YAYS)c2(ZAZS)
a3(XAXS)b3(YAYS)c3(ZAZS)
x,y为像点的框标坐标.x0,y0,f为影像的内方位元素.XS,YS,ZS为摄站点的物方空间坐标.XA,YA,ZA为物方点的物方空间坐标.ai,bi,ci为3个外方位角元素组成的9个方向余弦.
第四篇:摄影测量学VirtuoZo实习
摄影测量学VirtuoZo实习
一. 实习目的:
了解VirtuoZo NT系统的运行环境及软件模块的操作特点,了解实习工作流程,从而能对数字摄影测量实习有个整体概念。完成原始数字影像格式的转换。掌握创建/打开测区及测区参数文件的设置。掌握参数文件的数据录入。通过对模型定向的作业,了解数字影像立体模型的建立方法及全过程,并能较熟练地应用定向模块进行作业,满足定向的基本精度要求。掌握核线影像重采样,生成核线影像对。掌握匹配窗口及间隔的设置,运用匹配模块,完成影像匹配。掌握匹配后的基本编辑,能根据等视差曲线(立体观察)发现粗差,并对不可靠区域进行编辑,达到最基本的精度要求。掌握DEM格网间隔的正确设置,生成单模型的DEM。掌握正射影像分辨率的正确设置,制作单模型的数字正射影像。通过DEM及正射影像的显示,检查是否有粗差。掌握拼接区域的选定及确定拼接产品的路径。掌握DEM拼接及自动正射影像镶嵌。分析拼接精度。理解数据格式输出的意义。了解VirtuoZo NT系统的数据格式输出的具体操作。通过对实习成果的分析,了解数字产品的基本质量要求。总结实习中出现的问题以及实习成果的不足之处,并能分析其原因。
二、实习步骤:
1、数据准备:
在D盘准备好实习操作所需要的数据,包括images中的tif格式航飞图片,HammerIndex文件以及hammer.ctl控制点和rc30.cmr相机格式文件。2建立测区(打开测区):
新建一个测区,打开测区参数设置界面,分别进行测区参数的设置:主目录文件位置的确定,控制点文件,加密点,相机检校文件格式命名和输入,将DEM格网间隔设置为10。
3.设置相机参数
打开设置菜单下的相机参数设置,进行参数修改,从准备好的文件夹中输入rc30.cmr,在实际生产操作中,相机的参数是用户给定的。
4.设置控制点
打开设置菜单下的地面控制点,进行地面控制点的输入,引入hammer.ctl。5.引入影像
将文件中的影像资料通过设置——〉引入影像,并设置像素大小为0.045mm。
6.新建模型,设置模型左右影像及参数打开测区,根据处理的影像文件来进行命名,便于我们的识别,在此后面的操作都是以相片0-157和01-156进行操作,故输入模型名157_156(左相片名在左),进行模型参数的设置:分别在左影像和右影像中输入左右影像(vz格式); 在核线参数中选择水平核线(注意生成的产品目录文件所在的位置)。7.模型定向
内定向:点击处理目录下的模型定向下的内定向,将扫描坐标转换为像平面坐标。分别将框标进行移动,尽可能的使得十字丝在中心处,可适时查看当前的十字丝的中误差,不可一味追求中误差使得十字丝偏离中心。
相对定向:在内定向结束后,点击模型处理下面的自相对定向,右键选择自动相对定向,在定向结果中查看,删除中误差大于0.01的点。
绝对定向:在相对定向的基础上,对照给定的hmmerIndex网页文件选取控制点,找到控制点的大概位置后,在视图中进行粗调左右视图中控制点的位置,然后在右下角进行微调,使得控制点的中误差尽可能的小,左右视图中十字丝匹配。在当前视图下选择三个控制点之后,在进行一次自动相对定向,便可以将其余的控制点预测出来,把预测出来的控制点进行调整保存,退出。
内定向
绝对定向
8.生成核线影像
点击处理目录下的核线重采样进行核线采集。影像匹配:在核线采集之后,进行影像匹配。9.生成DEM 击产品下生成DEM中的DEM生成,在显示目录下的立体显示下的透视显示进行查看。
10.生成DOM,正射影像地图制作
点击产品目录下的生成正摄影像,在显示目录下的立体显示下的透视显示进行查看如下,可以看到图片中的房屋被当成地面进行来DEM格网生成,所以还需要进行DEM编辑,消除房屋上面的等高线影响。
11.DEM拼接
在系统主菜单中,选择菜单“镶嵌→设置”项,屏幕弹出拼接与镶嵌参数设置对话框。在系统主菜单中,选择“镶嵌→DEM拼接”项,进入DEM的拼接计算,屏幕弹出拼接进展显示条。当拼接完成后,将显示拼接中误差、总点数、误差分布统计及误差分布图。
三、实习总结
通过此次实习,了解了使用VirtuoZo 全数字摄影测量系统生产4D产品的过程,熟悉了VirtuoZo 全数字摄影测量系统的使用,加深了对相关知识的理解。
4D产品生产实习是一个综合性很强的实习,它是对所学摄影测量及相关专业的综合应用。该实习在数字摄影测量实习的基础上进行。
通过本次实习,了解到了VirtuoZo 全数字摄影测量系统的功能强大,在4d产品生产实习的过程中自动与半自动的快速生成功能。实习中需要注意:定义核线范围以将控制点划在作业区范围内为宜,但不能超控太多;其次应结合实际地形情况,如高山地或大比例城区,由于左右像片视差较大,就应适当将核线范围划大些。
单像空间后方交会程序
西南石油大学 土木工程与建筑学院 测绘工程 周凯强 学号:201308030143 输入文件形式如下:
C++源程序如下:
#include
//内方位元素 double m=39689;//估算比例尺
double B[4][5]={0.0},R[3][3],XG[6][1],AT[6][8],ATA[6][6],ATL[6][1];input(B,4,5);
//从文件中读取控制点的影像坐标和地面坐标,存入数组B double Xs=0.0, Ys=0.0, Zs=0.0,Q=0.0,W=0.0,K=0.0;
double X,Y,Z,L[8][1],A[8][6];
//确定未知数的出始值
for(int i=0;i<4;i++){Xs=Xs+B[i][2];
Ys=Ys+B[i][3];
Zs=Zs+B[i][4];} Xs=Xs/4;Ys=Ys/4;Zs=Zs/4+m*fk;int f=0;do//迭代计算
{f++;//组成旋转矩阵
R[0][0]=cos(Q)*cos(K)-sin(Q)*sin(W)*sin(K);
R[0][1]=-cos(Q)*sin(K)-sin(Q)*sin(W)*cos(K);
R[0][2]=-sin(Q)*cos(W);
R[1][0]=cos(W)*sin(K);
R[1][1]=cos(W)*cos(K);
R[1][2]=-sin(W);
R[2][0]=sin(Q)*cos(K)+cos(Q)*sin(W)*sin(K);
R[2][1]=-sin(Q)*sin(K)+cos(Q)*sin(W)*cos(K);
R[2][2]=cos(Q)*cos(W);
//计算系数阵和常数项
for(int i=0,k=0,j=0;i<=3;i++,k++,j++)
{
X=R[0][0]*(B[i][2]-Xs)+R[1][0]*(B[i][3]-Ys)+R[2][0]*(B[i][4]-Zs);
Y=R[0][1]*(B[i][2]-Xs)+R[1][1]*(B[i][3]-Ys)+R[2][1]*(B[i][4]-Zs);
Z=R[0][2]*(B[i][2]-Xs)+R[1][2]*(B[i][3]-Ys)+R[2][2]*(B[i][4]-Zs);
L[j][0]=B[i][0]-(x0-fk*X/Z);
L[j+1][0]=B[i][1]-(y0-fk*Y/Z);
j++;
A[k][0]=(R[0][0]*fk+R[0][2]*(B[i][0]-x0))/Z;
A[k][1]=(R[1][0]*fk+R[1][2]*(B[i][0]-x0))/Z;
A[k][2]=(R[2][0]*fk+R[2][2]*(B[i][0]-x0))/Z;A[k][3]=(B[i][1]-y0)*sin(W)-((B[i][0]-x0)*((B[i][0]-x0)*cos(K)-(B[i][1]-y0)*sin(K))/fk+fk*cos(K))*cos(W);A[k][4]=-fk*sin(K)-(B[i][0]-x0)*((B[i][0]-x0)*sin(K)+(B[i][1]-y0)*cos(K))/fk;
A[k][5]=B[i][1]-y0;
A[k+1][0]=(R[0][1]*fk+R[0][2]*(B[i][1]-y0))/Z;
A[k+1][1]=(R[1][1]*fk+R[1][2]*(B[i][1]-y0))/Z;
A[k+1][2]=(R[2][1]*fk+R[2][2]*(B[i][1]-y0))/Z;A[k+1][3]=-(B[i][0]-x0)*sin(W)-((B[i][1]-y0)*((B[i][0]-x0)*cos(K)-(B[i][1]-y0)*sin(K))/fk-fk*sin(K))*cos(W);A[k+1][4]=-fk*cos(K)-(B[i][1]-y0)*((B[i][0]-x0)*sin(K)+(B[i][1]-y0)*cos(K))/fk;
A[k+1][5]=-(B[i][0]-x0);
k++;} transpose(A,AT,6,8);multi(AT,A,ATA,6,8,6);inverse(ATA);multi(AT,L,ATL,6,8,1);multi(ATA,ATL,XG,6,6,1);Xs=Xs+XG[0][0];Ys=Ys+XG[1][0];Zs=Zs+XG[2][0];Q=Q+XG[3][0];W=W+XG[4][0];K=K+XG[5][0];}while(XG[3][0]>=6.0/206265.0||XG[4][0]>=6.0/206265.0||XG[5][0]>=6.0/206265.0);cout<<“迭代次数为:”< double AXG[8][1],V[8][1],VT[1][8],VTV[1][1],m0,D[6][6];multi(A,XG,AXG,8,6,1); for(i=0;i<8;i++) //计算改正数 V[i][0]=AXG[i][0]-L[i][0]; transpose(V,VT,1,8); multi(VT,V,VTV,1,8,1); m0=VTV[0][0]/2;for(i=0;i<6;i++) for(int j=0;j<6;j++) D[i][j]=m0*ATA[i][j];//屏幕输出误差方程系数阵、常数项、改正数 output(A,“误差方程系数阵A为:”,8,6);output(L,“常数项L为:”,8,1);output(XG,“改正数为:”,6,1);outFile.open(“aim.txt”,ios::app); //打开并添加aim.txt文件 outFile.precision(10);//以文件的形式输出像片外方位元素、旋转矩阵、方差阵 outFile<<“ 一、像片的外方位元素为:”< 二、旋转矩阵R为:”< outFile< outFile< 三、精度评定结果为:”< outFile< outFile< template for(i=0;i for(j=0;j mat2[j][i]=mat1[i][j]; return;} template int i,j,k;for(i=0;i {result[i][j]=0; for(k=0;k result[i][j]+=mat1[i][k]*mat2[k][j]; } } return;} template ifstream inFile;inFile.open(“控制点坐标.txt”);while(!inFile.eof()){for(int i=0;i for(int j=0;j inFile>>mat[i][j];} inFile.close();return;} template cout< cout< cout< double p; double q[n][12]; for(i=0;i for(j=0;j q[i][j]=c[i][j]; for(i=0;i for(j=n;j<12;j++) {if(i+6==j) q[i][j]=1; else q[i][j]=0;} for(h=k=0;k for(i=k+1;i {if(q[i][h]==0) continue; p=q[k][h]/q[i][h]; for(j=0;j<12;j++) { q[i][j]*=p; q[i][j]-=q[k][j]; } } for(h=k=n-1;k>0;k--,h--)// 消去对角线以上的数据 for(i=k-1;i>=0;i--){if(q[i][h]==0) continue; p=q[k][h]/q[i][h]; for(j=0;j<12;j++) {q[i][j]*=p; q[i][j]-=q[k][j];}} for(i=0;i q[i][j]*=p;} for(i=0;i c[i][j]=q[i][j+6];} 程序的结果输出如下:(包括文本输出结果和荧屏输出中间数据) 摄影测量学实习报告 为期两周的摄影测量学实习今天正式结束了,虽然两周时间并不长,但是对于我来说,学到的东西远不能用时间来衡量。在这两周里,我们完成了全数字摄影测量系统实习、数字影像分割程序编制、立体影像匹配程序编制等内容,这些东西让我们的两周很充实,很有意义。 其实刚开始时一直怀疑摄影测量学实习有什么意义,到了今天,我才发现这是有意义的。因为通过本次实习,我们可以将课堂理论与实践相结合,使我们深入掌握摄影测量学基本概念和原理,加强摄影测量学的基本技能训练,并且培养了我们的分析问题和解决实际问题的能力。通过使用数字摄影测量工作站,我们可以了解数字摄影测量的内定向、相对定向、绝对定向、测图过程及方法;通过开发数字影像分割程序和立体影像匹配程序,使自己掌握数字摄影测量基本方法与实现技术,为今后从事有关应用遥感技术应用和数字摄影测量打下坚实基础。所以,就算现在觉得没什么用,但是也为将来奠定了很好的基础。 正因为如此,在这两周中我们都很认真的在学习并且完成实习任务。其实说是两周,但时间真的更短,毕竟赶上了元旦假期,联欢晚会等一系列活动。所以如何在短暂的时间里,更出色的完成任务,是我们必须考虑的。记得实习动员的时候,老师花了很长时间又给我们讲了一次这次实习对我们的重要性,这很触动我们,毕竟老师的苦口婆心我们都看在眼里。不光如此,老师又耐心的把实习要求,实习任务,实习步骤讲解了一遍,让我们大致明白了这次实习从何入手,这 让本来很迷茫的我们瞬间找到了方向,也为我们接下来的工作提供了便利。动员结束的日子,我们便进入机房,正式开始了实习。 首先我们结束了全数字摄影测量系统,这款软件是我们从来没有接触过的,所以刚开始的时候很陌生,不知道怎么用,也不知道能用来做什么。还好,我们有老师的细致讲解,并且借助帮助向导可以解决我们很多问题。所以在这个实习中,我们没有遇到太多困难。让我印象深刻的是,我在做我们小组的绝对定向时,总是提示同名点数不够,就因为此,很难往下一步进行。后来在我们小组的讨论中,和老师的辅导后,我才得以解决这个困难。 第二周的时候,我们主要是利用Matlab进行程序的编写。因为之前的别的实习也要用到Matlab,所以对他已经不是很陌生了。但是当把Matlab和摄影测量的思路相结合的时候,还是出现了不少问题。毕竟摄影测量的原理也不是很容易理解的,加之需要利用计算机语言来实现程序就难上加难了。本来我想过放弃,因为编程实在是一件很麻烦的事。但在同组成员的鼓励下,以及老师的耐心讲解下,我还是坚持了下来,跟着我们小组一起商讨一起编写,虽然途中遇到了很多错误提示,遇到了很多无法实现的程序,但是我们都没有放弃,虚心的请教老师和同学,仔细的检查每一处错误,一一克服了这些问题。就这样,在磕磕绊绊中我们完成了立体影像匹配程序和立体影像匹配程序的编写。当看到最终执行出来的成果时,我们都很高兴,因为,这过程只有体会过的人才知道! 总的来说,这两周过得很难忘,毕竟这是这个学期最后一个实习,也是相当重要的一个实习。在这两周里,我们把平时课本上的知识又复习了一次,并且把它付诸于实践中。能把知识转化为技能是一个很好的过程。在这两周里,我体会到了老师的认真负责,如果没有老师的一遍一遍的耐心讲解,我们估计无法如此按时的完成任务。当然如果没有我们小组成员的通力合作,我也无法得到这么多财富。 感谢本学期的最后一次实习,因为在这过程中,我不光学到了知识,更体会了成长。这是多么宝贵的财富啊。摄影测量学实习真的是很难忘的两周。 王名洋测091第五篇:摄影测量学实习报告