Le portail francophone de la géomatique

Bonjour à tous,

Je me tourne vers vous aujourd'hui car je suis en panne sèche d'inspiration, de stratégies à mettre en place pour réaliser cette tâche complexe.

Mise en situation :

J'ai en ma possession une base de donnée reprenant environs 6.000.000 de points étendu sur tout un territoire, correspondant à autant d'observation de plante et d'animaux sur le terrain. A chacun de ses points sont liés tout un tas d'information (environ 40 colonnes) dans la table des attributs qui sont importantes à conserver pour le résultat final (dont les coordonnées géographique, des codes d'identification de l'espèces, les dates, etc.)

De l'autre côté, j'ai un couche qui reprend 45.000 polygones autour desquels j'ai réalisé 3 différents buffer (tampons) : buffer 500 m, buffer 1000 m et buffer 2000 m. Ces polygones représentent des zones dans lesquelles je souhaite savoir qu'elles ont été les observations réalisées.

Résultat recherché :

Je cherche à intégrer les 6.000.000 d'observations dans les polygones. Ces 6.000.000 d'observations sont séparés en 3 couches shapefile. Une couche correspondant aux observations que je souhaite inclure dans le buffer 500, une autre pour le buffer 1000 et une autre pour le buffer 2000.

Le résultat final serait de pouvoir, pour chacun des 45 000 polygones, pour dire par exemple pour le polygone n°4859:
- Dans une zone de 500m autour (buffer 500), on a observé les points n° 3, 85, 4569, 29394, 38888 et 983338
- Dans une zone de 1000m (buffer 2000) autour, on a observé etc, etc, etc
- IDEM zone 2000m (buffer 2000)

Ci-dessous, un exemple illustré avec 3 polygones du buffer 2000 et environ 400 points.
https://ibb.co/rt3yPXt

Avez-vous une idée de la stratégie à appliquer, comment est-ce que je pourrais réaliser ceci ?

Merci d'avance pour vos éventuelles réponses, Arthur