Le portail francophone de la géomatique

Salut,

Est ce que l'outil "rasterize (vector to raster)" de GDAL te donne le même résultat ?

Ouhlà...comme je ne connaissais que v.to.rast de GRASS, me voilà perdue avec GDAL..;

J'ai bien compris qu'il fallait que je choisisse quelle couche rasteriser, quel attribut affecter, le nom de la couche à la sortie, mais pour définir l'emprise sur laquelle je souhaite rasteriser, je suis perdue...

Comment configurer l'emprise sur laquelle je souhaite rasteriser? A quoi correspondent "horizontal" et "verticale"? Je souhaiterais tout simplement, comme dans GRASS, configurer l'emprise en fonction de ma couche du canevas (au final dans ce que je souhaite faire, il n'y a qu'une seule couche qui intervient...).

Quelqu'un saurait-il m'aider?

En vous remerciant par avance!

Bien cordialement,

Lucie D.

EDIT : je précise juste que je n'ai rien trouvé à ce sujet sur le web, ne tombant que sur des sites en anglais (aïe, va falloir que je révise mes cours...) et en plus sur...PYTHON, que je ne maitrise absolument pas...

Merci!

Dernière modification par Lucie D. (Tue 04 August 2015 08:41)

Super! Merci, ça fonctionne, et effectivement j'obtiens sensiblement la même chose sans mes deux grosses taches. Je me demande d'ailleurs pourquoi deux outils censés faire le même travail donnent des résultats différents...

En revanche, je remarque un léger décalage entre les pixels de l'ancienne couche (avec les taches) et la nouvelle (voir zoom ci-joint). De plus, tous les pixels ne sont pas positionnés de la même manière, mais cette différence est légère.

Dans mon travail, je me demande si cette différence risque de poser un souci puisque nous travaillons à une échelle plutôt fine (nous avons fixé la taille des pixels à 5mX5m).

En fait, je dois utiliser r.mapcalculator de GRASS pour compiler cette couche à d'autres couches correspondant à d'autres modes d'occupation du sol et ensuite comparer la couche finale avec celle d'un autre service avec lequel je travaille. La comparaison ne risque-t-elle pas d'être faussée à cause de cette couche rasterisée avec un autre outil que les autres?

Désolée si je m'exprime mal, pas évident à expliquer...

En tout cas, merci beaucoup pour votre aide!

Bien cordialement,

Lucie D.

EDIT : je précise que j'ai bien vérifié et mes couches ont bien la même projection!

Dernière modification par Lucie D. (Tue 04 August 2015 10:51)

Merci pour l'info.
Mais je suis désolée je vais encore devoir vous embêter...
Dans ma version de QGis, je trouve l'outil Rasterize (write over existing raster). Je suppose donc que c'est la même chose...?
Toutefois, quand j'utilise cet outil, il ne me demande que trois choses:
- Ma couche à rasteriser
- Mon attribut
- Mon raster emprise.

En réalité, depuis le début de mon affaire, j'utilise une couche raster que nos collègues "collaborateurs" nous ont envoyé. Elle me sert de guide depuis le début. Cette couche a été rasterisée avec un champ bien précis, et donc pas de 0 comme vous, mais plusieurs valeurs en fonction de l'occupation du sol. Rien de choquant jusqu'ici.

Quand je remplis ces trois champs dans la boite de dialogue GDAL, ça mouline un moment, puis la boite de dialogue disparait sans aucun message et je n'obtiens aucune nouvelle couche correspondant à ma couche rasterisée...

En gros il ne se passe rien...

Est-ce normal?

Merci!

Lucie D.

oui, c'est tout à fait normal. L'outil modifie les valeurs du raster emprise.
Si tu regardes les valeurs du raster qui t'as servi de guide, elles ont changé et doivent correspondre aux valeurs du champ que tu as rasterisé la où il y a de la donnée.

Attention donc a cet outil "rasterise over" car il altère le raster qui sert de guide... mieux vaut travailler sur une copie!
Le problème que tu vas avoir si ton raster emprise n'est pas vide c'est que seul les pixels compris dans la couche vecteur que tu rasterises vont être modifiés, les autres conserveront leur valeur d'origine. Tu vas donc te retrouver avec un mélange des deux couches.

L'idéal est donc de travailler sur un raster vide qu'on remplis avec les valeurs du champ de la couche vecteur.

Je réalise une étude qui me demande de croiser de nombreux rasters, un peu comme toi j'imagine. Ce que j'ai fais, c'est que j'ai crée un script qui permet de générer un raster vide sur l'emprise d'un polygone qui représente ma zone d'étude (création d'un raster rectangulaire 1000m plus grand que le polygone). C'est ce raster qui me sers ensuite de guide pour la création des différentes cartes que je souhaite croiser.

Les stagiaires n'ont pas de vacances

bon courage pour la suite de ton travail!

Bonjour!

Finalement, me voilà de retour car je rencontre un nouveau problème avec la rasterisation de mes couches via Gdal.

En effet, voilà comment j'ai procédé :

Je dois rasteriser 4 fois la même couche (toujours mes tronçons routiers+ferrés) suivant un champ différent (qui sont bien à chaque fois des entiers):

- code_nat
- comparaison
- ind
- syst

J'utilise donc Gdal rasterize (write over existing raster). Pour cela, j'ai besoin de l'emprise exacte que nos collaborateurs nous ont envoyée. Ainsi, j'utilise la calculatrice raster de GRASS r.mapcalculator sur cette couche emprise et j'indique la formule suivante : if(A<>0,0,0). J'obtiens donc bien mon emprise avec uniquement des 0. Je copie cette couche 4 fois (avec des noms différents) et ensuite je lance avec Gdal la rasterisation de ma couche vecteur. Mes 4 couches "raster emprise" sont donc modifiées.

Pour le code_nat, pas de souci, mes tronçons routiers ont bien pris les différentes valeurs de mon champ code_nat.

En revanche, pour comparaison, ind et syst, mes tronçons routiers ont pris une unique valeur qui ne correspond à rien (par exemple, pour comparaison, au lieu d'avoir 11000, j'avais un 92 alors que ce chiffre ne revient nulle part dans mon travail!!!

Alors comme je devais rendre mon travail assez rapidement, j'ai utilisé de nouveau la calculatrice raster r.mapcalculator, et, heureusement pour moi, mes codes comparaison, syst et ind n'étaient concernés que par des valeurs uniques! Donc dans r.mapcalculator, j'ai entré la formule suivante:


if(A<>0,"ma_valeur",0).

Et cela pour chacun de mes 3 champs restant.
Ça a bien fonctionné et j'ai utilisé r.stats pour vérifications plus contrôle visuel. Tout m'a semblé exact.

La première partie de mon travail est donc terminée et envoyée (les résultats étaient très satisfaisant à première vue), mais je vais devoir effectuer dans un deuxième temps la même manip, avec un code usage_bd qui lui comporte deux valeurs.
Je ne pourrai donc pas bidouiller avec r.mapcalculator!!!

Ou ai-je bien pu me tromper pour que ma couche raster ne prenne pas les bonnes valeurs?
Est-ce que mon bidouillage avec r.mapcalculator est risqué?

En vous remerciant par avance pour votre aide!!!

Bonne journée.

Bien cordialement,

Lucie D.

Salut,

Est ce que tu as une amplitude de valeurs supérieure à 256 à rasteriser?
Le problème peut venir de la. Si tu rasterises en Byte, le codage ne peut s'effectuer que sur 256 valeurs différentes... (0 --> 255)

voir par exemple: http://georezo.net/forum/viewtopic.php?pid=264210

Bonjour YoLecomte,

Et merci pour ta réponse.

Dans mes recherches j'étais effectivement tombée sur ce post mais j'étais loin de penser qu'il pouvait me concerner...
Dans mon code nat, j'ai 5 valeurs différentes : 209, 210, 211, 212, 213, valeurs toutes < à 256 effectivement.

Dans mes trois autres codes :

- comparaison : une seule valeur égale à 11 000 donc largement > 256!!!
- syst et ind : une seule valeur égale à 11 100 donc là encore largement > 256!!!

Pour mon code usage_bd, sur lequel il va falloir désormais que je travaille, j'ai deux valeurs : 411 et 412, qui sont bien elles aussi > 256, donc il y a de très fortes chances pour que le problème vienne bien de là.

En revanche, dans le lien que tu cites, j'ai cru comprendre qu'il fallait utiliser Python pour changer cette amplitude (aaahh ce python décidément il va vraiment falloir que je m'y mette!).

N'y a-t-il pas d'autre solution?

Si dans ma couche vecteur tronçons+ferrés je change les valeurs de mon code usage_bd, par exemple :

-pour 411 = 111
-pour 412 = 112

que je rasterise ensuite suivant ces codes, et qu'ensuite j'utilise r.mapcalculator en entrant la formule suivante :

if(A==111,411,if(A==112,412,A))

Penses-tu que ça en revienne au même?

Oui, je sais, j'aime bien les bidouillages, l'ennui c'est que parfois ça entraine des petites erreurs qui peuvent devenir énorme... =p

Encore merci pour ton aide!!!

Lucie D.

La ligne de commande concerne l'encodage des données raster.

Tu est limitée à 256 possibilités car un la valeur de tes pixels est codée sur 8 bits qui peuvent prendre la valeur soit 0 soit 1 (2^8 = 256).
mais d'autres systèmes d'encodage existe, par exemple int32 codes tes valeurs sur 32 bits, tu autorises donc une amplitude de plus de 2000 millions ((2^32)-1) de valeurs.

Pour l'outil rasterise (write over an existing layer), je ne sais pas trop comment ça fonctionne mais j'imagine que tu est limitée par l'encodage du raster que tu écrases donc si tu écris sur un raster en int32 ou même int16 tu devrais pouvoir obtenir les bonnes valeurs (c'est une supposition de ma part, il faudrait tester).

L'outil rasterise (vector to raster) te permet de choisir l'encodage des données que tu souhaites pour ton raster. Peut être il serait intéressant de passer par là pour ton traitement en choisissant Int16 comme paramètre Raster type.