DIBR程序代码，基于左右图

/*************************************** * Author : Ley Wang * Date : 2012-10-21 * Description : ***************************************/ #include "dibr.h" void DIBR(IplImage* grayimage,IplImage* Depth,IplImage* NewImg) { IplImage* mask = cvCreateImage(cvGetSize(grayimage),8,1); cvSetZero(mask); // caculate the shift table to look up. int table[256]; int knear = -128; // knear = 0 means everything is displayed behind the screen int kfar = 64; int xb = 6; // eye seperation distance 6cm int D = 300; // view distance 3m int Npix = 320; // standard definition display to reduce the parallax (maybe 720 better!!) for (int i = 0; i < 256; i++) { double A = i * ( knear/64.0 + kfar/16.0 ) / 255.0; double h = - xb * Npix * ( A-kfar/16.0) / D; table[i] = (int)(h/2); } int S = 26; // depth=0, maxmium shift. !!If Npix changes, this value will change!! int step = Depth->widthStep/sizeof(uchar); uchar* datadepth = (uchar*)Depth->imageData; uchar* data = (uchar*)grayimage->imageData; uchar* NewImgdata= (uchar*)NewImg->imageData; uchar* maskdata = (uchar*)mask->imageData; for (int k=0; k<256; k++) { for (int i=0; iheight;i++) { for (int j=Depth->width - 1; j>=0; j--) { int d = datadepth[i*step+j]; if (d==k) { int shift = table[d]; if (j-shift+Swidth-1) { NewImgdata[i*step+j-shift+S] = data[i*step+j]; maskdata[i*step+j-shift+S] = 1; } } } } } // Filling of disocclusions with Background Pixel Extrapolation /*for (int i=0; iheight; i++) { for (int j=0; jwidth; j++) { if (maskdata[i*step+j]==0) { if (j-7<0) NewImgdata[i*step+j] = data[i*step+j]; else NewImgdata[i*step+j] = (NewImgdata[i*step+j-3]+NewImgdata[i*step+j-4] +NewImgdata[i*step+j-5]+NewImgdata[i*step+j-6]) / 4; } } } // Smooth out the filled occlusions for (int i=2; iheight-2; i++) { for (int j=2; jwidth-2; j++) { if (maskdata[i*step+j]==0) { int sum=0; for (int k1=-2;k1<=2;k1++) { for (int k2=-2;k2<=2;k2++) { sum += (int)(NewImgdata[(i+k1)*step+j+k2]); } } NewImgdata[i*step+j] = sum / 25; } } }*/ }