geotrucs

Le format gdb est idéal tant qu’on reste dans le monde esri, il reste accessible dans les versions récentes de qgis (via gdal/ogr) moyennant certaines limitations. Par contre, impossible d’accéder aux données attributaires avec access, encore moins excel ou libre office.

Le format mdb est interopérable entre arcgis, qgis et access/libre office mais son gros défaut est que c’est un format obsolète qui n’est plus supporté par microsoft (son inventeur), qu’il comporte des limitations en termes de taille et qu’il sature rapidement la RAM en cas d’utilisation intensive.

Le problème, c’est donc de trouver un format moderne et largement diffusé qui combine les avantages de performance et d’ergonomie du format gdb, et l’interopérabilité du format mdb.

la solution: le geopackage

Le géopackage est un format de données géographique défini par l’OGC (donc libre et ouvert) basé sur le format de base de données sqlite.

Un fichier sqlite contient à la fois une base de données et le moteur de base de données qui permet d’exploiter la donnée, le tout totalement compatible SQL. Le format sqlite est très discret mais universellement répandu grâce à sa capacité à fonctionner sous n’importe quelle infrastructure matérielle et n’importe quel OS. On le trouve dans firefox, chrome, android... C’est donc le format interopérable par excellence.

Un geopackage c’est simplement une base sqlite qui contient une cartouche spatiale permettant de gérer des données géographiques vectorielles. La donnée attributaire attachée reste simplement au format sqlite classique (même principe que la donnée attributaire dans une mdb).

Pour plus de détails:

http://www.geopackage.org/

https://www.sqlite.org/index.html

https://en.wikipedia.org/wiki/GeoPackage

https://fr.wikipedia.org/wiki/SQLite

Dans la suite de cette note, on démontrera toute la chaîne d’opérations depuis la création jusqu’à l’édition de données dans un geopackage avec tous les outils à notre disposition: arcgis pro/desktop, qgis 3, access et libre office.

Prérequis

arcgis desktop/pro: aucun prérequis, les outils sont disponibles en standard

qgis: idem en version 3, interface presse-bouton

ms access: nécessite le driver sqlite 3.0 pour accéder à la donnée attributaire (http://www.ch-werner.de/sqliteodbc/)

libre office: idem, installation du même driver

Créer un geopackage

avec arcgis desktop/pro: data management tools/workspace/create SQLite Database. Bien spécifier le type spatial: GEOPACKAGE_1.2

avec qgis 3: couche>gestionnaire de sources de données>geopackage>nouveau (et plusieurs autres endroits dans l’interface)

avec access/libre office, on créera une base sqlite, mais sans cartouche spatiale

lire/écrire des données géo dans un geopackage

avec arcgis desktop/pro: le geopackage est un workspace comme n’importe quel autre. On peut enregistrer une couche, créer une nouvelle couche, utiliser le geoprocessing copy features exactement comme avec une gdb. Attention: l’écriture dans un gpkg nécessite une licence basic

avec qgis 3, c’est le même principe. Le geopackage est une source comme une autre, on la manipule exactement de la même manière

Réponse longue: voir ci-dessous

Space Shuttle Radar Topography Mission (SRTM)

Produit par la Nasa en 2000, la résolution est constante: 30 m sur toute la couverture (mondiale) depuis 2015. Le capteur est un radar inSAR (https://en.wikipedia.org/wiki/Interferometric_synthetic-aperture_radar). Ce qu’il faut savoir de ce capteur, c’est qu’il peine dans les zones de fort relief en bordure de scène: effets de distorsions ou zones complètement noires.

Téléchargement en dalles de 1 deg ici: http://earthexplorer.usgs.gov/ (sur inscription gratuite)

ASTER Global Digital Elevation Model

NASA - Japon: couverture mondiale avec résolution de base de 30 m. La technologie utilisée est une corrélation stéréoscopique d’images satellitaires dans le visible et le proche infrarouge. Une V2 sortie en 2011 corrige beaucoup de défauts de jeunesse bien que les images présentent toujours beaucoup de bruit. Ce bruit est gênant en zone de plaine, mais ASTER est considéré comme plus fiable que SRTM en zone accidentée.

Téléchargement en dalles de 1 deg ici: http://reverb.echo.nasa.gov/reverb/

JAXA’s Global ALOS 3D World

MNT mondial à 30 m dont la dernière version est sortie en mars 2017 comme une version dégradée de AW3D, un MNT mondial à 5 m (payant bien sûr). ALUS n’apporte pas d’amélioration de résolution par rapport à ses deux prédécesseurs, mais le gain de qualité est incontestable.

Au démarrage d’arcgis pro, on choisit un projet (nouveau ou en cours) et toute la session se déroulera dans ce projet. La doc officielle sur la notion de projet dans pro: https://pro.arcgis.com/fr/pro-app/help/projects/what-is-a-project.htm

En bref: qu’est-ce qui change par rapport aux mxd

Un mxd, comme un wor (mapinfo), contenait tous les éléments nécessaires à la production d’une carte papier:

une mise en page avec ses éléments graphiques

un, ou plusieurs blocs carte, qui s’insèrent dans cette mise en page

des appels à un ensemble de données (shp, jointures...) ainsi que des définitions de symbologies

Point de détail, mais important, un mxd était un fichier binaire qui n’était modifiable (en dehors d’arcgis) qu’au travers d’une api arcpy. Pour faire simple, un mxd = une carte papier. Dans une étude, on multiplie les mxd.

Un aprx (projet arcgis pro) c’est bien plus que ça et c’est conçu pour coller à la notion de projet sig avec des cartes papier, des cartes web, des métadonnées, des process, des tâches... Les éléments du projet sont regroupés dans le panneau “catalogue”:

mises en page (au pluriel)

cartes

boîtes à outils

bases de données (locales et connexions)

dossiers

styles

géocodeurs

tâches

De nouveaux éléments font leur apparition: des tâches (guides détaillés de process), des notifications, des boîtes à outils. L’objectif est de favoriser la collaboration, la mutualisation.

Techniquement, comment ça se présente

Un dossier projet (ici __nada__ qui est mon projet par défaut) contient une gdb, une boîte à outils, un fichier aprx (définition de projet), un dossier de log d’imports et un historique des géotraitements (.pyHistory).

Lorsqu’on ajoute une donnée distante dans le projet, on ajoute au projet une connexion à une gdb, à un dossier et ces connexions sont référencées. On peut ajouter des dossiers à l’arborescence de base du projet (ex: pdf, jpg, sources...).

On voit qu’on s’achemine vers une normalisation des dossiers étude (un de mes dadas !).

32 ou 64 bits?

64 bits !

En 2017 tous les processeurs sont multicoeurs. On choisira donc par défaut un installateur 64 bits, sauf contexte très particulier (machine virtuelle dotée d'un seul coeur et d'un OS 32 bits par exemple)

installateur indépendant ou réseau?

Installateur réseau !

Pour une installation simple sur un poste de travail on a spontanément tendance à choisir l'installateur indépendant. A mon avis, dans un cadre professionnel, il vaut mieux choisir l'installateur réseau qui fonctionne aussi très bien sur un poste indépendant et qui offre en plus des possibilités d'installation de dépendances, de logiciels complémentaires et surtout de mises à jour très intéressantes.

LTR ou latest ?

Lastest !

La version LTR actuelle est la 2.14, tandis que la dernière version stable est la 2.18.9. La LTR est réputée particulièrement stable et conseillée dans le cadre d'une utilisation professionnelle en production. Cependant le développement de qgis est rapide et de nombreuses nouvelles fonctionnalités sont régulièrement ajoutées dans la version à jour. Pour moi le bénéfice - risque penche largement en faveur de la dernière version, d'autant plus que les bugs rédhibitoires restent rares.

Comment installer qgis avec l'installateur réseau ?

télécharger l'installateur et le lancer, de préférence en tant qu'admin

choisir express desktop install

laisser tout par défaut et accepter les licences proposées

L'installateur installe qgis, grass avec son plugin, saga-gis et un certain nombre de dépendances, notamment gdal-fileGDB qui offre un acces rw aux géodatabases fichier d'esri.

Comment mettre à jour qgis avec l'installateur réseau ?

relancer l'installateur (dans le dossier SSGeo4W) et choisir l'installation avancée, toujours en tant qu'admin

Mais la superficie n'est pas le seul critère: une bande continue de 100 km par 3 m le long d'une autoroute ne fait pas un habitat utile pour un lynx.

Vient donc un deuxième critère qu'on pourrait définir comme la compacité: intuitivement un rapport périmètre/aire. Malheureusement celui-ci varie en fonction de la taille de l'objet: on préfèrera

q = 4piA/P²

qui est indépendant de la taille de l'objet:

pour un disque q = 1

pour un triangle équilatéral q =

pour un carré q = 1/4pi

pour un rectange q = 4piA/P²

q est en fait le quotient isopérimétrique: https://fr.wikipedia.org/wiki/Isop%C3%A9rim%C3%A9trie

Avec un couple de critères A et q, on peut éliminer facilement les polygones les moins compacts, soit ceux dont le rapport du périmètre à l'aire sont les plus éloignés du disque.

Mais ce quotient n'est pas très parlant. On peut en dériver un indicateur plus beaucoup plus explicite: la "largeur standard".

Assimilons le polygone à un rectangle et définissons le couple longueur (L) et largeur standards (W) tels que L et W vérifient

P = 2(L + W)

A = LW

La résolution du polynôme du second degré passe par la recherche du discriminant (Attention programme de math en seconde) :

Δ = 1/4P² - 4A

Rq: si Δ <= 0, le polygone n'est pas assimilable à un rectangle: il n'y a pas de couple L et W qui vérifient les propriétés d'aire et de périmètre énoncées plus haut. Dans ce cas

on en déduit

L = 1/4P - 1/2√Δ

et bien sûr

W = 1/2P-L

On des longueurs, beaucoup plus faciles à comprendre. Par exemple, on peut définir que la largeur standard limite est de 20 m pour une parcelle constructible, en plus du seuil de 800 m², ou que la largeur standard limite est de 5 km pour le territoire d'un lynx.