Compression des tif de la BD Parcellaire (noir et blanc)
- En DOS (windows)
- Script Bash (Linux)
Création d'un fond en relief pour la région Languedoc-Roussillon
- données disponibles
- Le script complet

Compression des tif de la BD Parcellaire (noir et blanc)

D'après cette discussion exhumée par David Michallet du Conservatoire des espaces naturels de l'Isère

En DOS (windows)

FOR %Q IN (*.tif) DO gdal_translate -co TFW=YES -co COMPRESS=CCITTFAX4 -co TILED=YES -co NBITS=1 -a_nodata 0 "%Q" "%~dpnQ_NEW.tif"

Script Bash (Linux)

#!/bin/bash
for tif in `ls -R *.tif`; do
	eval gdal_translate -co TFW=YES -co COMPRESS=CCITTFAX4 -co TILED=YES -co NBITS=1 -a_nodata 0 $tif $(echo $(basename $tif) | cut -d"." -f1 | tr '[A-Z]' '[a-z]')_compresse.tif
done

Création d'un fond en relief pour la région Languedoc-Roussillon

inspiré de l'article de Tim Sutton

données disponibles

les MNT de la BD TOPO en Lambert 93

1 - Fusionner les 5 départements

gdal_merge.py -o mnt_lr.tif -of GTiff DEPT48.asc DEPT30.asc DEPT34.asc DEPT11.asc DEPT66.asc

2 - Création d'une image en couleur du MNT

On utilisera ici le fichier de couleur “srtm” de grass dont voici le contenu :

-500 0 0 10
-300 0 0 20
-200 0 0 70
-100 0 0 130
-50 0 0 205
0 aqua
0.1   57 151 105
100  117 194 93
200  230 230 128
500  202 158 75
1000 214 187 98
2000 185 154 100
3000 220 220 220
5000 250 250 250
8850 255 255 255
nv white

gdaldem color-relief -of GTiff mnt_lr.tif /usr/lib/grass-6.3.0/etc/colors/srtm mnt_lr_couleur.tif

3 - Création d'un relief ombré

gdaldem hillshade -z 5 -of GTiff mnt_lr.tif mnt_lr_ombre.tif

4 - Récupération du script hsv_merge.py

wget http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/hsv_merge.py

5 - Permettre l’exécution du script

chmod 775 hsv_merge.py

6 - Fusion des deux images (ombre et couleur)

./hsv_merge.py mnt_lr_couleur.tif mnt_lr_ombre.tif mntlr_couleur_ombre.tif

7 - découpage selon limites régionales

création si nécessaire du contour régional

pgsql2shp -f region_lr.shp -h 192.168.1.230 -u mathieu sicen "SELECT 'lr' as nom, st_union(buffer(geometrie,1000)) from commune"

découpage à proprement parler

gdalwarp -s_srs EPSG:2154 -co compress=deflate -dstnodata 255 -cutline ./region_lr.shp mntlr_couleur_ombre.tif mnt_lr_final.tif

8 - Transformation en jpeg

gdal_translate -of jpeg mnt_lr_final.tif mnt_lr_final.jpeg

Le script complet

creation_mnt_couleur_ombre.sh

gdal_merge.py -o mnt_lr.tif -of GTiff DEPT48.asc DEPT30.asc DEPT34.asc DEPT11.asc DEPT66.asc 
gdaldem color-relief -of GTiff mnt_lr.tif ./srtm mnt_lr_couleur.tif
gdaldem hillshade -z 5 -of GTiff mnt_lr.tif mnt_lr_ombre.tif
./hsv_merge.py mnt_lr_couleur.tif mnt_lr_ombre.tif mntlr_couleur_ombre.tif
gdalwarp -s_srs EPSG:2154 -co compress=deflate -dstnodata 255 -cutline ./region_lr.shp mntlr_couleur_ombre.tif mnt_lr_final.tif
gdal_translate -of jpeg mnt_lr_final.tif mnt_lr_final.jpeg
rm -f mnt_lr_couleur.tif
rm -f mnt_lr_ombre.tif
rm -f mntlr_couleur_ombre.tif
rm -f mnt_lr_final.tif