数字图象处理边界区域生长（Matlab）

clear all image=imread('001.png'); figure(1),imshow(image); axis on; title('原始图像'); I=rgb2gray(image); figure(2),imshow(I); axis on; title('原始图像灰度图像'); I=double(uint8(I)); [M,N]=size(I); [y,x]=getpts; %获得区域生长起始点 x1=round(x); %横坐标取整 y1=round(y); %纵坐标取整 seed=I(x1,y1); %将生长起始点灰度值存入seed中 Y=zeros(M,N); %作一个全零与原图像等大的图像矩阵Y，作为输出图像矩阵 Y(x1,y1)=1; %将Y中与所取点相对应位置的点设置为白场 sum=seed; %储存符合区域生长条件的点的灰度值的和 suit=1; %储存符合区域生长条件的点的个数 count=1; %初始化, 记录每次判断一点周围八点符合条件的新点的数 threshold=30; %设置阈值 while count>0 s=0; %记录判断一点周围八点时，符合条件的新点的灰度值之和 count=0; for i=1:M for j=1:N if Y(i,j)==1 if (i-1)>0 & (i+1)<(M+1) & (j-1)>0 & (j+1)<(N+1) %判断此点是否为图像边界上的点 for u= -1:1 %判断点周围八点是否符合域值条件 for v= -1:1 %u,v为偏移量 if Y(i+u,j+v)==0 & abs(I(i+u,j+v)-seed)<=threshold%判断是否未存在于输出矩阵Y，并且为符合域值条件的点 Y(i+u,j+v)=1; %符合以上两条件就将其在Y中与之位置对应的点设置为白场 count=count+1; %满足条件的点的个数加1 s=s+I(i+u,j+v); %此点的灰度之加入s中 end end end end end end end suit=suit+count; %将count加入符合点数计数器中 sum=sum+s; %将s加入符合点的灰度值总合中 seed1=sum/suit; %计算新的灰度平均值 end figure imshow(Y); axis on; title('生成图像'); display('选取的种子点坐标及像素值'); X=x1,Y=y1 seed display('平均灰度值'); seed1 suit