近景攝影測量試驗(yàn)報(bào)告

1、課間實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)課程：近景攝影測量實(shí)驗(yàn)班級：學(xué)生姓名：學(xué) 號：指導(dǎo)教師：重慶交通大學(xué)測量與空間數(shù)據(jù)處理實(shí)驗(yàn)室 TOC o 1-5 h z HYPERLINK l bookmark2 o Current Document 實(shí)驗(yàn)一 近景單張像片空間后方交會(huì) 3 HYPERLINK l bookmark4 o Current Document 實(shí)驗(yàn)二 近景立體像對空間前方交會(huì) 8 HYPERLINK l bookmark6 o Current Document 實(shí)驗(yàn)三直接線性變換求L系數(shù)精確解13 HYPERLINK l bookmark8 o Current Document 實(shí)驗(yàn)四直接線性變換求

2、空間坐標(biāo)精確值 17實(shí)驗(yàn)一近景單張像片空間后方交會(huì)一、實(shí)驗(yàn)?zāi)康耐ㄟ^單張像片空間后方交會(huì)算法編程，掌握空間后方交會(huì)的基本原理和基 本步驟，為近景攝影測量解析處理方法打下理論基礎(chǔ)。二、實(shí)驗(yàn)內(nèi)容利用VB編程平臺(tái)，通過給定的單片像點(diǎn)的像點(diǎn)坐標(biāo)、內(nèi)方位元素和控制點(diǎn) 物方空間坐標(biāo)，計(jì)算出像片的外方位元素。三、實(shí)驗(yàn)步驟：1、編程流程圖2、程序代碼定義內(nèi)方位元素變量Dim x0, yo, f As Double定義外方位元素變量Dim Xs, Ys, Zs As DoubleDim Q, w, k As Double定義旋轉(zhuǎn)變量元素Dim a1, a2, a3 As DoubleDim b1, b2, b3

3、As DoubleDim c1, c2, c3 As Double定義A矩陣Dim a11, a12, a13, a14, a15, a16 As DoubleDim a21, a22, a23, a24, a25, a26 As Double定義像片共線方程變量Dim X8 As Double, Y8 As Double, Z8 As Double定義像點(diǎn)坐標(biāo)變量Dim x, y, z As Double定義像片L變量Dim L1 As Double, L2 As Double定義外方位角元素三角函數(shù)值Dim sk, sw, ck, cw As DoublePrivate Sub Comma

4、nd1_Click()x0 = 48.6401Txtx0 = x0y0 = 11.4765Txty0 = y0f = 4548.5949Txtf = fXs = 3391.658TxtXs = XsYs = -135.645TxtYs = YsZs = 97.25TxtZs = ZsTxtQ = -0 0 35 9.00 ”Txtw = 3 0 33 13.00 ”Txtk = -5 0 4 59.00 ”End SubPrivate Sub Command3_Click()計(jì)算外方位角元素的三角函數(shù)值sw = 0.061982457cw = 0.998077239sk = -0.088599

5、727ck = 0.996067311輸入像點(diǎn)坐標(biāo)x = 599.331y = 1894.198z = -4929.805輸入旋轉(zhuǎn)矩陣a1 = 0.995959: a2 = 0.089231: a3 = 0.010205b1 = -0.088434: b2 = 0.994152: b3 = -0.061983c1 = -0.015676: c2 = 0.06083: c3 = 0.998025計(jì)算像片共線方程系數(shù)X8 = a1 * (x - Xs) + b1 * (y - Ys) + c1 * (z - Zs)Y8 = a2 * (x - Xs) + b2 * (y - Ys) + c2 *

6、(z - Zs)Z8 = a3 * (x - Xs) + b3 * (y - Ys) + c3 * (z - Zs)計(jì)算矩陣Aa11 = (a1 * f + a3 * (x - x0) / Z8a12 = (b1 * f + b3 * (x - x0) / Z8a13 = (c1 * f + c3 * (x - x0) / Z8a14 = (y - y0) * sw - (x - x0) * (x-x0)* ck- (y- y0) * sk) / f + f * ck) * cwa15 = -f * sk - (x - x0) * (x - x0) *sk +(y -y0) *ck)/ fa1

7、6 = y - y0a21 = (a2 * f + a3 * (y - y0) / Z8a22 = (b2 * f + b3 * (y - y0) / Z8a23 = (c2 * f + c3 * (y - y0) / Z8a24 = -(x - x0) * sw - (y - y0) * (x - x0) * ck - (y - y0) * sk) / f - f * sk) *cwa25 = -f * ck - (y - y0) * (x - x0) * sk + (y - y0) * ck) / fa26 = -(x - x0)計(jì)算LL1 = x - (x0 - f * X8 / Z8)

8、L2 = y - (y0 - f * Y8 / Z8)Print all, a12, a13, a14, a15, a16Print a21, a22, a23, a24, a25, a26Print L1, L2Xs = 3391.658 - 0.768Ys = -135.645 - 0.016Zs = 97.25 - 0.819Q = -00 35 8.00 ”w = 3 33 9.00 ”Txtk = -5 0 4 59.00 ”MsgBox 改正后的 Xs 為：& Xs & Chr(13) & 改正后的 Ys 為：& Ys & Chr(13) & 改正后的 Zs 為：& Zs & Ch

9、r(13) & 改正后的小為：& Q & Chr(13) & 改正后的為：& w & Chr(13) & 改正后的k為：& Txtk, vbOKOnly,計(jì) 算結(jié)果End Sub四、實(shí)驗(yàn)結(jié)果 1、實(shí)驗(yàn)數(shù)據(jù)物方坐標(biāo):第一列是點(diǎn)號；單位：毫米第二列是X坐標(biāo)，朝南為正向；第三列是Y坐標(biāo)，朝上為正向；第四列為Z坐標(biāo)，朝前（西）為正向。像點(diǎn)坐標(biāo)：以毫米為單位點(diǎn)名，x, y4301,1105.312313,3367.97924243011617.7121883.870-4917.3084302,1111.483438,3203.99800943021615.5721687.514-4913.7504303

10、,1128.936009,2886.14637443031614.6401315.066-4916.7244304,1145.729428,2615.38530143041615.7841001.358-4916.6914305,1163.579521,2343.60358443051617.197690.040-4917.0704306,1179.483899,2076.67273843061615.933387.492-4917.2344307,1196.664025,1861.97188543071618.946145.491-4917.8074308,1216.048311,1565.

11、04608743081616.753-186.640-4918.5314309,1236.196347,1299.27541343091617.714-482.264-4918.5614310,1249.998721,1120.5242543101617.859-680.767-4918.799x0 = 48.6401mm , y0 = 11.4765mm , f = 4548.5949mm2、實(shí)驗(yàn)結(jié)果內(nèi)方位元素外方位元素xO:48, 6401Xs科匯65E-忘由1-0g 35 9. QOyO：11, 4765Ys I.13, 33H 13,00*f：|4548. 5949Zs97. 25k

12、I-5D 4* 59.A0*-8普睥33期T554S05KEHEMB0N6IWIF： SKI典T“U癢E-P?翦！|口也罐室印T一笫田髀要用筠l一中 |,網(wǎng)”制苒MBS 3fiZ32S35aEEGQT.甫哭碑景XQ；|48.6401XS|33&1. 658y0；1L47B5 Ys|玉-6如f.|4548.5949Zs|97725IBO 4時(shí)7淅掇“1a1 1S-529 哨T 軸 99外方位元素計(jì)艇HX中RIE目隙由（E）; 6.56ZSlEJ6SS5?tJm） : -2.1 T8133dLs的Z助附心力.1箱近雛學(xué)為：-1Oa.OA851的E百色妙為雨盼h班正聞舐為；-24279,n1 臉輸人

13、亮畢耨H五、心得體會(huì)在做此次試驗(yàn)之前，先對上學(xué)期所學(xué)的該知識點(diǎn)進(jìn)行了系統(tǒng)的復(fù)習(xí)，從整體上把握了單相空間后方交會(huì)原理及解算過程。 也理解了將控制點(diǎn)的地面點(diǎn)坐標(biāo)視為真值，把相應(yīng)的像點(diǎn)坐標(biāo)視為觀測值， 誤差方程，然后用矩陣表示，最后構(gòu)成慶 數(shù)。加入相應(yīng)的改正數(shù)，最后列出每個(gè)點(diǎn)的 誤差方程，從而求出外方位:元素的改正實(shí)驗(yàn)二近景立體像對空間前方交會(huì)、實(shí)驗(yàn)?zāi)康耐ㄟ^立體像對空間前方交會(huì)算法編程，掌握空間前方交會(huì)的基本原理和基 本步驟，為近景攝影測量解析處理方法打下理論基礎(chǔ)。二、實(shí)驗(yàn)內(nèi)容利用VB編程平臺(tái)，通過給定的立體像對像點(diǎn)的像點(diǎn)坐標(biāo)以及對應(yīng)的內(nèi)外方 位元素，實(shí)現(xiàn)待測控制點(diǎn)的物方空間坐標(biāo)的求解。三、實(shí)驗(yàn)步驟

14、：1、編程流程圖2、程序代碼Dim f#, x1#, y1#, x2#, y2#, i%, j%, u1#, u2#, v1#, v2#, w1#, w2#, fail#, kabl#, omg1#, fai2#, kab2#, omg2#Dim a12#, a13#, b11#, b12#, b13#, c11#, c12#, c13#, a21#, a22#, a23#, b21#, b22#, b23#, c21#, c22#, c23#Dim n1#, n2#, bu#, bv#, bw#Dim xs1#, xs2#, ys1#, ys2#, zs1#, zs2#Private Sub

15、 Command1_Click()fai1 = Val(Text1.Text): kab1 = Val(Text3.Text): omg1 = Val(Text2.Text)fai2 = Val(Text11.Text): kab2 = Val(Text10.Text): omg2= Val(Text7.Text)f = Val(Text13.Text)xs1 = Val(Text4.Text): xs2 = Val(Text9.Text): ys1 = Val(Text5.Text): ys2=Val(Text8.Text): zs1 = Val(Text6.Text): zs2 = Val

16、(Text12.Text)End SubPrivate Sub Command3_Click()x1 = Val(InputBox( y1 = Val(InputBox( x2 = Val(InputBox( y2 = Val(InputBox(輸入第一片坐標(biāo)x的值) 輸入第一片坐標(biāo)y的值) 輸入第二片坐標(biāo)x的值) 輸入第二片坐標(biāo)y的值)End SubPrivate Sub Command2_Click()Text1.Visible = FalseText2.Visible = FalseText3.Visible = FalseText4.Visible = FalseText5.Visib

17、le = FalseText6.Visible = FalseText7.Visible = FalseText8.Visible = FalseText9.Visible = FalseText10.Visible = FalseText11.Visible = FalseText12.Visible = FalseText13.Visible = FalseLabel1.Visible = FalseLabel2.Visible = FalseLabel3.Visible = FalseLabel4.Visible = FalseLabel5.Visible = FalseLabel6.V

18、isible = FalseLabel7.Visible = FalseLabel8.Visible = FalseLabel9.Visible = FalseLabel10.Visible = FalseLabelll.Visible = FalseLabel12.Visible = FalseLabel13.Visible = FalseLabel14.Visible = Falseall = Cos(fail) * Cos(kabl) - Sin(fail) * Sin(omgl) * Sin(kabl)a12 = -Cos(fai1) * Sin(kabl) - Sin(fail) *

19、 Sin(omgl) * Cos(kabl)a13 = -Sin(fai1) * Cos(omg1)b11 = Cos(omg1) * Sin(kab1)b12 = Cos(omg1) * Cos(kab1)b13 = -Sin(omg1)c11 = Sin(fai1) * Cos(kab1) + Cos(fai1) * Sin(omg1) * Sin(kab1)c12 = -Sin(fai1) * Sin(kab1) + Cos(fai1) * Sin(omg1) * Cos(kab1)c13 = Cos(fai1) * Cos(omg1)a21 = Cos(fai2) * Cos(kab2

20、) - Sin(fai2) * Sin(omg2) * Sin(kab2)a22 = -Cos(fai2) * Sin(kab2) - Sin(fai2) * Sin(omg2) * Cos(kab2)a23 = -Sin(fai2) * Cos(omg2)b21 = Cos(omg2) * Sin(kab2)b22 = Cos(omg2) * Cos(kab2)b23 = -Sin(omg2)c21 = Sin(fai2) * Cos(kab2) + Cos(fai2) * Sin(omg2) * Sin(kab2)c22 = -Sin(fai2) * Sin(kab2) + Cos(fai

21、2) * Sin(omg2) * Cos(kab2) c23 = Cos(fai2) * Cos(omg2)u1 = a11 * x1 + a12 * y1 - a13 * fv1 = b11 * x1 + b12 * y1 - b13 * fw1 = c11 * x1 + c12 * y1 - c13 * fu2 = a21 * x2 + a22 * y2 - a23 * fv2 = b21 * x2 + b22 * y2 - b23 * fw2 = c21 * x2 + c22 * y2 - c23 * fbu = xs2 - xs1bv = ys2 - ys1bw = zs2 - zs1

22、n1 = (bu * w2 - bw * u2) / (u1 * w2 - u2 * w1)n2 = (bu * w1 - bw * u1) / (u1 * w2 - u2 * w1)Dim xx#, yy#, yy1#, yy2#, zz#xx = Fix(xs1 + u1 * n1) * 1000 + 0.5) / 1000yy1 = Val(Text5.Text) + v1 * n1yy2 = Val(Text8.Text) + v2 * n2yy = Fix(yy1 + yy2) / 2) * 1000 + 0.5) / 1000zz = Fix(Val(Text6.Text) + w

23、1 * n1) * 1000 + 0.5) / 1000Forml.Print 地面坐標(biāo)為:Forml.Print xx=; xx, yy=; yy; zz=; zzEnd Sub四、實(shí)驗(yàn)結(jié)果1、實(shí)驗(yàn)數(shù)據(jù)相片1坐標(biāo)：點(diǎn)名，x, y4201,896.942591,3200.1324684202,903.373104,3012.7248834203,918.450557,2802.550964204,937.004622,2449.9082384205,957.037515,2086.3138984206,970.889881,1858.4535284207,979.044342,1714.819

24、8454208,989.950247,1549.5977224209,1001.556951,1384.7020554210,1015.53176,1189.988689相片2坐標(biāo)：點(diǎn)名，x, y4201,595.600431,3258.0235544202,601.34579,3061.8024614203,615.818572,2841.9570754204,633.003988,2472.9369064205,651.907837,2092.444594206,664.95803,1854.1427144207,672.855606,1703.7227394208,683.251761,

25、1530.7373154209,694.71749,1357.8616454210,708.184825,1154.038955第1張像片：Xs = 4310.470mm Ys = -145.342mm Zs = -295.236mm R =0.949693-0.086534-0.3009913 =0.1011730.9943090.0333630.296391-0.0621370.953043k = -4 58 25.00第2張像片：Xs = 4079.418mm Ys = -145.179mm Zs = -298.330mm R =0.973090-0.085100-0.2141360.0

26、961040.9945060.0414920.209429-0.0609550.975922k = -4 0 53 27.006=-12 0 6 42.00=3 29 40.00 x0 = 48.6401mm y0 = 11.4765mm , f = 4548.5949mm3、實(shí)驗(yàn)結(jié)果曰 Formll輸入主距和外方位元素電認(rèn)拓久匚通二3 forml地面坐標(biāo)為Qm）:1= -4030. 729 丫= 3261.342 17除 073地面坐標(biāo)為腦101-4002.09- 3373.6531.57.564弛面坐標(biāo)為口卜1= 3974 9643456.96227645.521地面坐標(biāo)為LQ.I- 39

27、09. 334 3SP9.S64Z-MS635.803地面坐:標(biāo)為）；1= 3327. ZY= 366Z 443Z=T024. 711地面坐標(biāo)為31= 3772. 301 V= 3612.JGSZ=-0G93. ES9地面坐標(biāo)為缶百：I- 3736. 329- 353L.036Z-9057. 26地面坐標(biāo)為（卜=3fi39.1L3= 3360 936Z=536b 颯地面坐標(biāo)為（“LI- 3669. 299- 3153.405MT 蝌.434地面坐:標(biāo)為（）： E-3M2,eO3 1=2602.03$ICU4生2五、心得體會(huì)通過這次編程實(shí)驗(yàn)，對于共線方程和立體像對前方交會(huì)公式的理解更為透徹 了，

28、對于方程中Bu矩陣和常數(shù)L矩陣的求解也更為熟悉。各個(gè)等式中每個(gè)字母 所代表的含義也更為清楚。前方交會(huì)的編程其實(shí)不算難，只是在計(jì)算中需要定義 的變量很多，很容易混淆，因此在定義變量時(shí)要加上注釋并定義得與變量所代表 的含義相近。由于在VB中計(jì)算矩陣比較困難，因此可以將矩陣中的每個(gè)元素定 義為一個(gè)相應(yīng)的變量，在VB計(jì)算出后再到Matlab中計(jì)算矩陣求出改正數(shù)，將所 求出的改正數(shù)在VB中輸入并與之前的近似坐標(biāo)計(jì)算再輸出最后物方坐標(biāo)。這樣 求解比在VB中計(jì)算矩陣要省時(shí)省力。所以在實(shí)驗(yàn)中，不妨多借助幾樣軟件來達(dá) 到高效運(yùn)算的效果。實(shí)驗(yàn)三直接線性變換求L系數(shù)精確解一、實(shí)驗(yàn)?zāi)康耐ㄟ^直接線性變換基本原理編程解算

29、 L系數(shù)精確值，掌握直接線性變換的 基本原理和基本步驟，為近景攝影測量解析處理方法打下理論基礎(chǔ)。二、實(shí)驗(yàn)內(nèi)容利用MATLA踹程平臺(tái)，通過給定的單片像點(diǎn)的像點(diǎn)坐標(biāo)和控制點(diǎn)物方空間坐標(biāo)，計(jì)算出像片的 三、實(shí)驗(yàn)步驟：L系數(shù)精確值。2、編程代碼矩陣AA =1.0e+007 *0 0.1788 0.2082 -0.54350.5448 0.6345 -1.65610 0.1796 0.1876 -0.54620.5176 0.5407 -1.57440 0.1823 0.1485 -0.55510.4660 0.3795 -1.41900 0.1851 0.1147 -0.56330.4226 0.26

30、19 -1.28590 0.1882 0.0803 -0.57210.3790 0.1617 -1.15240 0.1906 0.0457 -0.5800 TOC o 1-5 h z 0.0002 0.0002 -0.0005 0.00000000000 0.0002 0.0002 -0.0005 0.00000.00020.0002-0.00050.000000000000.00020.0002-0.00050.00000.00020.0001-0.00050.000000000000.00020.0001-0.00050.00000.00020.0001-0.00050.000000000

31、000.00020.0001-0.00050.00000.00020.0001-0.00050.00000000000 0.0002 0.0001 -0.0005 0.00000.0002 0.0000 -0.0005 0.0000000矩陣BB =1.0e+003 *1.10533.36801.11153.20401.12892.88611.14572.61541.16362.34361.1795L的近似值L =1.0e+003 *-0.00010.0000-0.0002-1.1128-0.00010.0000-0.0007-3.22460.0000-0.00000.0000矩陣MM =1.

32、0e+011 *-0.0000 -0.0000 0.0001 -0.00000000.0612 0.13080000 -0.0000 -0.0000 0.00010.1865 0.41560.0000 -0.0001 0.0004 -0.0000000.0463 0.12030000 -0.0001 -0.0001 0.00041.3334 0.36150.0000 -0.0000 0.0001 -0.0000000.0737 0.10190000 -0.0000 -0.0000 0.00010.2651 0.27110.0000 -0.0000 0.0001 -0.0000000.0917

33、0.08760000 -0.0000 -0.0000 0.00010.2094 0.20790.0000 -0.0000 0.0001 -0.0000000.1033 0.07450000 -0.0000 -0.0000 0.00010.2081 0.15580.0000 -0.0000 0.0001 -0.0000000.1019 0.06270000 -0.0000 -0.0000 0.00010.1794 0.1145矩陣Ww =1.0e+007 *0.12450.37930.94142.71370.14990.38330.18650.42580.21010.42320.2073-0.0

34、614 -0.0000-0.0152 -0.0159-0.4384 0.4579-0.0242 -0.0197-0.0619 -0.0504-0.0301 -0.0187-0.0688 -0.0426-0.0340 -0.0145-0.0684 -0.0292-0.0335 -0.0080-0.0590 -0.0141-0.0201 -0.02350-0.000000-0.000000-0.000000-0.000000-0.000000-0.0000五、實(shí)驗(yàn)結(jié)果分析1、實(shí)驗(yàn)數(shù)據(jù)像片編號 2-40-g11-k1-0012201,1358.55679,2940.5965532202,1375.1

35、04882,2714.0732622203,1396.439124,2359.6673272206,1454.452586,1516.0217162207,1467.812313,1322.0091122209,1493.982413,993.473865控制點(diǎn)空間坐標(biāo)2201 597.747 1785.157 -5930.7482202 600.377 1446.575 -5934.1462203 599.126 922.699 -5934.2242206 599.171 -299.231 -5933.7772207 597.520 -577.132 -5933.3442209 598.15

36、2 -1046.647 -5934.4452、實(shí)驗(yàn)結(jié)果LG =1.0e+003 *0.00070.0000-0.0008-5.16810.00550.00000.0004-6.6576-0.00000.0000-0.00000.0000五、實(shí)驗(yàn)心得通過本次實(shí)驗(yàn)我的體會(huì)很多，就比如在求解L的近似值，過程比較繁瑣，且 求解物方空間坐標(biāo)的精確值需要和法方程聯(lián)立才能得到結(jié)果。直接線性變換雖然 是直接建立了像點(diǎn)坐標(biāo)和物方空間坐標(biāo)的關(guān)系，但其編程求解過程很是繁瑣。矩 陣的運(yùn)算量非常的大。相比之下倒是共線方程更為快捷。實(shí)驗(yàn)四 直接線性變換求空間坐標(biāo)精確值、實(shí)驗(yàn)?zāi)康耐ㄟ^直接線性變換算法編程，掌握直接線性變換的

37、基本原理和基本步驟,為近景攝影測量解析處理方法打下理論基礎(chǔ)。二、實(shí)驗(yàn)內(nèi)容利用VB編程平臺(tái)，通過給定的立體像對的同名像點(diǎn)坐標(biāo)和L系數(shù)，解算待測點(diǎn)的物方空間坐標(biāo)。三、實(shí)驗(yàn)步驟：1、編程流程圖利用L矩陣求內(nèi)方位元素:結(jié)束利用L矩陣求外方位元素:利用立體像對像點(diǎn)坐標(biāo)和內(nèi)外方位元素求物方坐標(biāo):2、編程代碼(1)用VB編程計(jì)算內(nèi)方位元素和外方位元素角元素:,定義L的精確值變量Dim L1, L2, L3, L4, L5, L6, L7, L8, L9, L10, L11 As Double定義r3Dim r3 As Double,定義計(jì)算變量Dim A, B, C As Double,定義內(nèi)方位元素Dim

38、 x0, y0 As DoubleDim fx, fy As DoubleDim dS As DoubleDim dB As String,定義外方位角元素Dim Q, W, K As Double,定義計(jì)算變量Dim a3, c3, b3 As DoubleDim b1, b2 As DoublePrivate Sub Command1_Click()輸入L的精確值L1 = 0.691694694954L2 = 0.000608383183L3 = -0.825372309722L4 = -5168.07051610071L5 = 5.470559066537L6 = 0.008254963

39、815L7 = 0.441376855385L8 = -6657.62884269957L9 = -0.000695566593L10 = 0.000000547457L11 = -0.000023826626計(jì)算r3r3 = Sqr(L9 A 2 + L10 A 2 + L11 A 2),計(jì)算x0,y0 x0 = -(L1 * L9 + L2 * L10 + L3 * L11) / (L9 a 2 + L10 a 2 + L11 a 2)y0 = -(L5 * L9 + L6 * L10 + L7 * L11) / (L9 A 2 + L10 A 2 + L11 A 2)計(jì)算A,B,CA =

40、 r3 a 2 * (L1 a 2 + L2 a 2 + L3 a 2) - x0 a 2B = r3 a 2 * (L5 a 2 + L6 a 2 + L7 a 2) - y0 a 2C = r3 a 2 * (L1 * L5 + L2 * L6 + L3 * L7) - x0 * y0計(jì)算dB,dsdB = 90 dS = Sqr(A / B) - 1計(jì)算fx,fyfx = Cos(dB) * Sqr(A)fy =僅 / (1 + dS)計(jì)算 a3,c3,b3a3 = r3 * L9b3 = r3 * L10c2 = r3 * L11計(jì)算b1,b2b1 = (r3 * L2 + b3 *

41、x0) / fx + r3*(l6+l10*y0)*(1+ds)*cos(db) / fxb2 = (r3 * (L6 + L10 * y0) * (1 + dS) * Cos(dB) / fx計(jì)算外方位角元素Q = Atn(-a3 / c3)K = Atn(b1 / b2)Sin(W) = -b3MsgBox x0= & x0 & Chr(13) & y0= & y0 & Chr(13) & fx=0 & Chr(13) & fy=0, vbOKOnly,內(nèi)方位元素計(jì)算結(jié)果MsgBox 少=& Q & Chr(13) & = & W & Chr(13) & k = & K , vbOKOnl

42、y,外方位角元素計(jì)算結(jié)果End Sub(2)用MATLA球解外方位線元素:LI =1.0e+003 *0.0006916946&4954-5.163070518100706-0.0000006955665930,0000005083831830,0000082549638150l0000000005474570. 000825372309722。. 0004413768553850(3)用VB編程計(jì)算物方空間坐標(biāo):,定義兩個(gè)像點(diǎn)坐標(biāo)Dim x1, y1 As DoubleDim x2, y2 As Double,定義物方空間坐標(biāo)Dim X, Y, Z As Double,定義兩張像片的L系數(shù)D

43、im L1, L2, L3, L4, L5, L6, L7, L8, L9, L10, L11 As DoubleDim LL1, LL2, LL3, LL4, LL5, LL6, LL7, LL8, LL9, LL10, LL11 As Double,定義基本關(guān)系式分母Dim A1, A2, A3 As Double,定義內(nèi)方位元素Dim x0, y0 As Double,定義像點(diǎn)向徑Dim r1, r2 As Double,定義待定對稱徑向畸變系數(shù)Dim k1, k2 As Double,定義兩張像片的矩陣 WDim W1, W2 As DoubleDim W3, W4 As Double

44、定義矩陣V變量Dim V1, V2, V3, V4, V5, V6 As Double定義矩陣NDim N1, N2, N3, N4, N5, N6, N7, N8, N9, N10, N11, N12, N13, N14, N15, N16, N17, N18 As Double定義矩陣Q的變量Dim Q1, Q2, Q3, Q4, Q5, Q6 As DoublePrivate Sub Command1_Click()Dim i As IntegerDim j As IntegerDim h As IntegerDim k As IntegerDim u As IntegerTextl.Text = Val(Val(InputBox(Text2.Text = Val(Val(InputBox(Text3.Text = Val(Val(InputBox(Text4.Text = Val(Val(InputBox( 輸入第一張像片的x坐標(biāo)) 輸入第一張像片的y坐標(biāo)) 輸入第一張像片的x坐標(biāo)) 輸入第一張像片的y坐標(biāo))輸入L的精確值For i = 1 To 11L(i) = Val(Val(InputBox(輸入第一張像片的Li的精確值)NextFor j = 1 To 11LL(j) = Val