anbul 2004 
and flying 
ct field) is 
‘à to the 
nption that 
r exposure 
(4) 
je camera 
2» Ya Z5), 
position is 
International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Vol XXXV, Part B5. Istanbul 2004 
3.1.2 Initial Value of the Temporal GCPs 
The Church method is one of methods to calculate exterior 
orientation parameters, and exterior orientation parameters of a 
tilted photo are calculated using 3 GCPs which have accurate 
3D ground coordinates (Wolf, 1983). 
In Figure 4, images p;, p, and p; of GCPs P;, P, and P4. 
Angle @ , #,. 7 are 
yt 
cosy = —%— 
OR -OP, 
Ky 7 (Xo - Xi XX - X2) (Y - Y Ko - 12) (7.1) 
“(20-2120 7) 
  
Kp 
COS zz E 
OP, - OP, 
Kg - (Xo - X5 XXo — X3) (9 - Y Jos - 55) (7.2) 
* (Zo - ZaXZo - Z3) 
K 
COSE m mum em 
OR OR 
K, z( Xs za oza)2s09 2509 -A) (7.3) 
*(Zo - Z3 YZo- A) 
where: 
OR «(Xo - X «(5 - Wy -(Zo- AY 
OP, «(Xy - XY «(5 - Yo «Zo - Z2) 
OB, «(Xo - Xs) «(9 - 5 «(Zo - Z3 
On the other hand, cosa , cos, and cosy are expressed with 
respect to focal length and image coordinates using the 
following geometry equations: 
P 
Xp Xp2 + Yi V p? x 
  
  
  
  
  
  
cosa = (8.1) 
Op; : Op; 
2 
cos ß = Xp2'Xp3 * Yp2 Vp3 +f (8.2) 
Op» - Opi 
2 
X508 piu Vom iY pif : 
cos = ne / (8.3) 
Ops - Op, 
In equation (8), Op, , Op; , and Op; are length from the 
camera position to the respective image point, and cosa , 
cos//, and cosy in equation (7) are given by equation (9). 
Therefore, equation (7) consists of 3 unknowns (Xo, Yo, Zo). 
and equation (7) are solved with respect to unknowns using 
Taylor's theorem. 
However, X5, Y5, and Z4 are unknowns in this paper. Therefore, 
in order to obtain these unknowns, the Church method was used 
inversely under assumption Za coordinate equal 0, and solved 
unknowns are given as initial values for the temporal GCP of 
P. 
Similarly, initial values for the other temporal GCPs are 
obtained by repeating the same procedures. 
3.1.3 Initial Values of Orientation Parameters 
The initial values of exterior orientation parameters are given as 
calculated values from the space resection using initial values 
for the temporal GCPs, which were obtained by the above 
procedures. 
3.2 Calibration Procedures 
In the above procedure, initial values for the exterior orientation 
parameters and 3D coordinates for the 6 temporal GCPs are 
given. Therefore, unknown parameters for the interior/exterior 
orientation parameters and 3D coordinates for temporal GCPs 
are calculated by the resection using collinearity conditions and 
space distances. 
equation (9) shows collinearity conditions, and equation (10) 
shows distance conditions, and these equations for the both 
images are used simultaneously in calibration procedures. 
Consequently, 36 observation equations are acquired for 30 
unknown parameters, and unknown parameters are calibrated as 
the values that make the following function # minimum under 
the least-squares method (equation (11)) Furthermore, 3D 
coordinates for additional feature points are able to compute 
simultaneously. 
Collinearity condition: 
RE (x X f + kır?) 
ny HX — Xo myj(Y m + m3(Z = Zo) 
m3 (X = Xo)+ my, (Y = Yo) + mys (Z - Zg) 
G= (y — yg y + kır?) 
. M9] (X TE Xo) Ma y TE Yo nips (z = Zo) 
S ml X - Xo) my(Y - Yo) m3(Z - Zo) 
where: x, y are image coordinates of temporal GCPs, r is 
  
(9) 
  
distance from principal point to image point, k; is lens 
distortion parameter, and m;; are rotation matrixes. 
Observation condition for the space distance: 
  
D,- A - Xy «(C - Xy «(Z, -Zy (10) 
Function with weight: 
n=2 n=6 , \ n=2 n=6 
If SS P (ax? * Av + SS nn — min (11) 
i=l j=l i=l j=l 
where: Ax, Ay : residual for the image coordinates, AD : 
residual for space distance. 
Weight values in equation (11) should be carefully considered 
under statistically or measurement accuracy, and calibration 
results were investigated using various weights. However, 
calibration results don't show any significant improvements 
compared with various weight values. Therefore, equal weight 
( py = po =1) was adopted in this paper. 
4. PERFORMANCE EVALUATION 
In order to evaluate an accuracy of the IBIM, experiment was 
performed. Figure 5 shows test site, 20 circular points in Figure 
5 are temporal GCPs and check points. À stereo image was 
taken with 21.790m (left side), 21.957m flying height (right 
side) respectively and 8.44m base line under the fixed focal 
length. The 3D coordinates for all circular points and