gdal_proximity génère des zéros sur les gros fichiers raster

Question

J'ai un raster (fcc.tif) de forêt (valeur 1) et des pixels déboisés (valeur 0) sans valeurs de données (255) ayant une résolution de 30 m pour la Papouasie-Nouvelle-Guinée. C'est un grand raster avec 60459 x 40811 pixels.

J'essaie d'utiliser gdal_proximity.py (GDAL version 2.1.2) pour calculer la distance (en m) aux pixels déboisés. J'utilise la commande suivante:

gdal_proximity.py fcc.tif dist_defor.tif -values 0 -ot UInt32 -distunits GEO -use_input_nodata YES -nodata 4294967295 

Le calcul fonctionne bien pour la plupart de la région, mais produit des valeurs de zéros incorrectes sur une grande surface ayant une sorte de forme de disque (marqué en rouge sur la figure ci-dessous).

Distance to deforestation in Papua New Guinea

Dans ce domaine, les pixels déboisées sont présents, mais ils ne semblent pas être pris en compte pour le calcul de distances.

La trame fcc.tif est disponible here

Answer 1

Cela était dû à un débordement d'entier qui a eu lieu au cours du calcul de la distance sur les grands ensembles de données. Le problème a été résolu par https://trac.osgeo.org/gdal/ticket/7102. Merci aux développeurs et mainteneurs de GDAL qui ont résolu ce problème.