function [H,Xp]=hough(i,theta)

[H,Xp]=radon(i,theta);