摄影测量学实验报告
实验名称：空间后方交会
班级：测绘工程12-1
一．实验目的
掌握运用空间后方交会求解外方元素的过程。学会运用空间后方交会的原理，根据所给控制点的地面摄影测量坐标系坐标以及相应的像平面坐标系中的坐标，利用已知公式和计算机编程语言实现空间后方交会的过程，完成外方位元素的求解。
二．实验工具
个人电脑，visual basic 6.0
三．实验源代码
Private Sub Command1_Click()
Dim c(), kx(), ky(), l(), aa(), a11(), a12(), a13(), a14(), a15(), a16(), a21(), a22(), a23(), a24(), a25()
Dim a26(), v(), mm(5)
nn = Val(Text2.Text)
nnt = nn - 1
ReDim c(nnt, 5), kx(nnt), ky(nnt), l(2 * nnt + 1, 0), aa(2 * nnt + 1, 5), a11(nnt), a12(nnt), a13(nnt), a14(nnt), a15(nnt), a16(nnt), a21(nnt), a22(nnt), a23(nnt), a24(nnt), a25(nnt)
ReDim a26(nnt)
f = Val(Text1.Text)
fei = 0#
w = 0#
k = 0#
If Label2.Caption = "" Then
MsgBox "尚未读取文本", , "提示信息"
Else
a = Split(Label2.Caption, vbCrLf)
For i = 0 To nnt
hh = Split(a(i), " ")
For j = 0 To 5
c(i, j) = Val(hh(j))
Next j
Next i
m1 = (c(0, 1) - c(1, 1)) ^ 2 + (c(0, 2) - c(1, 2)) ^ 2 m2 = (c(0, 3) - c(1, 3)) ^ 2 + (c(0, 4) - c(1, 4)) ^ 2 m = Sqr(m2 / m1)
zs = m * f
xs = 0
ys = 0
For i = 0 To nnt
xs = xs + c(i, 3)
ys = xs + c(i, 4)
Next i
xs = xs / nn
ys = ys / nn
Do
a1 = Cos(fei) * Cos(k) - Sin(fei) * Sin(w) * Sin(k)
a2 = -Cos(fei) * Sin(k) - Sin(fei) * Sin(w) * Cos(k)
a3 = -Sin(fei) * Cos(w)
b1 = Cos(w) * Sin(k)
b2 = Cos(w) * Cos(k)
b3 = -Sin(w)
c1 = Sin(fei) * Cos(k) + Cos(fei) * Sin(w) * Sin(k)
c2 = -Sin(fei) * Sin(k) + Cos(fei) * Sin(w) * Cos(k)
c3 = Cos(fei) * Cos(w)
For i = 0 To nnt
kx(i) = -f * (a1 * (c(i, 3) - xs) + b1 * (c(i, 4) - ys) + c1 * (c(i, 5) - zs)) / (a3 * (c(i, 3) - xs) + b3 * (c(i, 4) - ys) + c3 * (c(i, 5) - zs))
ky(i) = -f * (a2 * (c(i, 3) - xs) + b2 * (c(i, 4) - ys) + c2 * (c(i, 5) - zs)) / (a3 * (c(i, 3) - xs) + b3 * (c(i, 4) - ys) + c3 * (c(i, 5) - zs))
Next i
For i = 0 To nnt
l(i * 2, 0) = c(i, 1) - kx(i)
l(i * 2 + 1, 0) = c(i, 2) - ky(i)
Next i
For i = 0 To nnt
zb = (a3 * (c(i, 3) - xs) + b3 * (c(i, 4) - ys) + c3 * (c(i, 5) - zs))
a11(i) = (a1 * f + a3 * c(i, 1)) / zb
a12(i) = (b1 * f + b3 * c(i, 1)) / zb
a13(i) = (c1 * f + c3 * c(i, 1)) / zb
a21(i) = (a2 * f + a3 * c(i, 2)) / zb
a22(i) = (b2 * f + b3 * c(i, 2)) / zb
a23(i) = (c2 * f + c3 * c(i, 2)) / zb
a14(i) = c(i, 2) * Sin(w) - (c(i, 1) * (c(i, 1) * Cos(k) - c(i, 2) * Sin(k)) / f + f * Cos(k)) * Cos(w)
a15(i) = -f * Sin(k) - c(i, 1) / f * (c(i, 1) * Sin(k) + c(i, 2) * Cos(k))
a16(i) = c(i, 2)
a24(i) = -c(i, 1) * Sin(w) - (c(i, 2) * (c(i, 1) * Cos(k) - c(i, 2) * Sin(k)) / f - f * Sin(k)) * Cos(w)
a25(i) = -f * Cos(k) - c(i, 2) / f * (c(i, 1) * Sin(k) + c(i, 2) * Cos(k))
a26(i) = -c(i, 1)
Next i
For i = 0 To nnt
aa(2 * i, 0) = a11(i)
aa(2 * i, 1) = a12(i)
aa(2 * i, 2) = a13(i)
aa(2 * i, 3) = a14(i)
aa(2 * i, 4) = a15(i)
aa(2 * i, 5) = a16(i)
aa(2 * i + 1, 0) = a21(i)
aa(2 * i + 1, 1) = a22(i)
aa(2 * i + 1, 2) = a23(i)
aa(2 * i + 1, 3) = a24(i)
aa(2 * i + 1, 4) = a25(i)
aa(2 * i + 1, 5) = a26(i)
Next i
ReDim at(5, 2 * nnt + 1), ata(5, 5), nata(5, 5), nataat(5, 2 * nnt + 1), detx(5, 0)
at = zhuanzhi(aa)
ata = cheng(at, aa)
nata = qiuni(ata)
nataat = cheng(nata, at)
detx = cheng(nataat, l)
xs = xs + detx(0, 0)
ys = ys + detx(1, 0)
zs = zs + detx(2, 0)
fei = fei + detx(3, 0)
w = w + detx(4, 0)
k = k + detx(5, 0)
Loop Until Abs(detx(0, 0)) < 0.000001 And Abs(detx(1, 0)) < 0.000001 And Abs(detx(2, 0)) < 0.000001 And Abs(detx(3, 0)) <
0.000001 And Abs(detx(4, 0)) < 0.000001 And Abs(detx(5, 0)) < 0.000001
ReDim v(2 * nnt + 1, 0), ax(2 * nnt + 1, 0)
ax = cheng(aa, detx)
v = jian(ax, l)
vv = 0
For i = 0 To 2 * nnt + 1
vv = v(i, 0) ^ 2 + vv
Next i
mo = Sqr(vv / (2 * nn - 6))
For i = 0 To 5
mm(i) = mo * Sqr(nata(i, i))
Next i
txtshow.Text = txtshow.Text & "外方元素为中误差为" & vbCrLf & "Xs " & xs & " " & mm(0) & vbCrLf & "Ys " & ys & " " & mm(1) & vbCrLf & "Zs " & zs & " " & mm(2) & vbCrLf & "Φ" & fei & " " & mm(3) & vbCrLf & "W " & w & " " & mm(4) & vbCrLf & "K " & k & " " & mm(5) & vbCrLf & vbCrLf
txtshow.Text = txtshow.Text & " 报告时间：" & Year(Date) & Format(Month(Date), "00") & Format(Day(Date), "00") & vbCrLf
txtshow.Text = txtshow.Text & " 审批人签字："
End If
End Sub
Private Sub Command2_Click()
CommonDialog1.Filter = "文档|*.txt"
CommonDialog1.ShowOpen
Label1.Caption = CommonDialog1.FileName
End Sub
Private Sub Command3_Click()
Label2.Caption = ""
Dim myline
If Label1.Caption <> "" Then
Open Label1.Caption For Input As #1
Do While Not EOF(1)
Line Input #1, myline
Label2.Caption = Label2.Caption + myline + vbCrLf
Loop
Close #1
End Sub
Private Sub Command4_Click()
CommonDialog1.DialogTitle = "保存文件" CommonDialog1.Filter = "文档|*.doc" CommonDialog1.Action = 2
If CommonDialog1.FileName <> "" Then
Open CommonDialog1.FileName For Output As #1
Print #1, txtshow.Text
Close #1
End If
End Sub
Private Sub Command5_Click()
txtshow.Text = ""
End Sub
Private Sub Command6_Click()
If MsgBox("确认退出程序？", 68, "提示信息") = 6 Then End
End If
End Sub
Public Function cheng(a(), b())
ReDim s(UBound(a, 1), UBound(b, 2))
If UBound(a, 2) <> UBound(b, 1) Then MsgBox "不符合乘法要求"
End
Else
For i = 0 To UBound(s, 1)
For j = 0 To UBound(s, 2)
s(i, j) = 0
For Q = 0 To UBound(a, 2)
s(i, j) = Val(s(i, j)) + Val(a(i, Q)) * Val(b(Q, j)) Next Q
Next j
Next i
cheng = s
End If
End Function
Public Function zhuanzhi(a())
Dim s
ReDim s(UBound(a, 2), UBound(a, 1))
For i = 0 To UBound(a, 1)
For j = 0 To UBound(a, 2)
s(j, i) = a(i, j)
Next j
Next i
zhuanzhi = s
End Function
Public Function qiuni(m())
Dim s()
If UBound(m, 1) <> UBound(m, 2) Then MsgBox ("求逆错误")
End
Else
n = UBound(m, 1)
Dim e(), m1, m2, m3
ReDim e(n, 2 * n + 1)
ReDim s(n, 2 * n + 1)
For i = 0 To n
For j = 0 To n
e(i, j) = m(i, j)
Next j
Next i
For i = 0 To n
For j = n + 1 To 2 * n + 1
e(i, j) = 0
Next j
Next i
For i = 0 To n
e(i, n + i + 1) = 1 Next i
For t = 0 To n
If e(t, t) = 0 Then
For i = t + 1 To n
For j = 0 To 2 * n + 1 s(i, j) = e(i, j)
e(i, j) = e(t, j)
e(t, j) = s(t, j)
Next j
If e(t, t) <> 0 Then GoTo daima1
End If
Next i
If e(t, t) = 0 Then MsgBox ("求逆错误") End
GoTo lastline
End If
daima1:
m1 = e(t, t)
For j = 0 To 2 * n + 1
e(t, j) = e(t, j) / m1
Next j
For i = t + 1 To n
m2 = e(i, t)
For j = 0 To 2 * n + 1
e(i, j) = e(i, j) - m2 * e(t, j) Next j
Next i
Next t
For t = 0 To n - 1
For i = t + 1 To n
m3 = e(t, i)
For j = 0 To 2 * n + 1
e(t, j) = e(t, j) - m3 * e(i, j) Next j
Next i
Next t
ReDim c(n, n)
For i = 0 To n
For j = 0 To n
c(i, j) = Round(e(i, j + n + 1), 2)
Next j
Next i
qiuni = c
End If
lastline:
End Function
Public Function jian(a(), b())
Dim s, t
ReDim s(UBound(a, 1), UBound(a, 2))
If UBound(a, 1) <> UBound(b, 1) Or UBound(a, 2) <> UBound(b, 2) Then
MsgBox "不符合减法要求"
End
Else
For i = 0 To UBound(a, 1)
For j = 0 To UBound(a, 2)
s(i, j) = Val(a(i, j)) - Val(b(i, j))
Next j
Next i
jian = s
End If
End Function
Private Sub Form_Load()
Form1.Image1.Top = 0
Form1.Image1.Left = 0
Form1.Image1.Width = Form1.Width Form1.Image1.Height = Form1.Height End Sub
Private Sub Form_Resize()
Form1.Image1.Top = 0
Form1.Image1.Left = 0
Form1.Image1.Width = Form1.Width Form1.Image1.Height = Form1.Height End Sub
四．实验框图
获取已知数据
量测控制点的像点坐标
确定未知数初始值
五．实验数据
1 -0.08615 -0.06899 36589.41 25273.3
2 2195.17
2 -0.05340 0.08221 37631.08 31324.51 728.69
3 -0.01478 -0.07663 39100.97 24934.98 2386.5
4 0.01046 0.06443 40426.54 30319.81 757.31 六．实验截图
七．实验心得
此次实验让我更加了解空间后方交会的计算过程，加强了对空间后方交会基本公式和误差方程式，法线方程式的记忆。需要注意的是，矩阵进行运算的时候应注意定义其维数和大小，要不然会出现错误；注意摄影机主距f，像点坐标和地面坐标之间的单位统一，并要自行计算比例尺m的大小；迭代的时候应注意将初值放在循环外，而且在限差比较时注意要加上绝对值再比较，而不是单纯的比较。通过本次的实验，我更加了解了如何对摄影测量算法进行编程，为以后的摄影
测量学的学习打下了基础。
空间后方交会编程实习报告 一实习目的 用程序设计语言（Visual C++或者C语言）编写一个完整的单片空间后方交会程序，通过对提供的试验数据进行计算，输出像片的外方位元素并评定精度。本实验的目的在于让学生深入理解单片空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。通过上机调试程序加强动手能力的培养，通过对实验结果的分析，增强学生综合运用所学知识解决实际问题的能力。 二实习内容 利用一定数量的地面控制点，根据共线条件方程求解像片外方位元素。 三实习数据 已知航摄仪的内方位元素：f k =153.24mm，x =y =0.0mm，摄影比例尺为1:50000； 4个地面控制点的地面坐标及其对应像点的像片坐标： 四实习原理 如果我们知道每幅影像的6个外方位元素，就能确定被摄物体与航摄影像的关系。因此，如何获取影像的外方位元素，一直是摄影测量工作者所探讨的问题。可采取的方法有：利用雷达、全球定位系统（GPS）、惯性导航系统（INS）以及星相摄影机来获取影像的外方位元素；也可以利用影像覆盖范围内一定数量的控制点的空间坐标与摄影坐标，根据共线条件方程，反求该影像的外方位元素，这种方法称为单幅影像的空间后方交会。 单像空间后方交会的基本思想是：以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的外方位元素Xs，Ys，Zs，t，w，k。 五实习流程 （1）获取已知数据。从摄影资料中查取影像比例尺1/m，平均摄影距离（航空摄影的航高、内方位元素x0，y0，f；获取控制点的空间坐标Xt，Yt，Zt。 （2）量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。 （3）确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值，在竖直航空摄影且地面控制点大体对称分布的情况下，可按如下方法确定初始值：
空间后方交会的解算 一． 空间后方交会的目的 摄影测量主要利用摄影的方法获取地面的信息，主要是是点位信息，属性信息，因此要对此进行空间定位和建模，并首先确定模型的参数，这就是空间后方交会的目的，用以求出模型外方位元素。 二． 空间后方交会的原理 空间后方交会的原理是共线方程。 共线方程是依据相似三角形原理给出的，其形式如下 111333222333()()() ()()() ()()()()()()A S A S A S A S A S A S A S A S A S A S A S A S a X X b Y Y c Z Z x f a X X a Y Y a Z Z a X X b Y Y c Z Z y f a X X a Y Y a Z Z -+-+-=--+-+--+-+-=--+-+- 上式成为中心投影的构线方程， 我们可以根据几个已知点，来计算方程的参数，一般需要六个方程，或者要三个点，为提高精度，可存在多余观测，然后利用最小二乘求其最小二乘解。 将公式利用泰勒公式线性化，取至一次项，得到其系数矩阵A ；引入改正数（残差）V ，则可将其写成矩阵形式： V AX L =- 其中 111333222333[,]()()()()()()()()()()()()()()T x y A S A S A S x A S A S A S A S A S A S y A S A S A S L l l a X X b Y Y c Z Z l x x x f a X X a Y Y a Z Z a X X b Y Y c Z Z l y y y f a X X a Y Y a Z Z =-+-+-=-=+-+-+--+-+-=-=+-+-+- 则1()T T X A A A L -= X 为外方位元素的近似改正数， 由于采用泰勒展开取至一次项，为减少误差，要将的出的值作为近似值进行迭代，知道小于规定的误差 三． 空间后方交会解算过程 1. 已知条件 近似垂直摄影
全站仪后方交会法步骤和 高程测量步骤 Revised final draft November 26, 2020
1、角度测量（angleobservation） （1）功能：可进行水平角、竖直角的测量。 （2）方法：与经纬仪相同，若要测出水平角∠AOB，则： 1）当精度要求不高时： 瞄准A点——置零（0SET）——瞄准B点，记下水平度盘HR的大小。 2）当精度要求高时：——可用测回法（methodofobservationset）。 操作步骤同用经纬仪操作一样，只是配置度盘时，按“置盘”（HSET）。 2、距离测量（distancemeasurement） PSM、PPM的设置——测距、测坐标、放样前。 1）棱镜常数（PSM）的设置。 一般：PRISM=0（原配棱镜），-30mm（国产棱镜） 2）大气改正数（PPM）（乘常数）的设置。 输入测量时的气温（TEMP）、气压（PRESS），或经计算后，输入PPM的值。 （1）功能：可测量平距HD、高差VD和斜距SD（全站仪镜点至棱镜镜点间高差及斜距） （2）方法：照准棱镜点，按“测量”（MEAS）。 3、坐标测量（coordinatemeasurement） （1）功能：可测量目标点的三维坐标（X，Y，H）。 （2）测量原理任意架仪器，先设置仪器高为0，棱镜高是多少就是多少，棱镜拿去直接放在已知点上测高差，测得的高差为棱镜头到仪器视线的高差，当然，有正有负了，然后拿出计算器用已
知点加上棱镜高，再加上或减去（因为有正有负）测得的高差就是仪器的视线高啊，因为仪器高为0，所以这个数字就是你的测站点高程，进测站点把它改成这个数字就行了，改完测站点了一般情况下都要打一下已知点复核一下。。。 若输入：方位角，测站坐标（，）；测得：水平角和平距。则有： 方位角： 坐标： 若输入：测站S高程，测得：仪器高i，棱镜高v，平距，竖直角，则有： 高程： （3）方法： 输入测站S（X，Y，H），仪器高i，棱镜高v——瞄准后视点B，将水平度盘读数设置为——瞄准目标棱镜点T，按“测量”，即可显示点T的三维坐标。 4、点位放样(Layout) （1）功能：根据设计的待放样点P的坐标，在实地标出P点的平面位置及填挖高度。 （2）放样原理 1）在大致位置立棱镜，测出当前位置的坐标。 2）将当前坐标与待放样点的坐标相比较，得距离差值dD和角度差dHR或纵向差值ΔX和横向差值ΔY。 3）根据显示的dD、dHR或ΔX、ΔY，逐渐找到放样点的位置。
单张相片后方交会
目录 ●作业任务 (3) ●解算原理 (3) ●具体过程 (4) ●算法描述及程序流程 (4) ●计算结果 (7) ●结果分析 (8) ●心得体会及建议 (8) ●参考文献 (9)
一，作业任务 已知摄影机主距f=153.24mm，四对点的像点坐标与相应地面坐标列入下表： 表1-1 计算近似垂直摄影情况下后方交会解。 二，解算原理 【关键词1】中心投影构像方程 在摄影测量学中，最重要的方程就是中心投影构像方程（图2-1）。这个方程 将地面点在地面摄影测量坐标系中的坐标（物方坐标）和地面点对应像点的像平 面坐标联系起来。在解析摄影测量与数字摄影测量中是极其有用的。在以后将要 学习到的双像摄影测量光束法、解析测图仪原理及数字影像纠正等都要用到该 式。 图2-1 在上述公式中：x和y分别为以像主点为原点的像点坐标，相应地面点坐标 为X,Y,Z，相片主距f以及外方位元素Xs，Ys，Zs，ψ，ω，κ。 而在此次作业中，就是已知四个地面控制点的坐标以及其对应的像点坐标， 通过间接平差原理来求解此张航片的外方位元素。 【关键词2】间接平差 在一个平差问题中，当所选的独立参数X的个数等于必要观测值t时，可将 每个观测值表达成这t个参数的函数，组成观测方程，然后依据最小二乘原理求 解，这种以观测方程为函数模型的平差方法，就是间接平差方法 间接平差的函数模型为： 随机模型为： 平差准则为：VtPV=min 【关键词3】单像空间后方交会 利用至少三个已知地面控制点的坐标A（Xa，Ya，Za）、B（Xb，Yb，Zb）、Z（Xc，
Yc，Zc），与其影像上对应的三个像点的影像坐标a（xa，ya）、b（xb，yb）、c（xc，yc），根据共线方程，反求该像点的外方位元素Xs，Ys，Zs，ψ，ω，κ。这种解算方法是以单张像片为基础，亦称单像空间后方交会。 在此次作业中，就是已知四个控制点在地面摄影测量坐标系中的坐标和对应的像点坐标。由此可以列出8个误差方程，存在两个多余观测数，则n=2。故可利用间接平差里，依据最小二乘法则，进行求解。由于共线条件方程是非线性函数模型，为了便于计算，需把非线性函数表达式用泰勒公式展开成现行形式，即“线性化”。而又因为仅取泰勒级数展开式的一次项，未知数的近似值改正是粗略的，所以必须计算采用逐渐趋近法，解求过程需要反复趋近，直至改正值小于限差为止。 三，具体过程 1，获取已知点数据：从摄影资料中查取像片比例尺1/m，平均航高，内方元素x0，y0，f；从外业测量成果中，获取控制点的地面测量坐标Xt，Yt，Zt，并转换成摄影测量坐标X,Y,Z。 2，量测控制点的像点坐标：将控制点标刺在像片上，利用立体坐标量测仪量测控制点的像框坐标，并经像点坐标改正，得到像点坐标x，y。 3，确定未知数的初始值：在竖直摄影测量情况下，角元素的初始值为0，及ψ=ω=κ=0； 线元素中，Zso =m*f+(Z[0]+Z[1]+Z[2]+Z[3])/4，Xso，Yso的取值可用四个角点上制点坐标的平均值，即：Xso=(X[0]+X[1]+X[2]+X[3])/4;Yso=(Y[0]+Y[1]+Y[2]+Y[3])/4；4，计算旋转矩阵R：利用角元素的近似值计算方向余弦，组成R阵。公式如下：R[0][0]=cos(ψ)*cos(k)-sin(ψ)*sin(w)*sin(k); R[0][1]=-cos(ψ)*sin(k)-sin(ψ)*sin(w)*cos(k); R[0][2]=-sin(ψ)*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(ψ)*cos(k)+cos(ψ)*sin(w)*sin(k); R[2][1]=-sin(ψ)*sin(k)+cos(ψ)*sin(w)*cos(k); R[2][2]=cos(ψ)*cos(w); 5，逐点计算像点坐标的近似值：利用未知数的近似值按共线方程计算控制点像点坐标的近似值（x）、（y）； 6，组成误差方程式：参照教材（5-8）式、（5-9b）式、（5-4）式逐点计算误差方程的系数阵和常数项。 7，组成法方程：计算法方程的系数矩阵与常数项。 8，解求外方位元素：根据法方程，按间接平差原理解求外方位元素改正值，并与相应的近似值求和，得到外方位元素的新的近似值。 9，检查计算是否收敛：将求得的外方位元素的改正值与规定的限差比较，小于限差则计算终止，否则用新的近似值重复第4至第8步骤计算，直至满足要求为止。 四，算法描述及程序流程。 算法描述（图4-1）：
单像空间后方交会和双像解析空间后方-前 方交会的算法程序实现 遥感科学与技术 摘要：如果已知每张像片的6个外方位元素，就能确定被摄物体与航摄像片的关系。因此，利用单像空间后方交会的方法，可以迅速的算出每张像片的6个外方位元素。而前方交会的计算，可以算出像片上点对应于地面点的三维坐标。基于这两点，利用计算机强大的运算能力，可以代替人脑快速的完成复杂的计算过程。 关键词：后方交会，前方交会，外方位元素，C++编程 0.引言： 单张像片空间后方交会是摄影测量基本问题之一,是由若干控制点及其相应像点坐标求解摄站参数(X S,Y S,ZS,ψ、ω、κ)。单像空间后方交会主要有三种方法:基于共线条件方程的平差解法、角锥法、基于直接线性变换的解法。而本文将介绍第一种方法，基于共线条件方程反求象片的外方位元素。 而空间前方交会先以单张像片为单位进行空间后方交会，分别求出两张像片的外方位元素，再根据待定点的一对像点坐标，用空间前方交会的方法求解待定点的地面坐标。可以说，这种求解地面点的坐标的方法是以单张像片空间后方交会为基础的，因此，单张像片空间后方交会成为解决这两个问题以及算法程序实现的关键。
1.单像空间后方交会的算法程序实现： （1）空间后方交会的基本原理：对于遥感影像,如何获取像片的外方位元素,一直是摄影测量工作者探讨的问题,其方法有:利用雷达(Radar)、全球定位系统(GPS)、惯性导航系统(I N S)以及星像摄影机来获取像片的外方位元素;也可以利用一定数量的地面控制点,根据共线方程,反求像片的外方位元素,这种方法称为单像空间后方交会(如图1所示)。 图中,地面坐标X i、Yi、Zi和对应的像点坐标x i、yi是已知的,外方位元素XS、Y S、ZS(摄站点坐标),ψ、ω、κ(像片姿态角)是待求的。 （2）空间后方交会数学模型：空间后方交会的数学模型是共线方程, 即中心投影的构像方程： 式中X、Y、Z是地面某点在地面摄影测量坐标系中的坐标,x,y是该地面点在像片上的构像点的像片坐标,对 于空间后方交会而言它们是已知的,还有主距f是已知的。而9个方向余弦a 1,a 2,a3;b1,b 2,b 3;c 1,c2,c 3是未知的，具体表达式可以取
全站仪放样，作为施工过程中一项重要环节，对技术员已上升为必须擅长的仪器操作内容。全站仪建站一般有两种方法，即极坐标法建站和后方交会法建站。现以尼康全站仪为例，讲述全站仪后方交会法建站、放样全过程。（其他品牌全站仪可参考进行） 一、建站 1.将仪器架于两已知点均可通视，且可完全看到放样目标点位置的高处。尽量保证视线夹角在60度左右，仪器架设高度适中，三脚架腿踩实，不可出现放样过程中架腿松动现象。（注意：整个放样过程中仪器附近不应有人来回走动，且放样人员应尽量站在一点不动，减少因人员走动导致仪器震动偏移。） 2.固定仪器，上下松动架腿大致调整圆水准器气泡基本居中，按下电源键开机，上下左右转动一下，按下“0”键，进入精平模式。 将水准管放于平行于两螺旋连线方向，关注屏幕上数值，“”过大，便同时向内或向外转动平行方向两螺旋至数值符合要求（一般数值处于5＂以内即可）；“”过大，便左转或右转垂直方向螺旋至数值符合要求。旋转60度，检查，若仍有些许偏差，再按上述调整。再旋转60度继续检查至完成。 3.按下“确定”键记录，按“建站”键进入建站模式，选择“后方交会法”按“确定”。①若全站仪内已有建站点坐标，可在“PT”栏输入点名（“MODE”键可切换数字与字母），按“确定”键自动跳出坐标，再输入棱镜高（本项目为1.35m和1.2m两种）；②若全站仪内无建站点坐标，于“PT”处按“确定”键进入坐标输入界面，XYZ
输完后，按“确定”回到界面，再输入仪器高。 CD数值暂时不输，按“确定”跳过进而记录，进入瞄准后视点1界面，视线内横竖丝卡住棱镜头“横竖尖头”（一般要求：竖向从镜杆底部瞄起，再翻转上去；横向以卡住两边尖为准），瞄准后，点击“测量1”（一般仪器内部设置“测量1”为棱镜模式且双频，“测量2”为免棱镜模式且单频，具体设置可内部调节变动）测量，待响两声后，在不转动仪器前按“确定”键记录，重复“PT”输入点坐标和棱镜高进行后视点2的瞄准，按“测量1”测量（若发现测量时后视瞄准有移动，再瞄准再按“测量1”测量）。 4.确定无误后，按“确定”键记录，自动开始计算建站误差，一般要求建站误差在5mm以内。（考虑仪器自身状态和其他情况，计算出结果有几种不正常情况：①建站误差过大，处理办法为按一次“ESC 键”返回测量后视点2，再次瞄准，测量，再计算，若还是很大，重新建站；②出现“输入第三个点”，处理办法为检查输入点坐标是否输入有误，确定无误，再次测量，若不行，重新建站） 建站误差符合要求后，按“确定”键记录，重新输入点名，其他可按“确定”或“”键跳过，最后“确定”键完成建站。 二、放样 点击“放样”键，按“确定”或“”键跳过界面，至下一个坐标输入界面，输入坐标，瞄准，“测量1”测量，按指示告知架镜人员左右前后移动至定点位置，通知定点。 一次“ESC”键返回，再按“确定”或“”键跳过界面，进入下
#include for(l=0;l 1、角度测量（angle observation） （1）功能：可进行水平角、竖直角的测量。 （2）方法：与经纬仪相同，若要测出水平角∠AOB ，则： 1）当精度要求不高时： 瞄准A 点——置零（0 SET ）——瞄准B 点，记下水平度盘HR 的大小。 2）当精度要求高时：——可用测回法（method of observation set ）。 操作步骤同用经纬仪操作一样，只是配置度盘时，按“置盘”（H SET ）。 2、距离测量（distance measurement ） PSM 、PPM 的设置——测距、测坐标、放样前。 1）棱镜常数（PSM ）的设置。 一般：PRISM=0 （原配棱镜），-30mm （国产棱镜） 2）大气改正数（PPM ）（乘常数）的设置。 输入测量时的气温（TEMP ）、气压（PRESS ），或经计算后，输入PPM 的值。（1）功能：可测量平距HD 、高差VD 和斜距SD （全站仪镜点至棱镜镜点间高差及斜距） （2）方法：照准棱镜点，按“测量”（MEAS ）。 3、坐标测量（coordinate measurement ） （1）功能：可测量目标点的三维坐标（X ，Y ，H ）。 （2）测量原理任意架仪器，先设置仪器高为0，棱镜高是多少就是多少，棱镜拿去直接放在已知点上测高差，测得的高差为棱镜头到仪器视线的高差，当然，有正有负了，然后拿出计算器用已知点加上棱镜高，再加上或减去（因为有正有负）测得的高差就是仪器的视线高啊，因为仪器高为0，所以这个数字就是你的测站点高程，进测站点把它改成这个数字就行了，改完测站点了一般情况下都要打一下已知点复核一下。。。 若输入：方位角，测站坐标（，）；测得：水平角和平距。则有： 方位角： 坐标： 若输入：测站S 高程，测得：仪器高i ，棱镜高v ，平距，竖直角，则有： 高程： （3）方法： 输入测站S （X ，Y ，H ），仪器高i ，棱镜高v ——瞄准后视点B ，将水平度盘读数设置为——瞄准目标棱镜点T ，按“测量”，即可显示点T 的三维坐标。 4、点位放样(Layout) （1）功能：根据设计的待放样点P 的坐标，在实地标出P 点的平面位置及填挖高度。（2）放样原理 1）在大致位置立棱镜，测出当前位置的坐标。 2）将当前坐标与待放样点的坐标相比较，得距离差值dD 和角度差dHR 或纵向差值Δ X 和横向差值Δ Y 。 3）根据显示的dD 、dHR 或ΔX 、ΔY ，逐渐找到放样点的位置。 5、程序测量（programs ） （1）数据采集(data collecting) （2）坐标放样(layout) （3）对边测量(MLM)、悬高测量(REM)、面积测量(AREA)、后方交会(RESECTION) 等。 程序运行环境为Visual Studio2010.运行前请先将坐标数据放在debug 下。 1.单像空间后方交会 C语言程序： #include 18．后方交会测量 后方交会通过对多个已知点的测量定出测站点的坐 标。 输入值或观测值 输出值 Ni.Ei.Zi：已知点的坐标值No.Eo.Zo：测站点的坐标值 Hi ：水平角观测值 Vi ：垂直角观测值 Di ：距离观测值 已知点(P1) 已知点(P2) 测站点(P0) 已知点(P4) 已知点(P3) BTS-800 通过观测2-10 已知点便可计算出测站点的坐标。 当观测的已知点超过 2 个,计算N、E 坐标时将采用最小二乘法进行平差，并给出平差结果的不确定度。而Z 坐标则通过计算平均值求取。因此，观测的已知点越多，计算所得的坐标精度也就越高。 后方交会测量也可在菜单模式下选取“后方交会”来进行。 使用“后方交会”，已知点输入应按顺时针顺序输入，否则计算结果可能不准确。 49 18.1 测量两已知点求取测站坐标 操作过程操作键 1．在测量模式第三页下按【后交】进入 显示 后方交会测量功能，显示已知点坐标 输入屏幕。 在菜单模式下选取“3.后方交会”也 可以进入后方交会测量 2.输入已知点1 的坐标，每输入一行数据 按【】，输入完成后，照准已知点 1 棱镜，按【测量】进行测量。 3.测量完成后，显示测量结果，并要求输 入已知点棱镜高。【后交】 【测量】 【输入测量已知点1】 N﹤m﹥: E﹤m﹥: Z﹤m﹥: 【后方交会】 S: 557.259m ZA: 97°31′05″ HAR: 351°15′06″ 连续测量模式需按【停止】停止测量。 4.按【确定】，进入已知点2 坐标输入及测 量。 重复2-3 完成已知点2 的输入及测量。 【确定】 50 { System; System.Collections.Generic; System.Linq; System.Text; namespace 单像空间后方交会 { class Program { static void Main( string [] args) for (j = 0; j < 5; j++) if (j < 3) "请输入第 {0} 个点的第 {1} 个地面坐标： ", i + 1, j + 1); double .Parse( Console .ReadLine()); "请输入第 {0} 个点的第 {1} 个像点 坐标： ", i + 1, j - 2); double .Parse( Console .ReadLine()); Console .WriteLine(); // 归算像点坐标 (i = 0; i < 4; i++) for (j = 3; j < 5; j++) if (j == 3) zuobiao[i, j] = zuobiao[i, j] - x0; else zuobiao[i, j] = zuobiao[i, j] - y0; // 计算和确定初值 double zs0 = m * f, xs0 = 0, ys0 = 0; for (i = 0; i < 4; i++) else using using using using x0 = y0 = int x0, y0, i, j; double f, m; Console .Write( " 请输入像片比例尺： "); double .Parse( Console .ReadLine()); Console .Write( " 请输入像片的内方位元素 x0:" ); // 均以毫米为单 位 int .Parse( Console .ReadLine()); Console .Write( " 请输入像片的内方位元素 y0:" ); int .Parse( Console .ReadLine()); Console .Write( " 请输入摄影机主距 f:" ); double .Parse( Console .ReadLine()); Console .WriteLine(); // 输入坐标数据 double [,] zuobiao = new double [4, 5]; (i = 0; i < 4; i++) for Console .Write( zuobiao[i, j] = Console .Write( zuobiao[i, j] = for 作业报告 空间后方交会 专业：测绘工程 班级：2008级(1)班姓名：陈闻亚 指导教师：陈强 2010 年 4 月16 日 1 作业任务------------------------------------------------------------------------------------ 3 2 作业思想 --------------------------------------------------------------------------------------- 3 3 作业条件及数据 -------------------------------------------------------------------- 3 4 作业过程--------------------------------------------------------------------------- 3 5 源程序----------------------------------------------------------------------------- 4 6 计算结果--------------------------------------------------------------------------- 17 7心得体会与建议----------------------------------------------------------------------------- 17 1 作业任务 计算近似垂直摄影情况下后方交会解。即利用摄影测量空间后方交会的方法，获取相片的6个外方位元素。限差为0.1。 2作业思想 利用摄影测量空间后方交会的方法求解。该方法的基本思想是利用至少三个一直地面控制点的坐标A（X A，Y A，Z A）、B（X B，Y B，Z B）C（X C，Y C，Z C），与其影像上对应的三个像点的影像坐标a（x a，y a）、b（x b，y b）、c（x c，y c），根据共线方程，反求该相片的外方位元素X S、Y S、Z S、φ、ω、κ。 3作业条件及数据 已知摄影机主距f=153.24mm，四对点的像点坐标与相应的地面坐标列入下表： 4作业过程 4．1 获取已知数据 相片比例尺1/m=1：10000，内方位元素f=153.24mm，x0，y0；获取控制点的地面测量坐标X t、Y t、Z t。 4．2 量测控制点的像点坐标： 本次作业中为已知。见表1。 后方交会在施工测量中的应用与分析(图文) 中得到了广泛的应用，该工程为一等工程，工程规模为大一型、设计洪水位标冷为100年一遇，抗震烈度为6度。该工程是一座以航运、发电和改善哈尔滨市水环境为主，同时具有交通、水产养殖和旅游等综合功能的低水头航电枢纽工程。 问题的提出：在大顶子山航电枢纽工程的施工控制网加密过程中，受到地理条件的限制，首级控制网点之间相互不通视或通视条件不好，为此笔者采用了后方交会的办法解决了施测过程中遇到的困难，在实际生产过程中取得了很好的效果。 一、观测方法与基本原理 结合现场实际情况，在首级控制网的基础上，布设了加密控制网。根据松花江大顶子山航电枢纽厂房、泄洪闸、船闸土建工程所处的施工部位，本着便于整体控制，易于保存的原则，以首级控制网为基础，在施工区周围布设了JK01、JK02、JK03、JK04四个加密点。这些加密点，分布均匀，通视条件好，地基稳定且不易被破坏，对整个施工区域可以进行全方位的观测。加密控制网布设原则以首级控制点为基础，并按二等的施测方安案做了一条闭合导线。 由于首级控制点江南SN01、SN02、SN03、02-1之间互不通视，江北SN04、SN05互不通视。受地形、通视条件的限制，采取边角后方交会的方法，加密了JK01点、JK02点，再由SN02-Jk01起算，复核JK02，在布网过程中，为了保证精度，在不同的测站使用不同仪器和由不同人 员观测，采取了增加多余观测、增加测回数、强制归心等措施，后视SN01 、SN02、02-1，使用徕卡TCR1800全站仪，观测9个测回，经过计算JK01点的误差为2.3mm，达到二等的精度要求。JK01与JK02、JK03 、JK04、SN02构成一条闭合导线。 精度指标严格执法《水利水电工程施工测量规范》（DL/T5173-2003）中二等控制网的技术要求。Mb1.0、Mp（5~7）mm（注：Mb：测量角中误差；Mp：平面控制网点的点位中误差）。 使用仪器及观测方法。使用仪器为瑞士徕佧TCR1800系列全站仪，新建控制点采用具有强制归心装置的混凝土观测墩，水平角观测采用测回法，施测9个测回，同测回盘左、盘右所得角值较差小于4，半测回归零差小于6，同方向各测回互差小6；2C值互差小于9，距离观测采用电磁波测距（往返测），并进行了温度和气压修正。 二、精度计算与分析 1、平面部分精度计算，边角后方交会法测量测站点的精度估算公式为： {[1+(sin2)/(K2-sin2)]m2s+[1+(cos2)/(K2-sin2)] 2(s2m2/2)} =2.3mm(5~7)mm 其中：Mp JK01为测站点JK01的点位中误差，单位为mm =270611.4722 K=363.9389273/363.9341726=1.000013065 ms=0.0012855m S=652.166462 M=0.0392 =206265 一． 实验目的： 掌握摄影测量空间后方交会的原理，利用计算机编程语言实现空间后方交会外方位元素的解算。 二． 仪器用具及已知数据文件： 计算机windows xp 系统，编程软件（VISUAL C++6.0），地面控制点在摄影测量坐标系中的坐标及其像点坐标文件shuju.txt 。 三． 实验内容： 单张影像的空间后方交会：利用已知地面控制点数据及相应像点坐标根据共线方程反求影像的外方位元素。 数学模型：共线条件方程式： )(3)(3)(3)(1)(1)(1Zs Z c Ys Y b Xs X a Zs Z c Ys Y b Xs X a f x -+-+--+-+--= )(3)(3)(3)(2)(2)(2Zs Z c Ys Y b Xs X a Zs Z c Ys Y b Xs X a f y -+-+--+-+--= 求解过程： （1）获取已知数据。从航摄资料中查取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄影测量坐标。 （2）量测控制点的像点坐标并做系统改正。 （3）确定未知数的初始值。在竖直摄影且地面控制点大致分布均匀的情况下，按如下方法确定初始值，即： n X X S ∑=0,n Y Y S ∑=0,n Z mf Z S ∑=0 φ =ω=κ=0 式中；m 为摄影比例尺分母；n 为控制点个数。 （4）用三个角元素的初始值，计算个方向余弦，组成旋转矩阵R 。 （5）逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面 坐标代入共线方程式，逐点计算像点坐标的近似值（x ）、（y ）。 （6）逐点计算误差方程式的系数和常数项，组成误差方程式。 （7）计算法方程的系数矩阵A A T 和常数项l A T ，组成法方程式。 （8）解法方程，求得外方位元素的改正数dXs ,S dY ，s dZ ，d φ,d ω,d κ。 （9）用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素 的新值。 空间后方交会程序 ————————————————————————————————作者：————————————————————————————————日期： ? 一． 实验目的: 掌握摄影测量空间后方交会的原理,利用计算机编程语言实现空间 后方交会外方位元素的解算。 二． 仪器用具及已知数据文件： 计算机wind ｏws xp 系统,编程软件(VI ＳＵA Ｌ C ＋+6．0)，地面控 制点在摄影测量坐标系中的坐标及其像点坐标文件shu ｊu.txt 。 三． 实验内容: 单张影像的空间后方交会:利用已知地面控制点数据及相应像点坐标根据 共线方程反求影像的外方位元素。 数学模型：共线条件方程式: ) (3)(3)(3) (1)(1)(1Zs Z c Ys Y b Xs X a Zs Z c Ys Y b Xs X a f x -+-+--+-+--= ) (3)(3)(3)(2)(2)(2Zs Z c Ys Y b Xs X a Zs Z c Ys Y b Xs X a f y -+-+--+-+--= 求解过程: （1)获取已知数据。从航摄资料中查取平均航高与摄影机主距;获取 控制点的地面测量坐标并转换为地面摄影测量坐标。 （2)量测控制点的像点坐标并做系统改正。 (３)确定未知数的初始值。在竖直摄影且地面控制点大致分布均匀 的情况下,按如下方法确定初始值,即: n X X S ∑=0,n Y Y S ∑=0,n Z mf Z S ∑=0 φ =ω=κ=0 式中;ｍ为摄影比例尺分母;ｎ为控制点个数。 （４)用三个角元素的初始值,计算个方向余弦，组成旋转矩阵R 。 (5)逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面坐标代入共 线方程式，逐点计算像点坐标的近似值（x ）、（y ）。 （6)逐点计算误差方程式的系数和常数项,组成误差方程式。 (7）计算法方程的系数矩阵A A T 和常数项l A T ,组成法方程式。 （8）解法方程，求得外方位元素的改正数dXs ,S dY ，s dZ ,d φ,ｄω，d κ。 （９）用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素的新值。 空间后交-前交程序设计 (实验报告) 姓名： 班级： 学号： 时间： 空间后交-前交程序设计 一、实验目的 用 C 、VB或MATLAB语言编写空间后方交会-空间前方交会程序 ⑴提交实习报告：程序框图、程序源代码、计算结果、体会 ⑵计算结果：像点坐标、地面坐标、单位权中误差、外方位元素及其精度 二、实验数据 f=150.000mm，x0=0，y0=0 三、实验思路 1.利用空间后方交会求左右像片的外方位元素 (1).获取m(于像片中选取两点，于地面摄影测量坐标系中选取同点，分别计算距离，距离比值即为m)，x,y,f,X,Y,Z (2).确定未知数初始值Xs,Ys,Zs,q,w,k (3).计算旋转矩阵R (4).逐点计算像点坐标的近似值（x），（y） (5).组成误差方程式 (6).组成法方程式 (7).解求外方位元素 (8).检查是否收敛，即将求得的外方位元素的改正数与规定限差比较，小于限差即终止；否则用新的近似值重复步骤（3）-（7） 2.利用求出的外方位元素进行空间前交，求出待定点地面坐标(1).用各自像片的角元素计算出左、右像片的方向余弦值，组成旋转矩阵R1,R2 (2).根据左、右像片的外方位元素，计算摄影基线分量Bx，By，Bz (3).计算像点的像空间辅助坐标（X1,Y1,Z1）和（X2,Y2,Z2） (4).计算点投影系数N1和N2 (5).计算未知点的地面摄影测量坐标 四、实验过程 ⑴程序框图 函数AandL %求间接平差时需要的系数 %%%已知 %a=像点坐标x,b=像点坐标y，f内方位元素主距 %φ=q,ψ=w,κ=k %像空间坐标系X,Y,Z %地面摄影测量坐标系Xs,Ys,Zs function [A1,L1,A2,L2]=AandL(a,b,f,q,w,k,X,Y,Z,Xs,Ys,Zs) %%%%%%%%%%%选择矩阵元素 a1=cos(q)*cos(k)-sin(q)*sin(w)*sin(k); a2=-cos(q)*sin(k)-sin(q)*sin(w)*cos(k); a3=-sin(q)*cos(w); b1=cos(w)*sin(k); b2=cos(w)*cos(k); b3=-sin(w); c1=sin(q)*cos(k)+cos(q)*sin(w)*sin(k); c2=-sin(q)*sin(k)+cos(q)*sin(w)*cos(k); c3=cos(q)*cos(w); %%%%%%%共线方程的分子分母 X_=a1*(X-Xs)+b1*(Y-Ys)+c1*(Z-Zs); Y_=a2*(X-Xs)+b2*(Y-Ys)+c2*(Z-Zs); Z_=a3*(X-Xs)+b3*(Y-Ys)+c3*(Z-Zs); %%%%%%%近似值 x=-f*X_/Z_; y=-f*Y_/Z_; %%%%%%%A组成L组成 a11=1/Z_*(a1*f+a3*x); a12=1/Z_*(b1*f+b3*x); a13=1/Z_*(c1*f+c3*x); a21=1/Z_*(a2*f+a3*y); a22=1/Z_*(b2*f+b3*y); a23=1/Z_*(c2*f+c3*y); a14=y*sin(w)-(x/f*(x*cos(k)-y*sin(k))+f*cos(k))*cos(w); a15=-f*sin(k)-x/f*(x*sin(k)+y*cos(k)); a16=y; a24=-x*sin(w)-(y/f*(x*cos(k)-y*sin(k))-f*sin(k))*cos(w); a25=-f*cos(k)-y/f*(x*sin(k)+y*cos(k)); a26=-x; lx=a-x; ly=b-y; %%%%%%%%%组成一个矩阵，并返回 A1=[a11,a12,a13,a14,a15,a16]; A2=[a21,a22,a23,a24,a25,a26]; L1=lx; L2=ly; 函数deg2dms 空间后方交会的直接解 空间后方交会，即由物方已知若干个控制点以及相应的像点坐标，解求摄站的坐标与影像的方位，这是一个摄影测量的基本问题。通常采用最小二乘解算，由于原始的观测值方程是非线性的，因此，一般空间后方交会必须已知方位元素的初值，且解算过程是个迭代解算过程。但是，在实时摄影测量的某些情况下，影像相对于物方坐标系的方位是任意的，且没有任何初值可供参考。这时常规的空间后方交会最小二乘算法就无法处理，而必须建立新的空间后方交会的直接解法。 直接解法的基本思想是将它分成两步：先求出三个已知点i P 到摄站S 的距离i S ；然后求出摄站S 的坐标和影像方位。 物方一已知点()i i i i ,Z ,Y X P 在影像上的成像()i i i ,y x p ，根据影像已知的内方位元素()0 ,y f,x 可求得从摄站()S S S S ,Z ,Y X 到已知点i P 的观测方向i ,βαi 。 () ??? ????-+-= -=2 020 tan tan x x f y y βf x x αi i i i i (1) 距离方程组可以写成如下形式： ?? ??? =+++=+++=+++020202312 1133123232 3322322122 2211221b x x x a x b x x x a x b x x x a x (2) 其中()j ;i ,,i,j S ,b a ij ij ij ij ≠===321cos ?。因此，解算摄站S 到三个 控制点的距离问题，被归结为解算一个三元二次联立方程组的问题。这个方程组的解算方法选用迭代法。 迭代计算公式可写成： 一、观测方法与基本原理 ?? 结合现场实际情况，在首级控制网的基础上，布设了加密控制网。根据松花江大顶子山航电枢纽厂房、泄洪闸、船闸土建工程所处的施工部位，本着便于整体控制，易于保存的原则，以首级控制网为基础，在施工区周围布设了JK01、JK02、JK03、JK04四个加密点。这些加密点，分布均匀，通视条件好，地基稳定且不易被破坏，对整个施工区域可以进行全方位的观测。加密控制网布设原则以首级控制点为基础，并按二等的施测方安案做了一条闭合导线。 ?? 由于首级控制点江南SN01、SN02、SN03、02-1之间互不通视，江北SN04、SN05互不通视。受地形、通视条件的限制，采取边角后方交会的方法，加密了JK01点、JK02点，再由SN02-Jk01起算，复核JK02，在布网过程中，为了保证精度，在不同的测站使用不同仪器和由不同人员观测，采取了增加多余观测、增加测回数、强制归心等措施，后视SN01 、SN02、02-1，使用徕卡TCR1800全站仪，观测9个测回，经过计算JK01点的误差为2.3mm，达到二等的精度要求。JK01与JK02、JK03 、JK04、SN02构成一条闭合导线。 ?? 精度指标严格执法《水利水电工程施工测量规范》（DL/T5173-2003）中二等控制网的技术要求。Mb<1.0”、Mp<（5~7）mm（注：Mb：测量角中误差；Mp：平面控制网点的点位中误差）。 ?? 使用仪器及观测方法。使用仪器为瑞士徕佧TCR1800系列全站仪，新建控制点采用具有强制归心装置的混凝土观测墩，水平角观测采用测回法，施测9个测回，同测回盘左、盘右所得角值较差小于4”，半测回归零差小于6”，同方向各测回互差小6”；2C值互差小于9”，距离观测采用电磁波测距（往返测），并进行了温度和气压修正。二、精度计算与分析 1、平面部分精度计算，边角后方交会法测量测站点的精度估算公式为： ?? {[1+(sin2β)/(K2-sin2β)]m2s+[1+(cos2β)/(K2-sin2β)] 2(s2m2β/ρ2)} ?? =±2.3mm<±(5~7)mm ?? 其中：Mp JK01为测站点JK01的点位中误差，单位为mm ?? β=27”06”11.4722” ?? K=363.9389273/363.9341726=1.000013065全站仪后方交会法步骤和高程测量步骤
摄影测量程序汇总(后方交会+前方交会+单模型光束法平差)
后方交会测量步骤
摄影测量学单像空间后方交会程序设计作业
作业4--空间后方交会
后方交会在施工测量中的应用与分析(图文)
空间后方交会程序
空间后方交会程序
空间后交-前交程序设计实验报告
空间后方交会的直接解
全站仪后方交会测量
相关主题
文本预览