2

OpenCVを使用して、C++で画像の緑色のコンポーネントのヒストグラムを作成したいと思います。次のコードはカラー画像では正常に機能していますが、画像をRGBコンポーネントに分割し、緑色のコンポーネントを使用してcalcHist関数を呼び出すと、次のエラーが発生します。

OpenCVエラー:histPrepareImages、ファイル/root/src/OpenCV-2.4.1/modules/imgproc/src/histogram.cppでアサーションが失敗しました(j <nimages)、148行目は「cv::Exception」のインスタンスをスローした後に呼び出されましたwhat():/root/src/OpenCV-2.4.1/modules/imgproc/src/histogram.cpp:148:エラー:(-215)j <関数histPrepareImagesのnimagesが中止されました(コアダンプ)

これが同じコードです。ヒストグラムを作成するために2枚の画像を撮りました。誰もがこの問題を解決するのに役立ちます。

#include <cv.h>
#include <highgui.h>

using namespace cv;

int main( int argc, char** argv )
{
    Mat src,src1, hsv, hsv1;
    if( argc != 3 || !(src=imread(argv[1], 1)).data || !(src=imread(argv[2], 1)).data)
        return -1;

    std::vector<cv::Mat> three_channels;
    cv::split(src,three_channels);

    std::vector<cv::Mat> three_channels1;
    cv::split(src1,three_channels1);


    //cvtColor(src, hsv, CV_BGR2HSV);
    //cvtColor(src1, hsv1, CV_BGR2HSV);

    // Quantize the hue to 30 levels
    // and the saturation to 32 levels
    int hbins = 30, sbins = 32;
    int histSize[] = {hbins, sbins};
    // hue varies from 0 to 179, see cvtColor
    float hranges[] = { 0, 180 };
    // saturation varies from 0 (black-gray-white) to
    // 255 (pure spectrum color)
    float sranges[] = { 0, 256 };
    const float* ranges[] = { hranges, sranges };
    MatND hist, hist1, difference;
    // we compute the histogram from the 0-th and 1-st channels
    int channels[] = {0, 1};

    calcHist( &three_channels[1], 1, channels, Mat(), // do not use mask
            hist, 2, histSize, ranges,
            true, // the histogram is uniform
            false );
    calcHist( &three_channels1[1], 1, channels, Mat(), // do not use mask
            hist1, 2, histSize, ranges,
            true, // the histogram is uniform
            false );
    double maxVal=0;
    minMaxLoc(hist, 0, &maxVal, 0, 0);
    minMaxLoc(hist1, 0, &maxVal, 0, 0);

    int scale = 10;
    Mat histImg = Mat::zeros(sbins*scale, hbins*10, CV_8UC3);
    Mat hist1Img = Mat::zeros(sbins*scale, hbins*10, CV_8UC3);
    Mat hist2Img = Mat::zeros(sbins*scale, hbins*10, CV_8UC3);

    double hist_diff =0;    
    hist_diff = compareHist(hist, hist1, CV_COMP_CORREL);

    absdiff(hist, hist1, difference);

    printf("\nHist Diff:   %f\n", hist_diff);

    for( int h = 0; h < hbins; h++ )
        for( int s = 0; s < sbins; s++ )
        {
            float binVal = hist.at<float>(h, s);
            int intensity = cvRound(binVal*255/maxVal);
            rectangle( histImg, Point(h*scale, s*scale),
                    Point( (h+1)*scale - 1, (s+1)*scale - 1),
                    Scalar::all(intensity),
                    CV_FILLED );
        }
    for( int h = 0; h < hbins; h++ )
        for( int s = 0; s < sbins; s++ )
        {
            float binVal = hist1.at<float>(h, s);
            int intensity = cvRound(binVal*255/maxVal);
            rectangle( hist1Img, Point(h*scale, s*scale),
                    Point( (h+1)*scale - 1, (s+1)*scale - 1),
                    Scalar::all(intensity),
                    CV_FILLED );
        }

    for( int h = 0; h < hbins; h++ )
        for( int s = 0; s < sbins; s++ )
        {
            float binVal = difference.at<float>(h, s);
            int intensity = cvRound(binVal*255/maxVal);
            rectangle( hist2Img, Point(h*scale, s*scale),
                    Point( (h+1)*scale - 1, (s+1)*scale - 1),
                    Scalar::all(intensity),
                    CV_FILLED );
        }


    namedWindow( "Source", 1 );
    imshow( "Source", src );

    namedWindow( "H-S Histogram", 1 );
    imshow( "H-S Histogram", histImg );

    namedWindow( "H-S Histogram1", 1 );
    imshow( "H-S Histogram1", hist1Img );

    namedWindow( "H-S Histogram2", 1 );
    imshow( "H-S Histogram2", hist2Img );
    waitKey();
}
4

1 に答える 1

4

分割したときに、チャネルが1つしかない画像から2つのチャネル(0と1)のヒストグラムを計算しようとしています。

コードの詳細は調べませんでしたが、分割を省略して、three_channels [1] /three_channels1[1]の代わりにsrc/src1をcalcHistに渡すことができると思います。channels = {1}

編集

コードで、に変更channels = {0,1}channels{0}ます。エラーは発生しません。に単一チャネルcalcHist()の画像を渡すので、チャネル0(唯一のチャネル)のみを使用する必要があります。three_channels [1]を入力画像として渡すことにより、入力画像の2番目のチャネルを実際に分析していることを確認できます。

または、次の手順を実行します。

    int channels[] = {1};

    calcHist( &src, 1, channels, Mat(), // do not use mask
        hist, 2, histSize, ranges,
        true, // the histogram is uniform
        false );

もうする必要はありcv::split(src,three_channels)ません。

これはコンパイルする2つのバージョンですが、実際には2Dヒストグラムではなく、緑色(1D)ヒストグラムを計算する必要があります。これがあなたの編集したコードです、それは(うまくいけば)あなたが望むことをします:


int main( int argc, char** argv )
{
    Mat src,src1;
    if( argc != 3 || !(src=imread(argv[1], 1)).data || !(src=imread(argv[2], 1)).data)
            return -1;



    // Quantize the green to 30 levels
    int greenbins = 30;
    int histSize[] = {greenbins};
    // green varies from 0 to 255 (pure spectrum color)
    float greenranges[] = { 0, 256 };
    const float* ranges[] = { greenranges };
    MatND hist, hist1, difference;
    // we compute the histogram from the 2nd channel (green, index is 1)
    int channels[] = {1};

    calcHist( &src, 1, channels, Mat(), // do not use mask
            hist, 1, histSize, ranges,
            true, // the histogram is uniform
            false );
    calcHist( &src1, 1, channels, Mat(), // do not use mask
            hist1, 1, histSize, ranges,
            true, // the histogram is uniform
            false );
    double maxVal1=0;
    double maxVal2 =0;
    minMaxLoc(hist, 0, &maxVal1, 0, 0);
    minMaxLoc(hist1, 0, &maxVal2, 0, 0);
    double maxVal = max(maxVal1, maxVal2);

    int scale = 10;
    int width = 50;
    Mat histImg = Mat::zeros(greenbins*scale, width, CV_8UC3);
    Mat hist1Img = Mat::zeros(greenbins*scale, width, CV_8UC3);
    Mat hist2Img = Mat::zeros(greenbins*scale, width, CV_8UC3);

    double hist_diff =0;    
    hist_diff = compareHist(hist, hist1, CV_COMP_CORREL);

    absdiff(hist, hist1, difference);

    printf("\nHist Diff:   %f\n", hist_diff);

    for( int h = 0; h<greenbins; ++h)
    {
        float binVal = hist.at<float>(h);
        int intensity = cvRound(binVal*255/maxVal);
        rectangle( histImg, Point(0, h*scale),
                    Point(width, (h+1)*scale),
                    Scalar::all(intensity),
                    CV_FILLED );
    }
    for( int h = 0; h<greenbins; ++h)
    {
        float binVal = hist1.at<float>(h);
        int intensity = cvRound(binVal*255/maxVal);
        rectangle( hist1Img, Point(0, h*scale),
                    Point(width, (h+1)*scale),
                    Scalar::all(intensity),
                    CV_FILLED );
    }

    for(int h = 0; h < greenbins; ++h)
    {
        float binVal = difference.at<float>(h);
        int intensity = cvRound(binVal*255/maxVal);
        rectangle( hist2Img, Point(0, h*scale),
                    Point(width, (h+1)*scale),
                    Scalar::all(intensity),
                    CV_FILLED );
    }

    imshow( "Source", src );
    imshow( "Source1", src1 );

    imshow( "src1 green Histogram", histImg );

    imshow( "src2 green Histogram", hist1Img );

    imshow( "diff green Histogram", hist2Img );
    waitKey();
}
于 2012-07-25T08:11:22.630 に答える