Matlab区域生长函数regiongrow.m源码

regiongrow.m源码：

function [g,NR,SI,TI]=regiongrow(f,S,T)
f=tofloat(f);
if numel(S)==1
    SI=f==S;
    S1=S;
else
    SI=bwmorph(3,'shrink',Inf);
    S1=f(SI);
end
TI=false(size(f));
for K=1:length(S1)
    seedvalue=S1(K);
    S=abs(f-seedvalue)<=T;
    TI=TI|S;
end
[g,NR]=bwlabel(imreconstruct(SI,TI));

上面需要用到的tofloat函数源码：

function [out, revertclass] = tofloat(in)
identity = @(x) x;
tosingle = @im2single;

table = {'uint8', tosingle, @im2uint8
    'uint16', tosingle, @im2uint16
    'int16', tosingle, @im2int16
    'logical', tosingle, @logical
    'double', identity, identity
    'single', identity, identity};

classIndex = find(strcmp(class(in), table(:, 1)));

if isempty(classIndex)
    error('Unsupported input image class.');
end

out = table{classIndex, 2}(in);

revertclass = table{classIndex, 3};