Page 299 -
P. 299

278                                                                6 Feature-based alignment























                Figure 6.3 A simple panograph consisting of three images automatically aligned with a translational model and
                then averaged together.


                                where x i is the consensus (average) position of feature i in the global coordinate frame.
                                (An alternative approach is to register each pair of overlapping images separately and then
                                compute a consensus location for each frame—see Exercise 6.2.)
                                   The above least squares problem is indeterminate (you can add a constant offset to all the
                                frame and point locations t j and x i ). To fix this, either pick one frame as being at the origin
                                or add a constraint to make the average frame offsets be 0.
                                   The formulas for adding rotation and scale transformations are straightforward and are
                                left as an exercise (Exercise 6.2). See if you can create some collages that you would be
                                happy to share with others on the Web.


                                6.1.3 Iterative algorithms

                                While linear least squares is the simplest method for estimating parameters, most problems in
                                computer vision do not have a simple linear relationship between the measurements and the
                                unknowns. In this case, the resulting problem is called non-linear least squares or non-linear
                                regression.
                                   Consider, for example, the problem of estimating a rigid Euclidean 2D transformation
                                (translation plus rotation) between two sets of points. If we parameterize this transformation
                                by the translation amount (t x ,t y ) and the rotation angle θ, as in Table 2.1, the Jacobian of
                                this transformation, given in Table 6.1, depends on the current value of θ. Notice how in
                                Table 6.1, we have re-parameterized the motion matrices so that they are always the identity
                                at the origin p =0, which makes it easier to initialize the motion parameters.
                                   To minimize the non-linear least squares problem, we iteratively find an update Δp to the
                                current parameter estimate p by minimizing

                                                                           2
                                      E NLS (Δp)=        f(x i ; p +Δp) − x                          (6.13)
                                                                         i
                                                      i
                                                                        2

                                                 ≈       J(x i ; p)Δp − r i                          (6.14)
                                                      i
   294   295   296   297   298   299   300   301   302   303   304