彩色图像区域生长代码

I=imread('C:\Users\XIAOKUNKUN\Desktop\图5.jpg');

imshow(I);

I=double(I);

[M,N,H]=size(I);

[y,x]=getpts; %获得区域生长起始点

x1=round(x); %横坐标取整

y1=round(y); %纵坐标取整

I(x1,y1,1)=255;

I(x1,y1,2)=0;

I(x1,y1,3)=0;

Y=zeros(M,N,H); %作一个全零与原图像等大的图像矩阵Y，作为输出图像矩阵suit=1; %储存符合区域生长条件的点的个数

count=1; %记录每次判断一点周围八点符合条件的新点的数目

while count>0

s=0; %记录判断一点周围八点时，符合条件的新点的灰度值之和

count=0;

for m=1:M

for n=1:N

r=I(m,n,1); %将此像素的R分量赋给r；

g=I(m,n,2); %将此像素的G分量赋给g；

b=I(m,n,3);

X=2-sqrt(r.^2+g.^2+b.^2)/sqrt(13619);

Y=(r*109+g*37+b*13)/sqrt(r^2+g^2+b.^2)/sqrt(13619);

d1=X-0.754;

d2=Y-0.96;

if (I(x1,y1,1)==255&&I(x1,y1,2)==0&&I(x1,y1,3)==0)

if (m-1)>0 && (m+1)<(M+1) && (n-1)>0 && (n+1)<(N+1) %判断此点是否为图像边界上的点

for u= -1:1 %判断点周围八点是否符合域值条件

for v= -1:1 %u,v为偏移量

if (d1>=0&&d2>=0)%判断是否未存在于输出矩阵Y，并且为符合域值条件的点I(m+u,n+v,1)=255;

I(m+u,n+v,2)=0;

I(m+u,n+v,3)=0;

count=count+1;

end

end

end

end

end

end

end

suit=suit+count; %将n加入符合点数计数器中

end

imshow(Y);