Quantcast
Channel: Forum SIG - Systèmes d'Information Géographique et Géomatique
Viewing all articles
Browse latest Browse all 8030

[Gdal] Python GetStatistics() : Ignorer le NoData

$
0
0
Bonjour,

je tente, sous Gdal-python , de calculer les statistiques et les pyramides
de mes images de la façon suivante :
Code:

        # open the image
        inDs = gdal.Open(dataset, GA_ReadOnly)

        for i in range(inDs.RasterCount):

            inBand = inDs.GetRasterBand(i+1)

            # flush data to disk, set the NoData value and calculate stats
            inBand.FlushCache()
            inBand.SetNoDataValue(-99)
            inBand.GetStatistics(0, 1)

        # build pyramids
        gdal.SetConfigOption('USE_RRD', 'YES')
        inDs.BuildOverviews(overviewlist=[2,4,8,16,32,64,128])

Mais le problème qui se pose, est qu'il prend en compte mon NoData (-99)
dans le calcul des stats. Par conséquent je ne vois pas grand chose à
l'affichage sur mon image. Juste du blanc et du noir.
Je précise que les valeurs sont correctes, c'est juste l'affichage qui fait
default.

Alors je voudrais savoir si vous avez une idée qui permettrai au "GetStatistics"
d'ignorer le NoData dans le calcul des stats ?!

Je précise également que je peux le faire avec un logiciel libre ou propriétaire
mais je voudrais le coder sur Python avec GDAL.

Merci d'avance pour vos éventuelles suggestions :D,

Viewing all articles
Browse latest Browse all 8030

Trending Articles