Discussion :

[KrokroAus]

Bonjour a tous,

Je travaille actuellement sur des images satellite calibrees en termes "d'indice de Vegetation" - i.e, mes images ont des valeurs correspondant ~ a un taux de vegetation. J'ai des images tous les 16 jours, sur environ 7 ans.

J'ai un programme me permettant de relever les valeurs sur ces images le long d'un profil.
Je souhaiterai representer les valeurs des differents profils obtenus sur un meme graphique colore. L'idee serait de representer les valeurs par des couleurs et de suprposer les differents profils les uns sur les autres pour pouvoir etudier la variabilite interannuelle de la vegetation.

J'ai mis en fichier joint une figure illustrant le type de graphique que je souhaiterai obtenir (expl : figure A = Profils des variations interannuelles de la vitesse le long d'un profil). Existe t'il une fonction Matlab permettant de faire ca?
Merci pour votre aide.

[Jerome Briot]

Pourquoi ne pas stocker les valeurs des différents profils dans une matrice (une ligne par profil) ? Tu pourrais ensuite utiliser SURF.

Si les profils n'ont pas tous le même nombre de valeurs, tu devras passer par une étape d'interpolation (avec INTERP1) pour obtenir un même nombre de valeurs.

[KrokroAus]

Merci pour l'idee! Je vais tester et je vous tient au courant!

[KrokroAus]

Bonjour,

J'ai teste le graphique avec la fonction surface, ca marche tres bien, merci!
Juste un petit soucis...
Mon graphique est sous la forme :

Code :

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

surface(distance,jour,Valeurprofils)

Les donnees "jour" sont des dates (espacees regulierement) sous le format jj/mm/aaaa. La fonction surface ne les reconnais pas ... comment puis-je les transformer en un format reconnaissable??

Merci beaucoup!

[Jerome Briot]

Tu peux utiliser DATENUM pour convertir les dates en format numérique et DATETICK pour remplacer les valeurs numériques par leur équivalent en chaines de caractères sur le graphique

[ol9245]

Avec la fonction image ça marche aussi :

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
% fabrication des matrices X et Y
x = 1:16:1600 ; % jours
y = 1:100 ; % km
[X, Y] = meshgrid(x, y) ;
 
% fabrication d'un jeu de données
Data = (200 ./ (100+Y)) .* sin(X.*(1/160*2*pi) + Y.*(1/100*2*pi)) ;
Data = Data - min(Data(:)) ;
Data = round(255 .* (Data ./ max(Data(:)))) ;
 
% affichage
image(Data) 
legend('show')
title('Evolution interannuelle du NDVI quilé joli')
xlabel('distance le long du transect (km)')
ylabel('date d''observation (jour)')