% lireimage3d
%
% Cette fonction lit l'image associée à un canal.
%
% Ligne de commande :
% ImageCanal=lireimage3d(Canal,Nom);
% ----> Canal=numéro du canal (entre 1 et 224)
% ----> Nom=nom de fichier (extension .img ou .rfl)
% ----> ImageCanal=image correspondant au canal choisi
% 
% Paramètres par défaut :
% ImageCanal=lireimage3d(Canal);
% ----> Nom='f960705t01p02_r05_sc01.c.img'
%
% ImageCanal=lireimage3d;
% ----> Canal=100
% ----> Nom='f960705t01p02_r05_sc01.c.img'

function ImageCanal=lireimage3d(Canal,Nom)
if nargin==0
    Canal=100; % de 1 à 224
    Nom='f960705t01p02_r05_sc01.c.img';
elseif nargin==1
    Nom='f960705t01p02_r05_sc01.c.img';        
end
Canal=Canal-1; % de 0 à 223
Colonnes=614;
Canaux=224;

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

ImageCanal=zeros(Lignes,Colonnes);

for lig=1:Lignes
    for col=1:Colonnes
        Pfort=fread(fid,1,'int8');
        Pfaible=fread(fid,1,'uint8');
        ImageCanal(lig,col)=256*Pfort+Pfaible;
        fseek(fid,2*(Canaux-1),'cof'); 
    end
end

figure;imagesc(ImageCanal)
colorbar

fclose(fid);