% liredensite3d
%
% Cette fonction lit le tableau des 224 canaux pour un
% point de l'image hyperspectrale.
%
% Ligne de commande :
% Densite=liredensite3d(lig,col,Nom);
% ----> lig=numéro de ligne
% ----> col=numéro de colonne
% ----> Nom=nom de fichier (extension .img ou .rfl)
% ----> Densite=tableau des 224 valeurs des canaux pour le point (lig,col)
% 
% Paramètres par défaut :
% Densite=liredensite3d(lig,col);
% ----> Nom='f960705t01p02_r05_sc01.c.img'

function Densite=liredensite3d(lig,col,Nom)
if nargin==2
    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);

Densite=zeros(Canaux,1);

fseek(fid,2*Canaux*(Colonnes*lig+col),'bof');

for Canal=1:Canaux
    Pfort=fread(fid,1,'int8');
    Pfaible=fread(fid,1,'uint8');
    Densite(Canal)=256*Pfort+Pfaible;
end

figure;plot(Densite)
grid on
zoom on

fclose(fid);