cvsnakeimage 的例子-

cvsnakeimage 的例子opencv例子里没有提供cvsnakeimage的使用方法，在此整理一个例子，可以形象的看看snake算法的结果，大致做法是：首先设定域值分割，把基本的轮廓找出来，见图中蓝色轮廓线，再将轮廓点传入cvSnakeImage函数，计算出绿色的snake轮廓线。其中参数alpha代表点相互靠拢的权值(0-1.0)，beta表示弯曲能量（越小越容易弯曲）(0-1.0)，gamma表示整体能量(0

opencv例子里没有提供cvsnakeimage的使用方法，在此整理一个例子，可以形象的看看snake算法的结果，大致做法是：

首先设定域值分割，把基本的轮廓找出来，见图中蓝色轮廓线，再将轮廓点传入cvSnakeImage函数，计算出绿色的snake轮廓线。

其中参数alpha代表点相互靠拢的权值(0-1.0)，beta表示弯曲能量（越小越容易弯曲）(0-1.0)，gamma表示整体能量(0-1.0)。其中参数我自己也不确定具体的范围，最好自己更改不同的范围试试.

//
TrainingTools.cpp:定义控制台应用程序的入口点。

//


#include
“
stdafx.h
“

#include
<
iostream
>

#include
<
string
.h
>

#include
<
cxcore.h
>

#include
<
cv.h
>

#include
<
highgui.h
>

#include
<
fstream
>

IplImage
*
image
=

0
;
//
原始图像


IplImage
*
image2
=

0
;
//
原始图像copy




using

namespace
std;

int
Thresholdness
=

141
;

int
ialpha
=

20
;

int
ibeta
=
20
;

int
igamma
=
20
;

void
onChange(
int
pos)
{



if
(image2)cvReleaseImage(
&
image2);

if
(image)cvReleaseImage(
&
image);

image2
=
cvLoadImage(
“
grey.bmp
“
,
1
);
//
显示图片


image
=
cvLoadImage(
“
grey.bmp
“
,
0
);

cvThreshold(image,image,Thresholdness,
255
,CV_THRESH_BINARY);
//
分割域值



CvMemStorage
*
storage
=
cvCreateMemStorage(
0
);
CvSeq
*
contours
=

0
;

cvFindContours(image,storage,
&
contours,
sizeof
(CvContour),
//
寻找初始化轮廓


CV_RETR_EXTERNAL,CV_CHAIN_APPROX_SIMPLE);

if
(
!
contours)
return
;

int
length
=
contours
->
total;

if
(length
<
10
)
return
;
CvPoint
*
point
=

new
CvPoint[length];
//
分配轮廓点



CvSeqReaderreader;
CvPointpt
=
cvPoint(
0
,
0
);;
CvSeq
*
contour2
=
contours;

cvStartReadSeq(contour2,
&
reader);

for
(
int
i
=

0
;i
<
length;i
++
)
{

CV_READ_SEQ_ELEM(pt,reader);
point[i]
=
pt;
}
cvReleaseMemStorage(
&
storage);

//
显示轮廓曲线



for
(
int
i
=
0
;i
<
length;i
++
)
{


int
j
=
(i
+
1
)
%
length;
cvLine(image2,point[i],point[j],CV_RGB(
0
,
0
,
255
),
1
,
8
,
0
);
}

float
alpha
=
ialpha
/
100.0f
;

float
beta
=
ibeta
/
100.0f
;

float
gamma
=
igamma
/
100.0f
;

CvSizesize;
size.width
=
3
;
size.height
=
3
;
CvTermCriteriacriteria;
criteria.type
=
CV_TERMCRIT_ITER;
criteria.max_iter
=
1000
;
criteria.epsilon
=
0.1
;
cvSnakeImage(image,point,length,
&
alpha,
&
beta,
&
gamma,CV_VALUE,size,criteria,
0
);

//
显示曲线



for
(
int
i
=
0
;i
<
length;i
++
)
{


int
j
=
(i
+
1
)
%
length;
cvLine(image2,point[i],point[j],CV_RGB(
0
,
255
,
0
),
1
,
8
,
0
);
}
delete[]point;

}

int
main(
int
argc,
char
*
argv[])
{

cvNamedWindow(
“
win1
“
,
0
);
cvCreateTrackbar(
“
Thd
“
,
“
win1
“
,
&
Thresholdness,
255
,onChange);
cvCreateTrackbar(
“
alpha
“
,
“
win1
“
,
&
ialpha,
100
,onChange);
cvCreateTrackbar(
“
beta
“
,
“
win1
“
,
&
ibeta,
100
,onChange);
cvCreateTrackbar(
“
gamma
“
,
“
win1
“
,
&
igamma,
100
,onChange);
cvResizeWindow(
“
win1
“
,
300
,
500
);
onChange(
0
);

for
(;;)
{


if
(cvWaitKey(
40
)
==
27
)
break
;
cvShowImage(
“
win1
“
,image2);
}


return

0
;
}