%直线上两点为(X1,Y1),(X2,Y2)
%直线方程为y=kx+A
%(m,n)为极坐标
function [m,n]=myhough(Point1,Point2)
if Point1(1)==Point2(1)
n=0;
m=Point1(1);
%plot(Point1,Point2);
else
K=(Point2(2)-Point1(2))/(Point2(1)-Point1(1));
A=Point1(2)-K*Point1(1);
n=atan(K);%弧度
n=n*180/pi;%角度
m=abs(A)/sqrt( K*K+1);
end
plot(Point1,Point2);
end