Le portail francophone de la géomatique

Bjr,

Pour compléter la réponse de Y:

Avec GRASS:
r.to.rast pour rasteriser votre vecteur
http://grass.itc.it/grass62/manuals/htm … .rast.html
r.mask pour limiter le raster à l'extent du vecteur
http://grass.itc.it/grass64/manuals/htm … .mask.html

Avec GDAL:
Découpe du raster en fonction d'un rectangle
gdal_warp -projwin xmin ymin xmax ymax
http://www.gdal.org/gdal_translate.html

Nicolab

Bonjour,

une petite précision au sujet de l'ecw.

Quand vous créez un fichier ecw, que ce soit le résultat d'une découpe ou d'un mosaicage ou autre, la donnée est compressée, c'est à dire subit un traitement de compression.

Quand le fichier de départ est déjà en ecw, le résultat de la découpe, si vous souhaitez le conserver en ecw va ainsi subir une seconde compression.

La compression ecw est une compression particulièrement adaptée aux données géographiques. Elle est néanmoins avec perte et les compressions sucessives peuvent finir par altérer la qualité.

Pour cela, il faut éviter de multiplier les compressions d'une part et quand on le fait utiliser le même taux objectif que celui utilisé lors de la compression précédente, ce qui évite au maximum la génération de bruit dans l'image.

Bref, si vous partez d'un ecw et que vous souhaitez obtenir un ecw à l'issue de votre découpe... veillez a minima à utiliser le même taux de compression que celui utilisé la première fois !

Salutations,

Agnès

Bonjour,

Il existe un outil dans l'extension sextante de gvSIG. Il s'agit de Crop grid with polygon layer.
Il faut ton raster soit un tif

Sylvain

+1 Rodolphe.
Pour ceux qui n'aime pas la ligne de commande, dans qgis c'est possible également => raster puis warp

La dernière fois que j'ai utilisé osgeo4w, le package par défaut gdal utilisait une version maintenant assez ancienne, dans la série 1.5. Par contre, il y a moyen d'installer en complément un package  - gdal-dev il me semble - qui utilise une version nettement plus récente qui doit disposer de l'option -cutline

Une raison possible est que les coordonnées de ton shapefile ne correspondent pas au système de projection de ton image source.  De plus, je trouve bizarre que ton fichier destination porte une extension .jpeg. gdalwarp ne peut pas écrire dans un fichier jpeg. Donc tu fichier interpol.jpeg doit être en fait un fichier TIFF ;-)
Je suis également surpris que tu aies un fichier .jpeg en entrée et que ça fonctionne sans paramètres -s_srs et -t_srs. Soit il y a un fichier .aux.xml à côté, soit ton fichier jpeg n'est pas un jpeg...

le sujet date un peu mais bon...            comment tu fais pour découper dans QGIS? il n'y a pas l'option "raster --> warp", enfin du moins je ne la trouve pas.

Aussi, j'ai essayé de découper un raster avec OrbisGIS. Par contre il y a un message d'erreur quand j'utilise "cropRaster" : il me dit que cette fonction n'existe pas ! faut-il installer un plugin pour que cela marche?

Enfin, j'ai essayé avec l'option de Sextante sur gvSIG, mais le résultat est en noir et blanc...

j'espère avoir des réponses! merci

Effectivement ça n'y est plus !
gdalwarp est utilisé par projection uniquement maintenant ...

Je te conseille de te reporter aux codes soulignés par Rodolphe

Edit: Le plugin a changé voilà tout. C'est pourquoi il n'y a plus de shapefile à renseigner

Dernière modification par Sylvain33 (Wed 04 August 2010 15:37)