Le portail francophone de la géomatique

Bonjour,

as-tu essayé d'enlever l'extension .py de la ligne de commande?

à+

je remonte des souvenirs (arf je n'ai pas mes notes avec moi, dsl et pas fwtools installé), il faut peut être essayer de se mettre dans un autre dossier de fwtools, il me semble que le gdal_merge.py est ailleurs ...?

Bonjour,

Je me souviens que j'ai eu le même problème il n'y a pas si longtemps.

Question=> D'où lancer vous la fonction gdal_merge ? de l'invité de commande DOS ou bien de l'invité de commande FWTOOLS ?

Pour ma part (je dois encore avoir des problèmes dans mes variables d'environnement ), la fonction marche très bien dans l'invité de commande de FWTOOLS mais pas dans l'inité de commande windows.
La syntaxe pour lancer la fonction est "gdal_merge" uniquement. En tapant uniquement ce mot, l'aide sur les paramètres de la fonction apparait (voir ci dessous).

---------------------------------------------------------------
C:\FWTools>gdal_merge
No input files selected.
Usage: gdal_merge.py [-o out_filename] [-of out_format] [-co NAME=VALUE]*
                     [-ps pixelsize_x pixelsize_y] [-separate] [-v] [-pct]
                     [-ul_lr ulx uly lrx lry] [-n nodata_value] [-init value]
                     [-ot datatype] [-createonly] input_files
                     [--help-general]

-----------------------------------------------------------------------

Voici mes variables d'environnement :

FWTOOLS_DIR: C:\FWTools
GDAL_DATA: C:\FWTools\data
GDAL_DRIVER_PATH: C:\FWTools\gdal_plugins
GEOTIFF_CSV: C:\FWTools\data
Path: C:\FWTools\bin;C:\Logiciel\ARCGIS93\Python;C:\Logiciel\ARCGIS93\Python\Lib
PROJ_LIB: C:\FWTools\proj_lib
PYTHONPATH : C:\Logiciel\ARCGIS93\Desktop\ArcGIS\Bin

! j'utilise la version de python fourni avec Arcgis et non celle de FWTOOLS.


En espérant que cela t'aide. Moi j'y avais passé pas mal de temps avant que ça fonctionne

Ahhh Merci ça marche, en effet avec la fenêtre de FWTools ça fonctionne.

Par contre nouveau problème ()) quand j'essaye d'assembler, il me dit qu'il ne peut écrire qu'en GTiff ou HFA. Or je voudrais l'assembler virtuellement, comme "HotShot" le décrit dans ce topic : http://www.forumsig.org/archive/index.php?t-6161.html (le message du 15/03/2006, 14h33)

Mais là, il me met:

"
ERROR 1: Writing through VRTSourcedRasterBand is not supported.
Traceback (most recent call last):
  File "C:\PROGRA~1\FWTOOL~1.2\bin\gdal_merge.py", line 443, in ?
    fi.copy_into( t_fh, band, band, nodata )
  File "C:\PROGRA~1\FWTOOL~1.2\bin\gdal_merge.py", line 240, in copy_into
    nodata_arg )
  File "C:\PROGRA~1\FWTOOL~1.2\bin\gdal_merge.py", line 60, in raster_copy
    data, t_xsize, t_ysize, t_band.DataType )
  File "C:\PROGRA~1\FWTOOL~1.2\pymod\gdal.py", line 869, in WriteRaster
    buf_string, buf_xsize, buf_ysize,buf_type)
TypeError: Writing through VRTSourcedRasterBand is not supported.
"

Ca te dit quelque chose ?
Sachant (je ne sais pas si ça a un lien) que je n'ai pas de dossier gdal_plugins, donc je ne sais pas définir le GDAL_DRIVER_PATH avec quel dossier. (j'ai bin, csharp, data, demo-data, html, include, lib, pics, proj_lib, pymod, python, ramps, tools, vb6, wms, xmlconfig)

salut,

le message d'erreur me fais penser à un plugin absent ...

étant un utilisateur de linux j'aurai un peu de mal à te donner la démarche exacte, mais je te conseillerais de regarder du coté de OSGeo4W (http://trac.osgeo.org/osgeo4w/).
Il s'agit d'une interface graphique windows qui permet d'installer facilement tou un tas d'outils sig libres initialement concus pour linux (contient entre autre GDAL/OGR, GRASS, MapServer, OpenEV, uDig, QGIS ), en gérant les dépendances et la config ... et qui propose normalement tous les plugins disponibles !

Ahhhhhh

Ca suffisait en fait, ce petit .bat. Et du coup j'ai essayé, et ça fonctionne. Donc un gros merci !


Bon dernier problème (promis !).
J'assemble virutellement mes photos, je reprojectionne virtuellement mes photos, je sors ma photo en ecw. L'image est correctement calée, et je n'ai plus les bords qui font tache.
MAIS, je ne saisis pas trop comment redécouper proprement...

(toujours par HotShot sur le lien)
"
2/ extraction de dalles reprojetées en lat/lon(wgs84) grâce à gdal_translate, selon un pas constant (à déterminer... en fonction de la résolution voulue des textures et leur dimension LxH). Par exemple au moyen d'un petit fichier batch qui appelle de façon répétée gdal_translate en faisant varier à chaque fois les Xmin/max et Ymin/max de la zone à extraire.
gdal_translate -of [format] -projwin [Xmin Ymax Xmax Ymin, ce sont les params qu'on fait varier par fichier batch] Reproj.vrt [Fichier.xxx, le fichier de sortie, avec un nom différent pour chaque dalle extraite]
"

OK il faut mettre les coordonnées des images que l'on veut sortir. Mais comment savoir quelles sont les coordonnées à mettre ?
Je ne trouve pas ces infos dans le fichier .vrt par exemple...
Je ne sais pas, il faut par exemple convertir les coordonnées (avec Circé par exemple ?) de la photo d'origine de chaque coin (haut gauche et bas droit) vers le nouveau système ?
Il y a une solution plus rapide ?

Dernière modification par Sam_Dalembert (Mon 10 August 2009 16:29)