close?all;clear;clc;

I=load（‘changjing.txt‘）;
figure;imshow（I[]）;

BW=edge（I‘sobel‘）;?
figure;imshow（BW）;title（‘sobel后‘）;?

%hough進行hough變換，H為hough矩陣，T是變換角度，R是變換半徑。
[HTR]=hough（BW）;
figure;imshow（H[]‘XData‘T‘YData‘R‘InitialMagnification‘‘fit‘）;title（‘houghpeak‘）;
xlabel（‘\theta‘）?ylabel（‘\rho‘）;
axis?on
axis?normal
hold?on
P=houghpeaks（H8）;?%計算極值
x=T（P（:2））;?%橫坐標?角度
y=R（P（:1））;?%縱坐標?半徑
plot（x