http://blog.csdn.net/zouxy09/article/details/8210176
很感谢这个博主,从他这真的获益颇多,非常感谢
下面这个网址上有很多关于real-time compressive tracking的内容
http://www4.comp.polyu.edu.hk/~cslzhang/CT/CT.htm
插一句,很早就关注zhang lei了,之前师兄师姐有看过他的文章,他的主页上有很多代码,对我们这种编程能力这么水的人真的是很需要的啊~~~
作者有提供了很完整的代码,只是自己更偏向于matlab,对C++代码不是太能了解,顺便把作者提供的两个用C和matlab混编的求积分图函数也给改成matlab形式了
function ii=integral_matlab(img)
[M,N]=size(img);
s=zeros(M,N);
ii=zeros(M,N);
s(1,:)=img(1,:);
for i=2:M
s(i,:)=img(i,:)+s(i-1,:);
end
ii(:,1)=s(:,1);
for j=2:N
ii(:,j)=s(:,j)+ii(:,j-1);
end
end
function sanmplesFtrVal = FtrVal_matlab(iH,sx,sy,px,py,pw,ph,pwt)
[len_F,len_R]=size(px);
len_S=size(sx,2);
sanmplesFtrVal=zeros(len_F,len_S);
temp=(px>0);
temp=sum(temp,2);
for i=1:len_F
for j=1:len_S
m=0;
x=sx(1,j);
y=sy(1,j);
for k=1:temp(i,1)
minJ=x-1+px(i,k);
maxJ=x-1+px(i,k)+pw(i,k)-1;
minI=y-1+py(i,k);
maxI=y-1+py(i,k)+ph(i,k)-1;
m=m+pwt(i,k)*(iH(minI,minJ)+iH(maxI,maxJ)-iH(minI,maxJ)-iH(maxI,minJ));
end
sanmplesFtrVal(i,j)=m;
end
end