程序代碼： (代碼標記  [code]...[/code]  )

%%%%%%%%%%%%%%%%%%%%%%%%定義祛除背景的函數%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function chayingfa = chayingfa(i);
figure;
imshow(i);
title('原圖顯示')
I=rgb2gray(i);%灰度轉化
figure;
imshow(I);
title('灰度圖顯示')
bw1=im2bw(I);
imshow(bw1);
title('原圖的二值化圖片')
s=strel('disk',5);%根據背景的復雜程度修改圓盤半徑
bg=imopen(I,s);
figure,imshow(bg);%對‘原始圖象’進行開操作得到圖象背景
title('背景圖顯示')
bw=imsubtract(I,bg);
figure,imshow(bw);%用原始圖象與背景圖象作減法
title('使用差影法后的圖片')
level=0.3;
bw2=im2bw(bw,level);
figure,imshow(bw2);%通過圖象的直方圖取得閾值，將圖象二值化
title('差影法后的二值化圖片')