摄影测量内定向课程设计.docx
- 文档编号:9621104
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:12
- 大小:179.97KB
摄影测量内定向课程设计.docx
《摄影测量内定向课程设计.docx》由会员分享,可在线阅读,更多相关《摄影测量内定向课程设计.docx(12页珍藏版)》请在冰豆网上搜索。
摄影测量内定向课程设计
数字摄影测量内定向课程设计
一.课程设计的要求,目的及任务
1.摄影测量学基础课程设计要求
1)课程设计的具体项目和内容完全,而且要符合给定的课程设计大纲的规定。
2)设计的依据正确,鲜明表示出设计人员的主张,想法,和意见,涉及到的公式应该使用正确,推理符合逻辑,使实施方案和实际VB编程过程及最终实现。
3)每个设计阶段都做到阶段设计需求,阶段设计步骤,阶段设计检验方案。
4)认真编写技术设计说明书,在文字语言表达方面做到精准确信,简约,清楚,平实。
2.摄影测量学基础课程设计目的
1)总结和检验《摄影测量学》基本知识的学习情况。
通过学习的数字影像解析基础,完成一个较完整的设计计算过程,深化已有知识,拓宽新的知识,掌握影片内定相的基本过程和公式,加深对所学理论的理解与应用。
2)综合运用已学习的《误差理论与测量平差》、《visualbasic程序语言设计教程》,《数字测图原理与方法》等课程知识,选择主要理论与计算完成设计,达到综合性训练的目的。
3)加深对影片变形对内定向的影响,内定向实施方案的步骤与条件的认识及达到熟练运用的目的
4)锻炼和力求提高综合运用已学的基础理论知识和专业知识来解决摄影测量程序设计问题的初步技能,全面分析考虑问题的思想方法、工作方法以及计算、绘图和编写设计文件的能力。
进行科技写作的锻炼,培养科技论文实际写作能力。
5)将《摄影测量学》涉及到的大量的、零散的、独立的观点与资料,按照设计任务通知书的要求,经过演绎与归纳、分析与综合、抽象与具体的逻辑组合,而完成技术方案设计任务,达到锻炼和提高我们的逻辑思维能力和创造性思维能力的目的。
6)力求设计方案能够得到最优实施,设计实现阶段有序进行,程序运行成果确实可靠。
3.摄影测量学基础课程设计任务
1)根据数字影像解析基础和实际影片情况选择公式,进行计算。
2)对公式中所求的6个参数进行平差处理,得到改正后的数据,然后对其原计算数据作出改正。
3)运用VB程序语言编写影片内定向的程序过程代码。
4)实现VB过程,求出框标点在像平面坐标系和框标点在像素坐标系之间的关系,即求出6个参数值。
二.课程设计准备
1.设计思想
内定向是摄影测量测图的第一步。
1)了解什么是影片的内定向。
内定向:
由于数字影像的像素坐标系(扫描坐标系)是建立在像素矩阵之上的,其坐标原点在矩阵的左上角,坐标系也与相平面直角坐标系轴系不平行,为此必须建立像素坐标系和像平面直角坐标系之间的关系,这一过程叫做数字影像的内定向。
2)掌握内定向的目的。
内定向的目的是将像片纠正到像片坐标,通常方法是像片的周边有一系列的框标点,通常有4个或8个,它们的像片坐标是事先经过严格校正过的,利用这些点构成一个仿射变换的模型(或多项式),把象素纠正到像片坐标系。
通过这一步基本上消除了像片因扫描、压平等因素导致的变形。
2.所需已知数据和方法内容
所需已知数据:
影像数据和相机参数文件(含有相机类型、框标点的理论坐标、物镜畸变差等信息)。
内定向时,必须准确量测框标的像素坐标,再根据框标点的理论坐标,用解析的方法计算内定向参数。
如图所示,框标坐标系与像素坐标系的位置关系:
方法如下:
1)在模拟测图仪上测图时,首先应该将像片的框标与像片盘上的框标重合。
2)在解析测图仪上测图时,在将像片安放在像片盘上之后,首先要观测像片框标之坐标,使像片坐标系统与像片车架坐标系统联系在一起。
3.相关公式
1)仿射变换公式公式:
6个参数,量测4个框标
x=a0+a1x’+a2y’
y=b0+b1x’+b2y’
2)双线性公式为:
8个参数,量测8个框标
x=a0+a1x’+a2y’+a3x’y’
y=b0+b1x’+b2y’+b3x’y’
3)正形变换公式:
4个参数,量测3个框标
x=a0+a1x’-a2y’
y=b0+a2x’-a1y’
*说明:
(x’,y’)为框标的像素坐标,(x,y)是以像主点为原点的像平面直角坐标系,ai,bi为6个变换参数。
本设计只采用第一个仿射变换公式。
4.使用仪器工具
VB或其他编程软件
数字摄影测量系统软件适普软件VirtuoZo
三.课程设计步骤
1.理论步骤
1)选取定向的数学模型;
2)读取标准框标点在框标坐标系中的坐标;
3)人工量取(或自动量取)各框标点的屏幕坐标;
4)逐框标点列误差方程式,并用最小二乘原理求解各内定向参数;
5)将屏幕坐标(扫描坐标)转为框标坐标系。
2.程序设计步骤
1.准确量测出框标的像点坐标,并输入到相应的位置;
2.量测出框标坐标的像素坐标(扫描坐标);
3.根据仿射变换公式,列出方程式;
4.计算方程式中的6个参数,求出像空间直角坐标系和像素坐标系的转换关系。
主要程序代码如下:
通用
Dimc()AsSingle
Dimd()AsSingle
Dime()AsSingle
Dimf()AsSingle
DimL()AsSingle
DimU()AsSingle
PrivateSubCommand1_Click()
Dimx1,x2,x3,y1,y2,y3,y4AsSingle
DimI1,I2,I3,I4,J1,J2,J3,J4AsSingle
Dima0AsDouble
Dima1AsDouble
Dima2AsDouble
Dimb0AsDouble
Dimb1AsDouble
Dimb2AsDouble
x1=Val(Text1.Text)
x2=Val(Text5.Text)
x3=Val(Text7.Text)
x4=Val(Text9.Text)
y1=Val(Text2.Text)
y2=Val(Text6.Text)
y3=Val(Text8.Text)
y4=Val(Text10.Text)
I1=Val(Text4.Text)
I2=Val(Text11.Text)
I3=Val(Text13.Text)
I4=Val(Text15.Text)
J1=Val(Text3.Text)
J2=Val(Text12.Text)
J3=Val(Text14.Text)
J4=Val(Text16.Text)
a00=(x1+x2+x3+x4)/4'ai,bi未平差的值,由仿射变换公式得到
b00=(y1+y2+y3+y4)/4
a11=(x1*I1+x2*I2+x3*I3+x4*I4)/4
a22=(x1*J1+x2*J2+x3*J3+x4*J4)/4
b11=(I1*y1+I2*y2+I3*y3+I4*y4)/4
b22=(y1*J1+y2*J2+y3*J3+y4*J4)/4
ReDimf(i,j)AsSingle'矩阵B
Fori=1To8
Forj=1To6
f11=1:
f12=I1:
f13=J1:
f14=0:
f15=0:
f16=0
f21=0:
f22=0:
f23=0:
f24=1:
f25=I1:
f26=J1
f31=1:
f32=I2:
f33=J2:
f34=0:
f35=0:
f36=0
f41=0:
f42=0:
f43=0:
f44=1:
f45=I2:
f46=J2
f51=1:
f52=I3:
f53=J3:
f54=0:
f55=0:
f56=0
f61=0:
f62=0:
f63=0:
f64=1:
f65=I3:
f66=J3
f71=1:
f72=I4:
f73=J4:
f74=0:
f75=0:
f76=0
f81=0:
f82=0:
f83=0:
f84=1:
f85=I4:
f86=J4
Next
Next
ReDime(1To6,1)AsSingle
e11=a0:
e21=a1:
e31=a2:
e41=b0:
e51=b1:
e61=b2
x=I1+I2+I3+I4
y=J1+J2+J3+J4
p=I1^2+I2^2+I3^2+I4^2
q=J1^2+J2^2+J3^2+J4^2
z=I1*J1+I2*J2+I3*J3+I4*J4
ReDimc(m,n)AsSingle'矩阵B*B的转置的矩阵
Form=1To6
Forn=1To6
c11=4:
c12=x:
c13=y:
c14=0:
c15=0:
c16=0
c21=x:
c22=p:
c23=z:
c24=0:
c25=0:
c26=0
c31=y:
c32=z:
c33=q:
c34=0:
c35=0:
c36=0
c41=0:
c42=0:
c43=0:
c44=4:
c45=x:
c46=y
c51=0:
c52=0:
c53=0:
c54=x:
c55=p:
c56=z
c61=0:
c62=0:
c63=0:
c64=y:
c65=z:
c66=q
Next
Next
ReDimd(v,w)AsSingle'矩阵B的转置矩阵
Forv=1To6
Forw=1To8
d11=1:
d12=0:
d13=1:
d14=0:
d15=1:
d16=0:
d17=1:
d18=0
d21=I1:
d22=0:
d23=I2:
d24=0:
d25=I3:
d26=0:
d27=I4:
d28=0
d31=J1:
d32=0:
d33=J2:
d34=0:
d35=J3:
d36=0:
d37=J4:
d38=0
d41=0:
d42=1:
d43=0:
d44=1:
d45=0:
d46=1:
d47=0:
d48=1
d51=0:
d52=I1:
d53=0:
d54=I2:
d55=0:
d56=I3:
d57=0:
d58=I4
d61=0:
d62=J1:
d63=0:
d64=J2:
d65=0:
d66=J3:
d67=0:
d68=J4
Next
Next
ReDimL(8,1)AsSingle'矩阵L
l11=x1:
l21=y1:
l31=x2:
l41=y2
l51=x3:
l61=y3:
l71=x4:
l81=y4
ReDimU(1To6,1)AsSingle
U=-(c)^-1*d*L'U表示改正数△ai,△bi
Ua0=U11:
Ua1=U21:
Ua2=U31
Ub0=U41:
Ub1=U51:
Ub2=U61
Text23.Text=U11
Text24.Text=U21
Text25.Text=U31
Text26.Text=U41
Text27.Text=U51
Text28.Text=U61
a0=a00+U11'平差后的六个参数值
a1=a11+U21
a2=a22+U31
b0=b00+U41
b1=b11+U51
b2=b22+U61
Text17.Text=a0
Text18.Text=a1
Text19.Text=a2
Text20.Text=b0
Text21.Text=b1
Text22.Text=b2
EndSub
PrivateSubCommand2_Click()
Text1.Text=""
Text2.Text=""
Text3.Text=""
Text4.Text=""
Text5.Text=""
Text6.Text=""
Text7.Text=""
Text8.Text=""
Text9.Text=""
Text10.Text=""
Text11.Text=""
Text12.Text=""
Text13.Text=""
Text14.Text=""
Text15.Text=""
Text16.Text=""
Text17.Text=""
Text18.Text=""
Text19.Text=""
Text20.Text=""
Text21.Text=""
Text22.Text=""
EndSub
PrivateSubCommand3_Click()
UnloadMe
EndSub
*说明:
公式中用(xi’,yi’)表示框标坐标在像素坐标中的坐标值,在VB程序中用(I,J)表示框标坐标在像素坐标中的坐标值,因为在VB编程中,“’”具有说明解说作用,不能参与运行,故采用(I,J)表示。
四.程序实现成果
经过标准值和程序计算值的比较,算得的误差为0。
五.心得体会
通过本次的课程设计,我更加熟悉了数字影像的内定向的内容,原理,步骤,过程,同时也更加熟悉了VB编程的过程,对VB编程更加熟练。
通过仿射变换公式,在最小二乘法的原则基础上,设计了解决航空影像内定向参数解算的程序,输入框标点的框标坐标和影像坐标,输出仿射变换的6个参数。
在做该设计的时候,也发现了一些问题:
1.输入框标在像空间坐标系的坐标和框标在像素坐标系中的坐标时,由于数据精确,位数多,故该类数据均采用long型数据来保证精度。
2.每次调试时,输入的数据量过大
3.若是用焦点,可以减少输入数据的时间
本次问题已经经过了研究,在以后的实习实验和设计中将加以改正。
利
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 摄影 测量 内定 课程设计
![提示](https://static.bdocx.com/images/bang_tan.gif)