function?[BVDV]=BVandDV（ImgTh）
%?Img=imread（‘LENA.bmp‘）;
%?Th=117;
%?
%?%rgb轉灰度
%?if?isrgb（Img）==1
%?????I_gray=rgb2gray（Img）;
%?else
%?????I_gray=Img;
%?end
%?Img=I_gray;

Itemp=Menlarge（Img33）;%擴大周圍
[mn]=size（Itemp）;
Bsum=0;
Bcount=0;
Dsum=0;
Dcount=0;
if?size（Th1）==1
????for?i=1:m
????????for?j=1:n
????????????if?Itemp（ij）>=Th
????????????????Rr（ij）=255;%前景
????????????else
????????????????Rr（ij）=0;%背景
????????????end
????????end
????end
????for?i=2:m-1
????????for?j=2:n-1
????????????if?Rr（ij）==0%背景
????????????????Bsum=Bsum+（（Itemp（ij）-Itemp（i-1j））^2+（Itemp（ij）-Itemp（i+1j））^2+（Itemp（ij）-Itemp（ij-1））^2+（Itemp（ij）-Itemp（ij+1））^2）/4;
????????????????Bcount=Bcount+1;
????????????else
????????????????Dsum=Dsum+（（Itemp（ij）-Itemp（i-1j））^2+（Itemp（ij）-Itemp（i+1j））^2+（Itemp（ij）-Itemp（ij-1））^2+（Itemp（ij）-Itemp（ij+1））^2）/4;
????????????????Dcount=Dcount+1;
????????????end
????????end
????end
end
if?size（Th1）==m-2
????for?i=2:m-1
????????for?j=2:n-1
????????????if?Th（i-1j-1）==0%背景
????????????????Bsum=Bsum+（（Itemp（ij）-Itemp（i-1j））^2+（Itemp（ij）-Itemp（i+1j））^2+（Itemp（ij）-Itemp（ij-1））^2+（Itemp（ij）-Itemp（ij+1））^2）/4;
????????????????Bcount=Bcount+1;
????????????else
????????????????Dsum=Dsum+（（Itemp（ij）-Itemp（i-1j））^2+（Itemp（ij）-Itemp（i+1j））^2+（Itemp（ij）-Itemp（ij-1））^2+（Itemp（ij）-Itemp（ij+1））^2）/4;
????????????????Dcount=Dcount+1;
????????????end
????????end
????end
end
BV=Bsum/Bcount;
DV=Dsum/Dcount;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1525??2013-09-12?20:37??Bian_Yuan_Zeng_Qiang\BVandDV.m

?????文件?????786486??2000-08-28?12:02??Bian_Yuan_Zeng_Qiang\LENA.BMP

?????文件???????1595??2013-09-12?20:36??Bian_Yuan_Zeng_Qiang\LOG_enhancement.m

?????文件????????628??2013-09-12?19:21??Bian_Yuan_Zeng_Qiang\Menlarge.m

?????文件????????808??2013-09-12?20:06??Bian_Yuan_Zeng_Qiang\Otsu.m

?????目錄??????????0??2013-09-12?20:53??Bian_Yuan_Zeng_Qiang

-----------?---------??----------?-----??----

???????????????791042????????????????????6