Newer
Older
__Occupations du sol et multi-usage de l’espace autour du Bassin d’Arcachon__
Baptiste Hautdidier
Support TD 'GIZC', 2018-2019
28-29/11/2018
Fichier `Données.zip` à télécharger via le [service Renater](https://filesender.renater.fr/?s=download&token=7f2fe7ae-5666-9287-25e2-fd67746ec82a).
## Exercice 1 - Construire l'emprise d'une zone d'étude sur la base de découpages administratifs
* _Shapes_ : `Communes33_40.shp`
* Tables : `epcicom2014.csv`
Nous allons nous intéresser au [pays](https://anpp.fr/payspetr/definition-pays-petr/) 'Bassin d'Arcachon Val de l'Eyre', un territoire dont le périmètre coincide notamment celui du SCOT (Schéma de Cohérence Territoriale) du même nom. On a ici affaire à une démarche de planification et d'urbanisme, qui interagit fortement avec les logiques de conservation (Trame verte et bleue...). La listes des communes rattachées à ce pays correspond à l'union des structures intercommunales [EPCI](http://www.insee.fr/fr/methodes/default.asp?page=definitions/etab-pub-cooper-intercom.htm) suivantes :
EPCI | N° SIREN
-------------------------------------------------|----------
CA Bassin d'Arcachon Sud-Pôle Atlantique (COBAS) | 243300563
CC du Val de l'Eyre | 243301405
CC du Bassin d'Arcachon Nord Atlantique (COBAN) | 243301504
1. Charger les géométries communales des départements des Landes et de la Gironde (fichier `Communes33_40.shp`, source [BD Carto](http://professionnels.ign.fr/bdcarto)).
2. Construire une requête sur le champ `SIREN_EPCI`, qui renvoie les communes du seul pays 'Bassin d'Arcachon Val de l'Eyre'.
(La requête est de la forme: `SELECT * FROM Communes33_40 WHERE "Champ"='valeur1' OR "Champ"='valeur1'` ...)
3. Fusionner les géométries sélectionnées, de manière à obtenir une nouvelle couche avec une seule entité, et enregistrer cette couche. Exporter préalablement la sélection. Puis utiliser `Géotraitement > Fusionner`, ou dans la _toolbox_ : `Outils de gestion des données > Généralisation > Fusionner`, en se basant sur un champ dont les valeurs sont connues comme étant identiques sur la zone.
## Exercice 2 - Analyser la construction des emprises des parcs naturels de la zone
* _Shapes_ : `Communes33_40.shp`, `pnm2014_07.shp`, `pnr2008.shp`,`Contour_PnrLG_2014_2026.shp`
Le bassin d'Arcachon est concerné par deux territoires de conservation construits sur la base d'une contractualisation avec des collectivités territoriales : un parc naturel marin (PNM) et un parc naturel régional (PNR). La délimitation de leurs emprises n'est pas un processus anodin. Le PNR a fait l'objet d'une révision de sa charte en 2013. Quant au PNM, il s'agit d'un dispositif
1. Reprendre les géométries communales (`Communes33_40.shp`, source BD Carto).
2. Ouvrir `CommunesPNR.xls` dans un tableur : que décrit le fichier ?
3. Joindre cette table avec les géométries communales. Attention aux formats des champs utilisés pour la jointure.
4. Fusionner en se basant sur les informations tabulaires, de manière à construire une carte avec une symbologie montrant l'évolution du zonage (le plus simple est d'utiliser deux champs pour la fusion). Charger `pnr2008.shp` (un export de l'[INPN](https://inpn.mnhn.fr/accueil/index) montrant l'emprise officielle )
5. Charger la couche `pnm2014_07`. Qu'observe-t-on au fond du bassin ? Charger `Contour_PnrLG_2014_2026.shp`, qui comme son nom l'indique est le contour du parc pour la période de la charte en vigueur, 2014-2026 (Voir l'article 7 dans [LégiFrance](http://www.legifrance.gouv.fr/affichTexte.do?cidTexte=JORFTEXT000025199018&dateTexte=&categorieLien=id)).
6. Se connecter en flux WMS à l'infrastructure de données spatiales régionale PIGMA. Ajouter la couche décrivant le fuseau de la LGV (choisir la couche `RFF - Projet LGV GPSO 1000m`). En complément, consulter le tracé définitif des LGV Bordeaux–Toulouse et Bordeaux– Espagne sur le [site du Grand Projet Ferroviaire du Sud-Ouest](http://www.gpso.fr/cartes.html). Quelles hypothèses peut-on faire ?
> ### Jointures sur codes INSEE
> En faisant la manipulation sur les codes INSEE, on a alors des codes numériques - mais qui doivent tout de même être stockés sous forme de chaines de caractères. Une manipulation sur le fichier Excel permet de forcer ArcGIS à traiter cette colonne au format `string`, il suffit d'insérer un ' devant le premier enregistrement. Alternativement, la jointure fonctionne sans problème sur les noms des communes.
> ### Pour la plateforme Pigma
> L'adresse pour les services WMS est : https://www.pigma.org/geoserver/wms
> Plusieurs guides sont disponibles sur leur site web. L'un concerne le [visualiseur](https://portail.pigma.org/outils-et-services/cartographie/visualiseurs-thematiques/), l'autre le [catalogue](https://portail.pigma.org/outils-et-services/catalogue/guides-dutilisation-du-catalogue/). Il existe aussi une [vidéo](http://vimeo.com/83304664).
<!-- .element height="50%" width="50%" -->
> Ajouter un serveur WMS, dans ArcCatalog*
## Exercice 3 - Confronter les zonages politico-administratifs avec la logique 'bassin versant'
* _Shapes_: `Communes33_40.shp`, `ZoneHydro33_40.shp`, `Secteurs33_40.shp`
1. Charger les couches hydro : _secteurs_ et _zones_ hydrographiques issues de la BD Carthage (Voir sur [EauFrance](http://www.glossaire.eaufrance.fr/concept/secteur-hydrographique) pour une définition).
2. Faire une sélection des secteurs intersectés par l'emprise créée pour l'exercice 1 (`Sélection > Sélectionner selon l'emplacement...`).
3. Exporter cette sélection sous une nouvelle couche.
4. Idem pour les zones hydro.
5. [_étape facultative dépendant du programme des séances précédentes avec S. Zaragosi_]. Charger les couches de vos bassins versants modélisés. Inspecter les différences avec les couches 'Carthage' en adaptant les symbologies. Croiser les géométries de vos BV avec celle de la zone d'étude (`Géotraitement > Union`).
6. Construire un tableau récapitulant les surfaces (en ha) des différents secteurs hydrographiques situés dans _et_ hors la zone d'étude: vérifier que les surfaces sont bien calculées dans un champ dédié, exporter la table attributaire et faire le traitement dans un tableur par un tableau croisé dynamique
> ### Le calcul des superficies
> C'est automatique lorsqu'on travaille en format _geodatabase_ `*.gdb` : un des intérêts de ce format ESRI.
> Mais cela requiert une manipulation manuelle avec l'outil `Calculer la géométrie` quand on part d'un _shapefile_.
## Exercice 4 - Caractériser l'occupation du sol d'un bassin versant
* _Shapes_: `CLC12_RALPC_RGF.shp`
* Symbologie: `CORINE_Land_Cover_06_Vecteur.lyr`
Corine Land Cover est la base de données pan-Européenne décrivant les changements d'occupation des sols entre els millésimes 1990, 200, 2006 et 2012. Voir [ici](http://www.statistiques.developpement-durable.gouv.fr/donnees-ligne/li/1825.html) pour une description du produit et [là](http://www.statistiques.developpement-durable.gouv.fr/donnees-ligne/t/nomenclature.html) pour le détail des postes de la nomenclature.
{: .shadow}*Les 44 postes de la nomenclature CORINE Land Cover*
1. Découper le _shapefile_ CORINE Land Cover par l'emprise constituée par les secteurs hydro (sauvegardée au 3.3), ou bien par celle des bassins versant modélisés. Réaliser le découpage en _un seul_ traitement (`Outils d'analyse > Extraire > Découper`)
2. Appliquer une symbologie avec le fichier `*.lyr` (le champ`CODE_06` du `*.lyr` correspond au `CODE_12` du `*.shp`, les nomenclatures n'ont pas changé entre les millésimes 2006 et 2012).
3. Choisir un des secteurs (ou un de vos bassins versants modélisés): l'utiliser pour découper la couche CLC.
4. Y calculer les surfaces en ha de chaque type d'occupation du sol (Cf. exercice 3.6) - mais cette fois-ci en réalisant une fusion préalable sur le champ de la nomenclature CLC (`Géotraitement > Fusionner`).
{: .shadow}*Ne pas oublier le traitement par lot, accessible par un clic-droit sur l'outil _ArcToolbox_.*
## Exercice 5 - Qualifier les transitions identifiées avec la base Corine Land Cover
* _Shapes_: `CHA06_R72_RGF.shp, CHA00_R72_RGF.shp, CHA12_RALPC_RGF.shp`
* Tables: `LCF lookup.xls`
LCF | Description
----|------------------------------------------------------
1 | _Urban land management_
2 | _Urban residential sprawl_
3 | _Sprawl of economic sites and infrastructures_
4 | _Agriculture internal conversions_
5 | _Conversion from forested & natural land to agriculture_
6 | _Withdrawal of farming_
7 | _Forests creation and management_
8 | _Water bodies creation and management_
9 | _Changes of Land Cover due to natural and multiple causes_
> Les 9 _Land Cover Flows_ de l'[agence européenne pour l'environnement](https://www.eea.europa.eu/fr) (AEE)
Introduite dans une [étude de l'AEE](http://www.eea.europa.eu/data-and-maps/data/land-cover-accounts-leac-based-on-corine-land-cover-changes-database-1990-2000), cette nomenclature propose une simplification des transitions possibles en neuf grands types de _Land Cover Flows_ (ou LCF).
### Consignes
1. Charger et découper les fichiers _shapes_ correspondant aux changements d'occupation du sol sur les périodes 1990-2000, 2000-2006 et 2006-2012 (CHA...).
2. Inspecter leurs tables attributaires. Combien de transitions différentes sont présentes dans les deux couches ?
Sachant que la nomenclature a 44 postes, quel est le nombre total de transitions possibles ?
3. Joindre aux trois couches CHA... la table `LCF lookup.xls` après l'avoir inspectée sous un tableur.
4. Récapituler les surfaces des différentes LCFs sur la zone d'étude pour les différentes périodes. Qu'en conclure ?
## Exercice 6 - Construire des taches urbaines à partir de différentes bases 'OccSol'
* _Shapes_: `CLC12_RALPC_RGF.shp`, `IPLI_1977.shp`, `LittoMOS_2006.shp`
Sur le littoral français métropolitain, deux efforts exhaustif de cartographie des occupations du sol sont à noter. Ils sont tous les deux téléchargeables en accès libre :
* l'_Inventaire permanent du LIttoral_ réalisé en 1977 (ou IPLI77, voir [ici](http://www.geolittoral.developpement-durable.gouv.fr/ipli-1977-r13.html))
* le _Mode d'occupation du sol sur le littoral_, basé sur des données de 2000 et 2006 (ou LittoMOS, voir sur le [visualiseur Cartélie](http://cartelie.application.developpement-durable.gouv.fr/cartelie/voir.do?carte=litto_MOS&service=CEREMA))
1. Reprendre le _shapefile_ de CLC 2012 découpé sur votre zone d'étude.
2. La nomenclature CORINE est à trois niveaux hiérarchiques. Créer un champ reprenant uniquement le 1<sup>er</sup> niveau. Utiliser par exemple une expression de type `MyStr = Left([MyField], XX)` dans l'outil `Calculer un champ`.
3. Sélectionner les 'territoires artificialisés', exporter sous la forme d'un nouveau _shape_.
4. Charger les couches `IPLI_1977.shp` et `LittoMOS_2006.shp`. En vérifier rapidement les nomenclatures : elles sont compatibles avec CORINE mais avec un niveau hiérarchique supplémentaire. (Noter que certaines entrées de la nomenclature sont moins développées que dans Corine Land cover).
5. Construire des couches 'territoires artificialisés' sur le même principe. Superposer les trois produits, les inspecter visuellement.
***
## Exercice 7 - Construire des taches urbaines à partir d'une couche décrivant le bâti
* _Shapes_: `Bati_Eyre_2004.shp, Bati_Eyre_2014.shp`
1. Charger les deux _shapefiles_, extraits d'une couche 'bâti' de deux millésimes de la BD Topo sur la communauté de commune du Val de l'Eyre. Inspecter leurs tables attributaires.
2. Les croiser par l'outil `Union`, de manière à construire une symbologie permettant de visualiser les évolutions. Zoomer sur un quartier. Inspecter la forme des transitions: que traduisent-elles ? (Voir p. 70 de la [documentation IGN](http://professionnels.ign.fr/sites/default/files/DC_BDTOPO_2-1.pdf))
3. En se basant sur le seuil retenu par l'INSEE pour définir ses _unités urbaines (voir la définition [ici](http://www.insee.fr/fr/methodes/default.asp?page=definitions/unite-urbaine.htm)), construire des taches urbaines sur les deux dates par dilatation-érosion: `tampon` positif de la valeur du seuil divisée par deux (avec union des résultats), puis `tampon` négatif sur le résultat.
4. Refaire la même chose avec un outil dédié de l'_Arctoolbox_ : `Cartographie > Généralisation > Délimiter les zones construites`(voir [l'aide d'ArcGIS](http://resources.arcgis.com/fr/help/main/10.2/index.html#/na/007000000047000000/)). Comparer visuellement avec la méthode précédente. Croiser ces deux nouvelles couches (en copiant les symbologies déjà créées). Comparer avec le produit de base, puis avec les taches urbaines issues des bases 'OccSol'.
## Exercice 8 - Analyser les évolutions de la réserve naturelle du Banc d'Arguin
* _Shapes_: `N_ENP_RNN_S_000.shp`
* _Décrets_: docs PDF
La réserve naturelle nationale (RNN) du Banc d'Arguin, située à l'embouchure du Bassin d'Arcachon, a un nouveau décret depuis 2017, qui modifie à la fois le réglement et l'emprise de l'aire protégée. Plus de renseignements sur les sites des [réserves naturelles de France](http://www.reserves-naturelles.org/banc-d-arguin) ou de la [SEPANSO](http://www.sepanso.org/reserves/arguin_p.php).
{: .shadow}*La réserve en 2015 (source SEPANSO)*
1. Inspecter les emprises de la RNN sur le shapefile (source INPN, mars 2017).
2. Consulter les documents des décrets de 1986 et 2017.
3. Construire un polygone de la nouvelle emprise de la RNN, en se basant sur les indications du décret de 2017.
4. Télécharger une couche décrivant les emprises des réserves naturelles nationales de métropole sur le site de l'[INPN](https://inpn.mnhn.fr/telechargement/cartes-et-information-geographique/ep/rnn). Comparer avec la création précédente !
***
[](https://zenodo.org/badge/latestdoi/113714660)
(https://doi.org/10.1016/j.landusepol.2016.04.026)
<script id="asp-embed-script" data-zindex="1000000" type="text/javascript" charset="utf-8" src="https://spark.adobe.com/page-embed.js"></script><a class="asp-embed-link" href="https://spark.adobe.com/page/rtjWHN5dLMVvH/" target="_blank"><img src="https://spark.adobe.com/page/rtjWHN5dLMVvH/embed.jpg?buster=1543425105742" alt="Le marais de Reysson" style="width:100%" border="0" /></a>