Discussion :

[sylvain62]

Bonjour,

Voilà ma question, est il possible de récupérer des couples de coordonnées sur une cartes qui ne comprend que les coordonnées des lignes de côtes? Pour être plus clair j'attache une image.

Je ne dispose que des traits de côtes (2 colonnes: longitudes et latitudes). J'aimerai disposer des couples de points, à une résolution données, de tous les points compris à l'intérieur des traits de côtes (les coordonnées terrestres en fait). Ou les coordonnées des océans de la même façon.

Je ne sais pas si c'est clair lol, si quelqu'un à envie de m'aider, je ré explique sans soucis.

Merci

[Jerome Briot]

Pour générer les couples de points, tu peux utiliser MESHGRID :

Code :

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

[x,y] = meshgrid(-200:200,-80:100);

Ensuite, utilise la fonction INPOLYGON pour déterminer la position de chaque couple par rapport aux contours. Le plus dur consiste à définir proprement chaque contour.

Dans un soucis d'optimisation, tu peux d'abord définir chaque contour et seulement ensuite définir la grille des couples réduites pour chaque contour. Le problème à résoudre sera donc plus petit.

[sylvain62]

OK Merci pour ta réponse.

Je regarde ce que INPOLYGON fait... Enfin j'essaie lol.

Pour les contours, je dispose d'une matrice toute faite (en fait héhé). Qui est assez bien définie. Donc en vérité il ne me reste plus qu'à faire un meshgrid et d'appliqer la fonction de dessus si j'ai bien compris.

[Jerome Briot]

Je ne suis pas sûr que ce soit si simple... tu dois définir des polygones distincts pour chacun des contours, et non pas une seule matrice.

En gros, les polygones doivent approximativement ressembler à ceux générés par la fonction CONTOUR.

Sinon, pourrais-tu attacher un fichier mat (zippé) contenant la matrice des contours ?

[sylvain62]

ha euh ok....

Heu je vais peut-être mettre un zip alors... forcément c'est gros je vais essayer de réduire le nombre de coordonnées, je ne travaille pas sut tout la terre.

[sylvain62]

Voilà wc. txt (pour world coast pas pour autre chose:lol:)
colonne 1= longitudes
colonne 2= latitudes

[sylvain62]

Ce coup-ci la pièce est jointe....

[Jerome Briot]

Et sur quelle partie du globe souhaites-tu travailler ?

[sylvain62]

Sur l'hémisphère Nord.


latitudes +25 à +85
longitudes -80 à +70.5

[Jerome Briot]

Les contours ne sont pas continus.
Il y a soit des raccords, soit des trous...

Est-ce normal ?

L'idéal serait bien entendu d'avoir un polygon fermé unique par limite de côte

[sylvain62]

Oui c'est normal.

Enfin c'est comme ça qu'on les a.

Le fichier brute est dispo sur un site international.

Il s'agit d'une suite de point, que je relie toujours.

Code :

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

plot(wc(:,1),wc(:,2),'-k','LineWidth',2);

Ca suffit pour ce que je veux faire.

Pour éviter que des points soient reliés quand il ne faut pas, il y a des couples de nan entre les coordonnées.

J'avais réduit mes coordonnées à la zone qui m'intéresse mais du coup, il y a avait pleins de lignes supplémentaires qui reliées des coordonnées qui ne l'étaient pas sur la map mondiale.
Normale car en réduisant le ficheir coordonnées de base, je supprime des trous "nan nan"

[sylvain62]

Bonjour,

Je reviens sur ce sujet après une pause.

Bon effectivement c'est moins simple que je le croyais au début. Mais j'ai fais qqchose de bien j'ai lu l'aide (mais en vrai) lol
Pas en diagonale comme je fais d'habitude
Donc je vais essayer de faire ça proprement.

Effectivement la 1ere étape consiste à bien définir les contours. Ce qui n'est pas vraiment le cas avec cette matrices des points de côtes.

Je vais m'y mettre pendant les vacances on verra bien ce que ça donnera.

Bonne St-Sylvestre à tous