Le portail francophone de la géomatique

Cf ma réponse à ton post précédent
http://georezo.net/forum/viewtopic.php?id=78124

Extrait >>>
Si j'ai bien compris tu as un raster de pentes avec certaine étendue et tu voudrais en extraire une étendue plus petite.
As-tu essayé de découper ton raster via le menu Raster > Extraction > Découper avec un mode de découpage par "Emprise". Tu peux alors tracer à la souris un rectangle d'emprise (qui apparait rouge) sur ton raster qui englobe ta zone d'études. Ensuite tu spécifies un fichier de sortie puis OK, ça devrait fonctionner sans prendre de temps. Évidemment comme ta zone d'études n'est sûrement pas un parfait rectangle tu auras des zones inutiles, mais ton raster de pentes sera déjà moins étendu.

Si tu veux absolument masquer ton raster de pentes par ta zone d'études, c'est à dire avoir un raster tout à 0 sauf à l'intérieur de ta zone d'études où tu veux tes valeurs de pente, tu peux choisir le mode de découpage "couche de masquage" (toujours dans le module d'extraction du menu raster) et tu sélectionnes la couche vectorielle limitant ta zone d'étude comme couche de masquage. Après comme tu l'as dit, il y a message d'erreur. Tu peux le contourner en faisant la manip suivante (du moins ça fonctionne chez moi) :
En bas de la fenêtre "Découper" il y a une zone de texte qui commence par "gdalwarp -q ..." Tu cliques sur le petit crayon jaune "Éditer" et la ligne de texte devient éditable. Tu changes alors le "-crop_to_cutline" par "-r cubic" et ça devrait fonctionner.
<<<

Bonjour,

Le problème peut se situer au niveau de la reprojection de la donnée que vous voulez couper...
La donnée que vous cherchez à découper doit être dans un SCR différent de celui du projet... Du coup quand vous définissez votre emprise de découpe celle-ci n'étant pas reprojeté QGIS (GDAL) ne trouve rien à découper...
Si c est le cas (SCR de la couche à découper différent de celui du projet) vous pourriez essayer d'attribuer au projet (propriétés du projet / onglet SCR) le même SCR que celui de la couche à découper... et recommencer le découpage...?

Renaud

et en passant par GRASS?

OK
je me suis attaqué au problème des SCR (cf http://qgis.spatialthoughts.com/2012/11 … shape.html)
il faut que le SCR du raster soit déjà défini et reconnu par QGIS (pas de message d alerte de QGIS demandant de spécifier le SCR de la couche, SCR non spécifié par les propriétés de la couche): pour cette étape soit faire "sauvegarder sous" le raster; soit menu raster > projection > assigner une projection.
En sortie un raster en geotiff sera généré avec un SCR correctement lu par qgis (et gdal tools)
ce raster là sera découpable par gdal sans message d erreur

Renaud

Dernière modification par Renaud_IDGEO (Tue 14 January 2014 16:54)

Bonjour,

Pour éviter la valeur à 0, tu peux essayer de décocher l'option "valeur nulle" puis spécifier que tu veux créer une bande de transparence.
En dehors de l'emprise, tu n'as aucune valeur associée, je ne sais pas si cela t'irait...

Hny29