%using?function
%AA=newrag（watnumregionXsizeYsize）;生成鄰接關(guān)系矩陣AA
%[XsizeYsizewat]=remove（XsizeYsizewat）;去掉wat中的邊界值
%BB=updatebb（TTBBnumregion）;用于更新BB矩陣的區(qū)域鄰接關(guān)系
resorce=imread（‘1000.jpg‘）;

hsvimage=rgb2hsv（resorce）;

[XsizeYsizeZsize]=size（resorce）;
figureimshow（resorce）title（‘原圖像‘）;
fr=resorce（::1）;
fg=resorce（::2）;?
fb=resorce（::3）;?
[x1map1]=rgb2ind（resorce8‘dither‘）;
figureimshow（x1map1）;title（‘加抖動‘）;
imwrite（x1map1‘dith.jpg‘‘jpg‘）;
map2=rgb2gray（map1）;
figureimshow（x1map2）title（‘灰度圖‘）;
imwrite（x1map2‘gray.jpg‘‘jpg‘）;
I=medfilt2（x1）;
%se=strel（‘disk‘1）;
%I=imerode（x1se）;
figureimshow（Imap2）title（‘中值濾波后‘）;
imwrite（Imap2‘fushi.jpg‘‘jpg‘）;
%gs=fspecial（‘gaussian‘）;
%I=imfilter（x1gs）;
%I=imfilter（Igs）;
hy?=?fspecial（‘sobel‘）;
hx?=?hy‘;
Iy?=?imfilter（double（I）?hy?‘replicate‘）;
Ix?=?imfilter（double（I）?hx?‘replicate‘）;
gradmag?=?sqrt（Ix.^2?+?Iy.^2）;
wat=watershed（gradmag）;
%wat=watershed（map2）;

figureimshow（wat）title（‘水線‘）;
imwrite（wat‘watershedline.bmp‘‘bmp‘）;
w=imread（‘watershedline.bmp‘）;
seg=imcomplement（w）;
resultfr=imadd（segfr）;
resultfg=imadd（segfg）;
resultfb=imadd（segfb）;
result=cat（3resultfrresultfgresultfb）;
figureimshow（result）title（‘水線分割結(jié)果‘）;
imwrite（result‘watershedConbine.jpg‘‘jpg‘）;
%水線分割完畢
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%computing?RAG?2005-5
%AA=newrag（watnumregionXsizeYsize）;
%1代表鄰接，0代表不鄰接
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
numregion?=?double（max（wat（:）））;
AA=newrag（watnumregionXsizeYsize）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%將原圖從RGB空間轉(zhuǎn)換到HSV空間
hsvimage=rgb2hsv（resorce）;
H=hsvimage（::1）;
S=hsvimage（::2）;
V=hsvimage（::3）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%計(jì)算水線分割后每個區(qū)域中象素的Hue值的和并存儲在數(shù)組HH中數(shù)組tmp中存儲的是每個區(qū)域包含的象素個數(shù).
%計(jì)算水線分割后每個區(qū)域中象素的Saturation值的和并存儲在數(shù)組SS中
%數(shù)組Hmean中存儲的是每個區(qū)域的Hue均值數(shù)組Smean中存儲的是每個區(qū)域的Saturation均值
HH?=?zeros（1numregion）;
SS?=?zeros（1numregion）;
tmp?=?zeros（1numregion）;
Hmean?=?zeros（1numregion）;
Smean?=?zeros（1numregion）;
????for?i=1:Xsize
????????for?j=1:Ysize
?????????????????k?=?wat（ij）;?
?????????????????if?（k?~=?0）
?????????????????????if?H（ij）==0
?????????????????????????HH（k）=?HH（k）+V（ij）;
?????????????????????else
?????????????????????????HH（k）=?HH（k）+H（ij）;?
?????????????????????end;
?????????????????????if?S（ij）==0
?????????????????????????SS（k）=?SS（k）+V（ij）;
?????????????????????else
????????????????????????SS（k）=?SS（k）+S（ij）;
????????????????????end;
????????????????????tmp（k）?=?tmp（k）?+?1;
????????????????end;????
????????end;
????end;
sumtmp=0;
????for?i=1:numregion
????????Hmean（i）=HH（i）/tmp（i）;
????????Smean（i）=SS（i）/tmp（i）;
????????sumtmp=sumtmp+tmp（i）;
????end;
Hmeanback=Hmean;?
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%矩陣BB存儲水線分割后每個區(qū)域的鄰接矩陣
%其中BB（i1）為每個區(qū)域的鄰接區(qū)域個數(shù)
%BB（inumregion+2）