Changeset 58

Timestamp:

06/14/13 17:31:23 (12 years ago)

Author:

vcheplygina

Message:

 

Location:

birds

Files:

• domilbirds.m (modified) (2 diffs)
• specwhiten.m (modified) (3 diffs)

birds/domilbirds.m

@@ -17 +17 @@
 figure(3); clf; showspec(imag(S), f,t);
 
+%%
 % smooth and threshold the spectrogram?
 G = fspecial('gaussian',[5 5],2); % window 5x5, sigma=2
@@ -33 +34 @@
 %%
 s2 = specwhiten(S,fmax); 
-figure(4); imagesc(t,f,s2);
+figure(5); imagesc(t,f,s2);
 axis xy; xlabel('time (s)'); ylabel('freq.');
+
+% smooth and threshold the spectrogram?
+G = fspecial('gaussian',[5 5],2); % window 5x5, sigma=2
+I = imfilter(abs(s2),G,'same');
+t = dd_threshold(I(:),0.95);
+mask = (I>t);
+mask(f<2000,:) = 0;
+
+figure(6); clf; imagesc(mask); axis xy;    %Ja, dit ziet er vergelijkbaar uit, iets minder kleine blobs aanwezig geloof ik? 
+% measure something:
+x = regionprops(bwlabel(mask),abs(s2));
+y = x([x.Area]>10)
+

birds/specwhiten.m

@@ -3 +3 @@
 function s2 = whiten(s,fmax)
 
-%Values output by spectrogram are already in [0, 1] range
+%Not sure about this, but the paper seems to assume this
+s = abs(s); 
+
+%Values output by spectrogram are already in [0, 1] range?
+
+
 
 %Square root of spectrogram (of the absolute value?)
-s1 = sqrt(abs(s)); 
-%s1 = sqrt(s);
+
+%s1 = sqrt(abs(s)); %Version 1
+s1 = s;             %Version for segmentation
 
 %For each frame, find a quantity similar to the energy"
@@ -21 +27 @@
 
 lowestix = find(energy<threshold);
-lowestenergy = energy(lowestix);
+%lowestenergy = energy(lowestix);
 
 
@@ -27 +33 @@
 s2 = s1;
 for f=1:Nf
-    P(f) = sqrt(eps + sum(s1(f,lowestix).^2));
+    %P(f) = sqrt(eps + sum(s1(f,lowestix).^2));
+    
+    P(f) = sum(s1(f,lowestix)) / length(lowestix);
     s2(f,:) = s1(f,:) ./P(f);
 end