Discussion :

[Abc0123]

Bonjour,

Je voudrais superposer une image obtenu grâce à un trimesh et une image jpg chargée depuis mon disque dur.

J'ai essayé le "hold on" en le mettant entre le trimesh et l'image, ainsi que l'inverse mais rien ne fonctionne, le résultat est plutôt bizarre et n'est pas celui escompté.

Y a-t-il un code spécial à appliquer ?

Merci

[Jerome Briot]

Si l'image doit être affichée dans le plan xy, tu peux faire simplement :

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
load mandrill
 
xp = rand(30,1)*500;
yp = rand(30,1)*480;
 
tri = delaunay(xp,yp);
 
figure
colormap(map)
 
imagesc(X)
hold on
trimesh(tri,xp,yp,'color','b','LineWidth',2)
 
axis equal tight

Si c'est dans un plan quelconque :

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
load mandrill
 
xp = rand(30,1)*500;
yp = rand(30,1)*480;
 
tri = delaunay(xp,yp);
 
xs = [1 1 ; size(X,2) size(X,2)];
ys = [1 size(X,1) ; 1 size(X,1)];
zs = zeros(2);
 
figure
colormap(map)
 
surface(xs,ys,zs,'facecolor','texturemap','cdata',rot90(X,-1));
hold on
trimesh(tri,xp,yp,'color','b','LineWidth',2)
 
axis equal tight

voir : http://www.developpez.net/forums/d13...s/#post7403447

[Abc0123]

Salut, merci de ta réponse.

Le premier code doit convenir. J'ai remplacé tes variables par les miennes mais j'ai un soucis dans le load. Pour charger mon image, le code

Code :

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

 load nom-de-mon-image

ne marche pas. Alors que je l'ai bien mise dans le workspace. Imload ne marche pas non plus.

De plus, à quoi correspond ton

Code :

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

imagesc(X)

?

Merci

[Jerome Briot]

Essaie ceci :

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
rgb = imread('lena.jpg','jpg');
 
xp = rand(30,1)*500;
yp = rand(30,1)*480;
 
tri = delaunay(xp,yp);
 
xs = [1 size(rgb,2) ; 1 size(rgb,2)];
ys = [size(rgb,1) size(rgb,1) ; 1 1];
zs = zeros(2);
 
figure
 
surface(xs,ys,zs,'facecolor','texturemap','cdata',rgb);
hold on
trimesh(tri,xp,yp,'color','b','LineWidth',2)
 
axis equal tight

[Abc0123]

Mon objet triangulé en 2D est sur le plan y = 12.

Du coup j'ai modifié ton code en conséquence, je pense avoir tout changé correctement.

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
 
rgb = imread('visage-homme.jpg','jpg');
 
xp = rand(30,1)*500;
zp = rand(30,1)*480;
 
tri = delaunay(xp,zp);
 
xs = [1 size(rgb,2) ; 1 size(rgb,2)];
zs = [size(rgb,1) size(rgb,1) ; 1 1];
ys = zeros(12);
 
figure
 
surface(xs,zs,ys,'facecolor','texturemap','cdata',rgb);
hold on
trimesh(tri,xp,zp,'color','b','LineWidth',12)

MATLAB me sort encore une petite erreur:

Error using surface
Matrix dimensions must agree.

Error in Implementation (line 52036)
surface(xs,zs,ys,'facecolor','texturemap','cdata',rgb);

Et rien ne s'affiche. En effet xs, ys et zs n'ont pas la même taille. De plus j'ai du mal à comprendre le début de ton code

[Jerome Briot]

Essaie ceci :

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
rgb = imread('lena.jpg','jpg');
 
xp = rand(30,1)*500;
yp = ones(30,1)*12;
zp = rand(30,1)*480;
 
tri = delaunay(xp,zp);
 
xs = [1 size(rgb,2) ; 1 size(rgb,2)];
ys = ones(2)*12;
zs = [size(rgb,1) size(rgb,1) ; 1 1];
 
figure(1)
 
surface(xs,ys,zs,'facecolor','texturemap','cdata',rgb);
hold on
patch('vertices',[xp(:) yp(:) zp(:)],'faces',tri,'facecolor','none','edgecolor','b','LineWidth',2)
 
grid on
view(3)

Les coordonnées xs, ys et zs définissent les quatre sommets du rectangle qui accueille l'image.

[Abc0123]

Ok c'est impec ça s'affiche merci.

Maintenant pour afficher ma figure triangulé affichée via trimesh et non des triangles au hasard,

Code :

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

1
2
3
4
5
6
7
8
9
x=Tab(:,1);
y=Tab(:,2);
z=Tab(:,3);
y(:) =12;
tri = delaunay(x,z);
figure
Bb = trimesh(tri,x,y,z);
axis equal
axis tight

Ce sont les "xp" et "zp" que je dois remplacer par mes "x" et "z" extraient du tableau j'imagine ?

[Jerome Briot]

Ce sont les "xp" et "zp" que je dois remplacer par mes "x" et "z" extraient du tableau j'imagine ?

Oui.

Utilise plutôt patch à la place de trimesh

[Abc0123]

Du coup mon code complet serait celui-ci :

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
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
x=Tab(:,1);
y=Tab(:,2);
z=Tab(:,3);
 
 
figure
tri = delaunay(Tab(:,1),Tab(:,2),Tab(:,3))
B = trimesh(tri,x,y,z);
axis equal
 
%% Projection sur un plan
 
y(:) =12;
tri = delaunay(x,z);
figure
Bb = trimesh(tri,x,y,z);
axis equal
axis tight
 
%% Superposition des images
 
rgb = imread('visage-homme.jpg','jpg');
 
xp = x;
yp = ones(30,1)*12;
zp = z;
 
tri = delaunay(xp,zp);
 
xs = [1 size(rgb,2) ; 1 size(rgb,2)];
ys = ones(2)*12;
zs = [size(rgb,1) size(rgb,1) ; 1 1];
 
figure(1)
 
surface(xs,ys,zs,'facecolor','texturemap','cdata',rgb);
hold on
patch('vertices',[xp(:) yp(:) zp(:)],'faces',tri,'facecolor','none','edgecolor','b','LineWidth',2)
 
grid on
view(3)

Mais du coup j'ai juste mon image .jpg qui s'affiche, il n'y a pas de superposition, ai-je oublier quelque chose ?

Merci beaucoup

[Jerome Briot]

Remplace :

Code :

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

yp = ones(30,1)*12;

par

Code :

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

yp = y;

[Abc0123]

Ça me donne ceci :



Le petit point bleu est la bonne chose à superposer cependant elle est toute petite. Sais-tu comment je pourrais arranger ceci ?

[Jerome Briot]

Je pensais que les dimensions du nuages de points correspondaient aux dimensions de l'image.

D'où l'utilisation de size dans les lignes suivantes :

Code :

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

1
2
3
4
 
xs = [1 size(rgb,2) ; 1 size(rgb,2)];
ys = ones(2)*12;
zs = [size(rgb,1) size(rgb,1) ; 1 1];

[Abc0123]

Je pensais pas que la taille importée.

Une idée de comment arranger ceci ?

Merci

[Abc0123]

Je voudrais que ce soit les images affichées qui se superposent, comme lorsque l'on utilise un "hold on", sans tenir compte de leur taille, je trouve pas où est le problème

[Invité]

Bonjour,

Si tu ne vois qu'un petit point bleu, c'est que tes coordonnées sont trop petites par rapport à la taille de l'image du visage...
Donc ajuste-les à la bonne dimension.

[Abc0123]

Impec c'est bon le problème est résolu, merci !