Daniel@0: function y = peaksegments(func,d,p,varargin)
Daniel@0: 
Daniel@0: y = cell(1,length(d));
Daniel@0: for h = 1:length(d)
Daniel@0:     if not(iscell(d{h})) % for histograms
Daniel@0:         p{h} = {p{h}(:,:,1)'};
Daniel@0:         d{h} = {d{h}'};
Daniel@0:     end
Daniel@0:     for i = 1:length(d{h})
Daniel@0:         %sum(var(p{h}{i},1,2))
Daniel@0:         %p{h}{i} = repmat(p{h}{i}(:,1,1),[1,size(p{h}{i},2),size(p{h}{i},3)]);
Daniel@0:         d{h}{i} = d{h}{i} + repmat(min(d{h}{i}),[size(d{h}{i},1),1,1]);
Daniel@0:         for j = 1:size(d{h}{i},3)
Daniel@0:             for k = 1:size(d{h}{i},2)
Daniel@0:                 yk = [];
Daniel@0:                 nb = 0;
Daniel@0:                 f = find(not(isnan(d{h}{i}(1:end,k,j))));
Daniel@0:                 while not(isempty(f))
Daniel@0:                     l = f(1);
Daniel@0:                     f = find(isnan(d{h}{i}(l:end,k,j)))+l-1;
Daniel@0:                     if isempty(f)
Daniel@0:                         f = size(d{h}{i},1)+1;
Daniel@0:                     end
Daniel@0:                     r = l:f(1)-1;
Daniel@0:                     if length(r)>1
Daniel@0:                         nb = nb+1;
Daniel@0:                         if nargin > 3
Daniel@0:                             for g = 1:length(varargin)
Daniel@0:                                 if g == 1
Daniel@0:                                     v{g} = repmat(varargin{g}{h}{i}{1,k,j}(nb),...
Daniel@0:                                               [length(r),1,1]);
Daniel@0:                                 else
Daniel@0:                                     v{g} = varargin{g}{h}{i}{1,k,j}(nb);
Daniel@0:                                 end
Daniel@0:                             end
Daniel@0:                         else
Daniel@0:                             v = {};
Daniel@0:                         end
Daniel@0:                         yk(nb,1) = func(abs(d{h}{i}(r,k,j)),...
Daniel@0:                                         p{h}{i}(r,k,1),...
Daniel@0:                                         v{:});
Daniel@0:                     end
Daniel@0:                     f = find(not(isnan(d{h}{i}(f(1):end,k,j))))+f(1)-1;
Daniel@0:                 end
Daniel@0:                 y{h}{i}{1,k,j} = yk;
Daniel@0:             end
Daniel@0:         end
Daniel@0:     end
Daniel@0: end