function lines = hough(img, Ntheta, Nrho, thr)

% hough
%
% Executa a transformada de Hough com backmapping para linhas
% retas, usando a parametrização rho = x * cos(theta) + y * sin(theta).
%
% Recebe como entrada uma imagem, as dimensões da matriz de
% acumuladores (Ntheta x Nrho) e o limiar thr para escolha dos
% picos na matriz de acumuladores. Devolve a uma imagem binária
% contendo as linhas retas detectadas.

Acc = backhough(img, Ntheta, Nrho);
lines = reverse(Acc, size(img), thr);

