Discussion :

[Marinedu26]

Bonjour,

Etant novice en Matlab, j'ai du mal à trouver la fonction adéquate. En fait, j'aimerais savoir s'il est possible de construire une grille sur un shapefile de points. La grille serait tracée en spécifiant le coin supérieur gauche et inférieur droit à l'aide de coordonnées x,y puis il faudrait que je puisse spécifier la taille de chaque cellule de la grille. Dans mon cas, 5x5m. Est-ce qu'il existerait une fonction "grid" dans laquellle je puisse spécifier ces paramètres ?

Merci pour votre aide.
Marine

[Jean Dumoncel]

Bonjour

Code :

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

help meshgrid

?

[Marinedu26]

Bonjour,

Merci pour votre réponse. Je suis arrivée à faire ce que je voulais grâce à votre message mais je n'arrive pas à afficher les barres horizontales. Je m'explique en faisant le plot(xx,yy), j'ai seulement les barres verticales et pas les barres horizontales pour avoir une grille.

Code :

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

1
2
3
4
x=[635000:5:636450];
y=[5004160:5:5006040];
[xx,yy]=meshgrid(x,y);
plot(xx,yy)

Merci pour votre aide
Marine

[Jean Dumoncel]

Pour l'affichage, tu peux utiliser surf ou mesh par exemple :

Code :

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

mesh(xx,yy,zeros(size(xx)))

[Jerome Briot]

MESHGRID est inutile ici :

Code :

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

1
2
3
4
5
x = 635000:5:636450;
y = 5004160:5:5006040;
 
figure
mesh(x,y,zeros(numel(y),numel(x)))

[Marinedu26]

Bonjour,

Merci pour votre réponse. Je préfèrerais que le plot soit en 2D et non en 3D. En fait, je veux juste pouvoir modéliser une grille en 2D.

Merci pour votre aide.
Marine

[Jean Dumoncel]

C'est juste une question de point de vue (au sens propre!) :

Code :

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

view(2)

[Marinedu26]

Merci pour la réponse. Mais les cellules de la grille ne font pas 5 m x 5 m. J'ai testé la fonction pcolor. Le resultat est plutôt correcte. Je vais voir s'il est possible d'enlever les couleurs pour avoir juste un fond blanc et un quadrillage noir.

Merci pour votre aide.
Marine

[Jean Dumoncel]

Mais les cellules de la grille ne font pas 5 m x 5 m.

normalement si puisque tu as définis x et y avec un pas de 5

Qu'est-ce qui te fais dire que ce n'est pas le cas?

[Marinedu26]

Bonjour,

Merci pour vos réponses. Ah oui désolé, j'ai parlé trop vite. Effectivement, j'ai bien des cellules 5 m x 5 m. J'aurais dû zoomer sur la figure.
A partir de mon plot de points et de ma grille, je souhaiterais aussi déterminer le nombre de points qui se trouvent dans chaque cellule de la grille. En explorant les fonctions de MATLAB, j'ai trouvé une fonction qui a l'air pas mal.

Code :

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

n = hist3(data);

avec data qui contient les coordonnées x et y de mes points.
Grâce à cela, j'obtiens le nombre de points dans chaque cellule de la grille. Le problème est que, par défaut, la grille est faite avec 10 x 10 bins cad une grille 10 x 10 (il est bien sûr possible de modifier le nombre de bins). Est-ce qu'il serait possible de modifier les paramètres de cette fonction de sorte d'avoir des cellules 5 m x 5 m et non en se basant sur le nombre de bins ?

J'ai vu qu'il y avait possibilité de faire ça dans l'aide:

Code :

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

hist3(...,param1,val1,param2,val2,...)

Mais je ne sais pas vraiment comment l'exploiter et si ça résoudrait mon problème.

Merci pour votre aide.
Marine

[Invité]

Et que dis-tu de:

hist3(X,'Edges',edges), where edges is a two-element cell array of numeric vectors with monotonically non-decreasing values, uses a two-dimensional grid of bins with edges at edges{1} in the first dimension and at edges{2} in the second. The (i, j)th bin includes the value X(k,:) if

edges{1}(i) <= X(k,1) < edges{1}(i+1)
edges{2}(j) <= X(k,2) < edges{2}(j+1)