Matlab寫的區(qū)域生長(zhǎng)圖像分割程序。 %區(qū)域生長(zhǎng)算法:region function LabelImage=region(image,seed,Threshold,maxv) %image:輸入圖像 %seed：種子點(diǎn)坐標(biāo)堆棧 %threshold：用鄰域近似生長(zhǎng)規(guī)則的閾值 %maxv：所有生長(zhǎng)的像素的范圍小于maxv % LabelImage:輸出的標(biāo)記圖像，其中每個(gè)像素所述區(qū)域標(biāo)記為rn [seedNum,tem]=size(seed);%seedNum為種子個(gè)數(shù) [Width,Height]=size(image); LabelImage=zeros(Width,Height); rn=0;%區(qū)域標(biāo)記號(hào)碼 for i=1:seedNum %從沒有被標(biāo)記的種子點(diǎn)開始進(jìn)行生長(zhǎng) if LabelImage(seed(i,1),seed(i,2))==0 rn=rn+1;% %對(duì)當(dāng)前生長(zhǎng)區(qū)域賦標(biāo)號(hào)值 LabelImage(seed(i,1),seed(i,2))=rn; % end stack(1,1)=seed(i,1);%將種子點(diǎn)壓入堆棧（堆棧用來在生長(zhǎng)過程中的數(shù)據(jù)坐標(biāo)） stack(1,2)=seed(i,2); Start=1;%定義堆棧起點(diǎn)和終點(diǎn) End=1; while(Start<=End) %當(dāng)前種子點(diǎn)坐標(biāo) CurrX=stack(Start,1); CurrY=stack(Start,2); %對(duì)當(dāng)前點(diǎn)的8鄰域進(jìn)行遍歷 for m=-1:1 for n=-1:1 % %判斷像素（CurrX，CurrY）是否在圖像內(nèi)部 % rule1=(CurrX+m)=1&(CurrY+n)=1; % %判斷像素（CurrX，CurrY）是否已經(jīng)處理過 % rule2=LabelImage(CurrX+m,CurrY+n)==0; % %判斷生長(zhǎng)條件是否滿足 % rule3=abs(double(image(CurrX,CurrY))-double(image(CurrX+m,CurrY+n)))<Threshold; % %條件組合 % rules=rule1&rule2&rule3; if (CurrX+m)=1&(CurrY+n)=1&LabelImage(CurrX+m,CurrY+n)==0&abs(double(image(CurrX,CurrY))-double(image(CurrX+m,CurrY+n)))<=Threshold&image(CurrX+m,CurrY+n)0 %堆棧的尾部指針后移一位 End=End+1; %像素（CurrX+m,CurrY+n)壓入堆棧 stack(End,1)=CurrX+m; stack(End,2)=CurrY+n; %把像素（CurrX,CurrY)設(shè)置成邏輯1 LabelImage(CurrX+m,CurrY+n)=rn; end end end %堆棧的尾部指針后移一位 Start=Start+1; end end