new2.m

自己编的阈值分割的源程序

clc;
clear all;
T=imread('tire.tif');
[m,n]=size(T);
T=double(T);
Seedsx=[];
Seedsy=[];
Test=zeros(size(T));
seedx=floor(m/2);
seedy=floor(n/2);
Seedsx(1)=seedx;
Seedsy(1)=seedy;
head=1;
tail=1;
delta=45;
h=waitbar(0,'请等待...');
while head==1||head~=tail
     if (seedx-1)>=1&&Test(seedx-1,seedy)==0&&abs(T(seedx,seedy)-T(seedx-1,seedy))<=delta
        tail=tail+1;
        Seedsx(tail)=seedx-1;
        Seedsy(tail)=seedy;
        Test(seedx-1,seedy)=1;
    end
    if (seedx+1)<=m&&Test(seedx+1,seedy)==0&&abs(T(seedx,seedy)-T(seedx+1,seedy))<=delta
        tail=tail+1;
        Seedsx(tail)=seedx+1;
        Seedsy(tail)=seedy;
        Test(seedx+1,seedy)=1;
    end 
    if (seedy-1)>=1&&Test(seedx,seedy-1)==0&&abs(T(seedx,seedy)-T(seedx,seedy-1))<=delta
        tail=tail+1;
        Seedsx(tail)=seedx;
        Seedsy(tail)=seedy-1;
         Test(seedx,seedy-1)=1;
    end
    if (seedy+1)<=n&&Test(seedx,seedy+1)==0&&abs(T(seedx,seedy)-T(seedx,seedy+1))<=delta
        tail=tail+1;
        Seedsx(tail)=seedx;
        Seedsy(tail)=seedy+1;
         Test(seedx,seedy+1)=1;
    end
    head=head+1;
    seedx=Seedsx(head);
    seedy=Seedsy(head);
     waitbar(head/(m*n));   
end
close(h);
B=zeros(size(T));
for i=1:head
        B(Seedsx(i),Seedsy(i))=255;        
end
T=uint8(B);
imshow(T);
title('区域生长后的图片');
save T