Le portail francophone de la géomatique

Bonjour à toutes et tous,

Mon besoin me paraît assez anodin, et je suis surpris de ne pas trouver d'outil pour y répondre facilement... Je me demande donc si j'ai bien cherché avant de poser ici la question... J'espère ne pas vous faire perdre votre temps.

Pour un raster de valeurs entières correspondant à des catégories, je souhaite :
1) lister toutes les catégories présentes dans le raster,
2) dénombrer le nombre de pixels par catégorie,
3) exporter ces valeurs dans un tableau csv ou xlsx.

J'aurai besoin de cette fonction dans des contextes variés, je souhaite donc un fonctionnement assez général. Je n'ai encore jamais fait, mais je suis prêt à coder en python s'il le faut.

J'ai bien vu la fonction de statistiques de zone, mais son usage supposerait que, pour chacune des valeurs présentes dans le raster, je commence par reclasser le raster en 1 (pixel == valeur) et 0 ( pixel != valeur), puis que j'applique cette fonction de statistique (pour obtenir la somme), puis que je traite l'ensemble des résultats pour les mettre dans un tableau de synthèse avec le bon code... Ce serait fort laborieux, il me semble. Alors avant de me lancer là dedans, je vais attendre vos retours...

D'avance, merci de votre aide, Jean-Baptiste