Discussion :

[jyjeanphi]

Bonjour, je dispose d'une image rectangulaire que j'aimerai déformer en coordonnées polaires (arc de 45°)
J'ai récemment programmé une fonction qui travaille sur les pixels de l'image et qui interpole pas à pas en polaire par rapport à la valeur du pixel en cartésien, mais ça prend un temps fou et c'est pas très facile à adapter quand on change de configuration d'image initiale.
Je pensais ce matin à travailler sur la matrice de mon image initiale (qui est en fait une matrice que je traçais avec imagesc), augmenter la résolution pour pouvoir gérer l'interpolation en polaire et retracer ligne par ligne et colonne par colonne ma matrice, mais ça va prendre un temps fou également.

Est-ce que qqun a déjà eu à réaliser ce type d'opération ? si oui, qu'a-t-il opté pour solution ?
N'existe-t-il pas une fonction sur matlab où on donne les équations de la surface d'arrivée et pour laquelle matlab conforme lui même l'image initiale à la nouvelle surface ? j'ai regardé l'aide sur conformal mapping mais ils n'expliquent pas comment définir la surface...

Merci beaucoup d'avance !

[Jerome Briot]

Tu pourrais attacher une image de ce que tu souhaites obtenir ?

[jyjeanphi]

Voici en gros ce que je voudrai effectuer

[jyjeanphi]

Voici ce que j'obtiens avec le travail pixel par pixel, mais je préfèrerai une méthode directe sur les matrices

[Jerome Briot]

Tu pourrais nous fournir les données de la dernière images pour pouvoir travailler sur les vraies valeurs ?

[jyjeanphi]

La matrice quand je l'enregistre en .mat fait 46 mo...
il n'y a pas une autre solution pour que je vous la transmette ?

[jyjeanphi]

Je vous fourni celle-ci qui contient les même informations que la précédente mais en valeur brute (pas de moyenne donc pas d'étalement comme l'image précédente)
Si on parvient à déformer celle-ci, on parviendra à déformer celle dont j'ai besoin

[jyjeanphi]

On garde l'orientation de la matrice (5128 lignes par 13 colonnes), les premières lignes vont se retrouver au sommet de mon triangle etc...

[Jerome Briot]

Peut être une idée :

Code :

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
load matrice
 
X = MatriceAmplitude_MOYENNE_complet;
 
figure
 
subplot(2,1,1)
imagesc(X);
 
R = 100;
 
[theta,r]=meshgrid(linspace(0,pi/4,size(X,2)),linspace(0,R,size(X,1)));
x = r.*cos(theta);
y = r.*sin(theta);
 
[xi,yi] = meshgrid(0:1:sin(pi/4)*R,0:1:R);
 
Xi = griddata(x,y,X,yi,xi);
 
subplot(2,1,2)
imagesc(Xi);
axis image

Il faudra peut être modifier l'argument method de la fonction GRIDDATA

Voici ce que ça donne avec la capture de l'image de ton message précédent :

[jyjeanphi]

Et ben dites donc ! chapeau bas pour l'efficacité !
je vais jouer avec les paramètres de manière à obtenir la résolution idéale.
En tout cas merci beaucoup !

[Jerome Briot]

Par simple curiosité, ça représente quoi cette image ?

[jyjeanphi]

La norme du courant marin obtenue par mesure radar.
J'ai testé d'autres jeux de paramètres, ça donne des qualités vraiment très bonnes, mais dieu que les temps de calculs sont longs...