Discussion :

[Ogrim]

Bonjour =)

J'ai joint en piece jointe deux cartes vu en Lambert.
Elles ont été généré par le code suivant :

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
    subplot(1,2,1)
    [lat1,lon1]=meshgrid(lon1,lat1);
	[xp,yp,alat,alon,xi,yi]=createcoord(resol);
	h=displaym(myworldlo('POpatch'));
    set(h,'facec',[.5 .5 .5])
    pcolorm(lon1,lat1,squeeze(w1(t,:,:)))
    maxiw1=max(max(max(squeeze(w1(t,:,:)))));
    caxis([0 (maxiw1)])
    title (['VAGUES carte des Hs pour la zone [',num2str(latitude(1)),',',num2str(latitude(2)),'] ; [',num2str(longitude(1)),',',num2str(longitude(2)),'] a la date t= ', num2str(t),' (', num2str(mm),'/',num2str(yy),')'])
    colorbar
 
    subplot(1,2,2)
    HS_anomalies=squeeze(w1(t,:,:))-squeeze(model(mm,:,:));
    [xp,yp,alat,alon,xi,yi]=createcoord(resol);
	h=displaym(myworldlo('POpatch'));
    set(h,'facec',[.5 .5 .5])
    pcolorm(lon1,lat1,HS_anomalies)
    maxiw1=max(max(squeeze(HS_anomalies(:,:))));
    caxis([0 (maxiw1)])
    title (['VAGUES anomalies de  vagues pour la zone [',num2str(latitude(1)),',',num2str(latitude(2)),'] ; [',num2str(longitude(1)),',',num2str(longitude(2)),'] a la date t= ', num2str(t),' (', num2str(mm),'/',num2str(yy),')'])
    colorbar

J'aimeria en fait represter la zone (des cartes intitulées 'zone') par un rectangle.
Autrement dit : que la carte soit la meme que celle ou est représenté tout le globe (image intitulée 'tout'), avec juste un rectangle pour situer la zone d'étude.

Je pense qu'il faut utiliser la fonction 'contourm' mais j'ai du mal à l'utiliser.

[Jerome Briot]

Je ne suis pas sûr de bien comprendre ta question...

Si le problème vient du changement d'échelle colorimétrique entre les deux groupes d'image, utilise la fonction CAXIS pour forcer les deux graphiques à utiliser la même échelle.

C'est ça ?

[Ogrim]

non pas tout à fait,
je voudrais juste faire apparaitre sur la première image le contour de la zone de la deuxième.
c'est à dire : tracer un rectangle correspondant à ma zone plutot que de supprimer tout ce qui est en dehors de la zone.
merci beaucoup =)

[Jerome Briot]

Pourquoi ne pas utiliser LINEM à la place de CONTOURM ?

[Ogrim]

oui effectivement, pourquoi pas ^^
je ne connaissais pas cette fonction,
je veux essayer.
En tout cas merci beaucoup =)