Annonce
Pour sécuriser votre compte sur les forums du GeoRezo, nous demandons de changer votre mot de passe.
Vous allez recevoir un message pour effectuer ce changement de mot de passe.
Merci de bien respecter les règles préconisées.
#1 Fri 19 February 2010 15:49
- Sam_Dalembert
- Participant actif
- Lieu: Mérignac
- Date d'inscription: 5 Aug 2009
- Messages: 121
[GDAL] Problème gdal_translate en ECW avec wktext
Bonjour à tous,
Je viens vous faire part d'un problème assez bizarre. C'est assez simple :
gdalwarp -of VRT -s_srs "+init=IGNF:LAMB2 +wktext" -t_srs "+init=IGNF:LAMB93" fichier_de_depart.ecw fichier_virtuel_reproj.vrt
gdal_translate -of ECW -co TARGET=90 -co LARGE_OK=YES -co PROJ=LMFRAN93 -co DATUM=RGF93 -projwin X Y X Y fichier_virtuel_reproj.vrt Extraction.ecw
Le fichier VRT est bien créé, mais gdal_translate n'arrive pas à créer le fichier ECW, ou plutôt si, il le crée, mais il reste à 1ko, dans la fenêtre Shell ça ne dépassé pas 0, et le processus gdal_translate ne dépasse pas les 15 000 ko (alors que ça monte normalement à 200 000, 500 000, ou même bien plus).
J'ai essayé d'enlever les options PROJ et DATUM, ça ne marche pas, j'ai essayé d'enlever le +wktext, ça fonctionne, mais le résultat est complètement faux.
Par contre, quand j'essaye de sortir le fichier en TIFF au lieu de le sortir en ECW, ça marche. Et ensuite la conversion en ECW fonctionne également.
Quelqu'un a une idée ? 
Hors ligne
#2 Fri 19 February 2010 16:09
- sylpingus
- Participant occasionnel
- Lieu: Aix en Provence
- Date d'inscription: 9 Jan 2006
- Messages: 34
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Bonjour,
je ne suis pas spécialiste dans le domaine mais j'ai déjà rencontré le problème sous MS Windows en utilisant FWTools. La même manipulation sur Linux compilé avait bien fonctionné. Après, de là à vous dire pourquoi... (mauvaise gestion du VRT sous Windows ?)
Par rapport à l'option +wktext, avez-vous la grille IGNF qui est installée correctement et qui est bien prise en compte ? (cf http://lambert93.ign.fr/index.php?id=28) ? Si le résultat est faux en enlevant ce paramètre, c'est tout à fait normal étant donné qu'il sert à prendre la grille de transformation officielle pour passer du RGF93 au NTF.
Bonne journée
Hors ligne
#3 Mon 22 February 2010 08:47
- Sam_Dalembert
- Participant actif
- Lieu: Mérignac
- Date d'inscription: 5 Aug 2009
- Messages: 121
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Oui désolé, je n'ai pas précisé, mais oui les résultats sont bons (j'ai la grille IGNF) avec wktext.
Je vais essayer de bidouiller le VRT, mais bon...
Hors ligne
#4 Mon 22 February 2010 23:18
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 166
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Hum, s'il s'agit d'un tout dernier FWTools 2.4.7 qui utilise proj 4.7.0, alors je pense qu'il pourrait s'agir d'une manifestation du problème suivant :
http://trac.osgeo.org/proj/ticket/63
Tu devrais essayer avec une version antérieure pour voir ?
Hors ligne
#5 Wed 24 February 2010 09:07
- Sam_Dalembert
- Participant actif
- Lieu: Mérignac
- Date d'inscription: 5 Aug 2009
- Messages: 121
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Désolé pour le temps de réponse, j'attendais d'avoir un accès à un ordi où les anciennes versions étaient installées.
C'est toi qui a posté le ticket sur trac.osgeo.org/proj/ ? 
En tous cas merci beaucoup et bien vu 
C'est exactement ça, j'ai testé avec FWTools 2.4.6, ça ne marchait pas, par contre, avec la version 2.4.2, ça fonctionne. J'ai pris le proj.dll de cette version, je l'ai mis à la place du proj.dll de la version 2.4.7, et ça fonctionne.
Par contre, si ce n'est pas trop demander, tu pourrais m'expliquer ce qu'il faut faire pour résoudre le problème ? Apparemment il y a un patch, mais je n'y comprends rien...:s
Ou est-ce que remplacer le proj.dll suffit ? Ca change quelque chose au niveau de la précision ou pas ?
Je vais continuer les vérifications, mais entre l'ECW converti à partir du TIFF (qui marche avec wktext) et l'ECW que je viens de créer (en remplaçant le dll), je ne vois pas de différence pour l'instant alors que sans le wtext, on voit une différence).
Dernière modification par Sam_Dalembert (Wed 24 February 2010 09:10)
Hors ligne
#6 Wed 24 February 2010 09:36
- Sam_Dalembert
- Participant actif
- Lieu: Mérignac
- Date d'inscription: 5 Aug 2009
- Messages: 121
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Par contre, y a encore un souci là. Les tests marchent avec le Lambert 93, mais quand j'essaye avec une zone CC, il ne reconnait tout simplement le code de la zone !!! (ex: LAMBCC48).
Je viens de tomber sur ce ticket http://trac.osgeo.org/proj/changeset/1440
Ca veut dire qu'on peut rajouter manuellement dans le fichier IGNF ? J'ai l'impression que c'est ça mais bon, je n'y connais pas grand chose donc je préfèrerais avoir l'avis de quelqu'un qui s'y connait justement...
EDIT: bizarre, les codes apparemment ne sont pas sous la forme LAMBCCXX mais RGF93CCXX
http://lambert93.ign.fr/index.php?id=28
Je vais tester si les résultats semblent bons ou non.
Mais dans ce cas, à quel problème correspond le ticket dont je parle ?
Dernière modification par Sam_Dalembert (Wed 24 February 2010 09:52)
Hors ligne
#7 Wed 24 February 2010 13:13
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 166
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Oui, c'est moi qui avait créé le ticket http://trac.osgeo.org/proj/ticket/63. Ca date d'environ 1 mois.
Plusieurs solutions :
1) Soit tu sais te débrouiller avec un compilateur et tu appliques le patch à proj 4.7.0, ou alors tu prends un snapshot du HEAD du subversion de proj.
2) Tu attends la prochaine release de proj. (Je n'en connais pas le planning)
3) Tu utilises effectivement une version antérieure de proj.dll (4.6.1 par ex). J'ai été un peu surpris que fwtools 2.4.6 ait aussi le problème mais c'est tout à fait possible qu'elle ait déjà proj 4.7.0. Faudrait que tu vérifies en lancant cs2cs et voir le numéro de version qu'il te dit (s'il te dit 4.7.1, faut comprendre 4.7.0 : il y a eu un petit typo apparemment à la génération de 4.7.0). En tout cas utiliser une version antérieure ne devrait pas poser de problème car l'interface de proj est stable (évidemment si tu utilises des fonctionnalités avec des anomalies corrigées dans des versions ultérieures, tu peux avoir des soucis...)
Quant au changeset http://trac.osgeo.org/proj/changeset/1440, il date de 2 ans, donc je pense qu'il ne faut pas se prendre la tête avec du moment que tu n'utilises pas de version proj trop vieille. Par contre, celui qui est intéressant, c'est celui qui suit http://trac.osgeo.org/proj/changeset/1444 où effectivement on voit que les codes LAMBCCXX ont été renommés en RGF93CCXX
Hors ligne
#8 Wed 24 February 2010 13:31
- Sam_Dalembert
- Participant actif
- Lieu: Mérignac
- Date d'inscription: 5 Aug 2009
- Messages: 121
Re: [GDAL] Problème gdal_translate en ECW avec wktext
Merci beaucoup pour ta réponse !
Avec le proj.dll de la version 2.4.7: Rel. 4.7.1, 23 September 2009
Avec le proj.dll de la version 2.4.2: Rel. 4.6.1, 21 August 2008
Par contre, je dois t'avouer qu'appliquer le patch, je ne sais pas du tout comment faire. Et "le HEAD du subversion de proj" c'est du chinois pour moi 
Concrètement, je reprojette des photos aériennes (en Lambert II Centre) en CC48 et en RGF93, je fais des assemblages, des extractions, mais ça s'arrête là. Donc utiliser la version 4.6.1 de proj ne changera rien non ?
Et bon, en faisant des tests visuels (je superpose le cadastre en vecteur), le résultat est impeccable (entre l'image d'origine en L2 et la transfo en CC48 ou L93). Et avec des tests de coordonnées, c'est pareil (je calcule avec Circé, puis avec gdaltransform, et les résultats sont identiques, contrairement aux codes EPSG où il y a une différence avec Circé).
Donc je suppose que c'est bon...
Dernière modification par Sam_Dalembert (Wed 24 February 2010 13:31)
Hors ligne
