hough变换检测直线，效果较好

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