Bonjour,
je souhaitais utiliser la FFT pour accélérer mes filtrages par convolution.
Je fais donc la procédure suivante :
- Transformation de la taille de mon image en puissance de 2 supérieure
- Padding de mon image (centrage) sur les nouvelles dimensions.
- FFT sur l'image.
- Padding et FFT sur mon masque de convolution (avec les nouvelles dimensions de l'image).
- Multiplication de mes deux résultats dans l'espace de Fourier.
- FFT inverse du résultat.
- Unpadding pour tout remettre dans mon image résultat.
Malheureusement, cela ne fonctionne pas pour des images non carrées ou non dyadiques :-(
J'ai testé deux algorithmes différents (Butterfly et Fast Hartley Transform de ImageJ), mais j'ai les même problèmes.
Y a t-il une solution pour gérer ce problème ?
Ou connaitriez vous des algorithmes qui gèrent correctement les images non dyadiques ?
Merci par avance
Partager