380 likes | 550 Views
A Local Adaptive Approach for Dens e Stereo Matching i n Architectural Scene Reconstruction. C . Stentoumis 1 , L. Grammatikopoulos 2 , I. Kalisperakis 2 , E. Petsa 2 , G. Karras 1 1. Laboratory of Photogrammetry, Department of Surveying,
E N D
A Local Adaptive Approach for Dense Stereo Matching in Architectural Scene Reconstruction C. Stentoumis1, L. Grammatikopoulos2, I. Kalisperakis2, E. Petsa2, G. Karras1 1. Laboratory of Photogrammetry, Department of Surveying, National Technical University of Athens, GR-15780 Athens, Greece 2. Laboratory of Photogrammetry, Department of Surveying, Technological Educational Institute of Athens, GR-12210 Athens, Greece 5th International Workshop on 3D Virtual Reconstruction and Visualization of Complex Architectures (3D-ARCH '2013), 25-26 February 2013
Outline • Introduction • Related Work • Proposed Algorithm • Experimental Results • Conclusion
Objective • Present a method • Combines pre-existing algorithms and novel considerations • With good sub-pixel accuracy
Related Work Stereo Matching : 5 14 (X,Y) (X-d,Y) Range: 0 - 16
Census on intensity principal derivatives • Census transformation based on gradients: • Less sensitive to radiometric differences and repetitive patterns Intensity Gradient
Census on Gradients Census transform window :
Hamming Distance • Left image • Right image Hamming Distance = 3 XOR
Comparison 2.5% less erroneous pixels • After aggregation step: Default census Census on gradients
Comparison • After aggregation step[13]: [13]Mei X., Sun X., Zhou M., Jiao S., Wang H., Zhang X., 2011. On building an accurate stereo matching system on graphics hardware. ICCV Workshop on GPU in Computer Vision Applications.
Absolute Difference on Image Color and Gradients AD ( color ) : AD ( Gradient ) :
Total Matching Cost • normalized by λ Census (gradient) AD (color) AD (gradient)
Census (gradient) AD (color) AD (gradient) Combined
[25]Zhang K., Lu J., Lafruit G., 2009. Cross-based local stereo matching using orthogonal integral images. IEEE Transactions on Circuits & Systems for Video Technology. Support Region • Cross-based support region[25]: • Threshold of cross-skeleton expansion:
[16]Stentoumis C., Grammatikopoulos L., Kalisperakis I., Karras G., 2012. Implementing an adaptive approach for dense stereo matching. International Archives of Photogrammetry, Remote Sensing & Spatial Information Sciences, Support Region • Threshold of cross-skeleton expansion[16]: • Lmax: largest semi-dimension of the window size • τmax : largest color dissimilarity between p and q 3X3 median filter P q lq
Aggregation Step • A. Normalized by the number of pixels in the support region • B. 3D Gaussian function is applied for smoothing the aggregated costs. • C. Winner-take-all
Comparison • After aggregation step[13]: • Run this step for 4 iterations to get stablecost values. • For iteration 1 and 3, aggregated horizontally and thenvertically. • For iteration 2 and 4, aggregated verticallyand then horizontally.
[13] method Proposed Method
Refinement • Left-right consistency check • Pixel p is characterized as valid (inlier) if the following constraint holds:
Refinement • Outlier cross-based filtering • The cross-based support regions provide a robust description of pixel neighborhoods • The median value of inliers in the support region is selected and attributed to the mismatched pixel.
Refinement • Occlusion / mismatch labeling • Remaining outliers are re-estimated • Mismatches: • The epipolar line of the mismatch pixel intersects with disparity function • Use median interpolation in a small patch around them • Occlusions • Use the second lowest disparity value in the neighborhood
Refinement • EpipolarLine Before After
Refinement • Sub-pixel estimation • Estimation at the sub-pixel level is made by interpolating a 2nd order curve to the cost volume C(d). • This curve is defined by the disparitiesof the preceding and following pixelsand their corresponding cost values • Choose minimum cost position through a closed form solution for the 3 curve points.
Refinement • Disparity map smoothing • Median filter is applied. The effect of overall post-processing refinement
Experimental Results • Evaluated on the Middlebury and EPFL multi-view datasets • Parameter values were kept constant for all tests.
Experimental Results • Middlebury evaluation Error Threshold = 1 Error Threshold = 0.75
[13] method Proposed Method
Experimental Results Threshold = 0.75 , % of wrong pixels
Experimental Results Herz-Jesu-K7 stereo pair
Experimental Results Herz-Jesu-K7 stereo pair