Annonce
#1 Thu 13 January 2011 17:53
[GDAL] du TIF à l'ECW: l'étendue n'est plus retrouvée...
Bonjour,
J'utilise FWTOOLS 2.4.7
Lorsque je convertis une couche TIF en ECW,
Code:
gdal_translate -of ECW image.tif image.ecw
l'étendue de l'ECW devient la suivante avec un départ des coord. à x=0,y=0:
Code:
gdalinfo image.ecw
Code:
Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 9988.0) Upper Right (12589.0, 0.0) Lower Right (12589.0, 9988.0) Center ( 6294.5, 4994.0)
La commande que je lance est la suivante:
gdal_translate -of ECW image.tif image.ecw
Code:
gdalinfo image.tif
Code:
Upper Left ( 1217693.926, 6204435.653) ( 9d19'22.91"E, 42d45'10.18"N) Lower Left ( 1217693.926, 6199441.653) ( 9d19'5.37"E, 42d42'29.05"N) Upper Right ( 1223988.426, 6204435.653) ( 9d23'58.47"E, 42d44'53.78"N) Lower Right ( 1223988.426, 6199441.653) ( 9d23'40.73"E, 42d42'12.66"N) Center ( 1220841.176, 6201938.653) ( 9d21'31.87"E, 42d43'41.44"N)
Le code complet de gdalinfo image.tif est:
Code:
Driver: GTiff/GeoTIFF
Files: temp/planche1_georef.tif
Size is 12589, 9988
Coordinate System is:
PROJCS["unnamed",
GEOGCS["unnamed",
DATUM["GRS_80",
SPHEROID["GRS 80",6378137,298.257222101]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",44],
PARAMETER["standard_parallel_2",49],
PARAMETER["latitude_of_origin",46.5],
PARAMETER["central_meridian",3],
PARAMETER["false_easting",700000],
PARAMETER["false_northing",6600000],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (1217693.925922197800000,6204435.653465110800000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 1217693.926, 6204435.653) ( 9d19'22.91"E, 42d45'10.18"N)
Lower Left ( 1217693.926, 6199441.653) ( 9d19'5.37"E, 42d42'29.05"N)
Upper Right ( 1223988.426, 6204435.653) ( 9d23'58.47"E, 42d44'53.78"N)
Lower Right ( 1223988.426, 6199441.653) ( 9d23'40.73"E, 42d42'12.66"N)
Center ( 1220841.176, 6201938.653) ( 9d21'31.87"E, 42d43'41.44"N)
Band 1 Block=12589x1 Type=Byte, ColorInterp=Red
Band 2 Block=12589x1 Type=Byte, ColorInterp=Green
Band 3 Block=12589x1 Type=Byte, ColorInterp=BlueQqn aurait une idée du problème?
Dernière modification par baobazz (Thu 13 January 2011 17:53)
geodata au cerema et petits billets en géomatique
Hors ligne
#2 Thu 13 January 2011 18:04
- jcr83
- Participant actif
- Date d'inscription: 17 Mar 2008
- Messages: 147
Re: [GDAL] du TIF à l'ECW: l'étendue n'est plus retrouvée...
Bonjour,
Il faut indiquer au driver ECW le nom du datum et le nom de la projection. Comme c'est du Lambert93, essayez:
Code:
gdal_translate -of ECW -co DATUM=RGF93 -co PROJ=LMFRAN93 image.tif image.ecw
Il vous manque aussi le paramètre TARGET, qui permet de fixer le taux de compression. En général, il faut le mettre autour de 90% ou 95%. Par défaut, il est de 75%.
Lisez la page http://www.gdal.org/frmt_ecw.html pour plus d'infos.
Hors ligne
#4 Thu 15 September 2011 16:10
- Guillaume D
- Participant occasionnel
- Lieu: Dijon
- Date d'inscription: 11 Jan 2007
- Messages: 14
- Site web
Re: [GDAL] du TIF à l'ECW: l'étendue n'est plus retrouvée...
J'ai été confronté au même problème lorsque je mettais :
Code:
for %f in (*.tif) do gdal_translate -of ECW %f %f.ecw
la conversion s'effectuait mais les cordonnées des .ecw n'étaient pas bonnes.
en rajoutant les options de projection, cela marche nickel.
Code:
for %f in (*.tif) do gdal_translate -of ECW -co "TARGET=90" [b]-co DATUM=RGF93 -co PROJ=LMFRAN93[/b] %f %f.ecw
Encore merci pour la qualité de ce forum.
Hors ligne
