Warszawa dn. 02.06.2008
ZAKŁAD FOTOGRAMETRII
INSTYTUT FOTOGRAMETRII I KARTOGRAFII
Temat ćwiczenia: Analityczne opracowanie pojedynczego stereogramu na stereokomparatorze.
Wykonali:
Mateusz Hut
Piotr Kurek
II rok, IV sem GiK
Cel ćwiczenia:
Celem ćwiczenia było zapoznanie się z obsługą stereokomparatora Carl'a Zeiss'a 1818 nr 222658 oraz nauczenie opracowywania zdjęć lotniczych w postaci stereogramu. Ogólna zasada pomiaru polegała na początkowym zestrojeniu instrumentu wzdłuż bazy, pomiarze punktów głównych zdjęcia, aby na końcu można było dobrze wykonać pomiar współrzędnych tłowych terenu gminy Grybów.
Stereogram składał się z diapozytywu lewego 11-8653 i prawego- 10-8652. Po wykonaniu pomiarów i manualnym wyliczeniu w dzienniku pomiarowym według podanego schematu przystąpiliśmy do policzenia współrzędnych na podstawie programu dostępnego na komputerach w pracowni.
Przebieg ćwiczenia:
Strojenie wzdłuż bazy- wyznaczenie współrzędnych instrumentalnych
Przeliczenie współrzędnych instrumentalnych punktów bazy do układu współrzędnych tłowych zdjęcia. Dane zawarte w dzienniku pomiarowym
Wyznaczenie współrzędnych punktów tłowych specjalnie nakłutych na diapozytywie.
Policzenie współrzędnych punktów nakłutych w układzie współrzędnych zdjęcia oraz zapisanie ich w pliku tekstowym który wykorzysta program.
100.26 46.46
1 0 0 -41.857 -0.015
2 41.28 0.02 0 0
1031 39.48 51.46 -3.564 51.528
1032 46.16 22.44 3.894 22.425
1033 49.14 -9.81 8.051 -9.838
1034 42.18 -28.20 1.153 -28.186
1035 51.98 -52.55 11.464 -52.390
1030 4.50 -58.08 -35.790 -57.956
1029 3.38 -34.30 -37.772 -34.227
1028 3.12 2.55 -39.015 2.238
1027 3.58 28.80 -40.007 28.854
1026 11.32 56.18 -32.143 56.298
Zestawienie faktycznych współrzędnych terenowych punktów zaobserwowanych na zdjęciach
1026 246265.25 630231.66 784.33
1027 247844.08 630753.81 834.03
1028 249358.70 630843.02 651.19
1029 251611.11 630924.08 542.92
1030 253144.08 630917.09 435.85
1031 246589.03 628591.15 730.77
1032 248293.24 628227.55 663.31
1033 250244.63 628050.32 493.59
1034 251359.01 628525.09 493.73
1035 252927.51 627947.34 431.13
Orientacja wzajemna zdjęć
ELEMENTY ORIENTACJI WZAJEMNEJ ZDJEC
na podstawie pliku - C:\Documents and Settings\dephos\Pulpit\studenci_2008\pm\mp.tlo
----------------------------------------------
Bx = 46.46000000
by = 0.01704384
bz = -0.00021286
Omega [g] = 0.41634850
Fi [g] = 0.09211064
Kappa [g] = -0.04235916
Raport: WYZNACZENIE ELEMENTÓW ORIENTACJI WZAJEMNEJ ZDJĘĆ
_________________________________________________________
Plik: C:\Documents and Settings\dephos\Pulpit\studenci_2008\pm\mp.tlo
Ck= 100.260
Bx wziete z pliku
Obserwacje:
1 0.000 0.000 -41.857 -0.015
2 41.280 0.020 0.000 0.000
1031 39.480 51.460 -3.564 51.528
1032 46.160 22.440 3.894 22.425
1033 49.140 -9.810 8.051 -9.838
1034 42.180 -28.200 1.153 -28.186
1035 51.980 -52.550 11.464 -52.390
1030 4.500 -58.080 -35.790 -57.956
1029 3.380 -34.300 -37.772 -34.227
1028 3.120 2.550 -39.015 2.238
1027 3.580 28.800 -40.007 28.854
1026 11.320 56.180 -32.143 56.298
Liczba obserwacji: 12
_________________________________________________________
________________________________________ i= 1
Paralaksy wejsciowe:
1 py= 0.017
2 py= 0.023
1031 py= -0.073
1032 py= 0.016
1033 py= 0.032
1034 py= -0.016
1035 py= -0.183
1030 py= -0.143
1029 py= -0.082
1028 py= 0.344
1027 py= -0.058
1026 py= -0.126
Obs nr. 1028 |py| max = 0.344
[vv]= 746.705532 m0= 10.328224
Tablica wariancyjno-kowariancyjna zwiekszona 10^6 razy:
dBy dBz dOmega dFi dKappa
0.3916934 0.0204124 0.1415710 0.0167849 0.0104259
0.0204124 0.0737650 0.0093513 -0.0781545 -0.0089984
0.1415710 0.0093513 0.0522985 0.0028363 0.0006397
0.0167849 -0.0781545 0.0028363 0.1790027 0.0159660
0.0104259 -0.0089984 0.0006397 0.0159660 0.0204870
Bx z obs. = 41.809 +/- 1.439
| Bx = 46.460 +/- 0.000
| dBy = 0.789 +/- 0.300
| dBz = -0.014 +/- 0.130
|
| dOmega = 0.41618882 +/- 0.01884562 g
| dFi = 0.09958844 +/- 0.03486547 g
| dKappa = -0.04203674 +/- 0.01179519 g
______________________________________________________
________________________________________ i= 2
Paralaksy:
1 py= -0.072
2 py= -0.025
1031 py= 0.026
1032 py= -0.015
1033 py= 0.004
1034 py= 0.006
1035 py= 0.018
1030 py= 0.024
1029 py= -0.088
1028 py= 0.256
1027 py= -0.100
1026 py= 0.009
Obs nr. 1028 |py| max = 0.256
[vv]= 746.684681 m0= 10.328080
Tablica wariancyjno-kowariancyjna zwiekszona 10^6 razy:
dBy dBz dOmega dFi dKappa
0.3960725 0.0218678 0.1421997 0.0201171 0.0107874
0.0218678 0.0729950 0.0110374 -0.0766060 -0.0081071
0.1421997 0.0110374 0.0523365 0.0011253 0.0004988
0.0201171 -0.0766060 0.0011253 0.1785697 0.0156586
0.0107874 -0.0081071 0.0004988 0.0156586 0.0204403
Bx z obs. = 41.653 +/- 1.483
| Bx = 46.460 +/- 0.000
| dBy = 0.003 +/- 0.302
| dBz = 0.004 +/- 0.130
|
| dOmega = 0.00015944 +/- 0.01885219 g
| dFi = -0.00747784 +/- 0.03482279 g
| dKappa = -0.00032245 +/- 0.01178157 g
______________________________________________________
________________________________________ i= 3
Paralaksy:
1 py= -0.075
2 py= -0.029
1031 py= 0.025
1032 py= -0.018
1033 py= 0.000
1034 py= 0.002
1035 py= 0.012
1030 py= 0.021
1029 py= -0.091
1028 py= 0.252
1027 py= -0.103
1026 py= 0.006
Obs nr. 1028 |py| max = 0.252
[vv]= 746.684681 m0= 10.328080
Tablica wariancyjno-kowariancyjna zwiekszona 10^6 razy:
dBy dBz dOmega dFi dKappa
0.3958186 0.0218356 0.1421518 0.0201350 0.0107870
0.0218356 0.0729566 0.0110319 -0.0765896 -0.0081024
0.1421518 0.0110319 0.0523356 0.0011232 0.0004989
0.0201350 -0.0765896 0.0011232 0.1785679 0.0156593
0.0107870 -0.0081024 0.0004989 0.0156593 0.0204405
Bx z obs. = 41.665 +/- 1.479
| Bx = 46.460 +/- 0.000
| dBy = -0.000 +/- 0.302
| dBz = -0.000 +/- 0.130
|
| dOmega = 0.00000025 +/- 0.01885203 g
| dFi = 0.00000004 +/- 0.03482261 g
| dKappa = 0.00000003 +/- 0.01178162 g
______________________________________________________
________________________________________ i= 4
Paralaksy:
1 py= -0.075
2 py= -0.029
1031 py= 0.025
1032 py= -0.018
1033 py= 0.000
1034 py= 0.002
1035 py= 0.012
1030 py= 0.021
1029 py= -0.091
1028 py= 0.252
1027 py= -0.103
1026 py= 0.006
Obs nr. 1028 |py| max = 0.252
[vv]= 746.684681 m0= 10.328080
Tablica wariancyjno-kowariancyjna zwiekszona 10^6 razy:
dBy dBz dOmega dFi dKappa
0.3958186 0.0218356 0.1421518 0.0201350 0.0107870
0.0218356 0.0729566 0.0110319 -0.0765896 -0.0081024
0.1421518 0.0110319 0.0523356 0.0011232 0.0004989
0.0201350 -0.0765896 0.0011232 0.1785679 0.0156593
0.0107870 -0.0081024 0.0004989 0.0156593 0.0204405
Bx z obs. = 41.665 +/- 1.479
| Bx = 46.460 +/- 0.000
| dBy = -0.000 +/- 0.302
| dBz = 0.000 +/- 0.130
|
| dOmega = 0.00000000 +/- 0.01885203 g
| dFi = -0.00000000 +/- 0.03482261 g
| dKappa = -0.00000000 +/- 0.01178162 g
______________________________________________________
Liczba - obserwacji: 12
- iteracji : 4
WYNIKI KOŃCOWE:
Bx z obs. = 41.665 +/- 1.479
| Bx = 46.460 +/- 0.000
| By = 0.792 +/- 0.302
| Bz = -0.010 +/- 0.130
|
| Omega = 0.41634850 +/- 0.01885203 g
| Fi = 0.09211064 +/- 0.03482261 g
| Kappa = -0.04235916 +/- 0.01178162 g
|
________________________________________________________
________________________________________________________
_______________________________________________________
1 py= -0.075
2 py= -0.029
1031 py= 0.025
1032 py= -0.018
1033 py= 0.000
1034 py= 0.002
1035 py= 0.012
1030 py= 0.021
1029 py= -0.091
1028 py= 0.252
1027 py= -0.103
1026 py= 0.006
_______________________________________________________
Obs nr. 1028 |py| max = 0.252
[vv]= 746.684681 m0= 10.328080 [mm]
Powyższe wyniki paralaks py są zadowalające na tyle, że nie trzeba poprawiać można kontynuować pracę na naszych danych.
Obliczenie współrzędnych modelowych
1 0.000 0.038 -111.728 py= -0.075
2 46.624 0.037 -113.239 py= -0.029
1031 42.803 55.778 -108.698 py= 0.025
1032 50.927 24.766 -110.614 py= -0.018
1033 55.761 -11.132 -113.769 py= 0.000
1034 47.916 -32.036 -113.895 py= 0.002
1035 59.832 -60.494 -115.405 py= 0.012
1030 5.187 -66.961 -115.573 py= 0.021
1029 3.821 -38.730 -113.343 py= -0.091
1028 3.454 2.697 -110.999 py= 0.252
1027 3.839 30.932 -107.504 py= -0.103
1026 12.188 60.483 -107.945 py= 0.006
Orientacja bezwzględna i obliczenie współrzędnych terenowych:
>>> ELEMENTY ORIENTACJI BEZWZGLEDNEJ
- Wyniki na podstawie: C:\Documents and Settings\dephos\Pulpit\studenci_2008\pm\mp
Wsp. terenowe (1).X,Y,Z (2).Y,X,Z : 1
Wsp. modelowe (1).x,y,z (2).y,x,z (3).x,y,-z (4).y,x,-z : 1
__________________________________________________________
| Xs : 249763.66400000
| Ys : 629501.11100000
| Zs : 606.08500000
| xs : -28.57280000
| ys : 3.46970000
| zs : 111.77450000
| skala : 54.13590091
| Omega : 0.56125984 g
| Fi : 199.70612996 g
| Kappa : 97.17169452 g
__________________________________________________________
blad om : 5.04368561
Odchylki :
N dX dY dZ dDp
1031 -3.271 2.782 -10.023 10.904
1032 3.371 3.533 13.466 14.324
1033 0.288 0.654 -1.389 1.562
1034 2.973 0.905 -6.365 7.084
1035 3.098 -1.170 2.172 3.960
1030 1.233 -2.305 -0.756 2.721
1029 -0.668 -0.776 -1.271 1.632
1028 -10.712 -1.541 -0.209 10.825
1027 2.332 -1.162 7.007 7.476
1026 1.356 -0.920 -2.630 3.099
Patrząc na powyższe wyniki można stwierdzić że błędy dx i dy mają dość zadowalającą dokładność jak na pierwszy taki pomiar w życiu. Kilkumetrowe odchyłki mogą być spowodowane niedoświadczeniem pomiarowych, błędami celowania na punkty, martwym ruchem śruby oraz błędami przypadkowymi. Widać również że niekiedy punkt pomiarowy był ustawiany za wysoko nad terenem. Być może jest to wina nie wytrenowania oka do stereoskopowego widzenia. Punktem który niewątpliwie wypacza pomiary jest punkt 1028, bowiem ponad 10 metrowa odchyłka px to stanowczo zbyt dużo, postanowiliśmy nie brać pod uwagę go przy wyrównywaniu. Następny plik zilustruje co się będzie działo po wyeliminowaniu punktu 1028 podczas wyrównywania.
>>> ELEMENTY ORIENTACJI BEZWZGLEDNEJ
- Wyniki na podstawie: C:\Documents and Settings\dephos\Pulpit\studenci_2008\pm\mp
Wsp. terenowe (1).X,Y,Z (2).Y,X,Z : 1
Wsp. modelowe (1).x,y,z (2).y,x,z (3).x,y,-z (4).y,x,-z : 1
__________________________________________________________
| Xs : 249808.66000000
| Ys : 629352.01000000
| Zs : 601.07333333
| xs : -31.36377778
| ys : 4.15488889
| zs : 111.86066667
| skala : 54.13263314
| Omega : -199.43915271 g
| Fi : 0.29827147 g
| Kappa : -102.81287978 g
__________________________________________________________
blad om : 4.74977621
Odchylki [m]:
N dX dY dZ dDp
1031 -4.471 1.795 -10.016 11.115
1032 2.362 2.930 13.508 14.023
1033 -0.559 0.502 -1.325 1.523
1034 2.078 1.052 -6.323 6.738
1035 2.438 -0.682 2.264 3.396
1030 -0.134 -1.585 -0.865 1.811
1029 -2.130 -0.419 -1.389 2.577
1027 0.681 -1.710 6.882 7.124
1026 -0.264 -1.884 -2.737 3.333
Współrzędne [m]:
N X Y Z
1 249506.864 631038.086 614.210
2 249619.101 628517.114 520.595
1031 246593.501 628589.355 740.786
1032 248290.878 628224.620 649.802
1033 250245.189 628049.818 494.915
1034 251356.932 628524.038 500.053
1035 252925.072 627948.023 428.866
1030 253144.214 630918.675 436.715
1029 251613.240 630924.499 544.309
1028 249370.991 630844.745 651.526
1027 247843.399 630755.520 827.148
1026 246265.514 630233.544 787.067
Są to wyniki ostateczne które SA sporządzone z bardzo zadowalającą dokładnością.
1