#1 Wed 22 April 2009 16:26
- martin.baussier
- Participant actif
- Lieu: Montpellier
- Date d'inscription: 22 Sep 2007
- Messages: 69
découpage de raster par du vecteur
Bonjour,
Je cherche un outil me permettant de découper une image raster au format Tif ou ecw par un contour vecteur.
Quelqu'un aurait il une idée ?
Merci d'avance pour tous renseignement allant dans cette direction...
Cordialement,
Martin
Martin Baussier
Institut national de l'origine et de la qualité (INAO)
Chargé de mission géomatique (Service Délimitation - BDDC)
Hors ligne
#2 Wed 22 April 2009 17:26
Re: découpage de raster par du vecteur
Bonjour,
GRASS permet de faire cela en quelques étapes, voire même GDAL :
1. rasterisation du vecteur (masque)
2. utilisation du masque comme couche de découpe
Si le contour vecteur est un rectangle, GDAL permet de découper un raster en définissant ce rectangle.
Y.
Yves Jacolin, bénévole de l'association GeoRezo.net, agit au nom et pour le compte de l'association - Partageons ce qui nous départage !! - GeoRezo vous aide ? Aidez GeoRezo !
Hors ligne
#3 Wed 22 April 2009 18:55
- Nicolab
- Membre
- Lieu: Osaka
- Date d'inscription: 9 Dec 2005
- Messages: 95
Re: découpage de raster par du vecteur
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
Hors ligne
#4 Thu 23 April 2009 11:35
- martin.baussier
- Participant actif
- Lieu: Montpellier
- Date d'inscription: 22 Sep 2007
- Messages: 69
Re: découpage de raster par du vecteur
Merci à tous les deux, je pense que je vais m'orienter vers GRASS.
Martin Baussier
Institut national de l'origine et de la qualité (INAO)
Chargé de mission géomatique (Service Délimitation - BDDC)
Hors ligne
#5 Thu 23 April 2009 11:57
- Arpn'ko
- Participant actif
- Date d'inscription: 5 Sep 2005
- Messages: 54
Re: découpage de raster par du vecteur
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
Hors ligne
#6 Thu 14 January 2010 10:17
- gwen.petit
- Participant occasionnel
- Lieu: Vannes
- Date d'inscription: 7 Jul 2006
- Messages: 33
Re: découpage de raster par du vecteur
Bonjour,
La discussion date un peu, mais je me permet de proposer une alternative à Grass (qui est très bon outil par ailleurs).
Avec OrbisGIS et la fonction Cropraster vous pouvez découper votre raster très facilement. Cette fonction se base sur l'enveloppe de la couche vecteur et découpe le raster en conséquence.
- L'avantage : pas besoin de convertir la couche vecteur en raster
- L'inconvénient : pour l'instant OrbisGIS (version 2.2) ne peut ouvrir que du Tiff. L'ecw fait partie des développements en cours.
Plus d'info : www.orbisgis.org et http://brehat.ec-nantes.fr/orbisgis_dow … aster.html
Cordialement,
Gwendall Petit
Equipe OrbisGIS
Hors ligne
#8 Thu 14 January 2010 13:29
- freesnowpetrel
- Participant assidu
- Lieu: Pretoria, Afrique du Sud
- Date d'inscription: 1 Mar 2007
- Messages: 222
Re: découpage de raster par du vecteur
Bonjour,
Il y a aussi GDAL avec l'utilitaire gdalwarp qui permet de decouper un raster par le contour d'un shapefile.
Code:
gdalwarp -cutline shapefile -r cubic inputrasterfile outputrasterfile
Il est d'abord recommander de decouper le raster par le rectangle d'emprise du shapefile.
Code:
gdal_translate -projwin westbound northbound eastbound southbound inputrasterfile outputrasterfile
Avantage de la ligne de commande on peux l'utiliser en batch.
Rodolphe
Dernière modification par freesnowpetrel (Thu 14 January 2010 13:30)
Hors ligne
#10 Thu 08 April 2010 20:47
- xavier78
- Participant actif
- Date d'inscription: 12 Oct 2006
- Messages: 148
Re: découpage de raster par du vecteur
Bonjour,
le -cutline n'existe plus dans la commande de gdalwarp? Du moins dans Osgeo4W...
C:\>gdalwarp --help
Usage: gdalwarp [--help-general] [--formats]
[-s_srs srs_def] [-t_srs srs_def] [-order n] ] [-tps] [-et err_threshold]
[-te xmin ymin xmax ymax] [-tr xres yres] [-ts width height]
[-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16]
[-srcnodata "value [value...]"] [-dstnodata "value [value...]"] -dstalpha
[-r resampling_method] [-wm memory_in_mb] [-multi] [-q]
[-of format] [-co "NAME=VALUE"]* srcfile* dstfile
Available resampling methods:
near (default), bilinear, cubic, cubicspline, lanczos.
Enfin voici comment j'ai fait:
C:\>gdalwarp -cutline "F:\Travail\GEOFLADept_FR_Corse_AV_L93\LIMITE_DEPARTEM
ENT.SHP" -r cubic "F:\Travail\IR\interpolationIDW.jpeg" "F:\Travail\
IR\essai.jpeg"
Même en enlevant les guillemets cela ne marche pas...
Merci
Xav'
Dernière modification par xavier78 (Thu 08 April 2010 20:57)
Hors ligne
#11 Thu 08 April 2010 22:56
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: découpage de raster par du vecteur
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
Hors ligne
#12 Thu 08 April 2010 23:43
- xavier78
- Participant actif
- Date d'inscription: 12 Oct 2006
- Messages: 148
Re: découpage de raster par du vecteur
J'ai abandonné l'idée d'utiliser gdal dans Osgeo4w car même en dev, cela ne marche pas.
Conséquence, j'ai fait une utilisation de Fwtools et voilà le résultat:
C:\Utilitaires\FWTools2.4.7>gdalwarp -cutline F:\Travail\essai2.shp -r cubic
F:\Travail\IR\interpolationIDW.jpeg F:\Travail\IR\interpol.jpeg
Creating output file that is 1654P x 920L.
Sauf que l'image obtenu fait un 1ko et est vide
Hors ligne
#13 Thu 08 April 2010 23:52
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: découpage de raster par du vecteur
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...
Hors ligne
#14 Fri 09 April 2010 00:17
- xavier78
- Participant actif
- Date d'inscription: 12 Oct 2006
- Messages: 148
Re: découpage de raster par du vecteur
Je vais regarder tout ça!
Ceci dit le jpeg et jpegw ont été réalisé par la méthode interpolation IDW de Qgis et cela n'a pas l'air de poser pb en entrée de Gdal dans Fwtools à partir de données en Lambert93 dans un projet en Lambert 93 et en ayant vérifier le descripteur qui indique bien du lambert93.
Après vérif, effectivement cela ne supporte pas le jpeg en sortie donc je l'ai en Geotiff
C:\Utilitaires\FWTools2.4.7>gdalwarp -of Gtiff -s_srs "+init=IGNF:LAMB93" -cutl
ine F:\Travail\essai3.shp -r cubic F:\Travail\IR\interpolationIDW.jpeg
F:\Travail\IR\interpol.tiff
Creating output file that is 1654P x 920L.
Ma dernière solution serait alors de voir pour mettre un tiff en entrée.
Vu que le s_srs est renseigné, le t_srs est obsolète (du moins en ligne de commandes c'est le cas pour ogr2ogr concernant les données vecteur)
Merci pour l'aide, je tiens le fil au courant
Hors ligne
#15 Wed 04 August 2010 14:34
- scout
- Participant occasionnel
- Lieu: Bordeaux
- Date d'inscription: 17 Dec 2007
- Messages: 31
Re: découpage de raster par du vecteur
+1 Rodolphe.
Pour ceux qui n'aime pas la ligne de commande, dans qgis c'est possible également => raster puis warp
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
Hors ligne
#16 Wed 04 August 2010 15:32
- scout
- Participant occasionnel
- Lieu: Bordeaux
- Date d'inscription: 17 Dec 2007
- Messages: 31
Re: découpage de raster par du vecteur
je réponds très partiellement à ma réponse :p (et comme ça peut servir à d'autres personnes...)
J'ai réussi à avoir la fonction "raster --> warp" en allant chercher le plugin sur http://www.faunalia.it/qgis/plugins.xml car elle était introuvable sur la récupération d'extension python de QGIS, je l'ai incorporé dans un dossier Plugin, dans Program Files. Une fois chargée, la voilà à l'écran! par contre, à partir de l'option "warp" ("projection" en francais je suppose), je ne vois pas comment découper un raster...
Hors ligne
#17 Wed 04 August 2010 15:36
Re: découpage de raster par du vecteur
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)
Hors ligne