280 likes | 385 Views
Space-time Surface Reconstruction using Incompressible Flow. Andrei Sharf Dan A. Alcantara Thomas Lewiner Chen Greif Alla Sheffer Nina Amenta Daniel Cohen-Or. 4D Data acquisition - setting. Synchronized static cameras 2-16 views Capture rates 10-30 fps.
E N D
Space-time Surface Reconstruction using Incompressible Flow Andrei Sharf Dan A. Alcantara Thomas Lewiner Chen Greif Alla Sheffer Nina Amenta Daniel Cohen-Or
4D Data acquisition - setting • Synchronized static cameras • 2-16 views • Capture rates 10-30 fps
4D Data acquisition - limitations • Persistent self occlusions • Low frame rate and resolution • Noise
time mass field 2D Motivation • Volume is incompressible across time • Explicitly modeling of mass field: • Compute inside/outside volume • Include physical assumptions • Object is watertight manifold
Contribution • Incompressible mass flow 4D reconstruction • Global system considers all frames simultaneously • Simple formulation on a grid • General un-constrained deformations time
Related work Guskov et al. 2003 • Marker based[Marschner et al. 2000; Guskov et al. 2003; White et al. 2007] • Template based[Allen et al. 2002; Anguelov et al. 2004; Zhang et al. 2004; Anguelov et al. 2005, Aguiar et al. 2008; Bradley et al. 2008] • Point correspondence and registration[Shinya 2004; Wang et al. 2005, Anuar and Guskov 2004, Pekelny and Gotsman 2008, Chang and Zwicker 2008; Li et al. 2008] • Surface based space-time[Wand et al. 2007; Mitra et al. 2007; Suessmuth et al. 2008] Zhang et al. 2004 Anuar et al. 2004 Li et al. 2008 Mitra et al. 2007
FLOW reconstruction • In: 3D point cloud frames • Out: watertight surface • Explicit volume modeling • 4D solid on a grid • Characteristic function
3D reconstruction techniques fail • Surface reconstruction of individual frames fails
known outside 0 known inside 1 unknown [0-1] 2D example
1D representation • Domain: space time grid • Material: characteristic function xti values mass amount at each cell • Flow: amount of material vti,j moving from cell xti to cell xt+1j xt+1i-1 xt+1i xt+1i+1 t+1 vti,i vti,i-1 vti,i+1 time xti t
Higher dimensions generalization • Regular 4D grid on top of 3D scan frames • Space-time adjacency relationships: 1D 2D 3D
FLOW physical constraints • Mass preservation: material in cell equals to material flowing into and out of the cell time
FLOW physical constraints • Spatial continuity: values spatially adjacent to be identical everywhere, except across boundaries space
FLOW Physical Constraints • Flow momentum: flow direction should be smooth across time
Constrained Minimization Problem Optimization: Constraints: Incompressibility constraints Boundary values
Challenges Sublinear exponent iterative reweighted least squares Huge matrices fine-tuned iterative solver Mass stability boundary constraints, clamping
Sublinear exponent Iteratively Reweighted Least Squares: from previous iteration small close to discontinuities converges with good init and few outliers time iteration
Huge data problem Problem size 20-200 frames x (28)3 grid resolution x 8 variables per cell in time Reduce initial number of unknowns Pre-assignment from visibility hull : inside/outside labels High resolution per-frame surface reconstruction Sharf et al. 07
Lagrange multipliers: Minimization problem
Matrix engineering Iterative, preconditioned with many eigenvalues 1fast convergence Augmenting approach:solve with CGMINRES solver with decreasing tolerance
Mass stability: clamping/back substitution : amount of mass at cell (i,t) : inside, outside clamp if then adjacent back-substitution reduces the system size time iteration
Results 3D+time • 20 frames at constant resolution • Solver converges in 100 iterations. • Time: 1 minute per-frame • 3.73 GHz CPU, memory requirements up to 4.5GB