Uniwersytet Warmińsko - Mazurski w Olsztynie 20.11.2012 rok
Wydział Geodezji i Gospodarki Przestrzennej
Katedra Fotogrametrii i Teledetekcji
SPRAWOZDANIE
ANALITYCZNA ANALIZA POJEDYNCZEGO
ZDJĘCIA LOTNICZEGO
Wykonała:
Martyna Stępkowska
Grupa 2
GiSzN, rok III
Spis Treści
1. Przestrzenny rysunek zdjęcia lotniczego.
2. Zawartość metryki kalibracji kamery.
3. Obliczenie położenia punktu I, N, Z wg otrzymanych danych i wykonanie rysunku zdjęcia.
4. Wykonanie obliczeń i analiza wpływu nachylenia ν na przesunięcie punktów na zdjęciu.
5. Wykonanie obliczeń i analiza wpływu deniwelacji terenu Δh na przesunięcie punktów na zdjęciu.
2. Zawartość metryki kalibracji.
Aktualna metryka jest warunkiem dopuszczenia kamery do wykonywania pomiarowych zdjęć lotniczych.
Metryka kalibracji kamery powinna zawierać następujące informacje:
nazwa i adres instytucji wykonującej kalibrację,
data kalibracji,
nazwa i numer fabryczny kamery,
numer fabryczny obiektywu,
kalibrowana odległość obrazowa stożka obiektywowego,
dystorsja radialna w mikrometrach wyrażona w funkcji promienia radialnego o początku w punkcie najlepszej symetrii (PPS) w interwale co 10 mm, wzdłuż każdej z czterech półprzekątnych zdjęcia,
odległość między znaczkami tłowymi wzdłuż boków i przekątnych, współrzędne w prostokątnym układzie współrzędnych,
położenie punktu głównego autokolimacji (PPA) i punktu głównego najlepszej symetrii (PPS) w układzie współrzędnych związanym ze znaczkami tłowymi,
zdolność rozdzielczą obiektywu radialną i tangencjalną, określoną przez producenta kamery.
Metryka kalibracji kamery może zawierać również inne, dodatkowe parametry, określone w procesie kalibracji, a mianowicie:
matematyczny model dystorsji radialnej i tangencjalnej wraz z wartościami parametrów występujących w tym modelu, wyznaczonych w procesie kalibracji,
określenie równoległości płaszczyzn filtra (filtrów) używanego wraz z kamerą,
wyniki kalibracji migawki zawierające efektywne czasy ekspozycji i sprawność świetlną migawki dla ciągu nominalnych ekspozycji,
płaskość płyty wypłaszczającej film,
fotograficzną zdolność rozdzielczą dla układu: obiektyw kamery wraz z filmem,
kalibrowaną odległość obrazową dla układu: obiektyw kamery wraz z filmem,
rodzaj filmu użytego w procesie kalibracji,
dystorsję radialną i tangencjalną dla układu: obiektyw kamery wraz z filmem,
płaskość płyty dociskowej poprzez analityczne opracowanie modelu przestrzennego zdjęć kolimatorów zarejestrowanych przez układ: obiektyw kamery wraz z filmem,
ekscentr źrenicy wejściowej obiektywu względem ramki tłowej (dla zastosowań dGPS).
3. Obliczenie położenia punktu I, N, Z wg otrzymanych danych i wykonanie rysunku zdjęcia.
I’ = tg$\frac{\mathbf{\nu}}{\mathbf{2}}$ ck
Ck1[mm] | 91 |
---|---|
ν1 | 1o17’ |
I | 1,02 |
Ck2[mm] | 301 |
---|---|
ν1 | 1o17’ |
I | 3,37 |
Ck1[mm] | 91 |
---|---|
ν2 | 4o17’ |
I | 3,40 |
Ck2[mm] | 301 |
---|---|
ν2 | 4o17’ |
I | 11,26 |
N’ = tgν ck
Ck1[mm] | 91 |
---|---|
ν1 | 1o17’ |
N | 2,04 |
Ck2[mm] | 301 |
---|---|
ν1 | 1o17’ |
N | 6,74 |
Ck1[mm] | 91 |
---|---|
ν2 | 4o17’ |
N | 6,82 |
Ck2[mm] | 301 |
---|---|
ν2 | 4o17’ |
N | 22,54 |
Z’ = ctgν ck
Ck1[mm] | 91 |
---|---|
ν1 | 1o17’ |
Z | 4062,11 |
Ck2[mm] | 301 |
---|---|
ν1 | 1o17’ |
Z | 13436,22 |
Ck1[mm] | 91 |
---|---|
ν2 | 4o17’ |
Z | 1214,99 |
Ck2[mm] | 301 |
---|---|
ν2 | 4o17’ |
Z | 4018,81 |
4. Wykonanie obliczeń i analiza wpływu nachylenia ν na przesunięcie punktów na zdjęciu.
δν = $\frac{{\mathbf{r}\mathbf{'}}^{\mathbf{2}}}{\mathbf{c}_{\mathbf{k}}}$ sinν sinγ
r'1 [mm] | 22 |
---|---|
r'2 [mm] | 107 |
Ck1 [mm] | 91 |
Ck2 [mm] | 301 |
ν1 | 1o17’ |
ν2 | 4o17’ |
γ1 [o] | 0,0 |
---|---|
γ2 [o] | 22,5 |
γ3 [o] | 45,0 |
γ4 [o] | 67,5 |
γ5 [o] | 90,0 |
δν obliczamy dla ck1, ν1, r'1 i zmiennych γ.
Ck1[mm] | 91 | 91 | 91 | 91 | 91 |
---|---|---|---|---|---|
ν1 | 1o17’ | 1o17’ | 1o17’ | 1o17’ | 1o17’ |
r'1 [mm] | 22 | 22 | 22 | 22 | 22 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 0,046 | 0,084 | 0,110 | 0,119 |
δν obliczamy dla ck1, ν2, r'1 i zmiennego γ.
Ck1[mm] | 91 | 91 | 91 | 91 | 91 |
---|---|---|---|---|---|
ν2 | 4o17’ | 4o17’ | 4o17’ | 4o17’ | 4o17’ |
r'1 [mm] | 22 | 22 | 22 | 22 | 22 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 0,152 | 0,281 | 0,367 | 0,397 |
δν obliczamy dla ck2, ν1, r'1 i zmiennego γ.
Ck2[mm] | 301 | 301 | 301 | 301 | 301 |
---|---|---|---|---|---|
ν1 | 1o17’ | 1o17’ | 1o17’ | 1o17’ | 1o17’ |
r'1 [mm] | 22 | 22 | 22 | 22 | 22 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 0,014 | 0,025 | 0,033 | 0,036 |
δν obliczamy dla ck2, ν2, r'1 i zmiennego γ.
Ck2[mm] | 301 | 301 | 301 | 301 | 301 |
---|---|---|---|---|---|
ν2 | 4o17’ | 4o17’ | 4o17’ | 4o17’ | 4o17’ |
r'1 [mm] | 22 | 22 | 22 | 22 | 22 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 0,046 | 0,085 | 0,111 | 0,120 |
δν obliczamy dla ck1, ν1, r'2 i zmiennego γ.
Ck1[mm] | 91 | 91 | 91 | 91 | 91 |
---|---|---|---|---|---|
ν1 | 1o17’ | 1o17’ | 1o17’ | 1o17’ | 1o17’ |
r'2 [mm] | 107 | 107 | 107 | 107 | 107 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 1,078 | 1,992 | 2,603 | 2,818 |
δν obliczamy dla ck1, ν2, r'2 i zmiennego γ.
Ck1[mm] | 91 | 91 | 91 | 91 | 91 |
---|---|---|---|---|---|
ν2 | 4o17’ | 4o17’ | 4o17’ | 4o17’ | 4o17’ |
r'2 [mm] | 107 | 107 | 107 | 107 | 107 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 3,596 | 6,645 | 8,682 | 9,397 |
δν obliczamy dla ck2, ν1, r'2 i zmiennego γ.
Ck2[mm] | 301 | 301 | 301 | 301 | 301 |
---|---|---|---|---|---|
ν1 | 1o17’ | 1o17’ | 1o17’ | 1o17’ | 1o17’ |
r'2 [mm] | 107 | 107 | 107 | 107 | 107 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 0,326 | 0,602 | 0,787 | 0,852 |
δν obliczamy dla ck2, ν2, r'2 i zmiennego γ.
Ck2[mm] | 301 | 301 | 301 | 301 | 301 |
---|---|---|---|---|---|
ν2 | 4o17’ | 4o17’ | 4o17’ | 4o17’ | 4o17’ |
r'2 [mm] | 107 | 107 | 107 | 107 | 107 |
γi [o] | 0,00 | 22,5 | 45,0 | 67,5 | 90,0 |
δν[mm] | 0,000 | 1,087 | 2,009 | 2,625 | 2,841 |
WNIOSKI
Przesunięcie δν zależy od wartości jakie przyjmują takie parametry jak: kąt v, stała kamery ck, promień radialny r’ oraz położenie punktu na zdjęciu γ. Biorąc pod uwagę parametr ck, wnioskujemy, że przy jego wzroście maleje wartość przesunięcia i odwrotnie.
Dla ν2=4o17’, r'2=107, γ2=22,5o stała:
- gdy ck1=91mm gdy to przesunięcie δh wynosi 3,596 mm;
- gdy ck2=301mm gdy to przesunięcie δh wynosi 1,087 mm.
Wartości pozostałych parametrów rosną proporcjonalnie do przesunięcia czyli:
Dla ν2=4o17’, ck2=301mm, γ3=45,0o:
- gdy r'1=22 to przesunięcie δh wynosi 0,085mm;
- gdy r'2=107 to przesunięcie δh wynosi 2,009mm.
Dla r'1=22’, ck2=301mm, γ3=45,0o:
- gdy ν1=1o17’ to przesunięcie δh wynosi 0,025mm;
- gdy ν2=4o17’ to przesunięcie δh wynosi 0,085mm.
Dla r'2=107’, ck2=301mm, ν2=4o17’:
- gdy γ3=45,0o to przesunięcie δh wynosi 2,009mm;
- gdy γ5=90,0o to przesunięcie δh wynosi 2,841mm.
Gdy mamy dane cK, v, r’ to przesunięcie będzie zależało od położenia punktu na zdjęciu czyli γ . Największą wartość przesunięcia uzyskamy, gdy γ będzie miała 90o zaś najmniejszą gdy 0o.
Biorąc pod uwagę kąt γ=0o to bez względu na inne parametry przesunięcie będzie zerowe.
Maksymalne przesunięcie o wartości 9,397 mm dla parametrów podanych w ćwiczeniu uzyskamy biorąc pod uwagę największy promień radialny równy107mm, najmniejsze cK = 91 mm oraz największy kąt v=4o17’o wartości 9,397 mm.
5. Analiza wpływu deniwelacji terenu Δh na przesunięcie punktów na zdjęciu
Δh1 [mm] | 3000,00 |
---|---|
Δh2 [mm] | 6000,00 |
Δh3 [mm] | 30000,00 |
Mz | 5000 |
r'1 [mm] | 22 |
r'2 [mm] | 107 |
δh = $\frac{\mathbf{r}^{\mathbf{'}}\mathbf{\text{Δh}}}{\mathbf{w}}$ w = Mz*ck
w1 = Mz*ck1=455000
w2 = Mz*ck2=1505000
r’1[mm] | 22 | 22 | 22 |
---|---|---|---|
w1[mm] | 455000 | 455000 | 455000 |
Δh[mm] | 3000 | 6000 | 30000 |
δh[mm] | 0,145 | 0,290 | 1,451 |
r’1[mm] | 22 | 22 | 22 |
---|---|---|---|
w2[mm] | 1505000 | 1505000 | 1505000 |
Δh[mm] | 3000 | 6000 | 30000 |
δh[mm] | 0,044 | 0,088 | 0,439 |
r’2 [mm] | 107 | 107 | 107 |
---|---|---|---|
w1[mm] | 455000 | 455000 | 455000 |
Δh[mm] | 3000 | 6000 | 30000 |
δh [mm] | 0,705 | 1,411 | 7,055 |
r’2 [mm] | 107 | 107 | 107 |
---|---|---|---|
w2 [mm] | 1505000 | 1505000 | 1505000 |
Δh [mm] | 3000 | 6000 | 30000 |
δh[mm] | 0,213 | 0,427 | 2,133 |
WNIOSKI
Przesunięcie δh zmniejsza się wraz ze wzrostem wysokości fotografowania.
Dla r’2=107mm i Δh=3m przy w2 =1505000mm (δh =0,044 mm) jest mniejsze przesunięcie niż przy w1=455000mm (δh =0,145mm). w2 > w1
Dla r’1=22mm i Δh=3m przy w2=1505000mm (δh =0,427 mm) jest mniejsze przesunięcie niż przy w1=455000mm (δh =1,411mm). w2 > w1
Na przesunięcie δh wpływa również deniwelacja terenu. Przesunięcie rośnie proporcjonalnie do Δh.
Dla r’1=22 i w1=455000mm :
- gdy Δh=3m to przesunięcie δh wynosi 0,145mm;
- gdy Δh=6m to przesunięcie δh wynosi 0,290mm;
- gdy Δh=30m to przesunięcie δh wynosi 0, 1,451mm.
Dla r’2=107 i w1=455000mm :
- gdy Δh=3m to przesunięcie δh wynosi 0,705mm;
- gdy Δh=6m to przesunięcie δh wynosi 1,411mm;
- gdy Δh=30m to przesunięcie δh wynosi 7,055mm.