I=imread（‘val1.tif‘）;
[mn]=size（I）;I=double（I）;
mark（11）=I（11）;mark（12）=1;
c=1;S=zeros（mn）;S（11）=c;
%c為聯(lián)通區(qū)域的數(shù)目；

%用標(biāo)記法進(jìn)行區(qū)域劃分；
%將當(dāng)前目標(biāo)依次與左前點、左上點、正上點、右上點比較，如果值相似，作相同標(biāo)記
%否則標(biāo)記+1為新的區(qū)域。
for?i=1:m
????for?j=1:n
????????if?i==1&j==1;
????????elseif?i==1&j>=2?%對第一行的元素進(jìn)行區(qū)域劃分
????????????if?abs（I（ij）-I（ij-1））<25
????????????????S（ij）=S（ij-1）;
????????????????a=S（ij）;
????????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
????????????????mark（a2）=mark（a2）+1;
????????????else
????????????????c=c+1;S（ij）=c;a=S（ij）;
????????????????mark（a1）=I（ij）;
????????????????mark（a2）=1;
????????????end
????????elseif?j==1&i>=2?%對第一列的元素進(jìn)行區(qū)域劃分
????????????if?abs（I（ij）-I（i-1j））<25
???????????????S（ij）=S（i-1j）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????elseif?abs（I（ij）-I（i-1j+1））<25
???????????????S（ij）=S（i-1j+1）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;?
????????????else
????????????????c=c+1;S（ij）=c;a=S（ij）;
????????????????mark（a1）=I（ij）;
????????????????mark（a2）=1;
????????????end
????????elseif?i>1&j==n?%對最右邊列進(jìn)行區(qū)域劃分
????????????if?abs（I（ij）-I（ij-1））<25
???????????????S（ij）=S（ij-1）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????elseif?abs（I（ij）-I（i-1j-1））<25
???????????????S（ij）=S（i-1j-1）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????elseif?abs（I（ij）-I（i-1j））<25
???????????????S（ij）=S（i-1j）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????else
????????????????c=c+1;S（ij）=c;a=S（ij）;
????????????????mark（a1）=I（ij）;
????????????????mark（a2）=1;
????????????end
????????else???????%?對中間元素進(jìn)行區(qū)域劃分
????????????if?abs（I（ij）-I（ij-1））<25
???????????????S（ij）=S（ij-1）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????elseif?abs（I（ij）-I（i-1j-1））<25
???????????????S（ij）=S（i-1j-1）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????elseif?abs（I（ij）-I（i-1j））<25
???????????????S（ij）=S（i-1j）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????elseif?abs（I（ij）-I（i-1j+1））<25
???????????????S（ij）=S（i-1j+1）;a=S（ij）;
???????????????mark（a1）=（I（ij）+mark（a1）*mark（a2））/（mark（a2）+1）;
???????????????mark（a2）=mark（a2）+1;
????????????else
????????????????c=c+1;S（ij）=c;a=S（ij）;
????????????????mark（a1）=I（ij）;
????????????????mark（a2）=1;
????????????end
????????end
????end
end

b=1;error=[00];
for?i=2:m??%尋找本屬

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

?????文件??????37680??2010-04-22?08:51??1val1.tif

?????文件???????5016??2010-04-22?14:52??regiongrow.m

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

????????????????42696????????????????????2