- Timestamp:
- 06/14/13 17:40:57 (12 years ago)
- Location:
- birds
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
birds/domilbirds.m
r58 r59 37 37 axis xy; xlabel('time (s)'); ylabel('freq.'); 38 38 39 % smooth and threshold the spectrogram?40 G = fspecial('gaussian',[5 5],2); % window 5x5, sigma=239 % do the same smoothing but change the threshold 40 G = fspecial('gaussian',[5 5],2); 41 41 I = imfilter(abs(s2),G,'same'); 42 42 t = dd_threshold(I(:),0.95); … … 44 44 mask(f<2000,:) = 0; 45 45 46 figure(6); clf; imagesc(mask); axis xy; % Ja, dit ziet er vergelijkbaar uit, iets minder kleine blobs aanwezig geloof ik?46 figure(6); clf; imagesc(mask); axis xy; %Looks pretty similar, but less small blobs 47 47 % measure something: 48 48 x = regionprops(bwlabel(mask),abs(s2)); -
birds/specwhiten.m
r58 r59 35 35 %P(f) = sqrt(eps + sum(s1(f,lowestix).^2)); 36 36 37 P(f) = sum(s1(f,lowestix)) / length(lowestix); 37 P(f) = sum(s1(f,lowestix)) / length(lowestix); %Version for segmentation 38 38 s2(f,:) = s1(f,:) ./P(f); 39 39 end
Note: See TracChangeset
for help on using the changeset viewer.