Discussion :

[ol9245]

Bonjour,

J'ai un relief défini par des points xyz. Je peux voir ce relief avec trisurf.
J'ai aussi une image orthophotographique de ce relief (une image du relief, vue parfaitement de dessus).

Je veux draper (overlay) l'image sur la surface 3D générée par trisurf.

une recherche rapide dans la doc et dans mathworks n'a rien donné. je fais donc appel à vous .

Pour préciser ma question, voici un bout de code : Les parties "image" et "relief" génèrent séparément l'image orthophotographique et le relief 3D. Je veux draper la première sur la seconde.

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
 
% Partie image :
Im = imread('D:\_Photo\Orthos\IMG_0816_ortho.tif') ; % lecture image
imagesc([-5 65],[-5 65],Im) ; % affichage image géoréférencée (les vraies coordonnées)
axis square
 
% Partie relief
X = rand(20, 1) * 70 - 5 ;
Y = rand(20, 1) * 70 - 5 ;
Z = rand(20, 1) * 10 ;
X = [X ; -5; -5; 65; 65] ;
Y = [Y ; -5; 65; -5; 65] ;
Z = [Z ; 0;0;0;0] ;
Tri = delaunay(X, Y) ;
trisurf(Tri, X,Y,Z) ;
axis equal

[ol9245]

Je me réponds à moi-même :

Code :

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

trisurf(Tri,X,Y,Z,C)

C est l'image à draper

[Jerome Briot]

Cela revient à faire, de manière cachée, l'exemple de la documentation : http://www.mathworks.com/access/help...f/surface.html