// 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");

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

// récupération de la largeur (m) et de la hauteur (n) de l'image
[m,n]=size(img);

// initialisation de la valeur de retour y à 0 partout
for i=1:m 
for j=1:n
        y(i,j)=0;
end;
end;

// On applique la FFT
for i=1:m // Ligne de l'image de retour
for j=1:n // Colonne de l'image de retour
        for k=1:m   // Somme sur les lignes
        for l=1:n   // Somme sur les colonnes
          // Calcul de la somme
          y(i,j)=y(i,j)+img(k,l)*exp(inve*2*%i*%pi*(((i-1)*(k-1)/m)+((j-1)*(l-1)/n)));
        end;
        end;
end;
end; 

z=real(y);  // On récupère la partie réelle
imshow(z);  // On affiche le résultat