Changeset 72 for birds/createmilbirds_19cl.m

Timestamp:

07/09/13 15:27:17 (11 years ago)

Author:

vcheplygina

Message:

 

File:

• birds/createmilbirds_19cl.m (modified) (2 diffs)

birds/createmilbirds_19cl.m

@@ -46 +46 @@
 bagid = [];
 for i=1:B
+    i
+   
         %load the signal;
         [signal,fs] = wavread(fullfile(dpath,'src_wavs',names{i}(2:end)));
@@ -59 +61 @@
 
         % run over blobs:
-        im = abs(S);
-        thisx = zeros(Nseg,7);
+        absim = abs(S);
+        realim = real(S);
+    imagim = imag(S);
+    
+    maskfeats = nan(Nseg,3);
+    absfeats = nan(Nseg, 7);
+    realfeats = nan(Nseg,7);
+    imagfeats = nan(Nseg,7);
+        
         for j=1:Nseg
+   
+        
                 ix = (bloblab==j);
+        
       % compute/add some blob-properties:
-                thisx(j,:) = [props(j).Area, props(j).Centroid, props(j).BoundingBox];
+          %     thisx(j,:) = [props(j).Area, props(j).Centroid, props(j).BoundingBox];
+      
+        pixtotal = sum(sum(ix));
+        pixheight = max(sum(ix,1));
+        pixwidth = max(sum(ix,2));
+
+        maskfeats(j,:) = [pixtotal pixheight pixwidth];
+
+
+        seg = absim(ix);
+        absfeats(j,1) = mean(seg); 
+        absfeats(j,2) = std(seg); 
+
+        absfeats(j,3) = quantile(seg(:),0);
+        absfeats(j,4) = quantile(seg(:),0.25);
+        absfeats(j,5) = quantile(seg(:),0.5);
+        absfeats(j,6) = quantile(seg(:),0.75);
+        absfeats(j,7) = quantile(seg(:),1);
+        
+           seg = realim(ix);
+        realfeats(j,1) = mean(seg); 
+        realfeats(j,2) = std(seg); 
+
+        realfeats(j,3) = quantile(seg(:),0);
+        realfeats(j,4) = quantile(seg(:),0.25);
+        realfeats(j,5) = quantile(seg(:),0.5);
+        realfeats(j,6) = quantile(seg(:),0.75);
+        realfeats(j,7) = quantile(seg(:),1);
+        
+           seg = imagim(ix);
+        imagfeats(j,1) = mean(seg); 
+        imagfeats(j,2) = std(seg); 
+
+        imagfeats(j,3) = quantile(seg(:),0);
+        imagfeats(j,4) = quantile(seg(:),0.25);
+        imagfeats(j,5) = quantile(seg(:),0.5);
+        imagfeats(j,6) = quantile(seg(:),0.75);
+        imagfeats(j,7) = quantile(seg(:),1);
+          
+        
                 % don't forget:
-                %bagid(end+1) = i;
-        
-        bagid = [bagid i];
+                bagid(end+1) = i;
+
         
         end
-        x{i} = thisx;
+        x{i} = [maskfeats absfeats realfeats imagfeats];
 
         %Get the labels right for the training bags: