Le portail francophone de la géomatique

Merci pour ta réponse rapide.
En fait les données proviennent déja de données vectorielles (fichiers de point). Je vois ce que tu veux dire, il faudrait créer des surfaces à partir des classe de données (une surface par classe de données) puis remplir ces surfaces avec des couleurs différentes pour chaque classe.
Mais je t'avouerais que je ne sais pas comment définir une surface en vecteur à patir de données ponctuelles (relier les points entre eux...)

As tu une piste ?

Syl

Hello,

Alors, plusieurs trucs :

1. A partir d'un fichier de points, tu interpoles une matrice de pixels dont la dimension est défini avec g.region (tape -? à la fin pour les détails). Le logiciel a donc "deviné" grâce à IDW les valeurs intermédiaires.
2. Ensuite avant de passer en vecteur comme le suggère Ludovic, tu peux classer ton image (visiblement on voit genre 6 classes) pour homogénéiser les pixels en régions de même classe. Pour ca tu peux utiliser r.reclass.
3. Une fois tes pixels regroupés en classe, tu peux passer tout en vecteur (r.to.vect).

Il doit y avoir de quoi lisser tout ca (à fouiller) mais voilà pour le principe.

Tu as un fichier sur lequel on peut faire des essais ?

Robin.

En fait je vais essayer d'être plus clair :
Voici ma démarche pour arriver à la carte jointe dans mon premier message :

1. J'ai un fichiers de points égalements réparti selon une trame de 20 km x 20 km avec pour chaque coordonnées X, Y une classe de donnée (fichier vecteur). je parle de classe de donnée car mes points sont déja classifiés à l'origine (donc je pense que je n'ais pas faire de r.reclass comme tu me l'a indiqué Robin)

2. je génère mes surfaces de couleur à partir de mon vecteur avec un nombre de point = 1 (donc pas d'interpolation) en utilisant v.surf.idw et j'obtiens ainsi mes "carrés" avec une couleur par classe

Ce résultat est correct mais en fait je souhaiterais lisser mes frontières entre classe sans avoir de carrés grossiers mais plutot  des frontières curviligne. Il ne s"agit donc pas de créer des données intermédiaire à mes classes, ce que fait l'interpolation (je ne veux que 6 couleurs sur ma carte), mais de lisser mes coordonnées spatiales.

En fait ce que tu m'as dit m'a bien aider Robin, J'étais parti sur l'interpolation mais en fait il s'agit plus de convertir mes données de points en surface, je pourrais alors utiliser les cartes thématiques. je continu mes essais. Merci beaucoup pour vos réponses
Ca fait plaisir de pouvoir être aidé

Sylvain

Rebonjour,

Voilà, j'ai trouvé, ca me turlupinait depuis hier, je savais bien que ca existait :

En gros, cette notice de r.to.vect explique que comme les vecteurs produits par un raster peuvent apparaître un peu carrés (blocky), la fonction peut lisser les bords (option -s). La notice recommande aussi de bien vérifier la reolsution de la region et de la donnée originelle pour estimer les erreurs introduites par le lissage.

La commande donnera ainsi quelque chose comme ça :

r.to.vect input=monfichierraster output=monfichiervecteur feature=area -s

J'espère que ca te dépannera ! J'ai l'impression qu'on touche au but

Robin.

Les carrés proviennent de la résolution assez basse de la région lors de l'interpolation.

Plusieurs options:

- augmenter la résolution de la région avant v.surf.idw avec g.region res=
- augmenter la résolution de la région après v.surf.idw et interpoler le raster résultant avec r.resample.*
- interpoler les points vectoriels originels avec v.surf.rst qui en jouant sur les paramètres "smooth" et "tension" qui permettent d'arriver à un résultat plus lisse (de nouveau: il vaut mieux augmenter la résolution avant)


Moritz