% liredensite3danim
% 
% Cette fonction affiche en continu le tableau des 224 canaux
% pour tous les points de l'image.
% 
% Ligne de commande :
% liredensite3danim(Nom)
% ----> Nom=nom de fichier (extension .img ou .rfl)
% Paramètres par défaut :
% liredensiteanim
% ----> Nom='f960705t01p02_r05_sc01.c.img'
% Pour stopper l'affichage, taper CTRL-C


function liredensite3danim(Nom)
if nargin==0
    Nom='f960705t01p02_r05_sc01.c.img';
end
Colonnes=614;
Canaux=224;

fid=fopen(Nom,'r');
fseek(fid,0,'eof');
taille=ftell(fid);
Lignes=taille/(2*Colonnes*Canaux);
fseek(fid,0,'bof');

Densite=zeros(Canaux,1);

figure

for lig=1:Lignes
    for col=1:Colonnes
        for Canal=1:Canaux
            Pfort=fread(fid,1,'int8');
            Pfaible=fread(fid,1,'uint8');
            Densite(Canal)=256*Pfort+Pfaible;
        end
        plot(Densite)
        axis([1 224 -1000 20000])
        grid on
        drawnow
    end
end
fclose(fid);