brede_vol_regident

 brede_vol_regident   - Region identification in a volume

       function T = brede_vol_regident(V, 'PropertyName',
          'PropertyValue') 
       
       Input:    V         Volume structure
       
       Property: info      [ {0} | An integer ] Amount of feedback/debug 
                           information 
                 neighbor  [ {6} | 10 ] Selects which
                           neighborhood rule is used to determine
                           the region 
 
       Output:   T         Labeled volume structure

       This function computes a list of connected regions (clusters) 
       by labelling all pixels by a region number, ie., so-called
       region identification (coloring or labelling or perhaps
       'clustering'). The labelling is performed on the V.volume
       field and no other field of V is used.

       The regions in T.volume is sorted so that the largest cluster
       has the label=1, the second largest label=2, ... The
       Non-labeled region is returned with zero.

       Example:
         V  = brede_read_analyze('spmF_0001')   % Load SPM result
         Vt = brede_vol_threshold(V, 26.71);    % Threshold
         Vl = brede_vol_regident(Vt);           % Region ident.
         Vl1 = Vl;                           
         Vl1.volume(Vl1.volume~=1) = 0;         % Only region 1
         figure, brede_ta3_frame                          
         brede_ta3_volsurf(Vl1, 'isoabsolute', 0.5);  % Plot

       See also BREDE, BREDE_VOL, BREDE_VOL_MAT2VOL,
                BREDE_VOL_GENCOORD, BREDE_VOL_VOL2LOC.

 $Id: brede_vol_regident.m,v 1.12 2008/04/23 14:42:55 fn Exp $