????clc
????close?all
????clear?all;???

????%radon變換
????t=strcat（‘Stest‘int2str（15）‘.jpg‘）;%讀入圖片名稱
????I=imread（t‘jpg‘）;
????I=imresize（I[128128]）;???%重置圖片大小
????I=mat2gray（I）;
????p1=figure（1）;
????theta=0:179;???????????????%radon變換。獲得從0到179每個幅角的?
????[Rx]=radon（Itheta）;
????[MN]=size（R）;
????imagesc（thetaxR）;????????%顯示radon變換結果
????colormap（hot）;?
????title（‘Radon?result‘）;?
????xlabel（‘theta（幅角）‘）;?
????ylabel（‘x（幅值）‘）;?
????colorbar
????
????p2=figure（2）;
????imshow（I）;
????%findpeaks函數只能對一維向量進行處理，因此對radon變換結果變形
????J=reshape（RM*N1）;????????????????????
????[PKSLOCS]=?findpeaks（J‘THRESHOLD‘6）;??%利用峰值檢測函數進行峰值檢測，可設置閾值閾值因圖而異
????L=size（?LOCS1）;?????%檢測到的峰值個數，即最終檢測出的直線條數
????d=zeros（size（LOCS））;?%存放檢測結果-距離
????thetar=zeros（size（LOCS））;%存放檢測的直線角度
?????for?i=1:L
?????????d（i）=mod（LOCS（i）185）-（M-1）/2;???????%實際是距離
?????????thetar（i）=?fix（LOCS（i）/M）;???????????%實際是度數
?????end
?????thetar
?????d
?????%將直線進行呈現，原理參考https://cn.mathworks.com/help/images/detect-lines-using-the-radon-transform.html?searchHighlight=radon%20transform%20peak&s_tid=doc_srchtitle
????for?i=1:L
?????????if?thetar（i）==0
?????????????x1=64+d（i）;
?????????????x2=x1;
???????