Author: Jules FLEURY

Correction différentielle des données GPS ArcPad

Introduction

Cette procédure traite de la correction différentielle des associés au fichier SSF, acquis avec des appareils Trimble Juno SB, équipés des logiciels ArcPad8 et GPSCorrect.

Etapes

0.Le fichier SSF ne contient que des positions et non de l’information attributaire. Les informations attributaires sont contenues dans les fichiers SHP ou les AXF.

Pour réaliser la correction différentielle, vous devrez corriger le fichier SSF et ensuite réaliser la synchronisation entre le fichier corrigé COR et le fichier SHP ou AXF.

1. Copier les données de votre PAD (ex. Trimble Juno SB) sur le poste local. Mettre le fichier .ssf et les shapefiles dans un même répertoire.

2. Pour effectuer la correction différentielle, la méthode est d’utiliser l’assistant de correction de GPS Pathfinder Office disponible par le menu “Outils/Correction différentielle”. Vous obtenez alors un fichier GPSCorrect.cor.

A noter dans cet assistant, l’importance d’inclure dans le fichier “Cor” à la fois les positions corrigées et non corrigées (paramètres disponibles dans l’assistant, à la fenêtre “Paramètres de correction”, en cliquant sur “Changer”, dans l’onglet “Sortir”). Ce paramètre permet d’obtenir une correspondance exact entre le fichier SSF et le fichier COR.

3. Avant de procéder à la synchonisation, assurez-vous que les fichiers SHP (DBF, SHX….ou AXF seul) se trouvent dans un répertoire spécifique, dans lequel se trouvent également les fichiers SSF et COR….

4. Aller ensuite dans le dossier “Outils/Autre/ShapeCorrect” afin de lancer l’assistant de synchronisation.

6. Dans “ShapeCorrect”, cliquer sur “Parcourir” afin de sélectionner les fichiers SHP (ou le fichier AXF). Il n’est pas nécessaire de sélectionner les fichiers SSF et/ou COR, mais ils se doivent d’être présents dans le meme répertoire.

7. C’est à vous de décider si vous souhaitez dans les fichiers résultants les positions corrigées et/ou non corrigées.

8. Cliquer sur “OK” afin de lancer la synchro. Vérifiez le résultat de la synchronisation.


	


	





	

		
		
Traiter les dalles SRTM avec GDAL
	


	
	
	

		
Cette procédure décrit la méthode pour télécharger et prétraiter les dalles SRTM v4.


Téléchargement :


se connecter sur le serveur CGIAR dans la partie SRTM Data Search and Download et choisir une ou plusieurs dalles de 5°x5°


on obtient une archive contenant 4 fichiers (ex. readme.txt ; srtm_41_05.hdr ; srtm_41_05.tif ; srtm_41_05.tfw)


Mosaïque :


Déplacer tous les fichiers dans un seul répertoire.

Puis lancer la commande :

gdal_merge.py -n -32768 -o out_filename input_files


- n -32768 : pour fixer les valeurs nodata

out_filename : nom du fichier en sorite

input_files : liste de fichiers en entrée


Plus d’options GDAL sur http://www.gdal.org/index.html


Découpage d’une zone


gdal_translate -projwin ulx uly llx lly) src_dataset dst_dataset

ulx uly llx lly sont les coordonnées projetées de la zone à découper, ex. en WGS84 sur la Crête 23 36 27 34

src_dataset : fichier en entrée

dst_dataset : fichier en sortie


Plus d’options GDAL sur http://www.gdal.org/index.html


Ombrage (ou estompage)


gdaldem hillshade -s 111120 input_dem output_hillshade


- s 111120 : ratio pour les unités verticales dans le cas où les coordonnées sont en degrés (e.g. Lat/Long WGS84) et les altitudes en mètres.


Plus d’options GDAL sur http://www.gdal.org/index.html

	


	





	

		
		
Dallage avec GDAL
	


	
	
	

		
Générer des tuiles avec gdal2tiles.py


création d’un tuilage avec cet utilitaire gdal ici


Description d’une procédure complète pour créer un tuilage et le mettre au format KML de Google Earth ici


Re-tuiler des tuiles existantes


utilitaire ici


exemple pour générer les tuiles en JPEG, compressées à un taux de 80, à partir d’une image ECW


gdal_retile.py -of JPEG -ps 2000 2000 -s_srs EPSG:27572 --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config JPEG_QUALITY_OVERVIEW 80 -tileIndex TI_ortho_2003.shp -tileIndexField TileName -targetDir /tmp/tileortho/ /mnt/srvstk0_geomor/Data_SIG/OHM_Mine/Fonds_Reference/Ortho_2003/Ortho_2003_L2e.ecw


voir la documentation sur le site GDAL (rubrique Utilities)

	


	





	

		
		
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