2012-12-11 21 views
4

我想测试背景减法的ViBe算法。目前我正在使用opencv库。我在opencv/samples/gpu/bgfg_segm.cpp和bgfg_vibe.cpp文件中找到了一个示例实现。这些文件在gpu模块下。现在我有一个没有GPU的系统。当我尝试运行代码时,它在第一帧的初始化时崩溃。任何人都可以告诉我如何解决这个问题?是否可以使用ViBe算法,在opencv中实现,用于没有GPU的系统?

在此先感谢。

回答

7

伪代码吸大的时间!这是非伪/版本。结果:有只尝试还是尽量不要

bgfg_vibe.hpp

#ifndef bgfg_vibe_hpp 
#define bgfg_vibe_hpp 
#include "opencv2/core/core.hpp" 
#include "opencv2/imgproc/imgproc.hpp" 
#include "opencv2/highgui/highgui.hpp" 

struct Model { 
    cv::Mat*** samples; 
    cv::Mat** fgch; 
    cv::Mat* fg; 
}; 

class bgfg_vibe 
{ 
#define rndSize 256 
    unsigned char ri; 
#define rdx ri++ 
public: 
    bgfg_vibe(); 
    int N,R,noMin,phi; 
    void init_model(cv::Mat& firstSample); 
    void setphi(int phi); 
    cv::Mat* fg(cv::Mat& frame); 
private: 
    bool initDone; 
    cv::RNG rnd; 
    Model* model; 
    void init(); 
    void fg1ch(cv::Mat& frame,cv::Mat** samples,cv::Mat* fg); 
    int rndp[rndSize],rndn[rndSize],rnd8[rndSize]; 
}; 

#endif 

bgfg_vibe.cpp

#include "bgfg_vibe.hpp" 

bgfg_vibe::bgfg_vibe():R(20),N(20),noMin(2),phi(0) 
{ 
    initDone=false; 
    rnd=cv::theRNG(); 
    ri=0; 
} 
void bgfg_vibe::init() 
{ 
    for(int i=0;i<rndSize;i++) 
    { 
     rndp[i]=rnd(phi); 
     rndn[i]=rnd(N); 
     rnd8[i]=rnd(8); 
    } 
} 
void bgfg_vibe::setphi(int phi) 
{ 
    this->phi=phi; 
    for(int i=0;i<rndSize;i++) 
    { 
     rndp[i]=rnd(phi); 
    } 
} 
void bgfg_vibe::init_model(cv::Mat& firstSample) 
{ 
    std::vector<cv::Mat> channels; 
    split(firstSample,channels); 
    if(!initDone) 
    { 
     init(); 
     initDone=true; 
    } 
    model=new Model; 
    model->fgch= new cv::Mat*[channels.size()]; 
    model->samples=new cv::Mat**[N]; 
    model->fg=new cv::Mat(cv::Size(firstSample.cols,firstSample.rows), CV_8UC1); 
    for(size_t s=0;s<channels.size();s++) 
    {  
     model->fgch[s]=new cv::Mat(cv::Size(firstSample.cols,firstSample.rows), CV_8UC1); 
     cv::Mat** samples= new cv::Mat*[N]; 
     for(int i=0;i<N;i++) 
     { 
      samples[i]= new cv::Mat(cv::Size(firstSample.cols,firstSample.rows), CV_8UC1); 
     } 
     for(int i=0;i<channels[s].rows;i++) 
     { 
      int ioff=channels[s].step.p[0]*i; 
      for(int j=0;j<channels[0].cols;j++) 
      { 
       for(int k=0;k<1;k++) 
       { 
        (samples[k]->data + ioff)[j]=channels[s].at<uchar>(i,j); 
       } 
       (model->fgch[s]->data + ioff)[j]=0; 

       if(s==0)(model->fg->data + ioff)[j]=0; 
      } 
     } 
     model->samples[s]=samples; 
    } 
} 
void bgfg_vibe::fg1ch(cv::Mat& frame,cv::Mat** samples,cv::Mat* fg) 
{ 
    int step=frame.step.p[0]; 
    for(int i=1;i<frame.rows-1;i++) 
    { 
     int ioff= step*i; 
     for(int j=1;j<frame.cols-1;j++) 
     { 
      int count =0,index=0; 
      while((count<noMin) && (index<N)) 
      { 
       int dist= (samples[index]->data + ioff)[j]-(frame.data + ioff)[j]; 
       if(dist<=R && dist>=-R) 
       { 
        count++; 
       } 
       index++; 
      } 
      if(count>=noMin) 
      { 
       ((fg->data + ioff))[j]=0; 
       int rand= rndp[rdx]; 
       if(rand==0) 
       { 
        rand= rndn[rdx]; 
        (samples[rand]->data + ioff)[j]=(frame.data + ioff)[j]; 
       } 
       rand= rndp[rdx]; 
       int nxoff=ioff; 
       if(rand==0) 
       { 
        int nx=i,ny=j; 
        int cases= rnd8[rdx]; 
        switch(cases) 
        { 
        case 0: 
         //nx--; 
         nxoff=ioff-step; 
         ny--; 
         break; 
        case 1: 
         //nx--; 
         nxoff=ioff-step; 
         ny; 
         break; 
        case 2: 
         //nx--; 
         nxoff=ioff-step; 
         ny++; 
         break; 
        case 3: 
         //nx++; 
         nxoff=ioff+step; 
         ny--; 
         break; 
        case 4: 
         //nx++; 
         nxoff=ioff+step; 
         ny; 
         break; 
        case 5: 
         //nx++; 
         nxoff=ioff+step; 
         ny++; 
         break; 
        case 6: 
         //nx; 
         ny--; 
         break; 
        case 7: 
         //nx; 
         ny++; 
         break; 
        } 
        rand= rndn[rdx]; 
        (samples[rand]->data + nxoff)[ny]=(frame.data + ioff)[j]; 
       } 
      }else 
      { 
       ((fg->data + ioff))[j]=255; 
      } 
     } 
    } 
} 
cv::Mat* bgfg_vibe::fg(cv::Mat& frame) 
{ 
    std::vector<cv::Mat> channels; 
    split(frame,channels); 
    for(size_t i=0;i<channels.size();i++) 
    { 
     fg1ch(channels[i],model->samples[i],model->fgch[i]);   
     if(i>0 && i<2) 
     { 
      bitwise_or(*model->fgch[i-1],*model->fgch[i],*model->fg); 
     } 
     if(i>=2) 
     { 
      bitwise_or(*model->fg,*model->fgch[i],*model->fg); 
     } 
    } 
    if(channels.size()==1) return model->fgch[0]; 
    return model->fg; 
} 

的main.cpp

#include "bgfg_vibe.hpp" 
using namespace cv; 

int main(int argc, char ** argv) 
{ 
    Mat frame; 
    bgfg_vibe bgfg; 
    VideoCapture cap("c:/toprocess/frame_%04d.jpg"); 
    cap >> frame; 
    bgfg.init_model(frame); 
    for(;;) 
    { 
     cap>>frame; 
     Mat fg = *bgfg.fg(frame); 
     imshow("fg",fg); 
     waitKey(1); 
    } 
    return 0; 
} 
+1

感谢您的回复。它已经这么久了。那么,我还没有测试ViBe的上述实现,但它似乎非常有帮助。至少,上面的代码是ViBe的非gpu版本,代表了纸张的伪代码。所以,有兴趣的人可以试试这个。 – sam18

+0

非常感谢您的代码@Zaw Lin。它完美的作品。但是,我注意到了一些我不确定的事情。在bgfg_vibe.cpp中应该是这一行: model-> samples = new cv :: Mat ** [N]; be: model-> samples = new cv :: Mat ** [3]; //用于3声道图像。 – tnq177

+1

嗨,N是正确的。对于多通道图像,每个通道都被视为灰色,实现只是为每个通道初始化单个模型。 –

0

我和你有同样的问题:想在OpenCV中使用没有GPU的Vibe。不幸的是,在当前版本的OpenCV(2.4.3)中没有非gpu实现。如果我错了,请告诉我。所以我基于论文中提供的伪代码实施了自己的氛围。这并不像我想的那么困难。

到目前为止我对此算法感到满意。总的来说这很好,但我认为还需要一些进一步的工作来克服幽灵问题。

+0

感谢您的回复。顺便说一句你的经验与伪代码?我的意思是,你能告诉我一些统计数据或者其他什么,比如分辨率和fps?我已经实现了MOG,对于640 X 480分辨率帧,我将达到30帧/秒,学习参数设置为0.您认为PBAS能够以自己的方式解决鬼影问题吗? – sam18

0

请注意,盛传是专利而且它的使用受许可费的限制。请前往vibeinmotion.com获取试用版和更多信息。我负责ViBE的业务开发,您可以直接联系我获取更多信息。问候。

+0

具有完整SQL异常的网站不存在任何信息。 – Horizon1710

相关问题