Le portail francophone de la géomatique

Merci Even pour ta réponse et bienvenu dans le forum GeoLibre

Y.

Cool pour l'amélioration notable des perfos.

Sinon, pour le problème avec GDALCreate(), je parierais sur un oubli d'un appel à GDALClose() à la fin de ton traitement. En effet, quand tu fais un GDALRasterIO(hBand, GF_Write, ...), les données ne sont pas directement écrites dans le fichier mais mises dans des blocs en cache mémoire. Au moment du GDALClose(), avant la destruction de l'objet C++ correspondant au Dataset, on procède à l'écriture des blocs modifiés encore en mémoire sur disque (cf GDALFlushCache()). De plus, pour certains drivers comme celui du GeoTIFF, des informations essentielles outre les blocs de données sont également écrites à la toute fin, ce qui fait qu'un oubli de GDALClose() peut faire en sorte que le fichier ne soit même pas lisible.

Si ce n'est pas ça la cause de ton problème, essaye d'isoler une séquence minimale où tu reproduis un problème et poste là ici.

Even

Tu essayes de faire quoi exactement ? parce que GDALCreate() et GDALCreateCopy() font des choses complètement différentes.

- GDALCreateCopy() duplique un dataset déjà existant en un nouveau, avec les mêmes données de pixel, géoréférencement, etc... Tu as juste besoin de faire un GDALClose() juste après et voilou.
- GDALCreate() crée un nouveau dataset complètement vierge. Tout reste à faire après : les appels GDALSetGeoTransform(), GDALSetProjection(), GDALRasterIO(hBand, GF_Write, etc...), ...

Ne te base pas sur le fait que GDAL alloue des buffers internes. A la limite il pourrait ne pas en allouer du tout et écrire/lire directement les fichiers sans cache mémoire. C'est un détail d'implémentation interne à GDAL. Après un appel à GDALCreate(), il est normal que GDALGetCacheUsed() te renvoie 0 car on n'a pas besoin d'allouer de mémoire alors que rien n'a été fait (it's not a bug, it's a feature ;-)). Justement l'intérêt de GDAL est de pouvoir manipuler des fichiers énoooormes avec très peu de RAM. Le cache interne à GDAL est juste là pour "lisser" un peu les perfos en lecture en écriture.

Sinon, ton extrait de code est correct, mais ne m'éclaire pas vraiment sur tes intentions. J'ai l'impression que tu as "zappé" la partie la plus intéressante entre tes 2 extraits de code. Je me demande ce que tu fais dans ton "cast". Passes-tu bien un tableau de "short" à GDALRasterIO() ?

Voilà un exemple de code qui écrit une ligne entière avec la valeur 0, la suivante avec 1, etc...

#include "gdal.h"

int main(int argc, char* argv[])
{
    GDALAllRegister();

    int nXSize = 20;
    int nYSize = 20;
    GDALDatasetH ds = GDALCreate(GDALGetDriverByName("GTiff"), "test.tif", nXSize, nYSize, 1, GDT_Int16, NULL);
    GDALRasterBandH hBand = GDALGetRasterBand(ds, 1);

    int i, j;
    short* data = (short*)malloc(sizeof(short) * nXSize);
    for(j=0;j<nYSize;j++)
    {
        for(i=0;i<nXSize;i++)
            data[i] = j;

        GDALRasterIO( hBand, GF_Write, 0, j, nXSize, 1,
                            data, nXSize, 1, GDT_Int16, 0, 0 );
    }
    free(data);

    GDALClose(ds);
    return 0;
}