1

http://smsoftdev-solutions.blogspot.com/2009/08/integral-histogram-for-fast-calculation.htmlから抽出した、このコードを使用しています。

現在、そのページのコードをopencv2 c ++実装にインポートしています。これは、これまでのところです。

主要

using namespace cv;
using namespace std;

//CONSTANTES
#define PI 3.142

//FUNCIONES
vector<Mat> extractIntegralHist(Mat in, const int _nbins);

int main(int argc, char *argv[])
{
//Definicion de variables
    Mat frame;
    vector<Mat> integrals(9);
    namedWindow("main_video",WINDOW_AUTOSIZE);
    namedWindow("temporal_window",WINDOW_AUTOSIZE);

//Inicializa la captura del video
    VideoCapture cap(0);
       if(!cap.isOpened()){
           return -1;
           cout << "No se puede capturar";
       }

//Loop principal
    while(1){
       cap >> frame;
       imshow("main_video",frame);
       waitKey(1);

       extractIntegralHist(frame, 9);

    }
    integrals.clear();
}

働き

//CALCULAR HISTOGRAMAS DE INTEGRALES
vector<Mat> extractIntegralHist(Mat in, const int _nbins){
    Mat in_gray;
    int x,y;
    float temp_gradient,temp_magnitude;
    namedWindow("temporal_window",WINDOW_AUTOSIZE);
//Convierte la imagen a escala de grises y normaliza
    cvtColor(in,in_gray,CV_RGB2GRAY);
    equalizeHist(in_gray,in_gray);
//Calcula la derivada de la imagen en gris en direccion x,y utilizando sobel
    Mat xsobel, ysobel;
    Sobel(in_gray, xsobel, CV_32FC1, 1, 0, 1);
    Sobel(in_gray, ysobel, CV_32FC1, 0, 1, 1);
    in_gray.release();
//Crea la matriz de 9 imagenes contenedoras -- REVISAR ESTE CODIGO
    vector<Mat> bins(_nbins);
    for (int i = 0; i < _nbins ; i++) {
        bins[9] = Mat::zeros(in.size(),CV_32F);
    }
//Crea la matriz donde se guardaran las integrales de la imagen
    vector<Mat> integrals(_nbins);
    for (int i = 0; i < _nbins ; i++) {
        integrals[9] = Mat(Size(in_gray.cols+1,in_gray.rows+1),CV_64F);
    }
//Calcula los contenedores
    for(y=0;y<in.rows;y++){

        float* xSobelPtr = (float*) (xsobel.row(y).data);
        float* ySobelPtr = (float*) (ysobel.row(y).data);

        float** binsRowPtrs = new float *[_nbins];
        for (int i = 0; i < _nbins; i++){
            binsRowPtrs[i] = (float*) (bins[i].row(y).data);
        }
    //Para cada pixel en la fila los valores de magnitud y orientacion son calculados
        for(x=0;x<in.cols;x++){
            if(xSobelPtr==0){
                temp_gradient = ((atan(ySobelPtr[x]/(xSobelPtr[x]+0.00001)))*(180/PI)+90);
            }
            else{
                temp_gradient = ((atan(ySobelPtr[x]/xSobelPtr[x]))+(180/PI)+90);
            }
            temp_magnitude = sqrt((xSobelPtr[x] * xSobelPtr[x]) + (ySobelPtr[x] *       ySobelPtr[x]));
            float binStep = 180/ /*_nbins*/9;
                for (int i=1 ; i<=/*_nbins*/9 ; i++){
                    if (temp_gradient <= binStep*i){
                    binsRowPtrs[i-1][x] = temp_magnitude;
                    break;
                    }
                }
        }
     }

     for (int i = 0; i<_nbins;i++){
        integral(bins[i],integrals[i]);
     }

     for (int i = 0; i <_nbins ; i++){
        bins[i].release();
     }

    return integrals;
}

プログラムはコンパイルされますが、実行すると停止し、デバッガーはSIGSEVエラーを表示し、61行目と38行目を表示します。

38-extractIntegralHist（frame、9）;

61-bins [9] = Mat :: zeros（in.size（）、CV_32F）;

すべてがうまくいくように思われるため、問題を特定できませんでした。コードを確認してください。

編集。アップデート。関数extractIntegralHistが実行されると、プログラムが停止しているように見えます。これは、（関数の）行61にあるとされています。これがメモリ処理の例外であることはわかっていますが、何が問題なのかがわかりません。さらに、デバッガーが何を言っているのかわかりません。

ありがとう。

4

1 に答える 1

0

この線

bins[9] = Mat::zeros(in.size(),CV_32F);

の 10 番目の (インデックスがゼロであることに注意してください) 要素にアクセスしようとしていますbins。_nbinsは の値で渡されるため9、これはベクトルの境界を超えています。また、 を使用.at()して要素にアクセスする代わりに を使用する[]と、この種のエラーは例外をスローしout_of_range、問題の診断により役立つはずです。

ベクトルの各要素cv::Matをゼロを含むに割り当てるには、ループを次のように変更することを検討してください。

for (int i = 0; i < _nbins ; i++) {
    bins[i] = Mat::zeros(in.size(),CV_32F);
}

binsまたは、より簡潔に言えば、ループを削除して、初期化を次のように置き換えます

vector<Mat> bins(_nbins, Mat::zeros(in.size(),CV_32F));

の初期化でも同じ問題があることに注意してくださいintegrals。

于 2013-03-14T18:26:44.477 に答える