Ilkka Korpela
MML:n kuvat saatiin TIF-muodossa, nelikanavaisina (RGBN), 8 tavua per kanava per pikseli.
Kuvausajankohta oli 29.6.2015 eli kuvat olivat viimeisintä LiDARia (2016) vanhempia, eivätkä siis vastanneet kevään 2018 tilannetta maastossa.
Kuvien numerot olivat (jono/kuva) - 4_96…4_99; 5_144…5_147 and 6_194..6_197.
Huom! Nämä kuvat olivat ortokuvien pohjamateriaalia Paikkatietoikkunassa and Kansalaisen karttapaikassa keväällä 2018. Olli Sirkiältä saatiin myös orientointitiedot.
Figure 1. Information from the image-index map. Note that the order was placed by Suomen Metsäkeskus.
The interior orientation parameters of a calibration certificate
UC kuvien sisäinen orientaatio huomioi pääpisteen paikan kuvatasolla (PPA) ja polttovälin. Koska kuvaa voi kiertää 90-asteen välein, tulee tämä kamerakoordinaatiston kierto huomioida orientointitiedostoja muodostettaessa.
Figure 2. The UCAM pixels are 6x6 microns, the principal point is not quite at the center of the image array but slightly off (120 microns) in the Y-coordinate direction. Focal length is 10.05 cm. Given the 5000 m flying altitude, the pixels are 5000 m * 0.000006 m / 0.1005m = 0.298 m in GSD (ground sampling density)
Figure 3. The camera coordinate system is defined in the image plane with 11310 columns and 17310 rows.
UltraCAMXp kamera koostuu useasta samanaikaisesti valotettavasta kamerasta, joista neljä kuvaa panktromaattista maisemaa ja neljä on eri aallonpituusalueille dedikoituja kameroita. Meillä oli käytössä ns. Level-3 8-bittiset pansharpenoidut RGBN kuvat eli edellisestä laskennallisesti muodostetut tarkat värikuvat (perspektiivikuvat). Kuvista muodostettiin kolmen RGB-kanavan väri- ja väärävärikuvat Matlab-koodilla, joka luki TIF-muotoa ja muutti kuvat raakamuotoon säilyttäen BIL-rakenteen.
imgdata=imread('c:\data\aerial_NLS\20151603A_04_0096.tif');
A = imgdata(:,:,1:4);
B = rot90(A);
fileID = fopen('c:\temp\20151603A_04_0096_rgb.raw','w');
fwrite(fileID,flip(B(:,:,1:3)));
fclose(fileID);
close();
fileID = fopen('c:\temp\20151603A_04_0096_cir.raw','w');
fwrite(fileID,flip(B(:,:,[4,1,2])));
fclose(fileID);
close();
Nämä raakakuvat luettiin Corel ohjelmaan ja tallennus
tehtiin RAW-muodossa BIP-rakenteisena (band interleaved by pixel), jota
käytetty KUVAMITT ohjelma tunnistaa. Normaalisti TIF-kuvat ovat sellaisenaan
käyttökelpoisia eli muunnos tehtiin vain KUVAMITT-ohjelman formaattivaatimusten
takia.
Kuvien ulkoisen orientoinnin tiedot saatiin UTM35/N2000 järjestelmässä. Lentokorkeus
oli 5 km. Yksi kuva kattoi 10310 x 17310 30 cm pikseliä eli 1606 hehtaaria.
The exterior orientation parameters in UTM35/ellipsoidal heights (angles in gradians):
Kuvatunnus X(Itä) Y(Poh) Z(Kork) Omega Phi Kappa OR.TIETOLÄHDE LASK.KOORD.
20151603A_04_0096 355246.891 6859708.574 5162.117 399.67416 0.09041 0.78509 ilmakolmiointi Euref (UTM)
20151603A_04_0097 356604.389 6859724.222 5165.344 399.71915 0.10823 399.95260 ilmakolmiointi Euref (UTM)
20151603A_04_0098 357960.990 6859728.221 5171.124 399.71859 0.08351 0.22373 ilmakolmiointi Euref (UTM)
20151603A_04_0099 359317.921 6859731.873 5170.126 399.76912 0.12310 0.11834 ilmakolmiointi Euref (UTM)
20151603A_05_0144 360674.284 6856345.204 5171.785 0.25271 399.90179 199.80479 ilmakolmiointi Euref (UTM)
20151603A_05_0145 359318.014 6856353.897 5173.860 0.29272 399.90303 199.64506 ilmakolmiointi Euref (UTM)
20151603A_05_0146 357960.063 6856364.442 5172.618 0.29536 399.88529 199.53146 ilmakolmiointi Euref (UTM)
20151603A_05_0147 356602.567 6856360.921 5177.188 0.30301 399.90453 200.34235 ilmakolmiointi Euref (UTM)
20151603A_06_0194 356604.452 6852996.673 5183.292 399.75666 0.08644 0.63406 ilmakolmiointi Euref (UTM)
20151603A_06_0195 357960.076 6853006.860 5181.681 399.69678 0.09283 0.34506 ilmakolmiointi Euref (UTM)
20151603A_06_0196 359316.693 6853009.591 5180.174 399.75101 0.09263 399.65863 ilmakolmiointi Euref (UTM)
20151603A_06_0197 360674.685 6853004.407 5185.060 399.75082 0.11048 399.80702 ilmakolmiointi Euref (UTM)
The exterior orientation parameters in KKJ2/N60
Saatujen kuvien (N=12) ulkoinen orientointi ratkaistiin keväällä 2018 myös KKJ2/N60 järjestelmässä kolmioimalla kuvat s.e. maastosta osoitettiin tukipisteitä KKJ2/N60 järjestelmässä. Näin kaikki ilmakuvat ilmakuva-aikasarjassa saatiin samaan järjestelmään, mikä helpottaa kuvien fotogrammetrista käsittelyä.
a |
b |
c |
|
d |
e |
omega |
phi |
kappa |
X0 |
Y0 |
Z0 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00535 |
0.000565 |
-1.51235 |
2513350.843 |
6859556.402 |
5161.881 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00461 |
0.000968 |
-1.52539 |
2514706.886 |
6859634.234 |
5165.849 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00482 |
0.001181 |
-1.52117 |
2516065.162 |
6859702.075 |
5172.129 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00393 |
0.001918 |
-1.52279 |
2517421.1 |
6859767.596 |
5171.04 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
0.003995 |
-0.00035 |
1.61391 |
2518912.782 |
6856444.546 |
5172.274 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
0.004718 |
-0.00109 |
1.61137 |
2517553.369 |
6856390.182 |
5174.913 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
0.004791 |
-0.00204 |
1.609583 |
2516192.726 |
6856338.035 |
5173.389 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
0.004755 |
-0.002 |
1.622234 |
2514835.368 |
6856272.651 |
5177.191 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00333 |
-9.3E-05 |
-1.51477 |
2515013.295 |
6852909.526 |
5182.481 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00442 |
0.0007 |
-1.51928 |
2516370.946 |
6852982.93 |
5181.992 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00374 |
0.001236 |
-1.5301 |
2517728.816 |
6853048.91 |
5181.004 |
166666.7 |
0 |
8655 |
0 |
166666.7 |
5675 |
-0.00375 |
0.002049 |
-1.52779 |
2519088.658 |
6853106.334 |
5185.324 |