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