opencv c++についての質問です

質問

int main (int argc, char **argv){
int i, c, counter;
int INIT_TIME = 100;
int w = 0, h = 0;

char *FileName = "動画名";
if((capture = cvCaptureFromFile(FileName)) == NULL){
printf("指定された動画がみつかりません\n");
return -1;
}

CvCapture *capture = 0;
IplImage *frame;
IplImage *av_img, *sgm_img;
IplImage *lower_img, *upper_img, *tmp_img;
IplImage *dst_img, *msk_img;
CvFont font;

capture = cvCreateFileCapture (FileName);

// (2)1フレームキャプチャし，キャプチャサイズを取得する．
frame = cvQueryFrame (capture);
w = frame->width;
h = frame->height;

// (3)作業用の領域を生成する
av_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_32F, 3);
sgm_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_32F, 3);
tmp_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_32F, 3);
lower_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_32F, 3);
upper_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_32F, 3);
dst_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_8U, 3);
msk_img = cvCreateImage (cvSize (w, h), IPL_DEPTH_8U, 1);

// (4)背景の輝度平均の初期値を計算する
printf ("Background statistics initialization start\n");
cvSetZero (av_img);
for (i = 0; i < INIT_TIME; i++) {
frame = cvQueryFrame (capture);
cvAcc (frame, av_img);
}
cvConvertScale (av_img, av_img, 1.0 / INIT_TIME);

// (5)背景の輝度振幅の初期値を計算する
cvSetZero (sgm_img);
for (i = 0; i < INIT_TIME; i++) {
frame = cvQueryFrame (capture);
cvConvert (frame, tmp_img);
cvSub (tmp_img, av_img, tmp_img);
cvPow (tmp_img, tmp_img, 2.0);
cvConvertScale (tmp_img, tmp_img, 2.0);
cvPow (tmp_img, tmp_img, 0.5);
cvAcc (tmp_img, sgm_img);
}
cvConvertScale (sgm_img, sgm_img, 1.0 / INIT_TIME);
printf ("Background statistics initialization finish\n");

// (6)表示用ウィンドウを生成する
cvInitFont (&font, CV_FONT_HERSHEY_COMPLEX, 0.7, 0.7);
cvNamedWindow ("Input", CV_WINDOW_AUTOSIZE);
cvNamedWindow ("Substraction", CV_WINDOW_AUTOSIZE);

// (7)取得画像から背景を分離するループ
counter = 0;
while (1) {
frame = cvQueryFrame (capture);
cvConvert (frame, tmp_img);

// (8)背景となりうる画素の輝度値の範囲をチェックする
cvSub (av_img, sgm_img, lower_img);
cvSubS (lower_img, cvScalarAll (Zeta), lower_img);
cvAdd (av_img, sgm_img, upper_img);
cvAddS (upper_img, cvScalarAll (Zeta), upper_img);
cvInRange (tmp_img, lower_img, upper_img, msk_img);

// (9)輝度振幅を再計算する
cvSub (tmp_img, av_img, tmp_img);
cvPow (tmp_img, tmp_img, 2.0);
cvConvertScale (tmp_img, tmp_img, 2.0);
cvPow (tmp_img, tmp_img, 0.5);

// (10)背景と判断された領域の背景の輝度平均と輝度振幅を更新する
cvRunningAvg (frame, av_img, B_PARAM, msk_img);
cvRunningAvg (tmp_img, sgm_img, B_PARAM, msk_img);

// (11)物体領域と判断された領域では輝度振幅のみを（背景領域よりも遅い速度で）更新する
cvNot (msk_img, msk_img);
cvRunningAvg (tmp_img, sgm_img, T_PARAM, msk_img);

// (12)物体領域のみを出力画像にコピーする（背景領域は黒）
cvSetZero (dst_img);
cvCopy (frame, dst_img, msk_img);

// メディアンフィルタでノイズを除去する
cvSmooth( dst_img, dst_img, CV_MEDIAN );

cvMoments( frame, &moment, 0 );
m_00 = cvGetSpatialMoment( &moment, 0, 0 );
m_10 = cvGetSpatialMoment( &moment, 1, 0 );
m_01 = cvGetSpatialMoment( &moment, 0, 1 );
gravityX = m_10 / m_00;
gravityY = m_01 / m_00;

if (frame->nChannels == 3 && cvGetImageCOI (frame) == 0)
cvSetImageCOI (frame, 1);

cvSetImageCOI (frame, 0);
cvCircle( dst_img, cvPoint( gravityX, gravityY ), CIRCLE_RADIUS,CV_RGB( 255, 0, 0 ), LINE_THICKNESS, LINE_TYPE, 0 );

// (13)処理結果を表示する
cvShowImage ("Input", frame);
cvShowImage ("Substraction", dst_img);
counter++;
c = cvWaitKey (10);
if (c == '\x1b')
break;
}

cvDestroyWindow ("Input");
cvDestroyWindow ("Substraction");

cvReleaseImage (&frame);
cvReleaseImage (&dst_img);
cvReleaseImage (&av_img);
cvReleaseImage (&sgm_img);
cvReleaseImage (&lower_img);
cvReleaseImage (&upper_img);
cvReleaseImage (&tmp_img);
cvReleaseImage (&msk_img);
return 0;
}
動画を指定して背景差分を行い、動画内の物体の重心を求めて円を描きたいのですがcvMoments()とcvCircle()関数を入れてコンパイルしたところウィンドウが出て動画が始まる前に終了してしまいます。
どのようにすれば正常に動くようになるのでしょうか？

回答No.1

kabe64

回答が付かないようなので、あまり詳しくはないですが。

cvMoments()で落ちてます。
cvMomentsは1チャンネルかCOIを持つ3チャンネル画像でしか使えないようなので、

if (frame->nChannels == 3 && cvGetImageCOI (frame) == 0)
cvSetImageCOI (frame, 1);

の部分をcvMoments()の前に移動すればとりあえず動くようになりますよ。
これで期待してる動作なのかは知りませんが。