Le portail francophone de la géomatique

Aucune erreur flagrante ne me saute aux yeux.

Plusieurs points à vérifier (désolé si ça parait stupide. c'est comme les questions des téléopérateurs d'un SAV) :
1) Ton dataset est-il bien ouvert avec GA_Update ?
2) Comme tu utilises le même objet GDALRasterBand en lecture et en écriture, j'en déduis que tu veux modifier ton dataset en live, et non pas en lire un pour en écriture un autre ? Sinon y a un bug...
3) Le type natif de band_mask est-il bien GDT_Float32 ? Sinon il y a un bug. De toute façon pour plus de sécurité, il est préférable que tu remplaces GDALGetRasterDataType(band_mask) par GDT_Float32. C'est toi qui impose le type du buffer lu/écrit en cohérence avec sa déclaration. Ca n'a a priori rien à voir avec le type natif retourné par le driver du dataset.

Tout dépend du driver utilisé.

Il me semble que le tiff supporte ce genre d'image (d'ailleurs, c'est le sujet de la section 3 de la TIFF 6.0 Specification : "Bilevel Images").

Mais peut-être que gdal ne supporte pas ce genre d'opération et qu'il faille utiliser une autre bibliothèque de manipulation d'image (gdal étant spécifique aux données géographiques).

OK, je pense qu'il faudrait que tu utilises 2 dataset :
* un ouvert en lecture uniquement (de type GDT_Float32)
* l'autre pour créer ton masque que tu crées à partir de rien avec GDALCreate(driver, filename, xsize, ysize, nbands, datatype, options ) et que tu écris. Le mieux est de prendre datatype = GDT_Byte. Si tu utilises GTiff comme driver, tu peux passer l'option NBITS=1 pour forcer la création d'une image sur 1 bit (par contre le buffer de manipulation vis à vis de GDAL reste du GByte. Le passage de 8bit->1bit est fait en interne dans le driver GTiff), avec éventuellement un type de compression. Par ex :

Quelque chose dans le goût de :

#include "cpl_string.h"

char** papszOptions = NULL;
papszOptions = CSLAddString(papszOptions, "NBITS=1");
papszOptions = CSLAddString(papszOptions, "COMPRESS=PACKBITS");
GDALDriverH hDriver = GDALGetDriverByName("GTiff");
GDALDatasetH hDstDS = GDALCreate(hDriver, "out.tif", XSize, YSize, 1, GDT_Byte, papszOptions);
GDALRasterBandH hDstBand = GDALGetRasterBand(hDstDS, 1);
CSLDestroy(papszOptions);

float *pafScanline;
pafScanline = (float*)CPLMalloc(sizeof(float)*XSize);
GByte* pabyScanline;
pabyScanline = (GByte*)CPLMalloc(sizeof(GByte)*XSize);

  for (int i = 0; i<YSize; i++)
  {
    GDALRasterIO(hSrcBand, GF_Read, 0, i, XSize, 1, pafScanline, XSize, 1, GDT_Float32, 0, 0);

    for(int k=0; k<XSize; k++)
    {      
      pabyScanline[k] =(pafScanline[k] != 0);
    }
    // Ecriture
    GDALRasterIO(hDstBand, GF_Write, 0, i, XSize, 1, pabyScanline, XSize, 1, GDT_Byte, 0, 0);
  }

  GDALClose (hDstDS);

La doc des options du format GeoTIFF : http://gdal.org/frmt_gtiff.html

Attention tous les visualisateurs de TIFF ne gèrent pas forcément toutes les différentes options disponibles dans le driver GDAL...

Edit : Ah je viens de voir que tu es intéressé par le format Imagine, mais bon c'est quasiment la même chose. Tu peux utiliser NBITS=1 mais et COMPRESSED=YES par exemple, cf http://gdal.org/frmt_hfa.html

Dernière modification par rouault (Tue 06 July 2010 18:22)

Tout fonctionne comme prévu.

Merci beaucoup !