Convex Formulation of the Upper Bound Approach with Noise in all Images

In this section, we propose a convex formulation of the principles sketched in section 7.2 that, compared to (144), accounts for noise in both the template and the input images. We can express this in terms of image-plane measurements. As in (164,166), our approach is formulated as an SOCP problem. However, contrary to (164,166), our approach is a point-wise method that does not require us to tune the relative influence of minimizing the reprojection error and maximizing the depths.

Noise in the Template Only

Let us first remark that the basic principles explained in section 7.2 can be formulated as SOCP problems. In this first formulation, the noise is only accounted for in the template image. The inextensibility constraint $ \Vert \mathbf{Q}_i - \mathbf{Q}_j \Vert \leq d_{ij} + \varepsilon_{\mathcal {T}}$ can be written:

$\displaystyle \Vert \mu_i \mathbf{u}_i - \mu_j \mathbf{u}_j \Vert \leq d_{ij} + \varepsilon_{\mathcal {T}}.$ (7.1)

Including the maximization of the depths, we obtain this SOCP problem:

 \displaystyle \max_{\mathbold{\mu}} & &...
...& \mu_i \geq 0 \qquad & & i \in \{ 1,\ldots,n_c \}
 \end{array}\end{displaymath} (7.2)

where  $ \mathbold{\mu}^\mathsf{T}= \begin{pmatrix}\mu_1 & \ldots & \mu_{n_c}\end{pmatrix}$, and  $ \mathcal {E}$ is a set of pairs of points to which the inextensibility constraints are applied.

Noise in Both the Template and the Input Images

Let us now suppose that the inaccuracies are expressed in terms of image-plane measurements. Suppose that points are measured in the image with a maximum error of  $ \varepsilon_{\mathcal {I}}$, i.e.

$\displaystyle \Vert \hat{\mathbf{q}}_i' - \mathbf{q}_i' \Vert \leq \varepsilon_{\mathcal {I}}, \qquad \forall i
 \in \{1, \ldots, n_c\}.$ (7.3)

Since we are searching for the true 3D position of the point  $ \mathbf{Q}_i$, we say that:

$\displaystyle \hat{\mathbf{q}}_i' = \frac{1}{\mathbf{p}_3^\mathsf{T}\mathbf{Q}_...
 \end{pmatrix}.$ (7.4)

Equation (7.3) can thus be rewritten:

$\displaystyle \left\Vert \frac{1}{\mathbf{p}_3^\mathsf{T}\mathbf{Q}_i} \begin{p...
 \end{pmatrix} - \mathbf{q}_i' \right\Vert \leq \varepsilon_{\mathcal {I}}.$ (7.5)

We finally add the inextensibility constraints and the maximization of the depths (which are given by $ \mathbf{p}_3^\mathsf{T}\mathbf{Q}_i$) and we obtain the following SOCP problem:

 \displaystyle \max_{\mathbf{Q}} & & \di...
...hbf{Q}_i \geq 0 & & \forall i \in \{1,\ldots,n_c\}
 \end{array}\end{displaymath} (7.6)

where  $ \mathbf{Q}$ is the concatenation of the 3D points  $ \mathbf{Q}_i$, for $ i \in \{1, \ldots, n_c\}$.

Contributions to Parametric Image Registration and 3D Surface Reconstruction (Ph.D. dissertation, November 2010) - Florent Brunet
Webpage generated on July 2011
PDF version (11 Mo)