Discussion :

[rouby11]

Bonjour,
j'aimerais appliquer la transformée de Fourier par rapport à la coordonnée de temps, c'est a dire le vecteur remplit par (pixel1.1 de im1;pixel1.1 de im2, pixel1.1 de im3,..., pixel1.1 de im N)
puis calculer la phase et l'afficher dans une courbe.

Code qt :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
void forwardDFT(QVector<double> s, const int N, QVector<double> a, QVector<double> b)
{
 
    for (int k = 0; k <= N-1; ++k) {
        a[k] = b[k] = 0;
        for (int x = 0; x < N; ++x) {
            a[k] += s[x] * cos(2 * M_PI / N * k * x);
            b[k] += s[x] * sin(2 * M_PI / N * k * x);
        }
        // normalization
        a[k] *= (k == 0 || k == N/2) ? 1. / N : 2. / N;
        b[k] *= 2. / N;
    }
}

Code qt :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
void PFE::fftMag(){
   QVector<double> xf(widget.imagelistWidget->count());
    QVector<double> sainef(widget.imagelistWidget->count());
   QVector<double> sig_in(widget.imagelistWidget->count());
   QVector<double> a(widget.imagelistWidget->count());
   QVector<double> b(widget.imagelistWidget->count());
    for(int k=0;k<liste.size();k++){
 
      img->load(liste[k]);
     *img = img->scaled(widget.imageLabel->width(), widget.imageLabel->height(),Qt::IgnoreAspectRatio,Qt::FastTransformation);
      QPixmap pixmap=QPixmap(QPixmap::fromImage(*img));
    QImage imagefft=pixmap.toImage();
 
 
QRgb pixelfft=imagefft.pixel(100,100);
 
          double  moyp= qGray(pixelfft);
        sig_in[k] = moyp;
    forwardDFT(sig_in, widget.imagelistWidget->count(), a, b);
for(int l=0;l<liste.size();++l)
         { 
     resfft=atan2(b[l],a[l]);
sainef[l]=resfft;
       xf[l]=(double)l;
 
   }
//creer la courbe
myCurve->setSamples(xf.data(),sainef.data(),xf.size());
myPlot->replot();
    }

la courbe qui s'affiche est nulle tous les points sont des zeros comment resoudre ce problème SVP.

[ToTo13]

Tu appliques la FFT sur un tableau 1D. Donc tu mets ce que tu souhaites dans ton tableau avant traitement. Si tu souhaites mettre tes valeurs temporelles, rempli donc le tableau comme tu le souhaites, puis tu applique la FFT.

[rouby11]

C'est ce que j'ai fait mais le tableau se remplit par des valeurs nulles je ne sais pas pourquoi peut être il manque des pointeurs?la déclaration est elle juste?

[ToTo13]

As tu essayé d'appliquer une librairie qui calcule la FFT, il y en a plein en C++ ?

[rouby11]

Non j'ai essayé de m'en sortir avec opencv mais ça s'applique sur toute l'image alors qu'en réalité je veux comme entrée de fft un vecteur de pixels en temps c'est à dire appliquer la fft sur une séquence d'images suivant z puis à la fin récupérer les images de phase et les stocker.est ce que c'est clair? Avec le code ci dessus j'ai appliqué la fft sur un seul vecteur puis je veux parcourir sur tous les pixels suivant z ?aidez moi svp

[rouby11]

Bonjour,
j'ai reussi a remplir un vecteur et afficher ses valeurs jusqu'a maintenant tout va bien mais lorsque j'ai essayé de generaliser et parcourir tous les pixels ça influe sur la performance du programme et comme si il entre dans une boucle infinie svp aider moi a optimiser mon code.

Code qt :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
void forwardDFT(QVector<double> &s, const int N, QVector<double> &a, QVector<double> &b)
{
 
    for (int k = 0; k <= N-1; ++k) {
        a[k] = b[k] = 0;
        for (int x = 0; x < N; ++x) {
            a[k] += s[x] * cos(2 * M_PI / N * k * x);
            b[k] += s[x] * sin(2 * M_PI / N * k * x);
        }
        // normalization
        a[k] *= (k == 0 || k == N/2) ? 1. / N : 2. / N;
        b[k] *= 2. / N;
    }
}

Code qt :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
 double resfft=0;
     img->load(liste[0]);
     *img = img->scaled(widget.imageLabel->width(), widget.imageLabel->height(),Qt::IgnoreAspectRatio,Qt::FastTransformation);
      QPixmap pixmapone=QPixmap(QPixmap::fromImage(*img));
    QImage one=pixmapone.toImage();
   int  larg=one.width();
     int haut=one.height();
     for(int i=0;i<haut;i++)
    {
 
         for(int j=0;j<larg;j++)
         {
 
     for(int k=0;k<liste.size();k++){
 
      img->load(liste[k]);
     *img = img->scaled(widget.imageLabel->width(), widget.imageLabel->height(),Qt::IgnoreAspectRatio,Qt::FastTransformation);
      QPixmap pixmap=QPixmap(QPixmap::fromImage(*img));
    QImage imagefft=pixmap.toImage();
 
 
QRgb pixelfft=imagefft.pixel(i,j);
 
          double  moyp= qGray(pixelfft);
        sig_in[k] = moyp;
 
    } 
 
 
    forwardDFT(sig_in, widget.imagelistWidget->count(), a, b);  // forward DFT sur un vecteur des pixels en temps
 
 
 
       for(int l=0;l<liste.size();++l)//recuperer les valeurs de phase et remplacer les pixels correspondants par ces valeurs
         { 
           img->load(liste[l]);
         *img = img->scaled(widget.imageLabel->width(), widget.imageLabel->height(),Qt::IgnoreAspectRatio,Qt::FastTransformation);
       pixmapfin=QPixmap(QPixmap::fromImage(*img));
    QImage imgfinal=pixmapfin.toImage();
     resfft=atan2(b[l],a[l]);
       imgfinal.setPixel(i,j,qRgb(resfft, resfft, resfft));
}
}
}
//enregistrer la nouvelle liste d'images dans un dossier"pas encore fait"