实验一数字内定向编程实验 2
实验二空间后方交会编程实验 5
实验三空间前方交会编程实验 12
实验四同名点量测及地面点的cad成图实验 19
六、实验成果 20
掌握数字内定向的原理，利用计算机编程语言实现数字内定向的过程，完成一立体像对两幅像片影像的数字内定向。
计算机，编程软件（vb）
立体像对两像片）上四个框标点的理论位置以及四个框标点的像素坐标，具体如表3-1。
表3-1 4.1 内定向参数求取
方法：利用航摄像片上的四个框标点的理论位置以及四个框标点的像素坐标为依据。
根据条件：框标点的像素坐标通过公式转换后应等于其理论坐标。
利用式（3-1）列出方程，通过最小二乘法计算内定向参数： h0、 h1、 h2、 k0、k1、 k2。
计算公式：
x=h0 + h1*i + h2*j
y=k0 + k1*i + k2*j （3-1）
求解过程：
列误差方程法化求参数
4.2 将像点扫描坐标转化为框标坐标
通过公式：
x=h0 + h1*i + h2*j
y=k0 + k1*i + k2*j
将测量的同名点的扫描坐标转换位像框标坐标系中坐标。
private sub read_click() 读入数据。
dim i%, j%, k%
open for input as #
for i = 1 to 4
line input #1, s
s = split(s, "
for j = 1 to 2
a(i, j + 1) =s(j - 1)
next j
line input #1, s
s = split(s, "
lx(i, 1) =s(0)
ly(i, 1) =s(1)
next i
for i = 1 to 4
a(i, 1) =1
next i
end sub
内定向求解。
private sub jisuan_click()
dim i%, j%, k%
for i = 1 to 3
for j = 1 to 4
at(i, j) =a(j, i)
next j
next i
矩阵求解。jzc at(),a(),3, 4, 3, b()
rect b(),3, bdn()
jzc bdn(),at(),3, 3, 4, bb()
jzc bb(),lx(),3, 4, 1, hh()
jzc bb(),ly(),3, 4, 1, kk()
for i = 1 to 3
text1 = text1 & format(hh(i, 1), 0.000000") vbcrlf
next i
text1 = text1 & vbcrlf
for i = 1 to 3
text1 = text1 & format(kk(i, 1), 0.000000") vbcrlf
next i
close #
end sub
坐标转换。private sub command1_click()
xx = val(text2)
yy = val(text3)
text4 = hh(1, 1) +xx * hh(2, 1) +yy * hh(3, 1)
text5 = kk(1, 1) +xx * kk(2, 1) +yy * kk(3, 1)
= format(text4, "0.000")
= format(text5, "0.000")
end sub
1504定向参数：
x方向：h0=-121.766062,h1=0.027976,h2=0.000523,精度=0.006349
y方向：k0=115.125954,k1=0.000516,k2=-0.027982,精度=0.001338
1505定向参数：
x方向：h0=-121.642018,h1=0.027976,h2=0.000607,精度=0.005566
y方向：k0=115.321828,k1=0.000600,k2=-0.027981,精度=0.002853
掌握空间后方交会的原理，根据所给控制点的地面物方坐标以及相应的像点在像平面坐标系中的坐标，利用计算机编程语言实现空间后方交会的过程，完成所给立体像对）中两张像片各自的6个外方位元素的解求。
计算机，编程软件（vb）
实验数据为一立体像对）的像片控制点数据，包括像点坐标及相应地面摄测坐标。具体如表4-1所示：
表4-1两张像片的内方位元素为：x0=y0=0mm,f=210.681mm。
利用已知地面控制点数据以及相应像点坐标，编写空间后方交会程序，求解所给像对）中两张像片各自的6个外方位元素。
4.1 数学模型：共线条件方程式
4.2 求解过程：
1）获取已知数据：
a. 像片比例尺1/m：
由像片上两点间距离与相应地面点间距之比求得。
b. 平均航高h：
h=m*f+σzi/n
其中，n为已知控制点数。
c. 内方位元素x0,y0,f：
d. 地面摄影测量坐标xtp,ytp,ztp
e. 控制点的像点坐标（表4-1）
2） 确定未知数初始值：
在竖直摄影情况下，角元素的初始值为0，即φ=ω0；线元素中，zs0=h,xs0,ys0的取值可用五个已知控制点坐标的平均值，即：xs0=σxtpi/5，ys0=σytpi/5。
3） 计算旋转矩阵r：
利用角元素值的近似值按下式计算方向余弦值，组成r阵：
a1 = cos(φ)cos(κ)sin(φ)sin(ω)sin(κ)
a2 = cos(φ)sin(κ)sin(φ)sin(ω)cos(κ)
a3 = sin(φ)cos(ω)
b1 = cos(ω)sin(κ)
b2 = cos(ω)cos(κ)
b3 = sin(ω)
c1 = sin(φ)cos(κ)cos(φ)sin(ω)sin(κ)
c2 = sin(φ)sin(κ)cos(φ)sin(ω)cos(κ)
c3 = cos(φ)cos(ω)
4） 逐点计算像点坐标的近似值： 利用未知数的近似值按共线方程式计算控制点像点坐标的近似值。
x），（y）
5） 组成误差方程式：
按下式逐点计算误差方程式的系数和常数项：
系数： zba = a3 * xa - xs) +b3 * ya - ys) +c3 * za - zs))
a11 = a1 * f + a3 * x) /zba
a12 = b1 * f + b3 * x) /zba
a13 = c1 * f + c3 * x) /zba
a21 = a2 * f + a3 * y) /zba
a22 = b2 * f + b3 * y) /zba
a23 = c2 * f + c3 * y) /zba
a14= y * sin(ω)x * x * cos(κ)y * sin(κ)f + f * cos(κ)cos(ω)
摄影测量实验报告
学号 班级 姓名 随着摄影测量和遥感技术的蓬勃发展，摄影测量越来越显得重要。摄影测量在实际生产中也逐步占据了重要位置。我更加认识到摄影测量学不仅要有扎实的理论基础，还要有熟练的软件操作能力。对于我们来说，如果只有理论知识，而实践操作为零，那也只是纸上谈兵。只知道操作，不懂理论知识，那也不能独立完成工...
摄影测量基础实验报告
一 像对的立体观察。一 实验目的。1 了解人造立体视觉的原理。2 掌握观察人造立体的四个条件。3 掌握人造立体观察方法，借助仪器进行立体观察。二 实验内容。借助立体镜观察正立体 反立体 零立体。三 实验原理。用双眼观察空间物体时，可以容易地判断物体的远近，得到景物的立体效应，这种现象称为人眼的天然立...
摄影测量基础实验报告
一 像对的立体观察。一 实验目的。1 了解人造立体视觉的原理。2 掌握观察人造立体的四个条件。3 掌握人造立体观察方法，借助仪器进行立体观察。二 实验内容。借助立体镜观察正立体 反立体 零立体。一 像对的立体观察。一 实验目的。1 了解人造立体视觉的原理。2 掌握观察人造立体的四个条件。3 掌握人造...