function p = vgg_PhogDescriptor(bh,L,bin)

p = [];
%level 0
for b=1:bin   
    p = [p;sum(sum(sum(bh==b)))];
end
        
cella = 1;
for l=1:L
    x = fix(size(bh,2)/(2^l));
    y = fix(size(bh,1)/(2^l));
    xx=0;
    yy=0;
    while xx+x<=size(bh,2)
        while yy +y <=size(bh,1) 
            bh_cella = [];
            
            bh_cella = bh(yy+1:yy+y,xx+1:xx+x,:);
            
            for b=1:bin
                p = [p;sum(sum(sum(bh_cella==b)))];
            end 
            yy = yy+y;
        end        
        cella = cella+1;
        yy = 0;
        xx = xx+x;
    end
end
if sum(p)~=0
    p = p/sum(p);
end