Changeset 53

Timestamp:

06/14/13 16:16:48 (11 years ago)

Author:

dtax

Message:

Iets meer.

File:

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

birds/domilbirds.m

@@ -2 +2 @@
 dpath = fullfile(mildatapath,'birds');
 
-fname = 'PC13_20090531_050000_0.wav';
+fname = 'PC13_20090531_050000_1.wav';
 %fname = 'PC8_20090531_050000_99.wav';
 [signal,fs] = wavread(fullfile(dpath,'wavs',fname));
@@ -9 +9 @@
 
 figure(2); clf; showspec(S,f,t);
-
 figure(3); clf; showspec(imag(S), f,t);
 
+% smooth and threshold the spectrogram?
+G = fspecial('gaussian',[5 5],2); % window 5x5, sigma=2
+I = imfilter(abs(S),G,'same');
+t = dd_threshold(I(:),0.8);
+mask = (I>t);
+mask(f<2000,:) = 0;
+
+figure(4); clf; imagesc(mask); axis xy;
+
+