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