Différentes étapes pour construire la carte :
*Etape 1 : Sélection des communes de la Région Urbaine de Grenoble *Etape 2 : Choix des parcelles correspondant à des espaces verts *Etape 3 : Sélection des parcelles choisies dans un tableur Excel *Etape 4 : Carroyage de la RUG *Etape 5 : Sélection des carreaux d’espaces verts *Etape 6 : Calcul de la densité d’espaces verts par commune *Etape 7 : Affichage *Etape 8 : Affichage des couleurs officielles des parcelles de la RUG
cf méthode de la partie “Diffusion”
La sélection des parcelles est faite à partir de la liste complète des différents types de paysages qui composent les parcelles.
La liste complète est visible ici. (Fichier Excel)
<u>Codes sélectionnés et libellés correspondants:</u>
*141 : Espaces verts urbains *213 : Rizières *221 : Vignobles *222 : Vergers et petits fruits *223 : Oliveraies *231 : Prairies *241 : Cultures annuelles associées aux cultures permanentes *243 : Surfaces essentiellement agricoles, interrompues par des espaces naturels importants *244 : Territoires agro-forestiers *311 : Forêts de feuillus *312 : Forêts de conifères *313 : Forêts mélangées *321 : Pelouses et pâturages naturels *322 : Landes et broussailles *323 : Végétation sclérophylle *324 : Forêt et végétation arbustive en mutation *333 : Végétation clairsemée
Utilisation du tableau des parcelles en ajoutant une 4e colonne contenant les 17 codes choisis.
Macro VBA (Visual Basic for Applications) pour la sélection des parcelles :
Sub Selection() Dim Nb As Integer Dim Nbb As Integer Dim Aj As Integer Dim Verif As Integer Dim Bool As Boolean Nbb = 2 Aj = 2 MsgBox ("Debut de la procedure") ‘parcours de toutes les parcelles For Nb = 2 To 6305 Bool = False ‘vérification de l’appartenance de la parcelle à celles choisies For Verif = 2 To 18 If Range("B" & Nb).Value = Range("D" & Verif).Value Then Bool = True End If Next Verif ‘Bool est vrai s’il faut retenir la parcelle If Bool Then ‘recopiage des infos de la parcelle dans les colonnes F,G,H Range("F" & Aj).Value = Range("A" & Nb) Range("G" & Aj).Value = Range("B" & Nb) Range("H" & Aj).Value = Range("C" & Nb) Aj = Aj + 1 End If Next Nb End Sub
Résultat obtenu :
Copie des 3 informations (code, libellé, superficie) des parcelles iséroises selectionnées dans les colonnes F,G,H du tableur Excel.
center|Sélection des parcelles
<center>Résultat de la sélection des parcelles</center>
<u>Problème :</u> Impossible de visualiser une carte de ces nouvelles informations sous QGIS sans le fichier .shp correspondant.
<u>Solution :</u> *Ouverture du fichier .xls sous QGIS
couche → ajouter une couche vecteur
*Jointure de celui-ci avec l’ancien fichier afin de récupérer le fichier .shp nécessaire
clic droit → propriétés → + → joindre la couche : fichier .xls ; joindre le champ : ID
300px|Ajouter une couche vecteur → 300px|proprietes → 600px|jointure
Découpage de la Région Urbaine de Grenoble en carreaux de 200x200m afin de traiter les données.
Création d’une grille rectangulaire englobant la RUG composée de carreaux de 200x200m via QGIS (Couche → Nouveau → Nouvelle couche shapefile) puis intersection spatiale avec les communes sélectionnées en étape 1 afin de conserver uniquement les carreaux qui touchent au moins une commune de la RUG
Vecteur → Outils de géotraitement → Intersection → Couche vectorielle de saisie : fichier grille ; Couche d’intersection : fichier “communes” de l’étape 1
400px|intersection avec la RUG
A partir de l’intersection entre les carreaux et les communes de la RUG, il est également possible d’associer à chaque carreau la commune à laquelle il appartient.
Précision : les carreaux à cheval sur plusieurs communes sont comptés plusieurs fois. Une fois pour chaque commune.
Le résultat de l’étape 4 est une liste de carreaux de 200x200m avec, pour chacun d’eux, les informations suivantes (entre autres) :
Les autres colonnes sont inutiles.
Il faut donc sélectionner, dans cette liste, les carreaux qui appartiennent (qui touchent) une des parcelles sélectionnées à l’étape 3.
Intersection spatiale de la liste des carreaux avec la liste des parcelles sélectionnées dans l’étape 3 (Vecteur → Outils de géotraitement → Intersection → Couche vectorielle de saisie : fichier carreaux de la RUG créé ci-dessus ; Couche d’intersection : fichier parcelles de l’étape 3).
<u>Résultat :</u> liste des carreaux d’espaces verts avec, pour chacun d’eux, l’information du nom de la commune à laquelle il appartient.
<center>Rouge : communes de la RUG</center> <center>Bleu : carreaux sélectionnés</center>
Le but de cette étape est d’obtenir la densité d’espaces verts pour chaque commune (cela permettra de les comparer avec les densités de population). Il faut donc :
Utlisation d’une 2e macro VBA sous Excel pour ce calcul :
Sub Densite() Dim Indice As Variant Dim IndiceId As Variant Dim Somme As Variant Dim Superficie As Variant Dim Densite As Variant Dim CommCour As Variant Somme = 0 N = 0 IndiceId = 2 Superficie = Range("S" & IndiceId).Value / 100 'en km² MsgBox ("Debut procedure") ‘parcours de tous les carreaux de la RUG For Indice = 2 To 110804 ‘Tant qu’il s’agit toujours de la meme parcelle, on augmente Somme If Range("I" & Indice).Value = Range("AC" & IndiceId).Value Then Somme = Somme + 40000 'en m² (un carreau a une surface de 40000m²) Else ‘Parcelle suivante : calcul de la densité pour la parcelle qui est ‘terminée, réinitialisation des variables et incrémentation des indices Densite = Somme / (Superficie * 100000) Range("AD" & IndiceId).Value = Somme Range("AE" & IndiceId).Value = Densite Range("AF" & IndiceId).Value = Superficie Superficie = Range("S" & Indice).Value / 100 'en km² IndiceId = IndiceId + 1 Somme = 0 End If Next Indice End Sub
L’affichage se fait exclusivement sous QGIS.
1. Affichage des densités d’espaces verts Coloration de chaque commune en dégradés de verts en fonction de la densité d’espaces verts.
Choix d’affichage : 10 catégories à effectifs égaux (autant de communes dans chaque catégorie)
clic droit → propriétés → Style → Gradué → Colonne : Densité d’espaces verts → Classes : 10 → Mode : Quantiles
600px|center|affichage des espaces verts
2. Affichage des densités de population Affichage d’un disque rouge dont la taille dépend de la densité de population dans la commune en question.
Choix d’affichage : la surface du disque est proportionnelle à la densité de population de la commune.
clic droit → propriétés → Diagrammes → Type de diagramme : camembert → Ajouter l’attribut “densité de population” en bas → Taille → taille fixe : non ; Mettre à l’échelle linéairement par rapport à l’attribut “densité de population”
600px|center|affichage de la population
3. Affichage du nom des communes Affichage du nom de certaines communes.
Choix d’affichage : Communes aux extrémités (Treminis, Charancieu, Chalon, Allevard) + communes à fortes densités, grandes superficies et noms courts.
clic droit → propriétés → Etiquettes → Etiqueter cette couche avec : Editer l’expression → Entrer par exemple : “CASE WHEN (POPULATION>2 AND SUPERFICIE>1000 AND length(NOM_COMM)<7 END”
600px|center|affichage du nom des communes
Résultats obtenus :
- Pour l'année 2000
Espaces verts + population en 2000|résultat 2000
- Pour l'année 2006
