Discussion :

[Programmeur_Aladdin]

Bonjour
Je suis entrain de faire un programme pour la registration des images astronomiques à l'aide de la Transformée de Fourier.
J'ai réussis à trouver la matrice de corrélation ainsi que sa transformée de Fourier inverse.
Maintenant, j'aimerais déterminer la position du pic dans cette transformée inverse en C++.
Avez-vous des pistes à me suggérer?
Merci d'avance

[pseudocode]

Juste pour voir si j'ai bien compris...

1. 2 signaux 2D (des images ?) dans le domaine spatial
2. Passage dans le domaine fréquentiel (FFT)
3. Calcul de la corrélation dans le domaine fréquentiel
4. Passage de la corrélation dans le domaine spatial (FFT-1)

Et maintenant tu cherches un pic (le maximum de corrélation) dans une surface 2D discrétisée.

J'ai bon ?

[Programmeur_Aladdin]

Exactement..
Est ce qu'il suffit de chercher le point max dans mon tableau qui represente l'image dans le domaine spatial?

[pseudocode]

Exactement..
Est ce qu'il suffit de chercher le point max dans mon tableau qui represente l'image dans le domaine spatial?

A priori: oui. Il faudra juste gérer le cas ou le "pic" est en fait un "plateau" => plusieurs occurrences de la valeur max.

[Programmeur_Aladdin]

Très bien.. Merci beaucoup pseudocode..
Maintenant, pour trouver bêtement mon pic (avant la gestion des cas d'exceptions), je fais tout simplement:

for(i=0; i<largeur*hauteur; i++)
{
if(InverseMatriceCorrelation[i+1]>InverseMatriceCorrelation[i])
{
*xPic = i+1;
*yPic = i/largeur;
}
}

Puisque l'image est compsée d'un tableau unidimensionnel, je divise par la largeur (x) pour que j'obtienne le (y) correspondant..

Qu'en pensez-vous? c'est correcte ou bien il existe une méthode plus facile?

[pseudocode]

Maintenant, pour trouver bêtement mon pic (avant la gestion des cas d'exceptions), je fais tout simplement:

for(i=0; i<largeur*hauteur; i++)
{
if(InverseMatriceCorrelation[i+1]>InverseMatriceCorrelation[i])
{
*xPic = i+1;
*yPic = i/largeur;
}
}

En quoi comparer systematiquement un element avec son suivant permet de trouver le maximum du tableau ?

Moi je verrai plutot quelquechose comme cela:

Code C :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
int xPic, yPic, vPic=0;
 
for(i=0; i<largeur*hauteur; i++)
{
	if(InverseMatriceCorrelation[i]>vPic)
	{
		xPic = i%largeur;
		yPic = i/largeur;						
		vPic = InverseMatriceCorrelation[i];
	}
}