Странный вывод при поиске энтропии кадров видео в opencv

#include <cv.h>
#include <highgui.h>
#include <iostream>
#include <cmath>
#include <cstdlib>
#include <fstream>

using namespace std;

typedef struct histBundle {
double rCh[256];
double gCh[256];
double bCh[256];
}bundleForHist;

bundleForHist getHistFromImage (IplImage* img, int numBins) {
float range[] = {
    0, 
    numBins
};
float *ranges[] = { 
    range 
};

bundleForHist bfh;

CvHistogram *hist = cvCreateHist (1, &numBins, CV_HIST_ARRAY, ranges, 1);
cvClearHist (hist);
IplImage* imgRed   = cvCreateImage(cvGetSize(img), 8, 1);
IplImage* imgGreen = cvCreateImage(cvGetSize(img), 8, 1);
IplImage* imgBlue  = cvCreateImage(cvGetSize(img), 8, 1);
cvSplit(img, imgBlue, imgGreen, imgRed, NULL);
cvCalcHist(&imgRed, hist, 0, 0);
for (int i = 0; i < numBins; ++i) {
    bfh.rCh[i] = cvQueryHistValue_1D(hist, i);
}
cvClearHist(hist);
cvCalcHist(&imgGreen, hist, 0, 0);
for (int i = 0; i < numBins; ++i) {
    bfh.gCh[i] = cvQueryHistValue_1D(hist, i);
}
cvClearHist(hist);
cvCalcHist(&imgBlue, hist, 0, 0);
for (int i = 0; i < numBins; ++i) {
    bfh.bCh[i] = cvQueryHistValue_1D(hist, i);
}
cvClearHist(hist);
return bfh;
}

int main (int argc, char** argv) {
int c;
IplImage* img = NULL;
int frame_number = 0;
CvCapture* capture = cvCaptureFromAVI ("Cricketc1.avi");
assert(capture);
int fps = ( int )cvGetCaptureProperty( capture, CV_CAP_PROP_FPS );
cvNamedWindow ("Video", 0);

while (1) {
    //IplImage * img = cvLoadImage("C:\\Users\\ANIMES~1\\Desktop\\bw.jpg");
    img = cvQueryFrame(capture);
    frame_number++;
    if (img) {
        cvShowImage("Video", img);
        int numBins = 256;
        bundleForHist bfh;
        bfh = getHistFromImage (img, numBins);
        double totalForR = 0;
        double totalForG = 0;
        double totalForB = 0;
        double probR[256];
        double probG[256];
        double probB[256];
        for (int i = 0; i < numBins-1; ++i) {
            totalForR += bfh.rCh[i];
            totalForG += bfh.gCh[i];
            totalForB += bfh.bCh[i];
        }
        double lengthHistogram = totalForR + totalForG + totalForB;
        for (int i = 0; i < 256; ++i) {
            probR[i] = bfh.rCh[i]/(double)lengthHistogram;
            probG[i] = bfh.gCh[i]/(double)lengthHistogram;
            probB[i] = bfh.bCh[i]/(double)lengthHistogram;
            //file << bfh.rCh[i] << "\t" << bfh.gCh[i] << "\t" << bfh.bCh[i] << "\t" << probR[i] << "\t" << probG[i] << "\t" << probB[i] << "\n";
        }

        double entropyR = 0.0;
        double entropyG = 0.0;
        double entropyB = 0.0;
        for (int i = 0; i < numBins; ++i) {
            entropyR += probR[i]*log(probR[i]);
            entropyG += probG[i]*log(probG[i]);
            entropyB += probB[i]*log(probB[i]);
        }
        cout << frame_number << "\t" << (-1.0)*(entropyR + entropyG + entropyB) << endl;
    }
    c = cvWaitKey(1000/fps);
    if (c == 27)
        break;
}
//cvReleaseImage(&img);
cvReleaseCapture(&capture);
cvDestroyWindow("Video");
return 0;
}

ВЫВОД:

.
.
254     -1.#IND
255     -1.#IND
256     -1.#IND
257     -1.#IND
258      5.5686
.
.

Сначала я нашел энтропию изображения, которая оказалась правильной. но почти 80% кадровой энтропии для видео получается -1.#IND.

Это видео ....скачать

Что может пойти не так?


person Animesh Pandey    schedule 17.12.2012    source источник
comment
2 секунды просмотра кода показывают, что вы делаете entropy = blah. Это не то, как рассчитывается энтропия, вам нужно сделать сумму, чтобы вы хотели entropy += blah.   -  person mmgp    schedule 17.12.2012
comment
Ой! Извините... Я сделал это, но все равно получаю результат как -1.#QNAN .....   -  person Animesh Pandey    schedule 17.12.2012
comment
Еще 2 секунды, глядя на него, и я вижу, что вы возвращаете локальную переменную в calcHistProb и обрабатываете ее как указатель позже в своей программе. Я действительно сомневаюсь, что вы скомпилируете это без предупреждений, и вы игнорируете предупреждения. Исправьте все эти языковые проблемы, не игнорируя предупреждения, и отредактируйте свой код.   -  person mmgp    schedule 17.12.2012
comment
@mmgp Я отредактировал код, и все предупреждения были устранены. Но все равно эта ошибка приходит .... !   -  person Animesh Pandey    schedule 18.12.2012
comment
ссылаясь на эту ссылку... Там говорят это незаконная операция ... но я не могу найти ... !!!!   -  person Animesh Pandey    schedule 18.12.2012
comment
Я проверил и обнаружил, что индекс 255 равен нулю для всех массивов с плавающей запятой... это вызывает ошибку.... но я не знаю, как это исправить...   -  person Animesh Pandey    schedule 18.12.2012
comment
Обратите внимание, что вы исправили свой исходный код ответом, но теперь этот вопрос не имеет смысла для тех, кто его посещает. Вам следует подумать о том, чтобы отредактировать свой пост, чтобы на самом деле включить неправильный код, чтобы люди могли понять смысл ответа.   -  person mmgp    schedule 19.12.2012


Ответы (1)


Возможно, prob[i] = 0 для некоторого i, поэтому вы вычисляете log(0), который не определен. Чтобы исправить это, вы просто отбрасываете такие «вероятности»:

for (int i = 0; i < numBins; ++i) {
    if (prob[i])
        entropy += prob[i]*log(prob[i]);
}

Что касается другой ошибки, которую вы обнаружили в отношении значения 0 в ячейке 255, это связано с указанным вами диапазоном. OpenCV рассматривает диапазон для соответствующей функции как [начало, конец), поэтому указание диапазона [0, 255) будет игнорировать конечное значение 255. Вам нужно сохранить как 0, так и 255, поэтому:

float range[] = {0, numBins};
person mmgp    schedule 17.12.2012
comment
Я решил эту проблему с нулевым значением для изображения, я спрашиваю, как тогда может быть ошибка для видео .... ???? - person Animesh Pandey; 18.12.2012
comment
Код не был должным образом исправлен в том, который обрабатывает видео. - person mmgp; 18.12.2012
comment
Смотрите, вы все еще делаете probR[i]*log(probR[i]); и т. д. в своем текущем коде. Вы не проверяете probR[i] == 0 перед применением log. - person mmgp; 19.12.2012
comment
ой! Я сделал это .... забыл опубликовать это в комментариях ... Теперь все в порядке ... ! Спасибо... - person Animesh Pandey; 19.12.2012
comment
Какие критерии я должен указать в цикле while, чтобы после видео оно автоматически вырывалось из цикла? Я не хочу, чтобы какой-либо cvWaitKey() выполнял эту работу! - person Animesh Pandey; 20.12.2012
comment
cvQueryFrame возвращает NULL после последнего кадра. - person mmgp; 20.12.2012
comment
Я использовал это, и это решило мою проблему ... но я заметил, что при использовании while (cvQueryFrame(capture)) количество моих кадров уменьшилось. Видео содержит 1250 кадров. При использовании cvQueryFrame частота кадров видео изменяется случайным образом, а цикл выполняется только для 625 кадров. Такое впечатление, что видео пыталось подстроиться под 625 кадров. Но когда я использую while(1), он показывает все кадры с одинаковой задержкой. Что может быть причиной ?? - person Animesh Pandey; 20.12.2012