压缩跟踪Compressive Tracking之matlab代码

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