Daniel@0: function i = intersect(a,b,thr)
Daniel@0: 
Daniel@0: if nargin<3
Daniel@0:     thr = .01;
Daniel@0: end
Daniel@0: 
Daniel@0: va = get(a,'PeakVal');
Daniel@0: pa = get(a,'PeakPos');
Daniel@0: ma = get(a,'PeakMode');
Daniel@0: if isa(a,'mirscalar')
Daniel@0:     xa = get(a,'FramePos');
Daniel@0: else
Daniel@0:     xa = get(a,'Pos');
Daniel@0: end
Daniel@0: vb = get(b,'PeakVal');
Daniel@0: pb = get(b,'PeakPos');
Daniel@0: mb = get(b,'PeakMode');
Daniel@0: if isa(b,'mirscalar')
Daniel@0:     xb = get(b,'FramePos');
Daniel@0: else
Daniel@0:     xb = get(b,'Pos');
Daniel@0: end
Daniel@0: 
Daniel@0: for j = 1:length(va)
Daniel@0:     for k = 1:length(va{j})
Daniel@0:         [nl nc np] = size(va{j}{k});
Daniel@0:         for c = 1:nc
Daniel@0:             for p = 1:np
Daniel@0:                 [pacp,ix] = sort(pa{j}{k}{1,c,p});
Daniel@0:                 vacp = va{j}{k}{1,c,p}(ix);
Daniel@0:                 macp = ma{j}{k}{1,c,p}(ix);
Daniel@0:                 [pbcp,ix] = sort(pb{j}{k}{1,c,p});
Daniel@0:                 xajk = xa{j}{k};
Daniel@0:                 xbjk = xb{j}{k};
Daniel@0:                 if isa(a,'mirscalar')
Daniel@0:                     xajk = mean(xajk);
Daniel@0:                 end
Daniel@0:                 if isa(b,'mirscalar')
Daniel@0:                     xbjk = mean(xbjk);
Daniel@0:                 end
Daniel@0:                 ia = 1;
Daniel@0:                 ib = 1;
Daniel@0:                 ii = 1;
Daniel@0:                 while ia <= length(pacp) && ib <= length(pbcp)
Daniel@0:                     if abs(xajk(pacp(ia))-xbjk(pbcp(ib)))<thr
Daniel@0:                         picp(ii) = pacp(ia);
Daniel@0:                         vicp(ii) = vacp(ia);
Daniel@0:                         micp(ii) = macp(ia);
Daniel@0:                         ia = ia+1;
Daniel@0:                         ib = ib+1;
Daniel@0:                         ii = ii+1;
Daniel@0:                     elseif xajk(pacp(ia))<xbjk(pbcp(ib))
Daniel@0:                         ia = ia+1;
Daniel@0:                     else
Daniel@0:                         ib = ib+1;
Daniel@0:                     end
Daniel@0:                 end
Daniel@0:             pa{j}{k}{1,c,p} = picp;
Daniel@0:             va{j}{k}{1,c,p} = vicp;
Daniel@0:             ma{j}{k}{1,c,p} = micp;
Daniel@0:             end
Daniel@0:         end
Daniel@0:     end
Daniel@0: end
Daniel@0: 
Daniel@0: i = set(a,'PeakPos',pa,'PeakVal',va,'PeakMode',ma);