diff mirex2012-matlab/computeCQT.m @ 2:8017dd4a650d

Add MIREX 2012 code

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/mirex2012-matlab/computeCQT.m	Wed Mar 19 09:09:23 2014 +0000
@@ -0,0 +1,30 @@
+function [intCQT] = computeCQT(filename)
+% Settings for computing CQT for music signals (by E. Benetos)
+
+% Load .wav file
+[x fs bits] = wavread(filename);
+if (size(x,2) == 2) y = mean(x')'; clear('x'); else y=x; clear('x'); end;
+if (fs ~= 44100) y = resample(y,44100,fs); end;
+y = 0.5*y/max(y);
+fs = 44100;
+
+
+% Compute CQT
+Xcqt = cqt(y,27.5,fs/3,60,fs,'q',0.80,'atomHopFactor',0.3,'thresh',0.0005,'win','hann');
+%Xcqt = cqt(y,27.5,fs/3,120,fs,'q',0.35,'atomHopFactor',0.3,'thresh',0.0005,'win','hann'); % old resolution
+absCQT = getCQT(Xcqt,'all','all');
+
+% Crop CQT to useful time regions
+emptyHops = Xcqt.intParams.firstcenter/Xcqt.intParams.atomHOP;
+maxDrop = emptyHops*2^(Xcqt.octaveNr-1)-emptyHops;
+droppedSamples = (maxDrop-1)*Xcqt.intParams.atomHOP + Xcqt.intParams.firstcenter;
+outputTimeVec = (1:size(absCQT,2))*Xcqt.intParams.atomHOP-Xcqt.intParams.preZeros+droppedSamples;
+
+lowerLim = find(outputTimeVec>0,1);
+upperLim = find(outputTimeVec>length(y),1)-1;
+
+%intCQT = absCQT(112:1200,lowerLim:upperLim); % old resolution
+intCQT = absCQT(56:600,lowerLim:upperLim);
+
+
+%figure; imagesc(imrotate(abs(intCQT),90));