[mjordan]

Ho visto ora che ti serve una routine di region growing...
Questo è quello che ho scritto tempo fa:

function [g, NR, SI, TI] = regiongrow(f, S, T)

f = double(f);

% Se S e' uno scalare, ottieni il "seme" dell'immagine.
if numel(S) == 1
         SI = f == S;
         S1 = S;
else
         SI = bwmorph(S, 'shrink', Inf);
         J = find(SI);
         S1 = f(J); % Questo e' l'array dei semi.
end

TI = false(size(f));

for K = 1:length(S1)
         seedvalue = S1(K);
         S = abs(f - seedvalue) <= T;
         TI = TI | S;
end

% Assegnamo un intero diverso ad ogni regione connessa.
[g, NR] = bwlabel(imreconstruct(SI, TI));

Tempo fa usavo Matlab 6.5 ma credo non dovresti avere problemi usando Matlab 7

Ciao