function his=hough(p,ib,sinf);
[r,c]=size(p);
his=zeros(1,r);
for io=1:c;
        pos=sinf*(ib+io-1);
        for jo=1:r;
                if (pos>=1) & (pos<=r)
                his(jo)= his(jo)+p(pos,io);
                end;
                pos=pos+1;
        end;
end;

