imageregiongrowing.m

希望大家给予支持,可以运行的呀,有空的话来学习学习

% 区域增长
clear;
clc;
bw=imread('source.bmp');
imshow(bw);
bw=double(bw);

%计算图像的大小
[ImageHeight ImageWidth]=size(bw);

%初始化标记
for i=1:ImageHeight
    for j=1:ImageWidth
        UnRegion(i,j)=0;
    end
end

%遍历图像的4邻域时使用的数组
nDx=[-1 0 1 0];
nDy=[0 1 0 -1];

% 设置种子点 SeedX, SeedY
SeedX=round(ImageWidth/2);
SeedY=round(ImageHeight/2);

%分配栈的空间GrowQueX，GrowQueY,把种子放入栈
GrowQueX(:,nEnd)=SeedX;
GrowQueY(:,nEnd)=SeedY;

% 定义的栈的起点和终点
nStart=1;
nEnd=1;

while(nStart<=nEnd)

    %当前点的坐标
    nCurrX=GrowQueX(nStart);
    nCurrY=GrowQueY(nStart);
    
    % 对当前点的4邻域的坐标，判断是否在图像内部
    
    for k=1:4
        % 对4邻域象素的坐标
        xx=nCurrX+nDx(k);
        yy=nCurrY+nDy(k);
        
        if (xx<ImageHeight)&(xx>1)&(yy<ImageWidth)&(yy>1)&(UnRegion(xx,yy)==0)&(abs(bw(nCurrX,nCurrY)-bw(xx,yy))<8)
            
            nEnd=nEnd+1;
            % 象素(xx,yy)压入栈
            
            GrowQueX(:,nEnd)=xx;
            GrowQueY(:,nEnd)=yy;
            UnRegion(xx,yy)=255;
            
        end

    end
    nStart=nStart+1;

end

figure;
imshow('result.bmp');
figure;
imshow(UnRegion);