diff toolboxes/MIRtoolbox1.3.2/MIRToolbox/combinepeaks.m @ 0:e9a9cd732c1e tip

first hg version after svn

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/toolboxes/MIRtoolbox1.3.2/MIRToolbox/combinepeaks.m	Tue Feb 10 15:05:51 2015 +0000
@@ -0,0 +1,40 @@
+function c = combinepeaks(p,v,thr)
+% dedicated function for (Klapuri, 99) that creates a curve made of burst
+% at position of peaks p and with amplitude related to peaks v.
+
+dp = get(p,'Data');
+dv = get(v,'Data');
+pp = get(p,'PeakPos');
+pv = get(v,'PeakPos');
+sr = get(v,'Sampling');
+l = length(dp);
+empty = cell(1,l);
+for i = 1:l
+    thr = round(thr*sr{i});
+    empty{i} = cell(1,length(dp{i}));
+    for h = 1:length(dp{i})
+        dih = zeros(size(dp{i}{h}));
+        for l = 1:size(pp{i}{h},3)
+            for k = 1:size(pp{i}{h},2)
+                j = 1;
+                ppkl = pp{i}{h}{1,k,l};
+                pvkl = pv{i}{h}{1,k,l};
+                while j < length(ppkl)
+                    if ppkl(j+1)-ppkl(j) < thr
+                        decreas = dv{i}{h}(pvkl(j+1),k,l) ...
+                                  < dv{i}{h}(pvkl(j),k,l);
+                        ppkl(j+decreas) = [];
+                        pvkl(j+decreas) = [];
+                    else
+                        j = j+1;
+                    end
+                end
+                dih(ppkl,k,l) = dv{i}{h}(pvkl,k,l);
+            end
+        end
+        dv{i}{h} = dih;
+    end
+end
+c = set(p,'Data',dv,'PeakPos',empty,'PeakVal',empty,...
+          'PeakPrecisePos',{},'PeakPreciseVal',{},'PeakMode',empty,...
+          'InterChunk',[]);
\ No newline at end of file