Discussion :

[habasque]

Salut à tous,

Je travaille sur une matrice de coordonnées X et Y.

Maintenant, je souhaiterais représenter graphiquement la densité de ces points sur une grille 1 degré de longitude * 1 degré de latitude ?

Quelqu'un aurait t-il une idée quant à la fonction à utiliser ?

Merci d'avance.

[vinc-mai]

Est-ce que hist3 ne suffirait-il pas?

PS: envoie-moi un pisco sour!

[habasque]

Hola,

J'aurais dû mettre ce post en résolu...
Quoique j'aimerais améliorer la figure en y ajoutant le trait de côte du Pérou.
En fait, il s'agit de données océanographiques...

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
 
% definition de la grille ou l'on va faire les moyennes
dlon = 1;
dlat = 1;
Xgrid = -90:dlon:-70;
Ygrid = -25:dlat:5;
 
% calcul de la densité de données par point de grille
distrib = NaN*ones(length(Xgrid),length(Ygrid));
for i=1:length(Xgrid)
    for j=1:length(Ygrid)
        indProfiles = find(X'<=Xgrid(i)+dlon/2 & X'>Xgrid(i)-dlon/2 & Y'<=Ygrid(j)+dlat/2 & Y'>Ygrid(j)-dlat/2); 
        distrib(i,j)=size(indProfiles,1);
    end
end
 
%Réalisation de la carte de densité
figure;    
    contourf(Xgrid,Ygrid,distrib',[0:1:10000])
    shading flat 
    axis([-90 -70 -20 0])
    caxis([0 1000])
    colorbar
    hold on
    grid

PS: pour le pisco sour, passes quand tu veux la licuadora est prête !