le?=?imread（‘wula.jpg‘）;
le?=?rgb2gray（le）;
le?=?im2double（le）;
%單一的閾值分割，展示不同閾值時的情形
subplot（131）;imshow（le[]）;title（‘原圖‘）;
bw1?=?im2bw（le0.2）;
bw2?=?im2bw（le0.9）;
subplot（132）;imshow（bw1[]）;title（‘閾值為0.2‘）;
subplot（133）;imshow（bw2[]）;title（‘閾值為0.9‘）;
%畫出直方圖，由直方圖判斷閾值