banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

#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

Yves
Membre du bureau
Lieu: Aix-les-Bains
Date d'inscription: 22 Mar 2006
Messages: 9869
Site web

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

 

#7 Thu 14 January 2010 12:29

Sylvain33
Participant assidu
Lieu: Mauritius
Date d'inscription: 13 Feb 2007
Messages: 320
Site web

Re: découpage de raster par du vecteur

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

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

 

#9 Thu 14 January 2010 13:42

Sylvain33
Participant assidu
Lieu: Mauritius
Date d'inscription: 13 Feb 2007
Messages: 320
Site web

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

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

Sylvain33 a écrit:

+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

Sylvain33
Participant assidu
Lieu: Mauritius
Date d'inscription: 13 Feb 2007
Messages: 320
Site web

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

 

Pied de page des forums

Powered by FluxBB