#1 Mon 24 September 2012 11:28
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
GDAL Coordonées absente fichier ECW
Bonjour et merci d'avance si vous pouvez m'aider.
Voila mon problème je souhaite extraire une dalle ecw d'un fichier vrt issu d'un assemblage de dalles IGN en Lamb93 que j'ai reprojeter en LambE.
j'utilise la commande suivante pour extraire ma dalle.
gdal_translate -of ECW -co TARGET=90 -co LARGE_OK=YES -co DATUM=NTF -co PROJ=LAMBE -a_srs EPSG:27572 -a_nodata 255 assemblage_LAMBE.vrt -projwin 847996 2465180 899707 2416140 Ortho_2009_VdL_LAMBE.ecw
je n'ai pas eu de problème sinon de mémoire mais que j'ai réglé.
j'obtiens un fichier de 3 Go en sortie.
Là ou j'ai un souci c'est que je n'ai pas de coordonnée affiché aux coins de mon image.
Lorsque que je fait un teste sur une surface plus petite mon fichier est nickel.
autre question je pensai que le bord noir de mon image serait devenu blanc et ce n'est pas le cas !
Voici le résultat de gdalinfo:
Driver: ECW/ERMapper Compressed Wavelets
Files: Ortho_2009_VdL_LAMBE.ecw
Ortho_2009_VdL_LAMBE.ecw.aux.xml
Size is 172370, 163467
Coordinate System is:
PROJCS["NTF (Paris) / Lambert zone II",
GEOGCS["NTF (Paris)",
DATUM["Nouvelle_Triangulation_Francaise_Paris",
SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
AUTHORITY["EPSG","7011"]],
TOWGS84[-168,-60,320,0,0,0,0],
AUTHORITY["EPSG","6807"]],
PRIMEM["Paris",2.33722917,
AUTHORITY["EPSG","8903"]],
UNIT["grad",0.01570796326794897,
AUTHORITY["EPSG","9105"]],
AUTHORITY["EPSG","4807"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Lambert_Conformal_Conic_1SP"],
PARAMETER["latitude_of_origin",52],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",0.99987742],
PARAMETER["false_easting",600000],
PARAMETER["false_northing",2200000],
AUTHORITY["EPSG","27572"],
AXIS["X",EAST],
AXIS["Y",NORTH]]
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0,163467.0)
Upper Right (172370.0, 0.0)
Lower Right (172370.0,163467.0)
Center (86185.0,81733.5)
Band 1 Block=172370x1 Type=Byte, ColorInterp=Red
NoData Value=255
Overviews: arbitrary
Band 2 Block=172370x1 Type=Byte, ColorInterp=Green
NoData Value=255
Overviews: arbitrary
Band 3 Block=172370x1 Type=Byte, ColorInterp=Blue
NoData Value=255
Overviews: arbitrary
Cordialement
Hors ligne
#2 Mon 24 September 2012 12:36
- Nicolas Granier
- Participant assidu
- Date d'inscription: 19 Apr 2007
- Messages: 271
Re: GDAL Coordonées absente fichier ECW
Salut,
-co PROJ=LAMBE
Sauf erreur de ma part LAMBE est le code de projection du registre IGNF, et ici il faut le code de projection d'ER-MAPPER qui est si je me souviens bien "LM2FRANC".
Essayez déjà en remplaçant ce code.
A+
Hors ligne
#3 Mon 24 September 2012 13:30
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
Re: GDAL Coordonées absente fichier ECW
Ok
je relance l’opération.
Merci
Hors ligne
#4 Mon 24 September 2012 17:32
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: GDAL Coordonées absente fichier ECW
Avec http://trac.osgeo.org/gdal/browser/trun … al_edit.py et GDAL >= 1.9, il n'y a pas besoin de relancer tout le processus de création pour modifier le géoréférencement. On peut modifier sur place pour un certain nombre de formats, en particulier l'ECW :
python gdal_edit.py ton.ecw -mo PROJ=LM2FRANC
et éventuellement -a_ullr avec les valeurs qui vont bien
Hors ligne
#5 Mon 24 September 2012 18:29
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
Re: GDAL Coordonées absente fichier ECW
Cool !
Merci
C'est exactement ce que je cherchai
J’essaie demain
Hors ligne
#6 Wed 26 September 2012 11:25
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
Re: GDAL Coordonées absente fichier ECW
Salut,
La manip avec le code de projection Er-mapper a marché.
Merci
Mais il faudrait que je fasse la même chose sur un fichier de 10 Go c'est pourquoi ta manip Rouault m'intéresse grandement.
je viens d'essayer ton patch (je ne sais pas si c'est comme ça que cela s’appel ?)
Voici le message que j'ai :
Z:\couches_sig\donnees_de_reference\orthophotographie\2009>python gdal_edit.py O
rtho_2009_SCOT_LAMBE.ecw -mo PROJ=LM2FRANC
Traceback (most recent call last):
File "gdal_edit.py", line 32, in ?
from osgeo import gdal
ImportError: No module named osgeo
J'utilise Fwtools 2.4.7
As-tu une idée de ce qui ne marche pas.
Merci d'avance
Hors ligne
#7 Wed 26 September 2012 14:31
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: GDAL Coordonées absente fichier ECW
Even a bien précisé :
GDAL >= 1.9
.
FwTools 1.7.2 est obsolète depuis plusieurs années. Il faut installer GDAL 1.9.1 .
Hors ligne
#8 Wed 26 September 2012 16:53
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
Re: GDAL Coordonées absente fichier ECW
Ok !
Mais faut savoir que Fwtools 2.4.7 est compilé avec un Gdal 1.7.2 !
Pas facile a trouver de premier abord.
Bon j'ai installé
python 3.2
GDAL 1.9.1
Les 3 Plugin ECW MrSID et Oracle
Et ...
j'ai toujour le meme message
Z:\couches_sig\donnees_de_reference\orthophotographie\2009>gdal_edit.py Ortho_20
09_SCOT_LAMBE.ecw -mo PROJ=LM2FRANC
Traceback (most recent call last):
File "Z:\couches_sig\donnees_de_reference\orthophotographie\2009\gdal_edit.py"
, line 32, in <module>
from osgeo import gdal
ImportError: No module named osgeo
???
Hors ligne
#9 Wed 26 September 2012 17:03
Re: GDAL Coordonées absente fichier ECW
Bonjour,
Ce script python doit se baser sur une version pas trop ancienne (genre gdal 1.8) par exemple car il y a eut des évolutions dans la manière d'importer les lib python.
Et ce message d'erreur m'y fait penser :
from osgeo import gdal
ImportError: No module named osgeo
Es tu sur sur que python-gdal a été mis à jour ?
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
#10 Mon 01 October 2012 12:27
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
Re: GDAL Coordonées absente fichier ECW
Bonjour,
Je ne sais pas si les mise a jour sont bien prise en compte.
Pour plus de sureté j'ai désinstallé Fwtool.
mais comment savoir si la liaison Python Gdal fonctionne bien ?
Je suis sous windows Xp.
Merci
Hors ligne
#11 Mon 01 October 2012 12:53
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: GDAL Coordonées absente fichier ECW
Essaie de taper directement la commande d'importation dans la console Python :
Code:
$ python Python 2.7.3 (default, Sep 25 2012, 12:34:17) [GCC 4.5.4] on linux2 Type "help", "copyright", "credits" or "license" for more information. >>> from osgeo import gdal >>>
S'il n'y a pas de message d'erreur, comme ci-dessus, c'est OK.
Dernière modification par jcr83 (Mon 01 October 2012 12:53)
Hors ligne
#12 Tue 02 October 2012 16:10
- Lathyrus
- Participant occasionnel
- Date d'inscription: 23 Jul 2010
- Messages: 14
Re: GDAL Coordonées absente fichier ECW
Merci à vous tous qui m'avez aidé.
la commande gdal_edit.py est génial !
J'ai réglé mon problème de liaison Python Gdal par une réinstalation des différents composants.
1-Python
2-Gdal core
3-Gdal 1.9.1 Py
Peut être que je n'avais pas respecter l'ordre.
Merci encore
Lathyrus
Hors ligne