superlets.m

function [duration,f_centroid]=superlets(Section,F_ds,fois)

Test=0; stop=0;
init_thresh=0.65;
srord   = [1, 3]; 
Clicks_bank= Section;                     % pick region of analysis
time_pulse= [0:1/F_ds:(1/F_ds)*(length(Clicks_bank)-1)]'; 
xS_slt=aslt(Clicks_bank, F_ds, fois, 3, srord, 0);              % calculate spectogram 


% figure;
% subplot(2,1,1);
%         imagesc(1e3*time_pulse, fois, xS_slt ); xlabel('t[ms]');
%         set(gca, 'ydir', 'normal');
%         colormap jet;

M_val_sel=max(max(xS_slt));                 % calculate the maximum value within the spectogram
slt_NORM=xS_slt./M_val_sel; 
[f_m,d_m]=find(slt_NORM==1);
 ref= [d_m f_m]'; 
while Test<0.05 && stop<100
    stop=stop+1;
    Binary=double(slt_NORM>init_thresh);                                    
    init_thresh=init_thresh-0.01;
    Binary(end,1:end)=zeros(1,size(Binary,2));     
    Binary(1,1:end)=zeros(1,size(Binary,2));
    Binary(1:end,end)=zeros(1,size(Binary,1));
    Binary(1:end,1)=zeros(1,size(Binary,1));
    Test=mean(mean(Binary));
end 
% subplot(2,1,2);
%         imagesc(1e3*time_pulse, fois, Binary ); xlabel('t[ms]');
%         set(gca, 'ydir', 'normal');
%         colormap jet;
%  

if Test>0.05
    m=iblobs(Binary);
    m_inds=find(m.area>150 & m.area<7000);
    M=m(m_inds);

    for b_ind=1:length(M)
      E_dist(b_ind)=norm(ref-M(:,b_ind).p);
    end
    roi=find(E_dist==min(E_dist));
    E_dist=[];

    % roi=find(M.vc==max(M.vc));

    duration=1e3*(M(roi).umax-M(roi).umin)/F_ds ;       % calculate the pulse duration [ms]
    f_centroid=1e-3*fois(f_m)   ;                       % calculate the dominant frequency  [khz]
    % f_centroid=1e-3*fois(round(M(roi).vc));           % calculate the frequency centroid [khz]

    if length(duration)>1
        duration=0.005;
        f_centroid=2;
    end
else
        duration=0;
        f_centroid=2;
end

end