function L = detection_dauphin(wavfile)
% 10, 100 : dauphin 
% 100, 600 : cachalot 1/3 s et 2 sec 

 outdir =/NAS3/LUCAS/txtcachalots
 
% relancé pour TARA 
% 25 mai 
 
    [s,Fs] = wavread(sprintf('%s%s','/NAS3/LUCAS/',wavfile));
 %   L=[ ];
    
    taille_clic=1; % taille moyenne d'un clic (en bins) pour 50 kHz fe ou 96kHz
    nb_clics=50; % nb de clics à garder par 10 sec % peut dépendre des saisons / activités
    % ICI moy 20 ms
%   taille_fen=160; % taille de la fenetre pour décimation, donc multiple de 50000 Hz

taille_fen=320; % pour 96 kHz 


    taille_var=11; % taille fenetre jitter apres decimation, environ 20 ms
    % tursiops adunctus ICI ( 20 - 60 - 280 ... ms )   
    lim_inf=100; % -> ICI min = lim_inf*taille_fen/Fs = environ 32 ms
    lim_sup=600; % -> ICI max = lim_sup*taille_fen/Fs = 320 ms % ROC peut augmenter avec fenetre plus longue apres une diminution, c'est un autre mode pour clic isolé, mais à valider en cas bruit chantier.

    for j=1:size(s,1)/(Fs*10) 
      % 0) DECOUPAGE FENETRE DE 10 SECONDES
      ps=s(1+(j-1)*Fs*10:j*Fs*10); 

      % 1) Passe haut robuste à une petite modulation
      pstk = tk_centre(ps); % 4 pts de + / tk

      % 2) EXTRACTION DES MAX LOCAUX 
      psml=(ps>[ps(1);ps(1:end-1)])&(ps>[ps(2:end);ps(end)]); 
      pstkml= pstk.*psml; % max piqué ET local : 2 pts de + / tk 
     
      % 3) CONSERVATION DES nb_clics MAX LOC LES PLUS GRANDS
      so=sort(pstkml); 
      pstkml=pstkml>so(end-nb_clics);

      % 4) DECIMATION : on ne garde que l'intensité dans la fenetre
      psd= sum(reshape(pstkml,taille_fen,length(ps)/taille_fen))>0;

      % 5) JITTER ICI de 20 ms
      psj=conv(double(psd),hamming(taille_var)); 

      % option centrage psj=psj(1+(taille_var-1)/2:end-(taille_var-1)/2);

      % 6) AUTOCORR
      x=xcorr(psj,psj);

      % 7) Pic local autocorr
      x = x((end+1)/2+lim_inf:(end+1)/2+lim_sup) / x((end+1)/2); % fenetre ICI, normee' par E 
      ml = (x>[x(1) x(1:end-1)])&(x>[x(2:end) x(end)]); % les max loc dans x
      indic(j) = max(x.* ml); % on annule les non maxloc, et on prend le max => 85.5 %auc

    %  presence(i,j) = indic(i,j) > seuil à choisir sur ROC 
    % fin 

    end;
    dlmwrite(sprintf('%s%s.txt',outdir,wavfile), indic);

end

