17.14. Premier exemple d’analyse

Note

Dans cette leçon, nous effectuerons une véritable analyse en utilisant uniquement la boîte à outils, afin que vous puissiez vous familiariser avec les éléments du module de traitements.

Maintenant que tout est configuré et que nous pouvons utiliser des algorithmes externes, nous avons un outil très puissant pour effectuer une analyse spatiale. Il est temps de réaliser un exercice plus important avec des données du monde réel.

Nous utiliserons l’ensemble de données bien connu que John Snow a utilisé en 1854, dans ses travaux révolutionnaires (https://en.wikipedia.org/wiki/John_Snow_%28physician%29), et nous obtiendrons des résultats intéressants. L’analyse de cet ensemble de données est assez évidente et il n’est pas nécessaire de recourir à des techniques SIG sophistiquées pour obtenir de bons résultats et conclusions, mais c’est un bon moyen de montrer comment ces problèmes spatiaux peuvent être analysés et résolus en utilisant différents outils de traitement.

Le jeu de données contient des fichiers Shape qui référencent les décès dus au choléra, les emplacements des pompes à eau ainsi qu’une carte raster issue d’OSM au format TIFF. Ouvrez le projet QGIS correspondant à cette leçon.

../../../_images/project2.png

La première chose à faire est de calculer le diagramme de Voronoï (alias polygones de Thiessen) de la couche des pompes, pour obtenir la zone d’influence de chaque pompe. L’algorithme du diagramme de Voronoï peut être utilisé pour cela.

../../../_images/voronoi.png

C’est assez facile, mais cela nous donne déjà des informations intéressantes.

../../../_images/voronoi2.png

Il est clair que la plupart des cas soient dans un des polygones.

Pour obtenir un résultat plus quantitatif, nous pouvons compter le nombre de décès dans chaque polygones. Comme chaque point représente un bâtiment où il y a eu un ou des décès, et que le nombre de décès est stocké dans un attribut, nous ne pouvons pas simplement compter les points. Nous avons besoin d’un nombre pondéré, donc nous utiliserons l’outil Compter les points dans les polygones (pondérés).

../../../_images/pointsinpoly.png

Le nouveau champ sera appelé DEATHS, et nous utiliserons le champ COUNT comme champ de pondération. La table de résultats reflète clairement que le nombre de décès dans le polygone correspondant à la première pompe est beaucoup plus grand que celui des autres.

../../../_images/pointsinpolytable.png

Une autre bonne façon de visualiser la dépendance de chaque point dans la couche Cholera_deaths avec un point dans la couche Pumps est de dessiner une ligne vers le plus proche. Cela peut être fait avec l’outil Distance au nœud le plus proche, et en utilisant la configuration montrée ci-après.

../../../_images/nearest.png

Le résultat ressemble à cela :

../../../_images/nearestresult.png

Bien que le nombre de lignes soit plus grand dans le cas de la pompe centrale, n’oubliez pas que cela ne représente par le nombre de décès, mais le nombre de lieux où des cas de choléra ont été détecté. C’est un paramètre représentatif, mais il ne considère pas que certains emplacements puissent représenter plus de cas qu’un autre.

Une couche de densité nous donnera aussi une vision très claire de ce qu’il se passe. Nous pouvons la créer avec l’algorithme Densité par noyau. En utilisant la couche Cholera_deaths et son champ COUNT comme champ de pondération, avec un rayon de 100 ainsi que l’emprise et la taille des cellule de la couche raster des rues, nous obtenons quelque chose comme ça.

../../../_images/density.png

Souvenez-vous que pour obtenir l’emprise de sortie, vous ne devez pas la taper. Cliquez sur le bouton sur le côté droit et sélectionnez Utiliser l’emprise de la couche/du canevas.

../../../_images/layerextent.png

Sélectionnez la couche raster des routes et son étendue sera automatiquement ajoutée au champ de texte. Vous devez faire de même avec la taille de cellule, en sélectionnant la taille de cellule de cette couche.

En combinant la couche des pompes, nous voyons qu’il y a une pompe qui se trouve clairement dans le point chaud où la densité maximale de cas de décès a été trouvée.