====== Compression des tif de la BD Parcellaire (noir et blanc) ======
//D'après [[http://georezo.net/forum/viewtopic.php?id=70210|cette discussion]] exhumée par David Michallet du [[http://avenir.38.free.fr/|Conservatoire des espaces naturels de l'Isère]]//
===== En DOS (windows) =====
<code dos>
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"
</code>
===== Script Bash (Linux) =====
<code bash>
#!/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
</code>
====== Création d'un fond en relief pour la région Languedoc-Roussillon ======
//inspiré de [[http://linfiniti.com/2010/12/a-workflow-for-creating-beautiful-relief-shaded-dems-using-gdal/|l'article de Tim Sutton]]//

{{:outils:gdal:mnt_lr_final.png?300|Le résultat obtenu}}
===== données disponibles =====
les MNT de la BD TOPO en Lambert 93
==== 1 - Fusionner les 5 départements ====
<code bash>
gdal_merge.py -o mnt_lr.tif -of GTiff DEPT48.asc DEPT30.asc DEPT34.asc DEPT11.asc DEPT66.asc 
</code>
==== 2 - Création d'une image en couleur du MNT ====
On utilisera ici le fichier de couleur "srtm" de grass dont voici le contenu :
<code>
-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
</code>
<code bash>gdaldem color-relief -of GTiff mnt_lr.tif /usr/lib/grass-6.3.0/etc/colors/srtm mnt_lr_couleur.tif</code>
==== 3 - Création d'un relief ombré====
<code bash>gdaldem hillshade -z 5 -of GTiff mnt_lr.tif mnt_lr_ombre.tif</code>
==== 4 - Récupération du script hsv_merge.py  ====
<code bash>wget http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/hsv_merge.py</code>
==== 5 - Permettre l’exécution du script ====
<code bash>chmod 775 hsv_merge.py</code>
==== 6 - Fusion des deux images (ombre et couleur) ====
<code bash>./hsv_merge.py mnt_lr_couleur.tif mnt_lr_ombre.tif mntlr_couleur_ombre.tif</code>
==== 7 - découpage selon limites régionales ====
=== création si nécessaire du contour régional ===
<code bash>pgsql2shp -f region_lr.shp -h 192.168.1.230 -u mathieu sicen "SELECT 'lr' as nom, st_union(buffer(geometrie,1000)) from commune"</code>
=== découpage à proprement parler ===
<code bash>gdalwarp -s_srs EPSG:2154 -co compress=deflate -dstnodata 255 -cutline ./region_lr.shp mntlr_couleur_ombre.tif mnt_lr_final.tif</code>
==== 8 - Transformation en jpeg ====
<code bash>gdal_translate -of jpeg mnt_lr_final.tif mnt_lr_final.jpeg</code>
===== Le script complet =====
<file sh 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
</file>