Suite

Rasteriser une couche de points avec plusieurs entités tombant à l'intérieur d'une cellule raster à l'aide de GDAL ?


Comment faire un raster avec des valeurs qui sont la somme de tous les points à l'intérieur de chaque cellule raster ?

La valeur de combustion est un attribut spécifique.

En utilisant gdal_rasterize ou gdal.RasterizeLayer(), le raster résultant n'a que la valeur d'un point aléatoire (ou probablement le premier ou le dernier) à l'intérieur d'une cellule raster. J'ai déjà essayé l'option gdal_rasterize ALL_TOUCHED, mais cela ne change pas le comportement.


L'ajout de 'MERGE_ALG=ADD' à la liste des options de gdal.RasterizeLayer() a fait l'affaire, grâce à kyle. Un exemple:

gdal.RasterizeLayer(target_ds, [1], src_layer, options = ['MERGE_ALG=ADD'])

Ceci est documenté dans http://www.gdal.org/gdal__alg_8h.html.


Voir la vidéo: Simple raster manipulation using GDAL (Octobre 2021).