%======================================main.m======================================
st = cputime;   %程序開始時的CPU時間 
%程序開始
s1;   %求圖像梯度
s2;   %圖像二值化
I = imread('pipe.jpg');   %重新讀入原始圖像的uint8型矩陣
s3;   %霍夫變換
s4;   %估算直徑
%程序結束
et = cputime - st    %程序運行所花的時間
%====================================main.m(end)====================================