// Augmentation du tampon
stacksize(10000000);

// Chargement de la bibliothèque SIP : gestion des images
exec("C:\Program Files\scilab-4.0\contrib\siptoolbox\loader.sce");

// Chargement de la fonction FFT sur une ligne
getf('C:\Documents and Settings\chatvache.info\Mes documents\_Jonathan_\FFT\fct_fft_rapide_optim.sci');

img=imread("C:\jpg.jpg"); // Lecture de l'image
img=normal(img);          // Normalisation de l'image
img=im2gray(img);         // Mise en niveau de gris
imgFinal=img;             // Image de retour
sensFFT=-1;               // Sens de la FFT : -1 directe, 1 indirecte

// Dimensions de l'image
LongueurImg=size(img,'c');
HauteurImg=size(img,'r');

// Tableau à 1 ligne
LigneImage=eye(LongueurImg);

//
// Traitement FFT par ligne
//

// On passe toutes les lignes : boucle sur le nombre de lignes
for h=1:HauteurImg

   // Récupération de la ligne à la hauteur h
   for l=1:LongueurImg
  	   LigneImage(l)=img(l,h);
   end;
  
  	// FFT sur la ligne
  	LigneImage=fct_fft_rapide_optim(LigneImage,-1,LongueurImg);
  	
  	// Recopie de la ligne obtenue dans l'image finale
  	for l=1:LongueurImg
  		  imgFinal(l,h)=LigneImage(l);
  	end;
	
end;

//
// Traitement FFT par colonne
//

// On passe toutes les colonnes : boucle sur le nombre de colonnes
for l=1:LongueurImg

   // Récupération de la ligne à la hauteur h
   for h=1:HauteurImg
  	   LigneImage(h)=imgFinal(l,h);
   end;
  
  	// FFT sur la colonne
  	LigneImage=fct_fft_rapide_optim(LigneImage,-1,HauteurImg);
  	
  	// Recopie de la colonne obtenue dans l'image finale
  	for h=1:HauteurImg
  		  imgFinal(l,h)=LigneImage(h);
  	end;
	
end;

imshow(imgFinal); // On affiche le résultat