Le portail francophone de la géomatique

Bonjour,

Tout d'abord, il ne suffit pas de modifier les coordonnées du point de calage indiqué dans l'entête du fichier ECW, comme cela vous a été dit. Il faut recalculer l'image entière.

Après, pour bosser uniquement au format ECW, la réponse précédente est tout a fait exacte. Par contre, si votre image n'est pas trop importante, vous avez la possibilité d'utiliser GDAL. Auquel cas, il faudra passer par un géoTiff qui ne peut excéder 4 Go (d'ou la problématique de taille de votre image), vous pourrez alors reprojeter votre image puis la reconvertir si vous le souhaitez en ECW ou en JPEG2000.

Pour se faire, je vous conseille d'utiliser l'ensemble d'outils FWtools dont une nouvelle version vient d'être compilée :
Pour windows : http://home.gdal.org/fwtools/FWTools131.exe
Pour linux : http://home.gdal.org/fwtools/FWTools-linux-1.3.1.tar.gz

La premiere étape sera de convertir votre ECW en Géotiff avec gdal_translate
-> gdal_translate -of GTIFF monfichier.ecw monfichier.tif
Puis de reprojeter votre image avec gdalwarp (utilisez les codes EPSG pour les projections)
-> gdalwarp -s_srs EPSG:27572 -t_srs EPSG:???? -of GTiff -rc monfichier.tif monfichier_reprojete.tif (a vous de chercher le bon code EPSG pour le lambert 3 sud)
Enfin de convertir en ECW ou JPEG2000 votre nouvelle image reprojeté avec gdal_translate
-> gdal_translate -of ECW monfichier_reprojete.tif monfichier_reprojete.ecw
-> gdal_translate -of JP2ECW monfichier_reprojete.tif monfichier_reprojete.jp2

Géomatiquement votre
Ludovic

Merci à tous pour ces réponses précieuses.

Cependant, j'ai un nouveau problème lorsque je tape la première ligne de commande :
gdal_translate -of GTIFF alea_induit.ecw alea_induit.tif

un message d'erreur apparaît et m'indique qu'il manque la dll sde91.
J'ai réinstallé FW Tools, sans succès...

Que dois-je faire ?

Gaëlle

oui, le .exe que j'ai téléchargé s'intitule "FWTools131.exe".

Au passage, une autre question (la dernière ?) : après le message d'erreur, le premier .tif a quand même été crée. Pour les étapes suivantes, le message s'est à nouveau affiché, mais les fichiers "nom_reprojete.tif" puis "nom_reprojete.ecw" ont également été générés.

Toutefois, comment puis-je me servir de cet ECW ? j'en ai besoin dans MapInfo. Si je ouvre le raster seul, MapInfo me demande si je souhaite caler l'image ! Où sont les coordonnées du raster ? dans un .twf ?

Je débute et m'y perds un peu...

D'avance merci

Gaëlle

Bonsoir,

j'ai retrouvé et placé à l'endroit approprié les dll manquantes. Donc, l'opération fonctionne parfaitement.
Ma dernière question reste la même : comment puis-je exploiter mon ecw, dans MapInfo par exemple ?

Encore merci

Gaëlle

Bonjour,
je souhaites reprojeter un ecw de NTF Paris France III à Lambert zone III
je ne connais pas du tout Gdal que je pense devoir utiliser
quelqu'un pourrait-il m'expliquer rapidement comment je dois m'y prendre ?

ce sont des "commandes à rentrer directement dans FWTools Shell ?

-> gdal_translate -of GTIFF monfichier.ecw monfichier.tif
-> gdalwarp -s_srs EPSG:CODE EPSG DU ECW INITIAL ??? -t_srs EPSG:???? -of GTiff -rc monfichier.tif monfichier_reprojete.tif

J'espère que quelqu'un m'aidera à me sortir de cette galère !

Merci
Yannick

Bonjour,
C'est ce que je pensais au tout début, je pensais que ça pouvait être du Lambert carto mais non !

Pas mal de paramètres (Arccatalog, propriétés) sont différents; par exemple :
unité angulaire (degrés pour LzIII, grad pour l'autre)
False Northing (3200000 au lieu de 200000 pour LzIII)
standard parallel 1 : 49 au lieu de 43.199
Lat of origin : 49 au lieu de 44
...
Et puis surtout, dans arcMap mon ecw en "NTF_Paris_France_III" ne se cale pas du tout sur ma zone où j'ai déjà des couches en Lambert zone III !!!

Si je ne me suis pas trompé le EPSG du "NTF_Paris_France_III" est le 27583
celui du Lambert zone III : 27573

Voilà où j'en suis.

Je ne sais pas du tout comment faire fonctionner Gdal.

HELP !!!

Dernière modification par yannick12 (Tue 12 June 2007 11:49)

Bonjour,

Il se peut que la dalle soit mal géoréférencé. Je dis ca sans aucune certitude. Mais le code 27573 et 27583 sont être les mêmes. En tout cas les codes 27572 et 57582 sont identique, et il est d'ailleurs recommandé de ne plus utiliser le code 275282.

Tu peux essayer de faire :

 gdalwarp -s_srs "EPSG:27583" -t_srs "EPSG:27573" -of GTiff -rc monfichier.tif monfichier_reprojete.tif

Y.

Re-bonjour,

j'ai repris mon problème à zero en "effaçant" la définition du système de coordonnées dans les propriétés du ecw; J'ai ensuite utilisé une couche en Lambert zone III à laquelle j'ai ajouté mon ecw.
Et là miracle !

En fait le système défini dans les propriétés n'était tout simplement pas le bon !
Je ne sais pas pourquoi. est-il possible qu'il y ait eu un problème du fait d'appelation différente entre arcgis 3.2 (utilisé par le bureau qui m'a fourni le ecw) et arcgis9.0 (celui que j'utilise) ?

Mon problème est donc résolu, le ecw était en fait bien en Lambert zone III !

Désolé de vous avoir dérangé pour rien, mais merci pour vos réponses !

PS : maintenant je ne ferai plus confiance aux propriétés des infos que l'on me fourniera et je revérifierai avant de me lancer dans des opérations fastidieuses !!! 

Dernière modification par yannick12 (Tue 12 June 2007 13:47)