Gestion des dalles raster avec GDAL

Ce document décrit différents scripts basés sur GDAL en Shell pour des dalles raster.
L’exemple d’utilisation est la récupération de données auprès de l’IGN, au format TIFF ou ECW

 

Afficher l’information sur une image raster

gdalinfo src_dataset

Affiche l’information sur un jeu de données raster. Permet notamment d’afficher la référence spatiale et l’étendue des données, le nombre de pixel, la taille du pixel. Documentation officielle gdalinfo

Affecter une référence spatiale à un ensemble de dalles

gdal_translate -a_srs srs_def src_dataset dst_dataset

Affecte le système de coordonnées « srs_def » à l’image en sortie. On utilisera la syntaxe EPSG:n (ex. EPSG :2154 pour du Lambert93 RGF93, 27572 pour du NTF Lambert II Etendu). ex. pour du NTF- Lambert II Etendu :

gdal_translate -a_srs EPSG:27572 im_src.tif im_dst.tif

Script pour un lot de dalles :

#!/bin/bash
for FILE in *.tif
do
BASE=`basename $FILE .tif`
NEWFILE=${BASE}_c.tif
gdal_translate -of GTiff -a_srs EPSG:2154 $FILE $NEWFILE
done

options de gdal_translate

  • compresser une image : l’option -co COMPRESS=DEFLATE permet de compresser l’image TIFF sans perte de qualité. L’option -co TILED=YES permet de spécifier que l’on fait un tuilage.
  • codage de l’image : l’option -co NBITS=1 permet de coder l’image en sortie à 1 bit. Documentation officielle gdal_translate

    Création d’une mosaique au format « Virtual Raster » VRT

    gdalbuildvrt mosaic.vrt *.tif

Crée un mosaïque « mosaic.vrt » à partir d’un ensemble d’une collection d’images tif. Ce format vrt est très léger et semblable aux algorithmes d’Er-Mapper, à savoir que les données d’origine sont conservées et inchangées. Le vrt est simplement un catalogue d’images. Toutes les images de la collection doivent avoir le même système de coordonnées. Pour éviter les zones noires là où il n’y a pas de données on rajoutera une option « -addalpha ». De cette manière, avec un éditeur d’image gérant le canal alpha (ex. ArcMap), on aura de la transparence dans les zones sans données source, et de l’opacité dans les zones avec données source. Ex. :

gdalbuildvrt -addalpha mosaic.vrt *.tif

Documentation officielle gdalbuildvrt

Construire des pyramides

Cette fonction permet d’accélérer l’affichage en créant des aperçus selon les niveaux de zoomes choisis. C’est l’équivalent de la création de pyramides dans ArcGIS.

gdaladdo -ro data_src 2 4 8 16 32

L’option « -ro » permet de créer un aperçu en tant qu’image externe et non inclus dans l’image tiff elle-même. Attention, si vous n’utilisez pas l’option -ro, les pyramides seront construites dans l’image d’origine et elle pourra être corrompue.
Pour compresser le plus possible en JPEG, on utilisera ce type d’options :

gdaladdo -ro --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config JPEG_QUALITY_OVERVIEW 80 data_src 2 4 8 16

Documentation gdaladdo

gdal_merge : mosaiquage

Mosaïque automatique d’un lot d’images. Toutes les images doivent être dans le même système de coordonnées et avoir le même nombre de bandes. Les images peuvent se chevaucher, et avoir des résolutions différentes. Dans les zones de chevauchement, les valeurs prendront celles de la dernière image.

gdal_merge.py [-o out_filename] [-of out_format] [-co NAME=VALUE]*
[-ps pixelsize_x pixelsize_y] [-tap] [-separate] [-v] [-pct]
[-ul_lr ulx uly lrx lry] [-n nodata_value] [-init "value [value...]"]
[-ot datatype] [-createonly] input_files

Exemple pour des dalles SRTM en Geotiff :

gdal_merge.py -o srtm_global.tif *.tif

Aide sur gdal_merge.py

Images SPOT 6-7 et Pléiades

 

Pour toute demande d’image Pléiades et SPOT 6-7 qui ne serait pas déjà accessible par la mutualisation Theia ou Geosud (voir ici), il est possible d’effectuer une demande ISIS selon la procédure décrite à ce lien :http://www.satelliteimageaccess.teledetection.fr/

L’accès à ces données est gratuit dans la limite d’une certaine emprise et payant au delà.

Dans le formulaire de demande, il faudra indiquer dans les deux premiers champs :

  • membre Geosud : OUI
  • référent Geosud : fleury-at-cerege.fr

Pour toute demande d’assistance, nous contacter sigeo@cerege.fr

Données dans Geosud et Theia

Geosud

Le CEREGE est membre de l’Equipex Geosud et du Pôle thématique Theia qui gère un ensemble des ressources satellitaires de ces structures, notamment un banque d’images Rapid Eye, SPOT, Sentinel, Pléiades ….

Le projet EQUIPEX GEOSUD vise à développer une infrastructure nationale de données satellitaires accessible gratuitement par la communauté scientifique et les acteurs publics. Site web Geosud

Tous les membres du CEREGE peuvent devenir membre de Geosud et ainsi avoir accès aux ressources. Pour s’inscrire c’est ici , en cas de question contacter le référent Geosud au CEREGE fleury-at-cerege.fr

Theia

Le pôle thématique surfaces continentales Theia (groupement de 9 institutions publiques françaises dont le CNRS) a pour vocation de faciliter l’usage des images issues de l’observation des surfaces continentales depuis l’espace. Theia met à disposition de la communauté scientifique nationale et internationale une vaste panoplie d’images à différentes échelles, de produits, de méthodes et de formation liés à l’observation des surfaces continentales en particulier depuis l’espace. Site web Theia

Accès aux données

L’accès au téléchargement pour les membres du CEREGE est décrit dans le document PDF suivant procedure_demande_Geosud_v01

La procédure de demande d’images Pléiades et SPOT 6-7 qui ne serait pas disponible au téléchargement est décrite par ailleurs dans l’article SPOT/Pléiades

Webservices IGN

Accès au webservices (WMS, WMTS, WFS, …) du Geoportail IGN

Les données IGN sont disponibles par webservices aux membres du CEREGE grâce à la licence Enseignement/Recherche (contrat de licence).

  1. documentation de base sur le site IGN
    ici
    et notamment celui pour QGIS, ceci
  2. Le serveur WMS à connecter dans votre SIG est de cette forme

on remplacera xxxx par la clé du contrat
Les codes d’accès et la clef sont à demander auprès du SIGéo

 

Serveurs WMS publics

L’accès aux données en WMS (ou WFS ou WCS) nécessite de disposer d’un client WMS (logiciel libre ou d’éditeur) : JUMP, Gv SIG, uDig, Quantum GIS, ArcGIS, Mapinfo (à partir de la version 7.8 et sup)… Après avoir lancé la commande permettant d’ouvrir une couche WMS, le logiciel demande de saisir une adresse WMS (ou WFS ou WFS).

Serveurs WMS
Organisme/Projet Portail WMS WFS WCS Description
IECDS ICEDS http://iceds.ge.ucl.ac.uk/cgi-bin/i… ? http://iceds.ge.ucl.ac.uk/cgi-bin/i… ? Données d’imagerie (principalement Europe et Afrique)
OneGeology OneGeology données géologiques sur le monde entier
BRGM Portail http://geoservices.brgm.fr/geologie, ou ancien WMS http://ogcpublic.brgm.fr/geologie http://geoservices.brgm.fr/geologie Géologie du BRGM
NASA OnEarth http://wms.jpl.nasa.gov/wms.cgi SRTM, Landsat, …Service ne semblant plus opérationnel
Forum SIG Forum SIG Liste de serveurs WMS
Eduterre Eduterre Liste de serveurs WMS
Georezo Georezo Liste de serveurs WMS
Demis Demis http://www2.demis.nl/worldmap/wms.asp ? Serveur WMS de la société Demis, des fonds de référence et notamment un beau fond hybride sur le monde.
wms-sites wms-sites Catalogue de WMS
Skylab Public OGC WMS server list Catalogue de WMS
Geosignal Geosignal http://wms.geosignal.fr/metropole ? ou http://www.geosignal.org/cgi-bin/wmsmap ? 7 niveaux de cartes raster, représentant des fonds cartographiques basiques , les plans des agglomérations de plus de 10000 habitants, les limites administratives, le réseau routier national et départemental.
CORINE LAND COVER CLC http://sd1878-2.sivit.org/geoserver/wms Occupation du sol, couverture européenne
GeoLittoral Geolittoral http://geolittoral.application.equi… ? Des informations relatives aux littoraux français et bien d’autres informations.
Gest’Eau Gest’eau http://gesteau.oieau.fr/service Des informations relatives aux SAGES notamment et bien d’autres informations.
Sandre Sandre Ouvrages : http://services.sandre.eaufrance.fr… La référence pour les données relatives à l’eau, il est possible d’accéder à un volume très important de référentiels et de bases géoréférencées.
Sandre Sandre Zonages : http://services.sandre.eaufrance.fr…
Sandre Sandre Stations : http://services.sandre.eaufrance.fr…
CRIGE-PACA CRIGE-PACA Geoservices WMS ortho : http://www.crige-paca.org… Un très bon serveur WMS notamment pour les orthos THR sur PACA

Données spatiales en ligne

Différents portails donnent accès à des données gratuites sur le monde. Données d’imagerie, fond de référence, MNT, données thématiques.
Ci-dessous une liste de portails.

Pour commencer un lien vers des données globales vectorielles, raster et thématiques assez complet : Global Datasets

Photos anciennes IGN

Ce lien explique tout : ici

NaturalEarth

Des données de référence raster et vecteur sur le site de NaturalEarth

Données du 1 : 10M au 1 : 110M

USGS Store

Des cartes sur le site de l’USGS Store

GLCF

GLCF : Catalogue de données Landsat, Aster, AVHRR, Land Cover, SRTM

Geogratis

Geogratis Des données principalement sur le Canada

Images SPOT

Catalogue d’images SPOT : SIRIUS (payantes)
Programme ISIS pour obtenir des données à tarif préférentiel dans le cadre d’une programme de recherche.
Des images SPOT d’exemple gratuites, pour chaque académie (nécessite de créer un compte). Voir ici

ResMap

Données d’imagerie gratuites (1984 à 2004). Voir le site

EROS

Earth Resources Observation and Science (EROS) Center. Voir le site
Catalogue de données de l’USGS.
Voir aussi le catalogue Glovis ou le EarthExplorer

SRTM

Le MNT SRTM version 4 est disponible sur CGIAR par dalles Geotiff de 5° sur 5°.
Site de référence du SRTM au JPL NASA.

Ionia Globcover

Carte d’occupation du sol 2008 mondiale à 300m de résolution par l’ESA. Site Globcover Ionia.
Des mosaïques MERIS sont également disponibles sur le site.

Liste issue du site Diva-gis

Diva-gis

Des données mondiales :

  • de référence (frontières, …)
  • altimétriques (MNT SRTM)
  • imagerie (Landsat)
  • climato
  • espèces
  • autres…

Vecteurs GPS pour ArcGIS

Note préliminaire

ToolBox ArcGis réalisée pour ArcGIS 9.x et 10.0. Non testé avec les versions ultérieures

Construction des vecteurs GPS et des ellipses d’incertitudes.

Boîte à outil ArcGIS

Nécessite en entrée une table contenant au minimum : “Longitude” “Latitude” “Deplacement Est” “Déplacement Nord” ” Incertitude Est” “Incertitude Nord”

INSTALLATION :

- Ajouter la boîte à outils correspondant à votre version dans ArcToolbox

UTILISATION :

- Suivre les étapes 0 à 5 (Si vous disposez déja d’un tableau avec coordonnées projetées, vous pouvez sauter l’étape 1)

- à partir de l’étape 2, utiliser le fichier projeté

Exemple de résultat :

JPEG - 35.5 ko
Vecteurs GPS

TELECHARGEMENT Vecteurs_GPS

BeachBalls pour ArcGIS

Note préliminaire

Macro créée en VB en 2011 avec ArcGis 9.x et 10.0. Non testé avec les versions récentes

Représentation des mécanismes au foyer dans ArcGis

Macro VBA Beach Balls
Macro à installer dans ArcMap

Mécanismes au foyer des séismes en beach balls

Archive à télécharger: Beach_Ball

Programmation VB et VBA dans ArcGIS

Note préliminaire

ceci est quasiment obsolète et je le conserve pour mes propres utilisations.

Un utilisateur qui voudrait scripter ArcGis devra plutôt utiliser Python.

Introduction

L’intégration de VBA dans la version 10 n’est plus assurée. Nous vous conseillons de vous orienter vers des scripts Python pour la Toolbox.

Le document ci-dessous est une introduction à la programmation pour les versions inférieures ou égales à la 9.3.1


Dans de nombreux cas, il est nécessaire de développer soi-même une application qui n’est pas implémentée dans le SIG :

  1. pour une fonction qui n’existe pas en standard ;
  2. pour une opération fastidieuse en standard ;
  3. pour une opération répétitive.

Pour cela, nous utilisons les outils de développement associés à ArcGIS. VBA intégré jusqu’à la version 9.3.1, VB si l’on dispose d’un compilateur. Pour la version 10 et supérieur d’ArcGIS, le support de VBA n’est plus assuré et il faudra préférer utiliser le langage Python, notamment via les scripts ArcToolbox.

Cette documentation est une introduction au développement VBA et VB pour ArcGis.

Doc pdf

documentation et exemples dans ce fichier DOC_CODE_ARCGIS_v0