Le portail francophone de la géomatique

Bonjour,

Je m'excuse par avance si la réponse à mon problème a déjà été donnée, mais je n'e l'ai pas trouvé sur le forum, ni par mes recherches sur Internet et dans les modules de Qgis.

Je dois réaliser des rasters pour construire des variables environnementales qui représentent soit une densité (ou moyenne) soit une surface totale d'habitats dans un rayon déterminé autour de chacun de mes pixels de 10 m de mùa zone d'étude (la Bretagne continentale, donc un Raster de taille assez importante).

J'ai une information de base, qui suivant les variables peut-être un raster de résolution variable ou un vecteur, contenant les habitats sur lesquels je travaille (par exemple les plans d'eau, les forets de tel ou tel type, les haies, les vergers, les zones urbaines...).

Je cherche donc une méthode assez simple et rapide (j'ai un paquet de rasters à construire) pour déterminer la surface de l'habitat déterminé recensé dans un rayon de 500 m / 2000 m / 5000 m suivant les variables considérées.

Je pensais trouver rapidement une méthode pour simplement dénombrer les pixels "positif" (codés 1) de mon habitat rasterisé dans un rayon donné autour de chaque cellules de mon masque (ma zone d'étude rasterisée à 10m), mais je ne trouve pas...

Par exemple, je souhaite connaître la "quantité" (nombre de pixels, ou surface totale ou ou densité ou moyenne... peut importe) de forets de feuillus recensée dans un rayon de 5000 m en tout "point" de ma zone d'étude. Je veux donc obtenir un raster de résolution 10 m où la valeur de chaque pixel est cette "quantité" de forets de feuillus dans les 5km environnants. Mon information de départ est un vecteur de polygones des forets de feuillus. Je l'ai rasterisé, et j'en suis à chercher une solution sans la trouver.

Il y a bien "Fragmentation" sous SAGA qui renvoie une couche de densité, mais qui génère aussi des sorties "obsolètes" pour moi (connectivité et fragmentation) qui alourdissent sérieusement les calculs inutilement.

Est-ce que l'un de vous aurait une méthode relativement simple à me proposer ?

Merci d'avance, Thomas