function [H] = my_hough(input_img)
    r = round(sqrt(power(size(input_img,1),2) + power(size(input_img,2),2)));
    H = zeros(r,180);
    for i1 = 1:size(input_img,1)
        for i2 = 1:size(input_img,2)
            if input_img(i1,i2) > 0
                for theta = 1:180
                     r = round(i1*cosd(theta) + i2*sind(theta));
                     if r > 0
                          H(r, theta) = H(r, theta)+1;
                     end
                end
            end
        end
    end    
end