Changeset 58 for birds/domilbirds.m

Timestamp:

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

Author:

vcheplygina

Message:

 

File:

• birds/domilbirds.m (modified) (2 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)
+