Analyses cartographiques avec R

L’objectif de cette note est de présenter un projet fil conducteur qui permette de dérouler un ensemble d’opérateurs SIG sous R, vecteur comme raster. Le projet fil conducteur retenu concerne la classification d’un ensemble de parcelles viticoles sur un petit bassin versant vis-à-vis du risque d’érosion ravinaire. Il s’appuie sur un schéma classique de croisement entre aléa et vulnérabilité par parcelle suivant un système de score (théorique et virtuel pour simplifier l’exercice). Nous négligerons dans l’approche développée beaucoup de questions et problèmes (effet de résolution spatiale des données, significativité physique de la démarche proposée, etc). Encore une fois elle n’est proposée ici qu’à des fins de simplification pédagogique. Plan du cours 1 Objectif 1.1 Aléa 1.2 Vulnérabilité 2 Description succincte des données 2.1 Site support : Les bassins versants de Roujan, du Bourdic et de la Peyne 2.2 MNTs 2.3 Images aériennes (NDVI) 2.4 Occupation du sol 2.5 Réseau hydrographique 2.6 Piérosité 2.7 Téléchargement et préparation des données 2.8 Chargement et installation des librairies 2.9 Table récapitulative des données calculées 2.10 Projection et étendue du projet 3 Import et analyse de données raster 3.1 Chargement et installation des librairies 3.2 Chargement du fichier Raster | raster 3.3 Découpage d’un raster | crop 3.4 Agrégation d’un raster | aggregate 3.5 Adaptation des résolutions entre deux rasters | resample 3.6 Reprojection d’un raster | projectRaster 3.7 Enregistrement d’un raster | writeRaster 3.8 Opérateurs sur données Raster | terrain 3.9 Opérateurs sur données Raster | focal 3.10 Import de données via des serveurs extérieurs | getData 3.11 Utilisation de plateforme tierce pour l’analyse raster | RSAGA 3.11.1 Création du repertoire de travail et export du fichier temporaire | dir.create 3.11.2 Pointage vers le logiciel SAGA | rsaga.env 3.11.3 Choix de la librairie SAGA | rsaga.get.libraries 3.11.4 Calcul et interfaçage avec SAGA | rsaga.geoprocessor 3.11.5 Récupération du résultat 4 Import et analyse de données vecteur 4.1 Chargement et installation des librairies 4.2 Chargement des données vecteur | readOGR 4.3 Calcul des aires de polygones | gArea 4.4 Découpe d’un Spatial Object crop […] 4.5 Calcul des centres de polygones coordinates 4.6 Accès aux données slot des vecteurs | slot 4.7 Manipulation de la table attributaire de données vectorielles | $ 4.8 Reprojection de données vectorielles | spTransform 4.9 Analyses vecteur: overlays | gUnion gDifference + - 4.10 Agrégation de données vectorielles | gUnaryUnion 4.11 Création d’une zone tampon | gBuffer 4.12 Distances entre données vectorielles | gDistance 4.13 Calcul de l’aire de chaque cellule d’un MNT | surfaceArea 4.14 Connectivité et topologies des données vectorielles | poly2nb 5 Intégration raster-vecteur 5.1 Chargement et installation des librairies 5.2 Chargement des données multicanaux RGB et IRC | merge stack aggregate crop 5.3 Extraction de données vecteur polygone-points | over 5.4 Rasterisation d’un fichier vecteur | rasterize 5.5 Polygonisation d’un raster | rasterToPolygons 5.6 Simplification des limites d’un SpatialPolygonsDataFrame | thinnedSpatialPoly 5.7 Extraction de données raster-polygone rasterisé | zonal 5.8 Extraction de données raster-vecteur | extract 5.9 Création d’une zone tampon | buffer 5.10 Interpolation spatiale | interpp 5.11 Calcul de la surface drainée et du Stream Power Index 5.11.1 Incrustation du réseau hydrographique dans le MNT 5.11.2 Bouchage des trous du MNT incrusté | rsaga.geoprocessor 5.11.3 Création du repertoire de travail et export du fichier temporaire | dir.create 5.11.4 Pointage vers le logiciel SAGA | rsaga.env 5.11.5 Choix de la librairie SAGA | rsaga.get.libraries 5.11.6 Calcul de la surface drainée | rsaga.geoprocessor 5.12 Evaluation du risqu