====== Projets_et_stages_2013-2014/EspacesVertsv1 ======

[[Fichier:fleche.jpg|100px|alt=Retour|link=http://im2ag-mwiki.e.ujf-grenoble.fr/mediawiki/index.php/Projets_et_stages_2013-2014#M.C3.A9thodologie_2]]

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


====Etape 1 : Sélection des communes de la Région Urbaine de Grenoble====
cf méthode de la partie “Diffusion”

====Etape 2 : Choix des parcelles correspondant à des espaces verts====
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 [[http://im2ag-mwiki.e.ujf-grenoble.fr/mediawiki/index.php/Projets_et_stages_2013-2014/nomenclature|ici]]. ([[https://drive.google.com/file/d/0BxpXy3QTc2FyNk5vTGhCelYxLXc/edit?usp=sharing|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

====Etape 3 : Sélection des parcelles choisies dans un tableur Excel====
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 :

 <nowiki>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 </nowiki>

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.

[[Fichier:EV01.png|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//

[[Fichier:10AjCoucheVect.png|300px|Ajouter une couche vecteur]] → [[Fichier:20Proprietes.png|300px|proprietes]] → [[Fichier:30Jointure.png|600px|jointure]]

====Etape 4 : Carroyage de la RUG====
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//

[[Fichier:40Intersection.png|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.

====Etape 5 : Sélection des carreaux d’espaces verts====
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 coordonnées de ce carreau
  * le nom de la commune à laquelle il appartient et toutes les informations relatives à cette commune

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.

[[Fichier:EV02.png|center|Carroyage]]

<center>Rouge : communes de la RUG</center>
<center>Bleu : carreaux sélectionnés</center>

====Etape 6 : Calcul de la densité d’espaces verts par commune====
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 :
  - Parcourir l’ensemble des carreaux
  - Pour chaque carreau parcouru, ajouter 40000m² (200*200) au décompte de la superficie d’espaces verts de la commune à laquelle il appartient
  - Puis, pour chaque commune, diviser la superficie d’espaces verts obtenue par la superficie de la commune.
  - Résultat : Densité (%) d’espaces verts pour chaque commune

Utlisation d’une 2e macro VBA sous Excel pour ce calcul :

 <nowiki>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 </nowiki>

====Etape 7 : Affichage====
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//

[[Fichier:50Graduation.png|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”//

[[Fichier:60Diagrammes.png|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” //

[[Fichier:70Etiquettes.png|600px|center|affichage du nom des communes]]

Résultats obtenus :

- Pour l'année 2000

[[Fichier:EV03.png|Espaces verts + population en 2000|résultat 2000]]


- Pour l'année 2006

[[Fichier:EV04.png|Espaces verts + population en 2006|résultat 2006]]