26262
%7B%22isLastPage%22%3Atrue%2C%22notes%22%3A%5B%5D%2C%22pid%22%3A%22t26262%22%2C%22tid%22%3A%2226262%22%2C%22mainForumsId%22%3A%5B%22134%22%5D%2C%22categoriesId%22%3A%5B%220%22%5D%2C%22tcId%22%3A%5B%5D%7D
%7B%22isEditMode%22%3Afalse%7D
向小俊同志学习，发一个基于Mean Shift的运动物体跟踪
当然没有小俊同志高深 。。。 = =
Mean Shift，业内人士称为Mean Sh!t，比较慢，剧烈运动的物体跟踪效果较差 好像是废话 = =
PDF截图两张 敏感部位
![viewer1.png]()
![viewer.png]()
“敏感部位”代码
[size=3]int MeanShift::CalcHisSp(unsigned char* frame, float* histogram){
int pxValue = 0;
for(int i = 0 ; i < HISLENGTH; i++)
histogram[i] = 0;
for(int j = max(0,cY - trackWinHeight / 2);j < min(cY + trackWinHeight / 2,imgH - 1); j ++){
for(int i = max(0,cX - TrackwWidth/ 2); i < min(cX + TrackwWidth/ 2,imgW - 1);i ++){
pixel r = frame[j * imgW * 3 + i * 3] / 16;
pixel g = frame[j * imgW * 3 + i * 3 + 1] / 16;
pixel b = frame[j * imgW * 3 + i * 3 + 2] / 16;
histogram[256 * (int)r + 16 * (int)g + (int)b] += 1;
pxValue ++;
}
}
for(int i = 0; i < HISLENGTH; i ++)
histogram[i] /= pxValue;
else
{
return 0;
}
return pxValue;
}
void MeanShift::MSPSp(unsigned char * frame){
float weights[HISLENGTH];
float newX = 0.0;
float newY = 0.0;
for (int i=0;i<HISLENGTH;i++)
{
if (currHistogram[i] >0.0 )
weights[i] = (float)tempHistogram[i]/(float)currHistogram[i];
else
weights[i] = 0;
}
float sumOfWeights = 0.0;
for(int j = max(0,cY - trackWinHeight / 2);j < min(cY + trackWinHeight / 2,imgH - 1); j ++){
for(int i = max(0,cX - TrackwWidth/ 2); i < min(cX + TrackwWidth/ 2,imgW - 1);i ++){
pixel r = frame[j * imgW * 3 + i * 3] / 16;
pixel g = frame[j * imgW * 3 + i * 3 + 1] / 16;
pixel b = frame[j * imgW * 3 + i * 3 + 2] / 16;
int ptr = (int)(256 * (int)r + 16 * (int)g + (int)b);
newX += (weights[ptr] * (float)i);
newY += (weights[ptr] * (float)j);
sumOfWeights += weights[ptr];
}
}
if(sumOfWeights != 0){
cX = int((newX/sumOfWeights) + 0.5);
cY = int((newY/sumOfWeights) + 0.5);
}
return;
}
void MeanShift::MSTP(unsigned char* frame,int frameNumber)
{
if(frameNumber == 0){
CalcHisSp(frame ,this->tempHistogram);
}
else
{
int stopThreshold = 10;
int iteratorCoumt = 0;
while(iteratorCoumt < stopThreshold){
CalcHisSp(frame ,this->currHistogram);
MSPSp(frame);
iteratorCoumt++;
}
}
}
[/size]
misuzu什么的最喜欢了
![11.jpg]()
Mean Shift，业内人士称为Mean Sh!t，比较慢，剧烈运动的物体跟踪效果较差 好像是废话 = =
PDF截图两张 敏感部位
“敏感部位”代码
[size=3]int MeanShift::CalcHisSp(unsigned char* frame, float* histogram){
int pxValue = 0;
for(int i = 0 ; i < HISLENGTH; i++)
histogram[i] = 0;
for(int j = max(0,cY - trackWinHeight / 2);j < min(cY + trackWinHeight / 2,imgH - 1); j ++){
for(int i = max(0,cX - TrackwWidth/ 2); i < min(cX + TrackwWidth/ 2,imgW - 1);i ++){
pixel r = frame[j * imgW * 3 + i * 3] / 16;
pixel g = frame[j * imgW * 3 + i * 3 + 1] / 16;
pixel b = frame[j * imgW * 3 + i * 3 + 2] / 16;
histogram[256 * (int)r + 16 * (int)g + (int)b] += 1;
pxValue ++;
}
}
for(int i = 0; i < HISLENGTH; i ++)
histogram[i] /= pxValue;
else
{
return 0;
}
return pxValue;
}
void MeanShift::MSPSp(unsigned char * frame){
float weights[HISLENGTH];
float newX = 0.0;
float newY = 0.0;
for (int i=0;i<HISLENGTH;i++)
{
if (currHistogram[i] >0.0 )
weights[i] = (float)tempHistogram[i]/(float)currHistogram[i];
else
weights[i] = 0;
}
float sumOfWeights = 0.0;
for(int j = max(0,cY - trackWinHeight / 2);j < min(cY + trackWinHeight / 2,imgH - 1); j ++){
for(int i = max(0,cX - TrackwWidth/ 2); i < min(cX + TrackwWidth/ 2,imgW - 1);i ++){
pixel r = frame[j * imgW * 3 + i * 3] / 16;
pixel g = frame[j * imgW * 3 + i * 3 + 1] / 16;
pixel b = frame[j * imgW * 3 + i * 3 + 2] / 16;
int ptr = (int)(256 * (int)r + 16 * (int)g + (int)b);
newX += (weights[ptr] * (float)i);
newY += (weights[ptr] * (float)j);
sumOfWeights += weights[ptr];
}
}
if(sumOfWeights != 0){
cX = int((newX/sumOfWeights) + 0.5);
cY = int((newY/sumOfWeights) + 0.5);
}
return;
}
void MeanShift::MSTP(unsigned char* frame,int frameNumber)
{
if(frameNumber == 0){
CalcHisSp(frame ,this->tempHistogram);
}
else
{
int stopThreshold = 10;
int iteratorCoumt = 0;
while(iteratorCoumt < stopThreshold){
CalcHisSp(frame ,this->currHistogram);
MSPSp(frame);
iteratorCoumt++;
}
}
}
[/size]
misuzu什么的最喜欢了
很遗憾PDF是英文的……
话说今天贴数居然能破60
话说今天贴数居然能破60
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
仔细看了遍pdf 似乎一知半解 又似乎什么都没懂。。。
大概还是后者的概率大些吧 =.=
大概还是后者的概率大些吧 =.=
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
引用第5楼ltl于2010-06-19 10:26发表的 :
代码不大美……
敢吐点有意义的槽么= =
晒下限自重
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
IEEE上DOWN的论文？那上面的论文不是一般的贵~~不过非常值的。。含金量非常高。。
CODEPROJECT上也有对象跟踪的文章。
CODEPROJECT上也有对象跟踪的文章。
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
鹰文不好
怎么IEEE那玩意上下载论文要钱么？他们发表论文的时候不是发到期刊上的么？
怎么IEEE那玩意上下载论文要钱么？他们发表论文的时候不是发到期刊上的么？
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
引用第8楼lovehongkong于2010-07-04 13:42发表的 :
鹰文不好
怎么IEEE那玩意上下载论文要钱么？他们发表论文的时候不是发到期刊上的么？
期刊也是要钱的。其实有某些帐号就free了，就是不要写个程序自动狂下就行了。
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
算法的代码通常是比较晦涩的，看不舒服很正常
非常鼓励发表算法的文章，特别是趣味算法。
非常鼓励发表算法的文章，特别是趣味算法。
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
我记得用这个在E6300配2GB双通道内存的机器上还不能做到实时。。
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
在下也想學習一下93大大有什麼資料可以讓我們學習一下?
引用
加载评论中，请稍候...
200字以内，仅用于支线交流，主线讨论请采用回复功能。
93°
学者 笔友
学者 笔友
200字以内，仅用于支线交流，主线讨论请采用回复功能。