Discussion :

[The-msx]

Salut .
Je cherche à implémenter la transformé de Fourrier (en C) .
J'ai préparé tout ce que je dois préparer ( les Structures , les variables ....tout est bien ) sauf que lors de la compilation le PC se bloque ! donc j'ai fais quelques modifications dans le code afin de localiser le problème . et j'ai trouver que mon compilateur se bloque lors du calcule de la transformé !!
Il prend énormément de temps pour calculer la transformer !et ceci est plus ou moins normale ! Car en fait pour calculer la transformé , j'ai fais 4 boucles POUR !

en fait pour une image 600*400 : j'ai fait comme ça :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
for(i=0;i<hauteur;i++)
      for(j=0;j<largeur;j++)
 
F[i][j]=0;// F est la transformé et Image est l'image initial 
for(x1=0;x1<hauteur;x1++)
      for(y1=0;y1<largeur;y1++)
 
F[i][j]= F[i][j] + Image [x1][y1]* cos(....)................;

donc la boucle contient 600*400*600*400 traitement

Ma question : Est ce qu'il y a autre méthode pour implémenter la transformé de Fourrier ou existe-il un outil pour accélérer le calcul ???


Merci d'avance

[millie]

Ma question : Est ce qu'il y a autre méthode pour implémenter la transformé de Fourrier ou existe-il un outil pour accélérer le calcul ???

Il ne vaut pas utiliser le calcul simple de la transformée de Fourier, sinon le temps de calcul va exploser. Il faut passer par une transformée de Fourier rapide (FFT). Si tu veux vraiment avoir les meilleurs performances, regarde du côté de la bibliothèque FFTW en C : http://humbert-florent.developpez.co...aitement/fftw/

[The-msx]

Il ne vaut pas utiliser le calcul simple de la transformée de Fourier, sinon le temps de calcul va exploser. Il faut passer par une transformée de Fourier rapide (FFT). Si tu veux vraiment avoir les meilleurs performances, regarde du côté de la bibliothèque FFTW en C : http://humbert-florent.developpez.co...aitement/fftw/

Salut
Merci bien millie pour votre réponse !


Je sais bien qu'il existe la FFTW , mais moi je dois écrire moi même la transformé !
Est ce que quelqu'un a une idée par exemple comment ils ont implémenter la FFTW ou comme j'ai dis une méthode pour l'implémentation?

[millie]

Salut
Merci bien millie pour votre réponse !


Je sais bien qu'il existe la FFTW , mais moi je dois écrire moi même la transformé !
Est ce que quelqu'un a une idée par exemple comment ils ont implémenter la FFTW ou comme j'ai dis une méthode pour l'implémentation?

Regarde sur le site de fftw, notamment : http://www.fftw.org/links.html

[The-msx]

Salut ,

J'ai utilisé la fftw (en C ) et j'ai obtenu la photo suivante ( à gauche l'image initial et à droite la nouvelle image après avoir appliquer la transformer de fourrier à la 1ére image :





Je veux vérifier si c'est bien déroulé ou pas?


PS : dans la lecture de l'image j'ai utiliser le format .bmp avec 24/bits (RVB), donc j'ai fais la transformé pour chaque canal de couleur .

[millie]

C'est quoi l'image de droite exactement ?

C'est le module de la transformée de Fourier ? (car ce sont des éléments complexes donc on ne peut pas a priori afficher directement l'image complexe)