function???[picbw_xzangle_xz]=myradon（temparea_gray）
%picbw_xz-校正后的圖片?angle_xz傾斜角度
[heightwidth]=size（temparea_gray）;
level?=?graythresh（temparea_gray）;
picbw=im2bw（temparea_graylevel）;
%figure（1）;imshow（picbw）;
sum_xz=zeros（21height）;
total_xz=zeros（121）;
max_total=0;
for??angle=-10:10
picbw_temp?=?imrotate（picbwangle‘crop‘）;
t=