Daniel@0: function c = combinepeaks(p,v,thr)
Daniel@0: % dedicated function for (Klapuri, 99) that creates a curve made of burst
Daniel@0: % at position of peaks p and with amplitude related to peaks v.
Daniel@0: 
Daniel@0: dp = get(p,'Data');
Daniel@0: dv = get(v,'Data');
Daniel@0: pp = get(p,'PeakPos');
Daniel@0: pv = get(v,'PeakPos');
Daniel@0: sr = get(v,'Sampling');
Daniel@0: l = length(dp);
Daniel@0: empty = cell(1,l);
Daniel@0: for i = 1:l
Daniel@0:     thr = round(thr*sr{i});
Daniel@0:     empty{i} = cell(1,length(dp{i}));
Daniel@0:     for h = 1:length(dp{i})
Daniel@0:         dih = zeros(size(dp{i}{h}));
Daniel@0:         for l = 1:size(pp{i}{h},3)
Daniel@0:             for k = 1:size(pp{i}{h},2)
Daniel@0:                 j = 1;
Daniel@0:                 ppkl = pp{i}{h}{1,k,l};
Daniel@0:                 pvkl = pv{i}{h}{1,k,l};
Daniel@0:                 while j < length(ppkl)
Daniel@0:                     if ppkl(j+1)-ppkl(j) < thr
Daniel@0:                         decreas = dv{i}{h}(pvkl(j+1),k,l) ...
Daniel@0:                                   < dv{i}{h}(pvkl(j),k,l);
Daniel@0:                         ppkl(j+decreas) = [];
Daniel@0:                         pvkl(j+decreas) = [];
Daniel@0:                     else
Daniel@0:                         j = j+1;
Daniel@0:                     end
Daniel@0:                 end
Daniel@0:                 dih(ppkl,k,l) = dv{i}{h}(pvkl,k,l);
Daniel@0:             end
Daniel@0:         end
Daniel@0:         dv{i}{h} = dih;
Daniel@0:     end
Daniel@0: end
Daniel@0: c = set(p,'Data',dv,'PeakPos',empty,'PeakVal',empty,...
Daniel@0:           'PeakPrecisePos',{},'PeakPreciseVal',{},'PeakMode',empty,...
Daniel@0:           'InterChunk',[]);