Le portail francophone de la géomatique

Petite question : lorsque que l'on veut reprojeter un raster possédant une table couleur, avec un ré échantillonnage autre que nearest, GDAL ne veut pas, il dit qu'il faut convertir l'image en 24Bit (et je veux le résultat en ECW).

Apparemment il y a deux solutions :
- soit on convertit avec pc2rgb le tif vers un tif en 24Bits, puis on reprojette en VRT et on sort en ECW
- soit on covertit le tif en ECW (avec -expand RGB pour garder la couleur), puis on reprojette en VRT et on sort en ECW

Les deux méthodes semblent fonctionner, aussi bien au niveau de la reprojection que de la qualité. J'ai juste une différence de taille de 5% (mais bon j'ai fait seulement un test).

Donc ma question, est-ce qu'il y a une méthode plus valable que l'autre ?
Sans tenir compte du temps des traitements (probablement plus rapide pour la première méthode que pour la seconde).

Dernière modification par Sam_Dalembert (Fri 16 April 2010 09:29)

C'est marqué
IAViewer, ou InterAtlas Viewer. Va voir sur leur site, je suis plus très sûr mais je crois qu'il est téléchargeable sur leur site.

Version 3.3 : http://up.sur-la-toile.com/iB2m

Ajout de 2 erreurs (#20 et #21).

A propos de #20: ok merci pour la précision. En fait le problème de couleurs -auquel expand RGB répond- je l'ai toujours rencontré en tentant de faire de l'ECW, donc j'avais jamais réfléchi au fait qu'on pouvait utiliser cette option avec du TIFF, je pensais que c'était un problème propre à l'ECW.

A propos de #21: oui oui, c'est marqué "Dans le cas où l’on essaye de créer un TIFF de plus de 4Go en utilisant une méthode de
compression (empêchant GDAL de calculer le poids final du fichier créé)"


Sinon oui, je sais que tu n'es pas fan de l'ECW, mais quel autre format aussi performant peut-on utiliser dans ce cas ? Parce que le JPEG2000, c'est bien bien compliqué pour quelqu'un qui ne s'y connait pas ou qui n'a pas le temps d'apprendre là-dessus...

Entre le      JPEG2000     (mais nécessite libjasper ? bien compliqué à installer), le JP2ECW    (mais même souci que l'ECW ou pas ? en tous cas nécessite le codec SDK), le JP2KAK (nécessite la Kakadu library), le JP2MrSID (nécessite le MrSID SDK) et le JP2OpenJPEG (nécessite la  OpenJPEG library), c'est bien bien compliqué de savoir ce que l'on peut faire avec chaque format et comment on peut arriver à utiliser chaque format...

Dernière modification par Sam_Dalembert (Mon 07 June 2010 10:40)

Pour #21, j'avais lu un peu trop vite ton doc effectivement !

Concernant le format JPEG2000, l'avantage par rapport à ECW est qu'il s'agit d'un format aux spécifications ouvertes (bien que probablement couvert par des brevets, mais c'est une autre histoire) et dont il existe différentes implémentations, libres (Jasper ou OpenJPEG) ou proprio ( ECW SDK, MrSID SDK, Kakadu), donc on est déjà dans une situation moins fermée qu'un fournisseur unique. GDAL n'est pas raciste puisqu'il permet d'offrir le support de JPEG2000 via une des 5 librairies sur lesquelles s'appuient les drivers que tu as cités. Il s'agit de 5 implémentations du même format, pas de 5 formats différents. Laquelle choisir ? C'est une histoire de goûts, sous, licence, contraintes ... Kakadu (non libre & payant) est pas mal réputée au niveau des perfos (jamais testé moi même). JP2ECW s'appuie sur le SDK ECW (donc si tu as le support de l'ECW dans GDAL, tu as aussi JP2ECW). Clairement les 2 implémentations libres ne sont pas au niveau que certaines implémentations commerciales qui bénéficient d'optimisations poussées via utilisation d'assembleur (instructions SSE, etc..) par ex. OpenJPEG est capable de lire de grosses images, à la différence de Jasper qui doit décompresser en mémoire toute l'image (impraticable sur de gros rasters). Par contre OpenJPEG est assez "expérimental" (le driver s'appuie sur une branche v2 qui n'a pas encore eu de version officielle), à la différence de Jasper qui lui ne connaît plus aucun développement apparent depuis plusieurs années.

Sur les bénéfices comparés ECW / JP2, quelques commentaires intéressants ici : http://georezo.net/forum/viewtopic.php?id=58391

Even tant que tu es dans le coin Je suis face à un autre type d'erreur, avec très très peu d'échos sur le net...

0...10...ERROR 2: CPLMalloc(): Out of memory allocating 148884 bytes.

This application has requested the Runtime to terminate it in an unusual way.
Please contact the application's support team for more information.

Cette erreur arrive quand j'essaye de sortir en ECW un assemblage d'ECW (donc plusieurs ECW assemblés avec gdalbuildvrt). A noter qu'hier, j'ai essayé à partir du même VRT, mais avec l'option GDAL_ECW_CACHE_MAXMEM, et là, il n'y a plus d'erreur, mais une fois arrivée à 0...10..., ça bloque. Je ne pense pas que ce soit une coïncidence...
Là je suis sur une machine à 2Go de RAM.
Je re-testerai quand je pourrai, mais il me semble que la semaine dernière sur une machine à 3Go, je n'ai pas eu cette erreur (mais j'ai eu l'erreur classique ERROR 1: IReadBlock failed at X offset 0, Y offset 2411, au bout de 40% je crois, je n'avais pas utilisé l'option GDAL_ECW_CACHE_MAXMEM).

Une idée du souci ?

Dernière modification par Sam_Dalembert (Tue 08 June 2010 10:21)

Arf ok merci.
A noter que le fichier est passé en TIFF, 92Go Illisible j'aurais pensé, mais non, QGIS arrive à le lire (en étant un peu patient), fort QGIS !

Sinon version 3.5: http://up.sur-la-toile.com/iB8w
pas de modif majeure, juste quelques précisions et corrections orthographiques

Je disais ça dans le sens où ce BIGTIFF (en effet ), je n'ai réussi à l'ouvrir qu'avec QGIS, les autres logiciels avec lesquels j'ai essayé n'ont pas réussi.

Sinon, si j'ai compris ta dernière phrase, je trouve que ça ferait doublon, la traduction de la doc GDAL d'Yves est très complète.
Mon tutoriel est très orienté reprojection raster, c'est un genre de document moins apte à être dans un wiki je trouve.
Et enfin il faudrait quelqu'un pour s'occuper de tout intégrer dans le wiki, ce qui serait très long je pense (remettre toutes les balises, les captures, etc.).

Ben c'est au même endroit...Tu vas sur le site de GDAL, tu cliques en haut sur Français, et là tu arrives sur les liens vers le tutoriel et vers la doc.

Dans la partie "Document", il y a à la fois la traduction de la doc officielle, et le tutoriel, je trouve ça très bien, les deux se complètent dans le sens où l'un décrit chaque commande de manière très technique, et l'autre décrit comment se servir de GDAL dans le domaine de la reprojection.
Ce n'est pas parce que l'un a été écrit en partie grâce à l'aide de l'autre que les deux vont ensemble...

Dernière modification par Sam_Dalembert (Fri 11 June 2010 08:39)

Juste une petite correction à apporter dans ce tuto suite aux différents tests que j'ai pu effectuer :

Dans le paragraphe II-3 - Problèmes de couleurs de fond/bords
Corriger le 2° ligne de l'exemple avec -wo init_dest=255 et non ini_dest (faute de frappe je présumes !!)
Et du coup, les bords seront blanc (255) et non noir. La phrase de conclusion est inversée.

Ensuite on pourrait étendre cet exemple en  construisant un VRT avec les bords blanc directement, sans construire un ECW dans la foulée. Avec cette ligne de commande et ces options, on obtient directement un VRT avec bords blanc :
gdalbuildvrt -vrtnodata "255" -input_file_list my_liste.txt assemblage.vrt pour du niveau de gris
gdalbuildvrt -vrtnodata "255 255 255" -input_file_list my_liste.txt assemblage.vrt pour du RGB.

Merci pour le manuel !!

Dernière modification par allfab (Mon 25 October 2010 16:18)

super ce tuto!

Dernière modification par baobazz (Fri 02 December 2011 09:12)

Merci beaucoup pour ce super tuto

Quelle chance vous avez !

J'ai pu le retrouver au fond de mon "historique" !

Y.

Une bière si jamais on se rencontre

Y.

Content de voir que le PDF sert toujours

J'en ai profité pour éditer le 1er message et enlever les liens pour ne laisser que la pièce jointe, plus visible