x=imread（‘hough.bmp‘）;
subplot（221）imshow（x）title（‘原圖象‘）;
x=double（x）;
[mn]=size（x）;
???????%???????zerocross檢測圖象???????%
X=edge（x‘zerocross‘）;??
subplot（222）imshow（X）title（‘zerocross檢測圖象‘）;
???????%??????????極坐標參數(shù)域圖象????%
X=double（X）;
M=499;N=456;
H=zeros（MN）;
i1=1;
for?i=1:m
????for?j=1:n
????????if?X（ij）==1
????????????for?t=2*pi/N:2*pi/N:2*pi
????????????????k=i*cos（t）+j*sin（t）;
????????????????H（round（k）+250round（t/2/pi*N））=H（round（k）+250round（t/2/pi*N））+1;
????????????end
????????end
????end
end
c=max（max（H））;
H1=H/c*255;
H1=uint8（H1）;
subplot（223）imshow（H1）title（‘極坐標參數(shù)圖象‘）;
?????%??????直線圖象?????????%
d=0;
H1=double（H1）;
for?i=1:499
????for?j=1:256
????????if?H1（ij）>=120
????????????d=d+1;
????????end
????end
end
B=zeros（2*d1）;
for?i=1:499
????for?j=1:256
????????if?H1（ij）>=120
????????????B（i11）=i;
????????????B（i1+11）=j;
????????????i1=i1+2;
????????end
????end
end
x1=double（x）;
for?l=1:2:d
????k=B（l1）-250;
????sita=B（l+11）*2*pi/256;
????x1=double（x1）;
????for?i=1:m
????????j=round（（k-i*cos（sita））/sin（sita））;
????????if?j>=1?&?j<=n
????????????x1（ij）=255;
????????end
????end
end
x1=uint8（x1）;
subplot（224）imshow（x1）title（‘直線圖象‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????43558??2009-06-06?23:45??實驗六\1.bmp

?????文件??????26464??2009-06-07?23:59??實驗六\hough.bmp

?????文件???????1270??2009-06-08?00:28??實驗六\hough.m

?????文件???????2865??2009-06-07?22:33??實驗六\Laplacian.m

?????文件???????2824??2009-06-07?22:22??實驗六\Sobel.m

?????文件?????155136??2009-06-08?21:52??實驗六\實驗報告.doc

?????文件????????128??2009-06-07?23:21??實驗六\實驗說明.txt

?????目錄??????????0??2009-07-01?20:01??實驗六

-----------?---------??----------?-----??----

???????????????232245????????????????????8