基于半全局和全局算法立体匹配与研究.docVIP

基于半全局和全局算法的立体匹配研究 摘要： 传统的基于像素点的匹配算法常常是算出初始匹配代价后直接采用贪心策略求取视差，虽然速度较快，但往往是局部最优的，以至精确度很低。针对这一问题，目前策略主要有：（1）半全局优化算法：扫描线算法和动态规划算法；（2）全局优化算法：置信度算法和图割算法。本文旨在通过详细讨论这四种算法原理本质，算法步骤与算法运行，从而深刻分析各自的优点与缺点，为进一步改进其不足，进而研究新的算法打下基础。 关键词：半全局优化，全局优化，扫描线，动态规划，置信度，图割 立体匹配介绍 图像的立体匹配即给定同一场景的两幅图像,寻找同一场景点投影到图像中的像素之间的对应关系。根据考虑的是基于像素点的还是基于区域块，可以分为基于像素点的匹配与基于区域的匹配。立体匹配算法通常是通过构建能量函数试图获得图像的某些全局性质,即全局能量最小化，但通常很难获得能量函数的全局最小化，鉴于此,很多学者更倾向于寻找局部小的求解.然而在一般情形下,局部小不能带来任何的全局性，所以匹配效果较差，准确率较低，基于像素点的匹配就是一种局部小的解，所以若想提高精度，研究的多是一种半全局或全局优化策略的区域匹配算法。 立体匹配的通常包括以下四步： 1 图像预处理Preprocessing —由于拍摄照片的时候难免会有传感器的噪声sensor noise）和光度的扭曲（photometric distortions）而这都会对视差的计算带来严重影响，常用的解决方法有，高斯拉普拉斯滤波 Laplacian of Gaussian LoG filtering [1]直方图均衡化Histogram Equalization/Matching ，中值滤波 Subtraction of mean values computed in the neighbours of each pixel [2]双边滤波Bilateral filtering [3]。 2 匹配代价计算Matching cost computation）—对匹配代价的计算通常有四种方法AD（1-1）、SAD（1-2）、SD（1-3）与SSD（1-4），计算公式，从而能得到元素的不同视差匹配代价所组成的初始视差空间。 1-1 1-2 1-3 1-4 3 视差的计算Disparity computation —真实 4）视差的优化 disparity refinement.）sub-pixel disparity Estimate ，即将原图像进行水平拉伸，再对行像素点进行模糊。 本文将详细研究的就是立体匹配的视差计算阶段，就是在初步计算出来的视差空间中进行半全局或全局优化，从而求出更好的视差值。 与基于像素点的匹配方法不同，基于半全局或全局的算法通常是将匹配问题转换为一个能量方程，然后通过求解该能量方程的最小值来求取视差值。能量方程通常具有以下的形式 1-5 其中Cdata dx 是数据项用来约束像素点在偏移前后的变化尽量小,V dx,dx-1 是光滑项，约束像素点在偏移前后与周围像素点的关系变化尽量小。不同的优化算法的数据项和光滑项的定义不同，本文采取的数据项计算方法为AD。而光滑项则根据不同的算法用不同的模型方法。 二．匹配算法与实现 2.1扫面线优化方法 SO Scanline Optimization 方法属于一种半全局能量最小化优化算法，比基于像素点匹配的算法有更好的适应性，通过最小化自身的能量方程来求得像素对应的视差，但是由于其只考虑一行的最优，并未考虑所有的像素点，所以本文称之为半全局能量最小化算法。 2.1.1算法原理： SO的能量方程的光滑项通常为： 2-1 2-2 其中opt_smooth表示图片中像素点x的光滑值，opt_grad_penalt用来对同纹理区域 灰度梯度值小于给定的阈值opt_grad_penalty 的区域进行惩罚的系数，从而不会在同纹理区域由于匹配代价相同而导致视差计算错误。SO是逐行考虑的算法，本质就是计算从每一行第一个像素点到当前像素点x在视差为dx时候的匹配代价，用sumcost x,dx 表示， 2-3 其中 2-4 而Cost dx 表示从前一个像素点视差为dx-1到当前像素点在视差为dx 时候的匹配代价。 2-5 当保存到本行最后一列元素x W时候，此时sumcost保存了x W处的所有视差对应的最小匹配代价，选取最小的sumcost，进行回溯，此时即得到本行最优视差路径。然后同理计算下一行元素。从而能得出所有像素点的视差值，而SO的最小化能量方程转化为： 2-6 2.1.2算法步骤： 初始化第一列元素 x 1 : 令sumcost x ; 1,2,3,…,dmax Cdata dx 1,2,3,…,dmax