Le portail francophone de la géomatique

Pour la reprojection à la volée, vois la Warp API : http://gdal.org/warptut.html. Il y a moyen de créer le "fichier" destination comme étant un Dataset MEM, donc sans passer par le disque dur.

Pour obtenir les 3 bandes à la fois, regarde du côté de GDALDatasetRasterIO. Le format ECW comprend justement une implémentation spécifique de l'API GDALDataset::IRasterIO, et non pas l'implémentation générique par défaut.

Concernant les mosaîques, de quoi parles tu exactement ? De mosaiques de plusieurs fichiers ECW ? Dans ce cas là, un fichier VRT permet de gérer ça de façon transparente. A noter qu'un des "modes" des fichiers VRT permet également la reprojection à la volée (cf VRTWarpedDataset, GDALAutoCreateWarpedVRT(),  GDALCreateWarpedVRT())

J'ai encore quelques questions concernant la reprojection à la volée, j'ai réussi à l'effectuer en utilisant le driver "MEM" (merci encore), par contre je me demandais s'il est possible de ne reprojeter qu'une seule partie de l'image car sinon les calculs sont très longs...
Première intuition, diminuer le contenu du dataset, par exemple en passant par un fichier VRT ne reprenant qu'une zone de l'image, mais il y a peut être plus simple là encore, par exemple modifier le dataset en lui-même ou en recréer un sans passer par le fichier VRT
(N'ayant pas trouvé de doc sur les classes spécifiques à l'ECW, j'essaye de décortiquer les fichiers cpp directement)

Une alternative pourrait être d'effectuer le reprojection une fois pour toute lors de l'ajout d'une nouvelle carte, mais le SDK qui est appelé par GDAL ne permet pas la compression de fichiers d'une taille > 500 mo, il y a cependant moyen de dépasser ce problème en ajoutant l'option LARGE_OK=YES lors du CreateCopy(), vu qu'il ne s'agit pas d'une compression d'un nouveau fichier ce n'est pas incorrect au niveau de l'éthique, je le précise car je l'ai découvert pendant la rédaction de ce message et ça pourrait également aider ceux qui ont le même problème

J'ai laissé tombé la reprojection à la volée et les fichiers VRT...

Par contre j'ai quelques soucis avec la fonction GDALDatasetRasterIO sur un fichier ecw, j'ai pourtant vérifié mes paramètres mais le résultat n'est pas une seule image en couleur, mais 9 images en gris (3 lignes et 3 colonnes)...

GDALDatasetRasterIO( hSrcDS, GF_Read    // dataset à lire/écrire et mode (GF_Read ou GF_Write)
            , start_x, start_y    // coordonnées du point de départ (haut gauche)
            , length_x, length_y    // taille de la vue    
            , ret            // buffer
            , window_x, window_y    // taille de la fenetre d'affichage (et donc du buffer)
            , eDT            // type (GDT_Byte pour ecw)
            , *nb_bands        // nombre de bandes
            , band_list            // tableau référençant les bandes : NULL <=> [1 2 .. *nb_bands]
            , 0, 0, 0        // écart entre pixels/lignes/bandes dans le buffer : 0 <=> valeurs selon le type (GDT_Byte)
            );

Je n'ai pas trouvé de doc sur la classe ECWDataset et je ne sais même pas quel .h ajouter pour pouvoir travailler dessus...

Je bloque un peu...

Merci d'avance pour toute aide

Merci beaucoup pour ton aide

Pour le Warping, j'ai justement repris le code des deux premiers exemples de http://gdal.org/warptut.html (le 2ème montre comment créer un nouveau fichier, et le 1er comment effectuer le warping une fois les deux fichiers prêts), je vais suivre ton conseil et analyser le fichier gdalwarp.cpp

Ok pour le buffer je vais me résoudre à le remanier, merci, je pensais que le résultats était directement exploitable car rien n'indiquait le contraire et je n'ai pas trouvé d'exemple de son utilisation, par contre l'image est répétée 9 fois : 3 fois 3, il doit donc rester une subtilité à éclaircir car j'ai aussi testé avec un fichier TIF et l'image résultante n'était divisée qu'en 3, j'espère qu'il n'y a pas de bug (EDIT : ça fonctionne :] )

Merci encore!

EDIT2 : J'ai finalement réussi à trouver le compromis que je cherchais en utilisant un dataset de type MEM que je peux créer sans avoir d'erreur avec la fonction Create() ce qui m'épargne le temps de la copie de CreateCopy(), je peux ensuite afficher ce dataset sans chargement supplémentaire par rapport au cas où j'effectue un affichage direct, pour le moment tout a l'air correct, merci encore une fois ^^

EDIT3 : Par contre je m'aperçois que c'est encore très lent pour des fichiers de très grande taille, car la reprojection est effectuée sur la totalité du dataset, je me demande donc s'il existe tout simplement un moyen de ne reprojeter qu'une partie de l'image du dataset (celle qu'on affiche dans la vue) même s'il faut refaire le calcul en modifiant la vue à chaque fois, je n'ai pas bien compris si cette possibilité existe (je précise que je n'arrive pas à reprojeter à partir d'un VRT, j'ai une erreur qui me dit que la bande n'est pas accessible en écriture, sur http://www.gdal.org/gdal_vrttut.html j'ai lu "The VRTRawRasterBand supports in place update of the raster, whereas the source based VRTRasterBand is always read-only.", et je n'ai pas trouvé d'info sur les WarpDataset)

Dernière modification par jza (Thu 11 June 2009 12:45)

Ok l'utilisation de VRT nécessite quand même de reprojeter toute l'image même si on souhaite l'afficher dans un miniature
Nous allons procéder sans la reprojection à la volée car finalement il semblerait que la solution la plus fiable niveau ressources et performances pour l'ensemble des utilisations possibles, soit de créer lun fichier VRT associé à chaque ECW lors de l'ajout de ce dernier et ce une fois pour toute à l'aide de la commande gdalwap
De cette façon le programme final ouvre ce VRT ce qui nous épargne les calculs de reprojection à chaque session d'utilisation par la suite

Je te remercie encore une fois pour ton aide Rouault

Je poste la solution trouvée par un collègue au cas où elle pourrait aider quelqu'un :

On prend un dataset duquel on extrait la vue souhaitée puis on crée un dataset temporaire (create() avec driver MEM) auquel on ajoute les données du tableau de pixels extrait
Puis on applique la reprojection uniquement sur ce dataset et de cette façon on contourne le problème des calculs sur la résolution totale et le résultat est très rapide et ne prend pas beaucoup de mémoire au final

En résumé, on utilise simplement un troisième dataset qui sert de source à la reprojection, méthode simple et logique :p

Dernière modification par jza (Tue 16 June 2009 16:34)