Changeset 61

Timestamp:

06/17/13 17:14:24 (12 years ago)

Author:

vcheplygina

Message:

 

Location:

birds

Files:

• domilbirds.m (modified) (2 diffs)
• gendatbirds_veronika.m (added)

birds/domilbirds.m

@@ -28 +28 @@
 % measure something:
 x = regionprops(bwlabel(mask),abs(S));
-y = x([x.Area]>10)
+y1 = x([x.Area]>10)
 
 
@@ -45 +45 @@
 
 figure(6); clf; imagesc(mask); axis xy;    %Looks pretty similar, but less small blobs 
-% measure something:
-x = regionprops(bwlabel(mask),abs(s2));
-y = x([x.Area]>10)
+
+labeledx = bwlabel(mask);
+Nseg = max(unique(labeledx)); %Segment labels, segment 0 is background
+
+absim = abs(S);
+realim = real(S);
+imagim = imag(S);
+
+
+maskfeats = nan(Nseg,3);
+absfeats = nan(Nseg, 7);
+
+
+for i=1:Nseg
+   
+    ix = (labeledx == i);
+    
+    pixtotal = sum(sum(ix));
+    pixheight = max(sum(ix,1));
+    pixwidth = max(sum(ix,2));
+
+    maskfeats(i,:) = [pixtotal pixheight pixwidth];
+    
+    
+    seg = absim(ix);
+    absfeats(i,1) = mean(seg); 
+    absfeats(i,2) = std(seg); 
+    
+    absfeats(i,3) = quantile(seg(:),0);
+    absfeats(i,4) = quantile(seg(:),0.25);
+    absfeats(i,5) = quantile(seg(:),0.5);
+    absfeats(i,6) = quantile(seg(:),0.75);
+    absfeats(i,7) = quantile(seg(:),1);
+    
+end
+
+mildata = [maskfeats absfeats];