Projet cartographie Foyers grippe avaire H5 en Belgique
Projet cartographie Foyers grippe avaire H5 en Belgique
Adapter le code R développé dans le projet précédent à un GeoJSON de la Belgique où les limites administratives sont les communes telles qu'elles étaient définies avant la fusion des communes
Permettre le paramétrage des distances retenues pour établir la zone de protection et la zone de surveillance selon les épizooties
Production d'une une carte HTML claire et interactive où le survol d’une commune affiche son nom et son code PSEUDONIS, le code postal n'étant pas une variable unique
Sources
Le jeu de données utilisé est AdminVector. Il correspond au jeu de données vectorielles des entités administratives et statistiques belges. Il est librement téléchargeable et disponible à partir de la page suivante :
https://www.geo.be/catalog/details/fb1e2993-2020-428c-9188-eb5f75e284b9?l=fr
La raison principale qui nous a orienté vers AdminVector est qu’il permet de récupérer les polygones des communes telles qu’elles étaient définies avant la fusion des communes.
En effet, via QGIS par exemple, nous constatons que le fichier adminvector_3812.gpkg permet d’ajouter plusieurs objets :
La couche MUNICIPALITY (565 entités) correspond aux communes actuelles. Depuis l’entrée en vigueur de la loi et de l’arrêté royal organisant la fusion (effective le 1er janvier 1977), la réorganisation a relégué au second plan les 2 664 entités préexistantes.
Or, ce sont ces entités historiques qui nous intéressent pour la mise en place d’une zone de protection et d’une zone de surveillance autour d’un foyer épizootique tout en s'appuyant sur des limites adminstratives précises.
Nous utilisons la couche MUNICIPALSECTION, qui contient les 2 664 entités qui correspondent aux territoires des communes autonomes au 1er janvier 1961.
La table des attributs est présentée à droite.
[zipcode] et [niscode_municipality] sont pas univoques : plusieurs communes peuvent avoir un code postal et/ou un niscode_municipality identiques. Seul l’attribut [.pseudonis] d’une section de commune est unique.
Selon la région linguistique à laquelle la commune est rattachée, son nom peut être en allemand [nameger], en français [namefre] ou néerlandais [namedut]. Une commune bilingue doit préserter les variables [namefre] et [namedut].
Nous procédons à une exportation de la couche [municipalsection] en format GeoJson.
Table des attributs
Données :
Communes (GeoJSON) : polygones communaux et création des codes LIBGEO et CODGEO en R pour s'harmoniser au code précédent.
Le libgeo n’existe pas tel quel . En Belgique, 3 langues coexistent : allemand, français et néerlandais. Les noms des communes sont donnés dans la langue officielle de la région linguistique. Dans une commune bilingue il faudra concaténer [namefre] ; « _ »;[namedut] lorsque les deux valeurs sont présentes.
Structure du GeoJSON après modifications
Foyers (Fichier Excel) : Commune, Latitude, Longitude, rayon_protection, rayon_ext
Note : pour le développement suivant, les coordonnées GPS des foyers H5 correspondent assez fidèlement aux centres des cercles de zones des zones de protection de la cartographie publiée sur le site de l'AFSCA
Nous avons opté pour un paramétrage possible du rayon de la zone de protection (ici 3 km) et la zone de surveillance (10 km) directement via le fichier de paramétrage Excel.
Méthodologie
1. Lecture et validation des sources.
2. Pour chaque foyer, génération vectorisée de points sur anneaux (angles × rayons) pour détecter les communes impactées.
3. Attribution de zones (Zone1 protection, Zone2 surveillance) et résolution des chevauchements par priorité.
4. Construction d’un jeu final utilisé pour l’affichage.
Visualisation
Carte Leaflet exportée en HTML autonome.
Polygones colorés par groupe.
Cercles tracés en polylines pour les rayons .
Tooltip au survol minimaliste : "Commune — Codegeo" pour garantir lisibilité.
Ici vu la présence de nombreux foyers sur une étendue restreinte, nous avons choisi de désactivation des marqueurs et leur détail avec popup au clic . Cette option peut être restaurée.
Application du Tooltip au survol minimaliste : "Libellé Commune - Codegeo - Statut Zone " aux communes survolées.
Livrables
• Fichier HTML interactif prêt à partager.
• CSV récapitulatif des communes et attributions