[42],Tommaselli and Lugnani(1988) [70], Heikkilä(1991) [32], Kubik(1992) [37], Petsaand Patias(1994) [50], Gülch(1995) [20], Wiman and Axelsson(1996) [76], Chen andShibasaki(1998) [12], Habib(1999) [26], Heuvel(1999) [34], Tommaselli(1999) [71],Vosselman and Veldhuis(1999)[73], Föstner(2000) [17], Smith and Park(2000) [65],Schenk(2002) [58], Tangelder et al.(2003) [68] and Parian and Gruen(2005) [49]. In1988 Mulawa and Mikhail[47] proved the feasibility of linear features for close-rangephotogrammetric applications such as the space intersection, the space resection, relativeorientation and absolute orientation. This was the first step to employ the linearfeature based method in the close-range photogrammetric applications. Mulawa[46]developed linear feature based methods for different sensors. While straight linearfeatures and conic sections can be represented as unique mathematical expressions,free-form lines in nature cannot be described by algebraic equations. Hence Mikhailand Weerawong[44] employed <strong>splines</strong> and polylines to express free-form lines as analyticalexpressions.Tommaselli and Tozzi[72] proposed that the accuracy of thestraight line parameter should be a sub-pixel <strong>with</strong> the representation of four degreesof freedom in an infinite line. Many researchers in photogrammetry have describedstraight lines as infinite lines using minimal representation to reduce unknown parameters.The main consideration of straight line expression is singularities. Habibet al.[29] performed the <strong>bundle</strong> <strong>block</strong> <strong>adjustment</strong> using 3D point set lying on controllinear features instead of traditional control points. EOPs were reconstructed hierarchicallyemploying automatic single photo resection (SPR). SPR was implementedby the relationship between image and object space points expressed as EOPs asfollowing equation (2.8).16
⎡⎢⎣x i − x py i − y p−f⎤⎡⎥⎦ = λR T ⎢(ω, φ, κ) ⎣X i − X oY i − Y oZ i − Z o⎤⎥⎦ (2.8)where λ scale, x i , y i the image coordinates, X i , Y i , Z i the object coordinates, x p , y p , finterior orientation parameters and X o , Y o , Z o , ω, φ, κ exterior orientation parameters.R T is the 3D orthogonal rotation matrix containing the three rotation angles ω, φ andκ.Linear features between image and object space were compared to calculate EOPs bythe modified iterated Hough transform. In the result, EOPs were robustly estimated<strong>with</strong>out the knowledge of the relationship between image and object space primitives.In addition, Habib et al.[28] demonstrated that straight lines in object spacecontributed for performing the <strong>bundle</strong> <strong>block</strong> <strong>adjustment</strong> <strong>with</strong> self-calibration. Theyproposed that a large number of GCPs were reduced in calibration because of linearfeatures having four collinearity equations from two end points defining each straightline. Their experiment proposed the usefulness of linear features in orientation notonly <strong>with</strong> frame aerial imagery, but also in <strong>bundle</strong> <strong>block</strong> <strong>adjustment</strong>.Habib et al.[25] summarized linear features extracted from mobile mapping system,GIS database and maps for various photogrammetric applications such as single photoresection, triangulation, digital camera calibration, image matching, 3D reconstruction,image to image registration and surface to surface registration. Matched linearfeature primitives are utilized in space intersection for reconstruction of object spacefeatures and linear features in the object space are used for control features in triangulationand digital camera calibration. Since linear feature extraction can meetsub-pixel accuracy across the direction of edge and linear features can be extracted17
- Page 2: c○ Copyright byWon Hee Lee2008
- Page 7 and 8: ACKNOWLEDGMENTSThanks be to God, my
- Page 10 and 11: 3. BUNDLE BLOCK ADJUSTMENTWITH 3D N
- Page 13 and 14: CHAPTER 1INTRODUCTION1.1 OverviewOn
- Page 15 and 16: y an intersection employing more th
- Page 17 and 18: similarity of geometric properties
- Page 19 and 20: straight linear features or formula
- Page 21 and 22: • Bundle block adjustment by the
- Page 23 and 24: Hessian. Interest point operators w
- Page 25 and 26: [60], Ebner and Ohlhof(1994) [16],
- Page 27: a complicated problem. The developm
- Page 31 and 32: x p = −f (X A + t · a − X C )r
- Page 33 and 34: surfaces and terrain models in 2D a
- Page 35 and 36: f(u) − e(u) = g(u)f(u) − e(u) =
- Page 37 and 38: Tankovich[69] used linear features
- Page 39 and 40: (a) 0th order continuity (b) 1st or
- Page 41 and 42: Cardinal splineA Cardinal spline is
- Page 43 and 44: 2.3.2 Fourier transformFourier seri
- Page 45 and 46: For other polyline expressions, Aya
- Page 47 and 48: Each segment of a natural cubic spl
- Page 49 and 50: ⎡⎢⎣2 11 4 11 4 1· · ·1 4 1
- Page 51 and 52: 3.2 Extended collinearity equation
- Page 53 and 54: R −1 = R T . The matrix R T (= R
- Page 55 and 56: dx p = M 1 dX C + M 2 dY C + M 3 dZ
- Page 57 and 58: In this research, the arc-length pa
- Page 59 and 60: =√∫ √√√ ()ti+1−f u′ (
- Page 61 and 62: This equation can be replaced with
- Page 63 and 64: order polynomial using Newton’s d
- Page 65 and 66: y collinearity equations, tangents
- Page 67 and 68: d tan(θ t ) = w′ (v ′ w − w
- Page 69 and 70: y each two points, which are four e
- Page 71 and 72: +M 14 db i3 + M 15 dc i0 + M 16 dc
- Page 73 and 74: collinearity model are described in
- Page 75 and 76: [ ] [ ] [ ]N11 N 12 ˆξ1 c1N12T =N
- Page 77 and 78: systematic errors in the image spac
- Page 79 and 80:
interval based on the normal distri
- Page 81 and 82:
1 ∂Φ2 ∂l= (X C + d 1 l − a i
- Page 83 and 84:
about splines, their relationships,
- Page 85 and 86:
cubic spline in the image and the o
- Page 87 and 88:
The redundancy budget of a tie poin
- Page 89 and 90:
of bundle block adjustment is requi
- Page 91 and 92:
ξ kiSP = [ da i0 da i1 da i2 da i3
- Page 93 and 94:
Spline location parametersImage 1 I
- Page 95 and 96:
Spline location parametersImage 1 I
- Page 97 and 98:
5.3 Recovery of EOPs and spline par
- Page 99 and 100:
Table 5.7 expressed the convergence
- Page 101 and 102:
Iteration with an incorrect spline
- Page 103 and 104:
Vertical aerial photographData 9 Ju
- Page 105 and 106:
All locations are assumed as on the
- Page 107 and 108:
of the Gauss-Markov model correspon
- Page 109 and 110:
estimation is obstacled by the corr
- Page 111 and 112:
Interior orientation defines a tran
- Page 113 and 114:
+ fu ( w2 31 (X i (t) − X C ) + s
- Page 115 and 116:
A.2 Derivation of arc-length parame
- Page 117 and 118:
+2f( t [1 + t 2) − 1 22s 12 (Y i
- Page 119 and 120:
+Du ′ ( t 1 + t 22)2r 11 t + Dv
- Page 121 and 122:
A 17 = t [2 − t 1 16 2 f(t 1) −
- Page 123 and 124:
1−u ′ w − w ′ u {w′ [s 21
- Page 125 and 126:
BIBLIOGRAPHY[1] Ackerman, F., and V
- Page 127 and 128:
[24] Haala, N., and G. Vosselman. 1
- Page 129 and 130:
[49] Parian, J.A., and A. Gruen. 20
- Page 131:
[73] Vosselman, G., and H. Veldhuis