<strong>Harita</strong> <strong>Dergisi</strong> Temmuz 2010 Sayı 144F.Yılmaztürk vd.Şekil 7. Test alanı hedef işaretlerinin resim koordinatlarını ölçme işlemi sonrası yazılım görüntüsü.c. Demet DengelemesiSelf kalibrasyonlu demet dengelemesininmatematik modeli, genişletilmiş kolineariteeşitliklerine dayanır.x − xy −ile,⎡R⎢⎢R⎢⎣R123y00Burada,+ ∆ x+ ∆ yR 1= − cRR 2= − cR⎤ ⎡X− X⎥ ⎢⎥= R⎢Y − Y⎥ ⎢⎦ Z − Z⎣003⎤⎥⎥⎥⎦03(2)x, y :Noktanın resim koordinatlarını,x 0 , y0, c :İzdüşümü merkezinin resim koordinatsistemindeki konumunu gösteren içyöneltme parametrelerini,X , Y,Z :Noktanın cisim koordinatlarını,0 0 0X , Y , Z :İzdüşümü merkezinin cisimkoordinatlarını,:Kameranın üç dönme açısıyla(ω, φ, κ) hesaplanan ortogonal dönmematrisini,∆ x, ∆y :Resim koordinatları için düzeltmeterimlerini gösterir.Ek parametrelerin bir setinin fonksiyonları olanresim koordinat düzeltme terimleri, doğadafiziksel olan kolinearite koşulunun idealgeometrisinden olan sapmaları modeller. Busapmalarının 4 ana kaynağı vardır. Bunlarsistematik radyal distorsiyon, teğetsel distorsiyon,görüntü düzleminin düz olmaması ve düzlemdekigörüntü distorsiyonudur. Herhangi bir görüntünoktasının teorik konumundan olan ötelenme bukaynakların her birinin etkisinin toplamı kadar olur(Fraser,vd., 1995, Habib vd., 2003).Geliştirilen yazılımda, iç yöneltme elemanları(x 0 , y 0 , c), mercek distorsiyon katsayıları (k 1 , k 2 ,k 3 ve p 1 , p 2 ), resim koordinat eksenleri arasındakiölçekleme farkı ve resim koordinat eksenlerininbirbirine dik olmamasından kaynaklanandüzeltme terimlerini içeren (b 1 , b 2 ) standart 10terimli fiziksel kalibrasyon modeli kullanılmıştır.(Fraser, 1997). Bu model,56
<strong>Harita</strong> <strong>Dergisi</strong> Temmuz 2010 Sayı 144Sayısal Kameraların Self-Kalibrasyonu için Bir Yazılım Tasarımı∆x= −x+ ( r02∆y= −y06−+ 2x−+ yr k(3)3xc2) p+ 2 p112∆c+ xr k + xr k1+ 2 py 2 4∆c+ yr k1+ yr k2c2x y + ( r242x y + b x + b+ 2y12+ xr) p262ky3düzeltmeleri,∆ X 0 ,..., ∆κ :Dış yöneltme elemanlarınınyaklaşık değerlerine getirilecekdüzeltmeleri,AP i:Ek parametreleri,2σ 0:Öncül varyans faktörü,Σ :Gözlemlerin varyans-kovaryansmatrisini gösterir.ile,r =x2+ y2; x = x − x ; y = y − y00Katsayılar matrisi (A), eşitlik 3’ünbilinmeyenlere göre kısmi türevlerini içerir.Parametrelerin kestirimi en küçük karelerçözümüyle verilir.eşitlikleri ile ifade edilmektedir.T −1Txˆ= ( A PA)A Pl(5)Demet dengelemesinde, resim üzerindeölçülen her bir nokta için iki kolinearite eşitliğiyazılır. Cisim noktalarının karşılıklı bütün resimnoktaları için yazılan eşitliklerin toplamı,bilinmeyen parametrelerin belirlenmesinimümkün kılan bir denklem sistemiyle sonuçlanır.Eşitlikler bilinmeyenlere göre lineer olmayan birdenklem sistemidir ve en küçük kareleryöntemiyle çözüm için 1. Derece Taylor serisikullanılarak lineer hale getirilmektedir. Kontrolnoktalarının ölçülen cisim koordinatları ve ekparametreler de gözlem denklemleri formundamodele dâhil edilir. Ek parametrelerin sayısı veseçimi, standart sapmaları kontrol edilerekuygulamanın gereksinimlerine bağlı olarakdeğiştirilebilmektedir. Gözlem denklemleri sistemiGauss-Markov modelinde aşağıdaki gibi formüleedilir:− e = Ax − lP = σile,x T20Σ−1(4)[ ∆X∆Y,∆Z,∆X, ∆Y, ∆Z, ∆ω,∆ϕ,∆ , ∆ ... ]= κ, 0 0 0AP1Burada,e :Gözlem hataları vektörünü,A :Katsayılar matrisini,x :Bilinmeyen parametreler vektörünü,l :Gözlem vektörünü,P :Ağırlık matrisini,∆ X , ∆Y,∆Z:Bir noktanın cisim koordinatlarınınyaklaşık değerlerine getirilecekLineer olmayan demet dengelemesi problemilineer problemlerin bir dizisi olarak çözülür. Heriterasyonda bir düzeltme vektörü kestirilir ve budeğerler yaklaşık değerlere eklenerekparametreler güncellenir (Mikhail, vd., 2001;Wolf, vd., 2000). Gözlemlerin düzeltmeleri vesoncul varyans faktörü aşağıdaki eşitliklerkullanılarak hesaplanır.v = Axˆ− l2σ 0PV= ;rv Tr = n − u(6)Burada, r serbestlik derecesini, n gözlemdenklemlerinin sayısını ve u da bilinmeyenlerinsayısını göstermektedir.Yazılımda iterasyonu durdurma kriteri olarak,ardışık iterasyonlar arasında birim ağırlıklıölçünün ortalama hata değerleri ( σ 0 ) arasındakifark için girilen bir sınır değer kullanılmıştır. Busınır değerle birlikte maksimum iterasyon sayısı,resim koordinatlarını ölçme prezisyonu ve kabahata ölçü testinde kullanılacak yönteminbelirlenmesinden oluşan dengeleme kontrolelemanları dengeleme öncesi demetdengelemesi diyalog penceresinde tanımlanır(Şekil 8). Dengeleme sonrası sonuçları, yine budiyalog penceresi üzerinde özet olarakgösterilirken dengeleme sonuçlarının ayrıntılıgösterimi bir metin dosyası biçiminde bilgisayarınproje klasörü içine yazdırılır.57