Hassan Afzal
email: hassan.afzal@uni.lu

Djamila Aouada

Michel Antunes

David Fofi

Bruno Mirbach

Björn Ottersten

Bi-objective Framework for Sensor Fusion in RGB-D Multi-View Systems: Applications in Calibration

Complete and textured 3D reconstruction of dynamic scenes has been facilitated by mapped RGB and depth information acquired by RGB-D cameras based multi-view systems. One of the most critical steps in such multi-view systems is to determine the relative poses of all cameras via a process known as extrinsic calibration. In this work, we propose a sensor fusion framework based on a weighted bi-objective optimization for refinement of extrinsic calibration tailored for RGB-D multi-view systems. The weighted bi-objective cost function, which makes use of 2D information from RGB images and 3D information from depth images, is analytically derived via the Maximum Likelihood (ML) method. The weighting factor appears as a function of noise in 2D and 3D measurements and takes into account the affect of residual errors on the optimization. We propose an iterative scheme to estimate noise variances in 2D and 3D measurements, for simultaneously computing the weighting factor together with the camera poses. An extensive quantitative and qualitative evaluation of the proposed approach shows improved calibration performance as compared to refinement schemes which use only 2D or 3D measurement information.

Introduction

RGB-D cameras provide simultaneous image and range data of the environment, offering enhanced sensing capabilities when compared to using single sensor modality. A non-exhaustive list of applications includes 3D telepresence systems [START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF], creation of viewpoint free 3D videos [START_REF] Kuster | FreeCam: A Hybrid Camera System for Interactive Free-Viewpoint Video[END_REF], simultaneous localization and mapping [START_REF] Henry | RGB-D mapping: Using kinect-style depth cameras for dense 3D modeling of indoor environments[END_REF], or the acquisition of textured 3D surface models of static and dynamic scenes [START_REF] Newcombe | KinectFusion: Real-time Dense Surface Mapping and Tracking[END_REF][START_REF] Kainz | OmniKinect: real-time dense volumetric data acquisition and applications[END_REF][START_REF] Sturm | CopyMe3D: Scanning and printing persons in 3D[END_REF][START_REF] Dou | Temporally enhanced 3D capture of room-sized dynamic scenes with commodity depth cameras[END_REF].

The acquisition of complete and textured 3D models of scenes required in domains such as security and surveillance, health, and entertainment, can be accomplished by using two different approaches. The first consists of using a single moving RGB-D camera with its location constantly being tracked [START_REF] Henry | RGB-D mapping: Using kinect-style depth cameras for dense 3D modeling of indoor environments[END_REF][START_REF] Newcombe | KinectFusion: Real-time Dense Surface Mapping and Tracking[END_REF][START_REF] Sturm | CopyMe3D: Scanning and printing persons in 3D[END_REF]. This solution is simple and attractive, however, has the drawback of not allowing to fully reconstruct dynamic scenes at each time-step. As an example, Dou et al. [START_REF] Dou | 3D Scanning Deformable Objects With a Single RGBD Sensor[END_REF] presented a 3D scanning system for deformable objects using a single Kinect sensor. The scanning results are promising, remark that, however, the obtained 3D reconstructions are rigid, even if the objects in the environment were constantly moving and deforming. This issue can be solved by using multiple fixed RGB-D cameras covering the entire scene [START_REF] Kainz | OmniKinect: real-time dense volumetric data acquisition and applications[END_REF][START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF][START_REF] Dou | Temporally enhanced 3D capture of room-sized dynamic scenes with commodity depth cameras[END_REF][START_REF] Palasek | A flexible calibration method of multiple Kinects for 3D human reconstruction[END_REF]. In this case, the relative poses of all RGB-D cameras are required for aligning the partial 3D reconstructions. The 1 arXiv:1905.09939v1 [cs.CV] 23 May 2019 problem of estimating the relative poses of cameras in a multi-view system is known as extrinsic calibration.

Most of the works for extrinsic calibration of RGB-D multi-view systems rely on well established 2D camera based calibration routines [START_REF] Zhang | Flexible camera calibration by viewing a plane from unknown orientations[END_REF]10] and pose refinement procedures, e.g. Bundle Adjustment (BA) [START_REF] Triggs | Bundle Adjustment -A Modern Synthesis[END_REF][START_REF] Amplianitis | Calibration of a Multiple Stereo and Rgb-D Camera System for 3d Human Tracking[END_REF][START_REF] Dou | Temporally enhanced 3D capture of room-sized dynamic scenes with commodity depth cameras[END_REF], using 2D feature points extracted from the RGB images [START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF][START_REF] Berger | Markerless Motion Capture using multiple Color-Depth Sensors[END_REF][START_REF] Berger | The capturing of turbulent gas flows using multiple Kinects[END_REF]. The 3D information from the depth sensor has mainly been used in subsequent refinement steps using, e.g., the Iterative Closest Point (ICP) algorithm [START_REF] Segal | Generalized-ICP[END_REF][START_REF] Penelle | Geometrical 3D reconstruction using real-time RGB-D cameras[END_REF][START_REF] Yang | Multi-Kinect scene reconstruction: Calibration and depth inconsistencies[END_REF]. In this regard, the following question arises: how to optimally use both sources of complementary information.

In this paper, we investigate a strategy for RGB-D sensor fusion for the extrinsic calibration of multiple cameras. Instead of using 2D data from RGB images and 3D data from depth images independently, we propose a weighted bi-objective optimization scheme. We analytically derive a Least Squares (LS) based cost function, via the Maximum Likelihood (ML) method, that optimally combines the BA based 2D cost function with the ICP based 3D cost function. The sensor fusion is achieved by using a weighting factor that depends on two types of noise, one contaminating the 2D feature locations in the RGB images, and the second one contaminating the 3D point positions provided by the depth sensor. The experiments suggest that using the proposed joint cost for relative pose refinement provides more accurate results than the refinement schemes using 2D and 3D information separately.

In the absence of information regarding noise levels in the 2D and 3D feature points we propose an iterative scheme which simultaneously estimates the noise along with the estimation of calibration parameters. The proposed scheme is completely automated requiring no manual intervention and no heuristic parameter setting. The quantitative and qualitative experiments show that the proposed scheme is able to perform sensor fusion for accurate camera calibration without any prior information about noise characteristics.

The present work extends and consolidates our previous work called BAICP+ [START_REF] Afzal | RGB-D Multi-view System Calibration for Full 3D Scene Reconstruction[END_REF], which experimentally showed that in many cases, using a heuristically constructed weighted bi-objective refinement approach that combines 2D and 3D information provides better results than refinement approaches based on cost functions using only 2D or 3D information.

Contributions

• We present a sensor fusion framework based on weighted bi-objective optimization for refinement of extrinsic calibration of an RGB-D multi-view system. We derive an analytic expression for the weighting factor, in the bi-objective optimization, in terms of noise in measurements of RGB and depth sensors.

• We propose an iterative scheme for extrinsic calibration of an RGB-D multi-view system, which alternates between camera pose estimation, and the computation of measurement noise levels.

• We perform a thorough experimental evaluation on synthetic and real data, and show that fusing the RGB-D information using the proposed bi-objective optimization provides superior results when compared to refinement schemes that only use 2D or 3D feature information.

Article Overview

We start by giving a brief overview of state-of-the-art methods for extrinsic calibration used in RGB-D multi-view systems, as well as of bi-objective pose estimation approaches in Section 2. In Section 3, the extrinsic calibration problem is formally presented. Section 4 gives a brief introduction of BA and ICP together with our previous work i.e., BAICP+ [START_REF] Afzal | RGB-D Multi-view System Calibration for Full 3D Scene Reconstruction[END_REF], in which BA and ICP are heuristically combined. Section 5 analytically derives the expression for the weighted bi-objective cost function for refinement of extrinsic calibration parameters. Section 6 presents an automated iterative approach for camera pose estimation and estimation of the measurement noise parameters. In Section 7 and Section 8, we analyze and illustrate the benefits of the proposed approach via extensive experiments using synthetic and real data respectively. This is followed by a conclusion in Section 9.

Related Work

In this section we review state-of-the-art techniques for extrinsic calibration in an RGB-D multiview system with a focus on the modality of data used. We also briefly overview the sensor fusion approaches, based on bi-objective optimization, for solving the pair-wise pose estimation problem.

A considerable amount of research has been, and is still being, carried out in the domain of RGB-D cameras based multi-view systems [START_REF] Berger | A State of the Art Report on Multiple RGB-D Sensor Research and on Publicly Available RGB-D Datasets[END_REF]. We are interested in analyzing the extrinsic calibration method used in such multi-view systems. Extrinsic calibration in a multi-view system requires information about the same points (feature points), in 3D space, to be acquired in different views. This information can be extracted from 2D or 3D camera acquisitions of objects with, e.g., known textural and/or geometrical properties [START_REF] Berger | Markerless Motion Capture using multiple Color-Depth Sensors[END_REF][START_REF] Kainz | OmniKinect: real-time dense volumetric data acquisition and applications[END_REF][START_REF] Alexiadis | Reconstruction for 3D immersive virtual environments[END_REF][START_REF] Miller | Unsupervised extrinsic calibration of depth sensors in dynamic scenes[END_REF][START_REF] Palasek | A flexible calibration method of multiple Kinects for 3D human reconstruction[END_REF]. A planar checkerboard pattern first proposed by Zhang [START_REF] Zhang | Flexible camera calibration by viewing a plane from unknown orientations[END_REF] is one of the most widely used objects for this purpose. Assuming correct mapping of RGB and depth sensors, Maimone and Fuchs [START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF], and Yang et al. [START_REF] Yang | Multi-Kinect scene reconstruction: Calibration and depth inconsistencies[END_REF] make use of such an object, with corners extracted from RGB and IR images respectively, to calibrate their systems with the help of Bougouet's calibration toolbox [10]. Berger et al. [START_REF] Berger | The capturing of turbulent gas flows using multiple Kinects[END_REF][START_REF] Berger | Markerless Motion Capture using multiple Color-Depth Sensors[END_REF], on the other hand, try to improve accuracy by calibrating all RGB and depth sensors together by extracting 2D feature points from both RGB and IR images with the help of a special checkerboard pattern consisting of diffuse and mirroring patches.

A major drawback of 2D only calibration approaches is their inability to tackle noise specific to depth sensors which causes problems in alignment of 3D data from multiple cameras. Various methods try to tackle this problem with the help of an explicit depth correction step based on comparing known and measured depths for each camera, e.g., [START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF][START_REF] Kainz | OmniKinect: real-time dense volumetric data acquisition and applications[END_REF][START_REF] Nakazawa | Dynamic scene reconstruction using asynchronous multiple Kinects[END_REF][START_REF] Dou | Temporally enhanced 3D capture of room-sized dynamic scenes with commodity depth cameras[END_REF] etc. There are other methods, such as [START_REF] Yang | Multi-Kinect scene reconstruction: Calibration and depth inconsistencies[END_REF], which add a final refinement step based on the ICP algorithm, using 3D data only, which tries to mitigate the pose misalignment due to depth specific sensor noise. Penelle et al. [START_REF] Penelle | Geometrical 3D reconstruction using real-time RGB-D cameras[END_REF] propose to use only 3D points, corresponding to 2D feature points (from checkerboard) extracted from RGB images, in the ICP based calibration scheme where initial estimates are provided via the RANSAC algorithm. Nakazawa et al. [START_REF] Nakazawa | Dynamic scene reconstruction using asynchronous multiple Kinects[END_REF], on the other hand, use 2D feature points in the BA algorithm for pose refinement but perform corrections using alignment of 3D data as well. Miller et al. [START_REF] Miller | Unsupervised extrinsic calibration of depth sensors in dynamic scenes[END_REF] propose to use 3D information related to foreground objects to obtain an occlusion aware energy minimization scheme for auto-calibration. Deng et al. [START_REF] Teng | Registration of Multiple RGBD Cameras via Local Rigid Transformations[END_REF], use 3D feature points observed at different locations in the scene to compute a smooth field of rigid transformation instead of a single rigid transformation to align the measurements of two RGB-D cameras.

Dou and Fuchs [START_REF] Dou | Temporally enhanced 3D capture of room-sized dynamic scenes with commodity depth cameras[END_REF], in their work on multi-view 3D reconstruction, proposed to combine 2D and 3D information in a weighted bi-objective optimization scheme derived from their previous work on pair-wise pose tracking for mono-view 3D reconstruction [START_REF] Dou | Exploring High-Level Plane Primitives for Indoor 3D Reconstruction with a Hand-held RGB-D Camera[END_REF]. They propose to use matching feature points extracted via Scale-Invariant Feature Transform (SIFT) from RGB images with matching planes extracted from 3D/depth images in a weighted bi-objective BA scheme. The weighting factor is selected empirically for all experiments. A similar approach is proposed by Henry et al. [START_REF] Henry | RGB-D mapping: Using kinect-style depth cameras for dense 3D modeling of indoor environments[END_REF], using a global ICP scheme to align 2D visual feature points and 3D/depth measurements from multiple views but the weights are, again, selected empirically. Tykkala et al. [START_REF] Tykkala | Direct Iterative Closest Point for real-time visual odometry[END_REF] use what they call an image based direct ICP approach for pairwise pose estimation. They propose to compute the weighting factor via a heuristic measure using ratio of the median intensity and the depth values of selected points. Michot et al. [START_REF] Michot | Bi-Objective Bundle Adjustment With Application to Multi-Sensor SLAM[END_REF] propose to use a weighted bi-objective BA scheme for the multi-sensor Simultaneous Localization and Matching (SLAM) problem. They discuss the dependence of the weighting factor on the ratio of the noise variance for each sensor's measurement and formulate their bi-objective optimization by using a Mean Squared Error (MSE) based cost function from individual sensors. They investigate three methods for automatic weight computation namely L-Curve, L-Tangent Norm and cross validation with experiments showing that the L-Curve based method performs better than the others.

This brief survey shows that most work done on extrinsic calibration of RGB-D multi-view systems is based on calibration schemes which use 2D or 3D information independently. These methods do not take into account the relative accuracy of 2D and 3D measurement in a systematic manner to, e.g., give more importance to less noisy measurements.

In this work, we perform sensor fusion by formulating a weighted bi-objective optimization scheme based on 2D and 3D cost functions for performing global refinement of poses in RGB-D multi-view systems. We draw upon work in the Robotics domain [START_REF] Michot | Bi-Objective Bundle Adjustment With Application to Multi-Sensor SLAM[END_REF][START_REF] Tykkala | Direct Iterative Closest Point for real-time visual odometry[END_REF][START_REF] Han | Bi-objective optimization for robust RGB-D visual odometry[END_REF][START_REF] Dou | Exploring High-Level Plane Primitives for Indoor 3D Reconstruction with a Hand-held RGB-D Camera[END_REF], where several sensor modalities are exploited for pair-wise pose estimation. The proposed bi-objective optimiza-tion uses cost functions from both the BA and ICP algorithms, based on 2D and 3D measurements, in a single unified cost function. The key to combining the two cost functions is the information about noise in the 2D and 3D measurements which is reflected in the weighting factor. In the absence of this information we propose to use a simple approach which iteratively estimates the noise parameters given the current camera poses and vice versa. We show the validity of the proposed approach by achieving improved results in the experiments performed under different conditions.

Extrinsic Calibration Problem of Multiple RGB-D Sensors

Notation: The following notation will be adopted estimates of p and T l , respectively.

I n : identity matrix of dimension n × n. 0 n : null vector of dimension n × 1.
In this section, we will formulate the extrinsic calibration problem for an RGB-D multi-view system. Let us consider a multi-view system composed of N , intrinsically calibrated, RGB-D cameras with intersecting FOVs, as shown in Fig. 1. Every RGB-D camera l, with l = 1, • • • , N , acquires an RGB image C l and a 3D vertex map V l , with associated known matrix of intrinsic parameters K l .

In order to correctly align the partial 3D reconstructions {V l }, where l = 1, • • • , N , acquired by N RGB-D cameras, it is necessary to accurately estimate their positions with respect to a global reference frame, referred to as world and denoted by w, as shown in Fig. 1. Each camera's relative position with respect to w is defined by:

T l = R l t l 0 3 1 , (1)
where T l ∈ SE(3) represents the rigid transformation, from camera l to w. The matrix R l is rotation matrix in SO(3) and t l ∈ R 3 is translation vector. Therefore the same point p ∈ R 3 in w viewed by camera l as p l and by cameras k as p k can be related to the cameras' reference frames as follows:

R l p l + t l = R k p k + t k . (2)
Similarly, for a given point x ∈ R 3 in w, its projection on each camera's image plane results in 2D pixel coordinates q l , such that:

q l = ψ (K l , T l , x) , ∀l, (3)
where ψ(.) is world to image plane projection function.

Let us assume that all cameras are of resolution M . The color pixel positions in C l may be represented by the points q

i l ∈ [q 1 l , • • • , q M l ] where i ∈ {1, • • • , M }. Similarly, the 3D coordinates in V l may be represented by points p k l ∈ [p 1 l , • • • , p M l ] where k ∈ {1, • • • , M }.
The problem at hand may therefore be stated as follows. Given N RGB-D cameras in a multiview system with acquired RGB images

{C 1 , • • • , C N } and 3D vertex maps {V 1 , • • • , V N },
we assume knowledge of H ≤ M matching points in each camera's RGB image plane referred to as 2D feature points and denoted as

[q 1 l , • • • , q H l ]
. Similarly, we assume knowledge of J ≤ M matching 3D points in each camera's 3D vertex map called 3D feature points and denoted as

[p 1 l , • • • , p J l ]. Moreover we assume knowledge of each camera's intrinsic parameters, K = [K 1 , • • • , K N ]
. Using this information, we want to find the estimates of the parameters

T = [T 1 , • • • , T N ].

Background

In this section, we introduce the two pose refinement algorithms namely BA and ICP, from which we derive the 2D and 3D cost functions used in our proposed bi-objective cost function. We also introduce BAICP+ [START_REF] Afzal | RGB-D Multi-view System Calibration for Full 3D Scene Reconstruction[END_REF], in which we proposed a bi-objective pose refinement scheme, but, with heuristically defined weight and scaling parameters.

Bundle Adjustment (BA)

For refinement of extrinsic calibration parameters using 2D feature points only, we use a cost function from the Bundle Adjustment (BA) algorithm [START_REF] Triggs | Bundle Adjustment -A Modern Synthesis[END_REF]. It has been the method of choice for problems related to multi-view 3D reconstruction and pose refinement based on 2D feature points extracted from RGB images [START_REF] Triggs | Bundle Adjustment -A Modern Synthesis[END_REF]. Bundle Adjustment (BA) requires an initial estimate of the pose parameters. Moreover, it also requires an estimate of 3D points i.e., [

x 1 , • • • , x H ], corresponding to available 2D feature points [q 1 l , • • • , q H l ]
. These estimates are then refined by computing the error of projection of estimate of each 3D point x h , h = 1, • • • , H, corresponding to the 2D feature point q h l to camera l via:

a h l (S h l ) = q h l -ψ K l , T l , x h , (4)
where a h l (S h l ) ∈ R 2 and S h l = T l , x h1 . Therefore, the total BA cost to be minimized for the refinement of estimates of each camera's pose parameters together with the estimates of 3D points corresponding to 2D feature points is given as:

V BA (S) = N l=1 tr(A l (S l )A l (S l )), (5)
where

S = (T, X), S l = (T l , X), X = [x 1 , • • • , x H ] and A l (S l ) = [a 1 l (S 1 l ), • • • , a H l (S H l )].

Iterative Closest Point (ICP)

For refinement of the extrinsic calibration parameters using 3D feature points only, we use the cost function from the global Iterative Closest Points (ICP) algorithm [START_REF] Besl | A method for registration of 3-d shapes[END_REF]. ICP algorithm has been the de facto solution for pose refinement problems when only 3D feature points are available [START_REF] Nüchter | Study of parameterizations for the rigid body transformations of the scan registration problem[END_REF][START_REF] Besl | A method for registration of 3-d shapes[END_REF]. ICP algorithm also uses initial estimates of the pose parameters and minimizes the Euclidean distance between corresponding 3D feature points from different views, such that:

b j l,k (T l , T k ) = (R l p j l + t l ) -(R k p j k + t k ), (6)
where

b j l,k (T l , T k ) ∈ R 3 and j ∈ [1, • • • , J].
Therefore, the total ICP cost to be minimized for refinement of each camera's pose parameters is given as:

V ICP (T) = 1<l,k<N l =k tr(B l,k (T l , T k )B l,k (T l , T k )), (7)
where

B l,k (T l , T k ) = [b 1 l,k (T l , T k ), • • • , b J l,k (T l , T k )].

BAICP+

In our previous work [START_REF] Afzal | RGB-D Multi-view System Calibration for Full 3D Scene Reconstruction[END_REF], we used the information provided by RGB-D cameras in a bi-objective optimization scheme for extrinsic calibration refinement. It combines the BA and ICP cost functions based on 2D and 3D feature points respectively resulting in the cost function for BAICP+:

V BAICP (S) = (1 -c) a V ICP (T) + sc b V BA (S), (8)
where c ∈ [0, 1] is the weighting factor and s = ( avgDepth avgF ocal ) 2 is a heuristic scaling factor, based on the ratio of the average depth of points in x versus average focal length across all views. The parameters a and b denote the total number of 3D point correspondences and 2D feature points points across all views.

Bi-Objective Extrinsic Calibration

In this section, we present the bi-objective optimization for refinement of the extrinsic calibration parameters in an RGB-D multi-view system. We use cost functions defined in the previous section which use 2D and 3D feature points extracted from RGB images and vertex maps, respectively.

In this work, we propose to formally analyze and derive an expression for the cost function, based on ML estimations, of the bi-objective optimization taking into account the noise affecting both 2D and 3D measurement/feature points. We assume the presence of independent additive Gaussian noise in each coordinate of the 3D feature points such that:

pj l ∼ N p j l , σ 2 3D I 3 , (9)
where pj l is the noisy 3D point and p j l is the noise free point. Similarly for 2D feature points we have:

qh l ∼ N (q h l , σ 2 2D I 2 ), (10)
where qh l is the noisy 2D point and q h l is the noise free point. This means that we have to use the noisy 2D and 3D feature points to estimate the pose parameters. This leads to redefining the 3D error function b j l,m (T l , T m ), given in [START_REF] Berger | A State of the Art Report on Multiple RGB-D Sensor Research and on Publicly Available RGB-D Datasets[END_REF], such that it computes the error between noisy points pj l and pm l projected to w, from camera l and m, using the pose parameters T l and T m , respectively. Similarly the 2D error function a h l (S h l ), given in (4) where S h l = (T l , x h ), is redefined such that it computes the 2D error between back projection of the estimated 3D point x h to camera l, using T l and K l , and the corresponding noisy 2D feature point qh l . Now, we can define the distribution the 3D error b j l,m (T l , T m ) is drawn from by considering the noise free 3D points p j l and p j m in (2) such that [START_REF] Segal | Generalized-ICP[END_REF]:

b j l,m (T l , T m ) ∼ N (R l p j l + t l ) -(R m p j m + t m ), R l σ 2 3D I 3 R l + R m σ 2 3D I 3 R m = N (0 3 , 2σ 2 3D I 3 ). (11)
Similarly, considering the noise free 2D measurements in (3), we have a h l (S h l ) ∼ N (0 2 , σ 2 2D I 2 ). It is clear from [START_REF] Dou | Temporally enhanced 3D capture of room-sized dynamic scenes with commodity depth cameras[END_REF] that since b j l,m (T l , T m ), which is based on the ICP algorithm, uses two noisy 3D feature points, hence, the variance of the corresponding distribution is two times the variance of noise in each 3D feature point. This is in contrast to the variance of distribution corresponding to a h l (S h l ), which is based on the BA algorithm and uses only one noisy 2D feature point [START_REF] Michot | Bi-Objective Bundle Adjustment With Application to Multi-Sensor SLAM[END_REF]. Using b j l,m (T l , T m ) and a h l (S h l ), we want to find the likelihood cost function, maximum of which gives the Maximum Likelihood Estimate (MLE) of the parameters S = (T, x). Since the MLE with Gaussian model is equivalent to the Least Squares Estimate (LSE) [START_REF] Prince | Computer Vision: Models, Learning, and Inference[END_REF], we can directly get:

Ŝ = arg min S 1<l,m<N l =m 1 2σ 2 3D tr B l,m (T l , T m )B l,m (T l , T m ) + N l=1 1 σ 2 2D tr (A l (S l )A l (S l )) . (12)
Therefore, the total cost to be minimized is:

V (S) = V ICP (T) + wV BA (S), (13)
where w =

2σ 2 3D σ 2 2D
is the weighting factor. The cost function in [START_REF] Dou | Exploring High-Level Plane Primitives for Indoor 3D Reconstruction with a Hand-held RGB-D Camera[END_REF] optimally combines information from RGB and depth sensors, to be used in the pose refinement scheme, by taking into account the noise levels in the 2D and 3D points. It formally defines the the relationship of measurement noise in the 2D and 3D feature points with the weighting factor w. In case the assumption of noise with same variances affecting all 2D and 3D points respectively, does not hold and information about the noise variances affecting each point is available, it can be incorporated in the proposed framework. Moreover, the use of the ICP based cost also allows the use of all the 3D points acquired by each sensor (with the help of nearest neighbor correspondence) in the optimization scheme when only 2D feature points are available.

The cost function ( 13) is a non-linear function of the parameters S and we resort to numerical search methods [START_REF] Myung | Tutorial on maximum likelihood estimation[END_REF] to optimize the criterion. Please refer to Appendix 11 for further discussion.

Weighting Factor Estimation

In this section we discuss the automatic and simultaneous estimation of the weighting factor w in [START_REF] Dou | 3D Scanning Deformable Objects With a Single RGBD Sensor[END_REF], together with the camera poses in the absence of information regarding noise affecting both the 2D and 3D measurements. We propose an approach which alternates between camera pose estimation and estimation of the 2D and the 3D noise variances to arrive at a suitable solution.

In the previous section, the estimates of camera pose parameters and 3D points in w corresponding to 2D feature points, were computed based on known 2D and 3D feature points and the noise affecting them. We assumed the presence of Gaussian noise with zero mean and variances of σ 2 2D and σ 2 3D in 2D and 3D measurements, respectively. These parameters, in turn, define the weighting factor w which is instrumental in constructing the sensor fusion framework by optimally combining the 2D and 3D cost functions to estimate the camera poses. In real-world scenarios, however, information about the noise affecting one or both sensor measurements is often unavailable. This makes the computation of a correct w difficult. As mentioned in Section 2, researchers have tried to estimate the optimal weighting factor, for their proposed bi-objective schemes, for solving mainly the pair-wise pose estimation problem. The commonly used used methods range from using simple heuristic measures such as in the case of [START_REF] Tykkala | Direct Iterative Closest Point for real-time visual odometry[END_REF] to more complex methods, based on analysis of trade-off between residuals of two cost functions and based on learning via cross-validation, such as in the case of [START_REF] Michot | Bi-Objective Bundle Adjustment With Application to Multi-Sensor SLAM[END_REF].

In this work, we propose to use a simple method for automatic estimation of the weighting factor w which finds its basis in finding the MLE of noise variances, σ 2 2D and σ 2 3D , using the 2D and 3D feature points together with the current estimates of camera poses and 3D points in Ŝ. The MLE of the variance σ 2 3D is given as [START_REF] Prince | Computer Vision: Models, Learning, and Inference[END_REF]:

σ2 3D = 1<l,m<N l =m tr(B l,m (T l , T m )B l,m (T l , T m )) 2a , ( 14
)
where a is the total number of 3D feature correspondences across all views. Similarly, the MLE of the variance σ 2 2D is computed via:

σ2 2D = N l=1 tr(A l (S l )A l (S l )) b , ( 15
)
where b is the total number of 2D feature points found across all views. We follow an iterative approach whereby using 2D and 3D feature points and an initial estimate Ŝ, the MLE estimates of noise variances and hence of w are obtained via ( 14) and ( 15). This initial estimate of w is then used to find an updated estimate of S using (13) via non-linear optimization which, in turn, is used to update the estimate of w. This process is repeated for a fixed number of iterations until the estimates of S and w converge.

Experiments with Synthetic Data

In this section, we carry out a quantitative performance analysis of the proposed bi-objective refinement with a known and an unknown weighting factor.

Evaluation Methodology and Parameters

We use V-REP [1] to simulate 2 and 4 cameras based RGB-D multi-view systems, with overlapping FOVs, as shown in Fig. 2 and Fig. 1, respectively. In both cases, the global reference frame w lies in camera C1. We simulate a scene containing several objects such as chairs, a table, sofas etc. The acquired noise-free data, in the form of RGB and depth images, is assumed to be perfectly mapped in each camera's RGB sensor's reference frame with known intrinsics. After data acquisition, random points, visible to all cameras, are extracted as feature points in both RGB and depth images as shown in Fig. 3 (points on the floor are discarded). Features extracted from depth maps are converted to the corresponding 3D points via known intrinsics.

In the next step, noise is added to the extracted 2D and 3D feature points. We assume the presence of independent Gaussian noise in each coordinate of position of 2D feature points with zero mean and standard deviation σ 2D similar to [START_REF] Zhang | Extrinsic calibration of a camera and laser range finder[END_REF]. The value of σ 2D is varied between 0.2 to 1.8 pixels with a step size of 0.4 pixels. Depth sensor measurements in RGB-D cameras suffer from different types of systematic and non-systematic errors as investigated in [START_REF] Smisek | 3D with Kinect[END_REF][START_REF] Fankhauser | Kinect v2 for Mobile Robot Navigation: Evaluation and Modeling[END_REF]. For our scheme we propose to counter, beforehand, the systematic errors in depth measurements of each camera via a correction step, based on comparing known and measured depths [START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF][START_REF] Nakazawa | Dynamic scene reconstruction using asynchronous multiple Kinects[END_REF]. Therefore, for all remaining errors we assume the presence of additive independent Gaussian noise in each -Initial pose obtained using a DLT like approach (2D feature points and corresponding 3D points are used) [10,[START_REF] Hartley | Multiple View Geometry in Computer Vision[END_REF]; ICP -refinement of Init. using Iterative Closest Point (only 3D feature points are used); Ourrefinement of Init. using our bi-objective optimization with known w (2D and 3D feature points are used); OurAuto -refinement of Init. using our bi-objective optimization with unknown w (2D and 3D feature points are used); BA -refinement of Init. using Bundle adjustment (only 2D feature points are used). noise is added to the data, being the variance of the 3D noise fixed (σ 3D = 18mm), and the 2D noise σ 2D is varied between 0.2 and 1.8 pixels (horizontal axes). coordinate of 3D feature points in each view with zero mean and standard deviation σ 3D . The value of σ 3D is varied between 6 to 30 mm with a step size of 6 mm to keep it in the range of errors computed in [START_REF] Fankhauser | Kinect v2 for Mobile Robot Navigation: Evaluation and Modeling[END_REF]. We test the performance of the proposed scheme under various conditions by varying the number of cameras and their positions as shown in Fig. 2 and Fig. ??, by varying the noise magnitude in 2D and 3D feature points as explained above, and by varying the number of 2D and 3D feature points. For each configuration, 50 noise realizations are generated. For each noise realization, 2D feature points and their corresponding noisy 3D measurements from vertex maps are used to initialize the pose estimates via a Direct Linear Transform (DLT) based approach [10,[START_REF] Hartley | Multiple View Geometry in Computer Vision[END_REF]. Using the initial pose estimates, optimization is carried out via the proposed scheme, with known noise parameters as explained in Section 5, and with unknown noise parameters using the automatic iterative estimation scheme as explained in Section 6 (required 3 iterations to converge in most cases). Furthermore, optimization is also carried out via ICP algorithm using 3D feature points only, and via BA algorithm using 2D feature points only.

Accuracy of the estimated poses is computed by comparison with the ground truth poses as done in [START_REF] Zhang | Extrinsic calibration of a camera and laser range finder[END_REF]. Two measures of accuracy are computed. First is the angular magnitude of residual rotation computed via RT l R l , and second is the relative translation error which is computed via

tl -t l t l
. The results of 50 realizations showing the accuracy, of each initialization and of each refinement approach, for each configuration are plotted by using the function boxplot in MATLAB as shown in Fig. 4 -12. The horizontal line inside each box marks the median, the edges mark the 25th and the 75th percentiles, the whisker edges show most extreme data points with outliers plotted separately as red crosses.

The implementation of the proposed bi-objective optimization scheme and ICP is based on the non-linear optimization via Levenberg Marquardt (LM) algorithm [START_REF] Moré | The Levenberg-Marquardt algorithm: Implementation and theory[END_REF], while the implementation of BA is based on a sparse variant of the LM algorithm called Sparse Bundle Adjustment (SBA) [START_REF] Rabaud | Vincent's Structure from Motion Toolbox[END_REF][START_REF] Triggs | Bundle Adjustment -A Modern Synthesis[END_REF].

System Composed of Two Sensors

This section compares the performance of the proposed bi-objective optimization scheme, with known and unknown weighting factor, ICP and BA for refinement of camera pose parameters in a two camera setup shown in Fig. 2. The pose of camera C2 with respect to camera C1 is estimated. After initialization, pose refinement is carried out using the four refinement methods and results are plotted in Fig. 45678.

Varying Noise Levels:

In this experiment, the extrinsic calibration is carried out using 100 2D feature points and 100 3D feature points. Fig. 4 shows the error distribution for fixed 3D noise and varying 2D noise, while Fig. 5 shows the distribution in case the 2D noise is kept fixed, and the 3D noise is varied.

As expected, the accuracy of the extrinsic calibration decreases with increasing noise levels. Also, all pose refinement approaches are able to improve the initial pose estimates, explained by the fact that only inlier data points are generated (no wrong matching feature points are included). A careful analysis of the results shows that our bi-objective optimization scheme with known w, which uses simultaneously the 2D and 3D data, provides better pose estimations when compared to ICP and BA, where only 3D feature points and 2D feature points are used, respectively. The weighting factor based on the noise variance information in [START_REF] Dou | 3D Scanning Deformable Objects With a Single RGBD Sensor[END_REF] automatically gives prominence to more reliable data, decreasing the impact of the other sensor modality. Moreover, it shows that our proposed automatic iterative estimation scheme used in the absence of information regarding noise parameters, and hence unknown w, is robust and also more accurate when compared to BA and ICP, and in most cases nearly as accurate as the method with known w.

Varying Data Ratio:

In this experiment, the extrinsic calibration is carried out using fixed noise variance (σ 2D = 1pix, σ 3D = 18mm). Fig. 6 shows the error distribution for a fixed number of 3D points and a varying number of 2D points, while Fig. 7 shows the distribution in case the 2D points are kept fixed, and the number of 3D points is varied. Since the initial poses are obtained by using 2D feature points and their corresponding 3D points, the initialization varies in Fig. 6 as number of 2D feature points vary but stays approximately the same in Fig. 7 as the number of 2D feature points remain fixed. The conclusions drawn in the previous section regarding improved accuracy of the proposed approaches hold, and these results show that the proposed scheme generalizes for different ratios between the number of 2D and 3D points. Increasing the number of data points of one of the sensor modalities always improves the extrinsic calibration accuracy for the algorithms using those modalities. Moreover, the results in Fig. 4, Fig. 5, Fig. 6, and Fig. 7 show the increased robustness of the proposed approach and ICP to bad initialization as compared to BA.

System Composed of Four Sensors

This section compares performance of the proposed bi-objective optimization scheme with ICP and BA for refinement of camera pose parameters in a four camera setup shown in Fig. ??. The poses of cameras C2, C3 and C4 are aligned with camera C1. After initialization, pose refinement is carried out using the four refinement methods and results are plotted.

Varying Noise Levels:

In this experiment, the extrinsic calibration is carried out using 100 2D feature points and 100 3D feature points. Fig. 8 shows the mean error distribution for computed poses of all cameras, for fixed 3D noise and varying 2D noise. Fig. 9 shows the mean distribution in case the 2D noise is fixed. These results again show the improved performance of the proposed approaches due to the use of both 2D and 3D information together, with the help of correct weighting factor. The performance of all methods gets affected as the noise in 2D and 3D data increases. These results also show improvement in performance of all methods as compared to the multi-view system composed of two cameras due to increased number of 2D and 3D points available. Moreover, these results show that the proposed scheme generalizes for different numbers of cameras used in the multi-view system. We also notice an interesting behavior where in some cases the proposed automatic iterative scheme based on alternative computation of camera poses and w gives better results compared to the scheme with known w. Apart from increase in the number of measurements per feature point, a reason for this can be that for the case of known w we are assuming that for all the 2D and 3D feature points the variances of noise affecting them are the same and constant; but depending on a particular realization, the noise will be a bit higher or lower than the fixed value. Therefore the automatic procedure which tries to compute the variances directly from the noisy data is, in many cases, better able to capture the noise characteristics. For BA, Fig. 9 shows a decrease in its performance as the 3D noise increases. The reason being that apart from its dependence on the initial camera poses, the initial guess of the 3D points corresponding to 2D feature points also gets worse due to increased 3D noise.

In Fig. 10, we compare the mean error distribution with error distributions of individual cameras for the single case of 2D and 3D noise variance (σ 2D = 1pix, σ 3D = 18mm). These results show that while the initial guess for camera C3 is comparatively worse, the performance of optimization schemes is comparable across all views.

Varying Number of Points:

In this experiment, the extrinsic calibration is carried out using a fixed noise variance (σ 2D = 1pix, σ 3D = 18mm). Fig. 11 shows the mean error distribution for a fixed number of 3D points and a varying number of 2D points. Fig. 12, on the other hand, shows the mean distribution in case the 2D points are kept fixed, and the number of 3D points are varied. Here, again, the conclusions drawn in the previous sections hold, while also showing that increasing the number of data points of one of the sensor modalities always improves the extrinsic calibration accuracy for the methods using those modalities.

Experiments with Real Data

In this section, we carry out a qualitative performance analysis of the proposed bi-objective refinement scheme using a real setup. Our setup consists of 4 Asus Xtion Pro Live cameras [START_REF]Xtion PRO[END_REF] with their positions shown in Fig. 13. Each camera acquires an RGB image and a depth image which is mapped to the RGB image. The first step is to perform intrinsic calibration to find the intrinsic and distortion parameters for each camera. For this purpose, we use the method proposed by Zhang [START_REF] Zhang | Flexible camera calibration by viewing a plane from unknown orientations[END_REF] which uses 2D corners extracted from RGB images of a checkerboard pattern viewed at different poses to compute these parameters [10]. As mentioned before, the measurements of these RGB-D cameras suffer from inherent depth bias. Therefore, we perform a depth bias correction procedure, similar to the one used in [START_REF] Nakazawa | Dynamic scene reconstruction using asynchronous multiple Kinects[END_REF], for each camera separately. This procedure requires placing the camera at known distances away from an object (a plane in our case). Using known and measured depth values, we estimate the coefficients of a polynomial which computes the depth correction as a function of measured depth value. These coefficients are unique to each camera and, hence, are used to correct the depth measurements acquired by that camera. After intrinsic calibration and depth bias correction, the next step is to perform the extrinsic calibration using the proposed bi-objective scheme. We first need to extract matching 2D and 3D feature points using RGB and depth images acquired by all 4 cameras. We again use different views of a (two-sided) planar checkerboard pattern as shown in Fig. 13 and extract matching corners from RGB images to be used as 2D feature points and use the corresponding depth values from depth images to get the 3D feature points. The 3D feature points are filtered via a plane detection approach based on RANSAC algorithm to remove outliers if any exist. For this experiment, only 59 2D and 3D feature points were used. The initial pose estimates are generated in the same manner as explained in Section 7, via a Direct Linear Transform (DLT) based approach [10,[START_REF] Hartley | Multiple View Geometry in Computer Vision[END_REF]. These initial poses are then refined via the proposed iterative pose estimation and weight estimation approach explained in Section 6, BA and ICP. Once the refined poses are obtained, they can be used to produce full, textured, 3D reconstructions using data acquired by all 4 cameras as shown in Fig. 14. A qualitative comparison of 3D reconstructions obtained via different calibration methods is shown in Fig. 15. It can be seen that the partial reconstructions are better aligned using the proposed method, which means that the quality of the extrinsic calibration is superior when compared to the other approaches. Note that we are only showing the alignment of the partial point clouds, and no Fig. 13. Multi-view system consisting of 4 Asus Xtion Pro Live Cameras C1, C2, C3 and C4 mounted on a ceiling lift. This system is used to acquire measurements of a real scene. A two-sided planar checkerboard calibration pattern used to extract feature points is also shown.

post-processing step such as smoothing or meshing are applied. We chose to do so to better assess, visually, the accuracy of the extrinsic calibration.

Conclusion

In this paper we have proposed a framework for RGB and depth sensor fusion based on biobjective optimization, for refinement of extrinsic calibration in RGB-D multi-view systems. Our bi-objective optimization scheme makes use of a cost function from the BA algorithm for 2D feature points extracted from RGB images and a cost function from the ICP algorithm for 3D feature points extracted from depth images. We analytically derive an expression for the weighted biobjective cost function. It also analytically relates the weighing factor to the noise in the 2D and 3D measurements, thus making the cost function free of any parameter that needs to be tuned. In case the information regarding measurement noise in 2D and 3D data is not available, we propose an iterative scheme which alternates between estimation of noise parameters assuming known poses, and estimation of camera poses assuming known noise parameters. Thus, it enables us to automatically compute the correct weighting factor when information about measurement noise is not available. A thorough investigation of the performance of the proposed approach for both synthetic and real data showed improved accuracy compared to refinement schemes which only use 2D or 3D information, and comparative performance of proposed approaches with known and unknown noise parameters. These experiments also showed the invariance of the proposed approach under various conditions which include varying the number and position of cameras, varying the 2D and 3D noise and varying the number of the 2D and 3D feature points.

Acknowledgment

This work was supported by the National Research Fund (FNR), Luxembourg, under the CORE project C11/BM/1204105/FAVE/Ottersten. Rearranging [START_REF] Kuster | FreeCam: A Hybrid Camera System for Interactive Free-Viewpoint Video[END_REF], we get the parameter update rule as:

∆S = -(J T B J B + wJ T A J A ) -1 (J T B B + wJ T A A). (20)
We can also rearrange [START_REF] Maimone | Encumbrance-free telepresence system with real-time 3D capture and display using commodity depth cameras[END_REF] according to Levenberg-Marquardt LM [START_REF] Moré | The Levenberg-Marquardt algorithm: Implementation and theory[END_REF] algorithm get the parameter update rule as:

(( 1 2σ 2 3D J T B J B + 1 σ 2 2D J T A J A )+ λdiag( 1 2σ 2 3D J T B J B + 1 σ 2 2D J T A J A ))∆S = -( 1 2σ 2 3D J T B B + 1 σ 2 2D J T A A), ( 21
)
where λ is the damping factor.

Fig. 1 .

 1 Fig. 1. RGB-D Multi-View System with full scene 3D reconstruction in a simulated setup. (a) RGB-D Multi-View System (4 cameras) with field of view (FOV) of each camera. The highlighted region represents overlapping FOVs of all cameras. The global reference frame w is aligned with camera C1. (b) Steps required for Full 3D Scene Reconstruction using an RGB-D Multi-View System. Each camera acquires a RGB image and a depth image, which are used to estimate the relative pose of each camera with respect to w. After extrinsic calibration, estimated poses are used to put all acquisitions in w to get complete reconstruction.

Fig. 2 .

 2 Fig. 2. RGB-D Multi-View System (2 cameras) with field of view (FOV) of each camera. The highlighted region represents overlapping FOVs of all cameras. The global reference frame w is aligned with camera C1.

Fig. 3 .

 3 Fig.3. Features extracted from RGB and depth images of camera C1 in the multi-view system composed of 2 cameras as shown in Fig.2. The extracted feature points are also visible to camera C2.

Fig. 4 .

 4 Fig.4. Error distribution of pose estimates for camera C2 in a two camera setup. 100 2D and 100 3D feature points are used. The following methods are compared: Init. -Initial pose obtained using a DLT like approach (2D feature points and corresponding 3D points are used)[10,[START_REF] Hartley | Multiple View Geometry in Computer Vision[END_REF]; ICP -refinement of Init. using Iterative Closest Point (only 3D feature points are used); Ourrefinement of Init. using our bi-objective optimization with known w (2D and 3D feature points are used); OurAuto -refinement of Init. using our bi-objective optimization with unknown w (2D and 3D feature points are used); BA -refinement of Init. using Bundle adjustment (only 2D feature points are used). noise is added to the data, being the variance of the 3D noise fixed (σ 3D = 18mm), and the 2D noise σ 2D is varied between 0.2 and 1.8 pixels (horizontal axes).

Fig. 5 .

 5 Fig. 5. Error distribution of pose estimates for camera C2 in a two camera setup. 100 2D and 100 3D feature points are used. Gaussian noise is added to the data, being the variance of the 2D noise fixed (σ 2D = 1pix), and the 3D noise σ 3D is varied between 8mm and 30mm (horizontal axes).

Fig. 6 .

 6 Fig. 6. Error distribution of pose estimates for camera C2 in a two camera setup. Gaussian noise is added to the data (σ 2D = 1pix, σ 3D = 18mm), 250 3D feature points and a varying number of 2D feature points (horizontal axes) is considered.

Fig. 7 .

 7 Fig. 7. Error distribution of pose estimates for camera C2 in a two camera setup. Gaussian noise is added to the data (σ 2D = 1pix),σ 3D = 18mm), 250 2D feature points and a varying number of 3D feature points (horizontal axes) is considered.

Fig. 8 .

 8 Fig. 8. Mean error distribution, of pose estimates for cameras C2, C3 and C4. in a four camera setup. 100 2D and 100 3D feature points are used. Gaussian noise is added to the data, being the variance of the 3D noise fixed (σ 3D = 18mm), and the 2D noise σ 2D is varied between 0.2 and 1.8 pixels (horizontal axes).

Fig. 9 .

 9 Fig. 9. Mean Error distribution, of pose estimates for cameras C2, C3 and C4. in a four camera setup. 100 2D and 100 3D feature points are used. Gaussian noise is added to the data, being the variance of the 2D noise fixed (σ 2D = 1pix), and the 3D noise σ 3D is varied between 8mm and 30mm (horizontal axes).

Fig. 10 .

 10 Fig. 10. Comparison of error distributions, of the extrinsic calibration of a four camera setup, using 100 2D and 100 3D feature points. The results are based on mean error distribution and error distribution for camera C2, camera C3 and camera C4. Gaussian noise is added to the data, being the variance of the both 2D noise and 3D noise fixed (σ 2D = 1pix, σ 3D = 18mm).

Fig. 11 .

 11 Fig. 11. Mean error distribution, of pose estimates for cameras C2, C3 and C4. in a four camera setup. Gaussian noise is added to the data (σ 2D = 1pix), σ 3D = 18mm), 250 3D feature points and a varying number of 2D feature points (horizontal axes) is considered.

Fig. 12 .

 12 Fig. 12. Mean error distribution, of pose estimates for cameras C2, C3 and C4. in a four camera setup. Gaussian noise is added to the data (σ 2D = 1pix), σ 3D = 18mm), 250 2D feature points and a varying number of 3D feature points (horizontal axes) is considered.

Fig. 14 .

 14 Fig. 14. 3D reconstruction of a human using a real scene acquired from the multi-view system shown in Fig. 13. Acquisition: Each of the 4 cameras acquire an RGB image and a depth image. 3D Reconstruction: Point clouds based 3D reconstruction using pose estimates refined by the proposed bi-objective scheme with the help of automated weighting.

Fig. 15 .

 15 Fig.[START_REF] Han | Bi-objective optimization for robust RGB-D visual odometry[END_REF]. Comparison of 3D reconstructions of a human using a real scene as shown in Fig.13, via different calibration methods namely Init., ICP, BA and OurAuto described in Fig.5. The acquisitions from cameras C1, C2, C3 and C4 are assigned the colors red, green, blue and magenta, respectively. Misalignments are highlighted via black boxes. Top Row shows side view of the 3D reconstruction and misalignment of views in the results of Init. and BA can be seen clearly, while Bottom Row shows the frontal view and misalignment of views in the results of Init., ICP and BA are visible. It can also be seen that OurAuto gives better results compared to the other methods.

  . Subscripts indicate camera or reference frame indexes and superscripts indicate point indexes:

	w :	world reference frame.
	A :	matrix.
	p :	vector.
	l, M :	scalars.
	tr(A) : trace of matrix A.
	A :	transpose of matrix A.

ψ(•), : w to image plane projection. p, Tl :

BA can also refine the estimate of intrinsics K l if required

(a) Mean Rotation Error (b) Mean Relative Translation Error

Appendix

Non-linear Optimization for Proposed Bi-Objective Framework

Due to the non-linear dependence of cost function in [START_REF] Dou | 3D Scanning Deformable Objects With a Single RGBD Sensor[END_REF] on parameters in S, the MLE Ŝ is to be computed via a numerical scheme based on non-linear optimization. In this scheme at every iteration a small change is introduced in the current set of parameters leading to comparatively improved performance or lower residual [START_REF] Myung | Tutorial on maximum likelihood estimation[END_REF]. First step in this scheme is to linearize b j l,m (T l , T m ) and a h l (S h l ) about current estimate Ŝ assuming very small error ∆S using Taylor expansion to get:

and:

where

and J a h l are Jacobians of b j l,m (T l , T m ) and a h l (S h l ), with respect to S, respectively. Replacing ( 16) and ( 17) in [START_REF] Dou | 3D Scanning Deformable Objects With a Single RGBD Sensor[END_REF] and concatenating b j l,m , a h l , ∆S and corresponding Jacobians we have:

After that we take the derivative of V (S) with respect to S and equate it to zero to get: