Annonce
L'appel à participation est ouvert jusqu'au 19 janvier 2025!
#1 Tue 12 April 2022 18:53
commandes GDAL pour retailler une image
Salut à toutes et à tous,
j'espère que c'est le bon forum, bien que je cherche à lancer une commande dans OSGeo4W Shell
J'ai une image tif fabriquée et géoréférencée à partir d'un pdf
Code:
gdal_translate --config GDAL_PDF_DPI 600 --config GDAL_PDF_BANDS 4 --config GDAL_CACHEMAX 1024 -co NUM_THREADS=ALL_CPUS -co COMPRESS=DEFLATE -co ZLEVEL=9 -co PREDICTOR=2 -co TILED=YES -outsize 4100 0 -a_srs EPSG:3943 -gcp 0 0 1988120.32 2237341.37 -gcp 4100 0 1988171.38 2237090.26 -gcp 0 2870 1987944.54 2237305.63 -mo CUTLINE="POLYGON((1987944.54 2237305.63,1987995.6 2237054.52,1988171.38 2237090.26,1988120.32 2237341.37,1987944.54 2237305.63))" image.pdf image.tif"
Je voudrais détourer cette image, car elle contient une marge blanche qui pose problème
je pensais même pouvoir clipper en précisant seulement la couleur, mais ça n'a pas l'air d'exister.
J'ai trouvé des pistes là:
https://trac.osgeo.org/gdal/ticket/1599
sauf que gdal n'existe pas dans OSGeo4W
Puis ça:
https://gdal.org/programs/gdalwarp.html
du coup j'ai essayé ça mais ça marche pas car il n'accepte pas de cutline autre que sous forme de fichier (super la cohérence enter les commandes)
Code:
gdalwarp -overwrite -to SRC_METHOD=NO_GEOTRANSFORM -to DST_METHOD=NO_GEOTRANSFORM -cutline CUTLINE="POLYGON((1988003.5 2237266.59,1988034.64 2237113.48,1988112.42 2237129.29,1988081.29 2237282.41))" image.tif imageclip.tif
j'aimerai mieux pouvoir eviter de créer un fichier csv à part.
le rectangle de découpe, je peux le connaitre sous forme de coordonnées lambert, mais je peux aussi l'exprimer en pixel.
si vous avez une idée de syntaxe simple,
merci ...
a+
gégé
Hors ligne
#2 Tue 12 April 2022 19:22
Re: commandes GDAL pour retailler une image
Un essai en me résignant à crer un fichier csv:
mon csv:
Code:
id,WKT 1,"POLYGON((0 0, 0 10, 10 10 , 10 0, 0 0))"
je lance:
Code:
gdalwarp -overwrite -to SRC_METHOD=NO_GEOTRANSFORM -to DST_METHOD=NO_GEOTRANSFORM -cutline "E:\ACLSP\dwg14\$Temporaire\png.csv" "E:\ACLSP\dwg14\$Temporaire\png.png" "E:\ACLSP\dwg14\$Temporaire\pngA2crop.tif"
le résultat:
Code:
Using band 4 of source image as alpha. Creating output file that is 938P x 656L.
J'arrive bien à découper un petit carré, mais il me le place en haut à gauche d'une image d'une taille 938x656,
je ne comprend pas d'où ça vient ...
Hors ligne
#3 Wed 13 April 2022 11:16
- tumasgiu
- Membre
- Lieu: Ajaccio
- Date d'inscription: 5 Jul 2010
- Messages: 1160
Re: commandes GDAL pour retailler une image
Salut,
votre tiff est probablement enregistré en TOPLEFT (origine du repère de l'image en haut à gauche).
Sinon, Qgis est capable de faire exactement ce que vous voulez et je crois que les versions récente du logiciel de retouche d'image Gimp conserve les métadonnées de géoreferencement pour le tiff !
Dernière modification par tumasgiu (Wed 13 April 2022 11:29)
Hors ligne
#4 Wed 13 April 2022 14:02
Re: commandes GDAL pour retailler une image
Bonjour,
sauf que gdal n'existe pas dans OSGeo4W
Vous avez essayé en choisissant l'option "Advanced install" ?
Dernière modification par jmarsac (Wed 13 April 2022 14:03)
Jean-Marie
Azimut
Hors ligne
#5 Wed 13 April 2022 16:25
Re: commandes GDAL pour retailler une image
Bonjour,
merci pour vos réponses,
@tumasgiu
oui, je pense que tous les géotif ont une origine haut gauche ?
ça ne me dit pas pourquoi la découpe selon un polygone agrandi la taille du canevas
Qgis le fait, c'est sur, mais bon but c'est justement de produire des images géoréférencées à partir d'Autocad, pour les utiliser sous qgis ...
@jmarsac
oui, c'est peut être la solution, ajouter des commandes gdal originale.
je vais voir ça.
Dernière modification par gegematic (Thu 14 April 2022 11:35)
Hors ligne
#6 Wed 13 April 2022 18:17
Re: commandes GDAL pour retailler une image
Salut,
en fait, a force de tâtonner et de faire des essais, j'avais perdu le crop-to-cutline
ça marche beaucoup mieux comme ça:
a écrire dans fichier cutline.csv:
Code:
id,WKT 1,"POLYGON((0 0, 0 100, 100 100 , 100 0, 0 0))"
commande a taper dans le shell 0sgeo4w:
Code:
gdalwarp -overwrite -to SRC_METHOD=NO_GEOTRANSFORM -to DST_METHOD=NO_GEOTRANSFORM -cutline "cutline.csv" -crop_to_cutline "source.tif" "dest.tif"
(j’obtiens une image 100x100 pixels depuis le coin haut gauche.)
a+
gégé
Dernière modification par gegematic (Tue 26 April 2022 13:31)
Hors ligne
