%本程序運用hough變換檢測出原圖像中最長的直線
%由于采用多次for循環，程序運行速度較慢
clear?all;
yuantu=imread（‘yan.jpg‘‘jpg‘）;
yuantu1=rgb2gray（yuantu）;%二值化
figureimshow（yuantu1）;
title（‘原圖像‘）;
pn=imnoise（yuantu1‘salt?&?pepper‘0.02）;
pt=imnoise（pn‘gaussian‘0.02）;
figureimshow（pt）;
title（‘加入椒鹽噪聲和高斯噪聲‘）;
[mn]=size（pt）;
p_max=round（sqrt（m^2+n^2））;
AC_matrix=zeros（2*p_max180）;%直線標準式共線點集生成正弦曲線
for?i=1:m
??for?j=1:n
????if?yuantu1（ij）==0
????????????x=j-1;
????????????y=m-i;
???????for?ii=1:2*p_max
????????for?a=1:180
??????????p=ii-p_max;
????????if??abs（p-x*cos（a*pi/180）-y*sin（a*pi/180））<1;
??????????