clear?all
rgb1?=?imread（‘Road1.tif‘）;
RGB=rgb2gray（rgb1）;
%對圖像進行邊緣檢測
BW?=?edge（RGB‘sobel‘0.1）;
imshow（BW）;
title（‘Figure1‘）;
%對檢測結(jié)果進行Hough變換
[HTR]?=?hough（BW‘ThetaResolution‘1‘RhoResolution‘10）;
imshow（H[]‘XData‘T‘YData‘R‘InitialMagnification‘‘fit‘）;
xlabel（‘\theta‘）ylabel（‘r‘）;
axis?onaxis?normalhold?on;
%對Hough變換中的密集交點進行檢測
P=?houg