Je travaille sur un modèle informatique de l'abondance de pollinisateurs sauvages à travers un paysage. Le modèle lui-même est complet et je suis maintenant aux prises avec une étape de post-traitement.
J'ai mon raster d'approvisionnement en pollinisateur GDAL qui ressemble à ceci (des couleurs plus claires signifient une fréquentation plus élevée des pollinisateurs au pixel près):
Et j'ai un fichier de formes OGR de points représentant des exemples d'emplacements dans le paysage:
J'essaie d'analyser les pixels situés sous ces points, mais pour ce faire, je dois pouvoir extraire la valeur d'un pixel sous un point.
Est-il possible d'extraire la valeur d'un pixel sous un point en utilisant uniquement OGR et GDAL via Python? Je préférerais éviter de lire l'intégralité du raster en mémoire ReadAsArray()
car mes rasters en sortie sont très très gros (trop gros pour être stockés dans la mémoire).
J'ai remarqué ce message , qui est similaire, mais nécessite un appel de ligne de commande.