Recette Hillshading sur opensnowmap, pour noter çà quelque part et si
certains son d'humeur joueuse.
Rien de particulier, si ce n'est ce réglage là avec imagemagick, qui ne
marche pas sur des tifff énormes, bien sûr :
composite -compose Multiply hillshade.tif slopes.tif composed.tif
mogrify -level 28%x70% composed.tif
Yves
gdalbuildvrt -srcnodata -32768 DEM.vrt tiles/*.tif
gdal_translate -of GTiff -co "BIGTIFF=YES" -co "TILED=YES" DEM.vrt DEM.tif
gdaldem hillshade -of GTiff -co "BIGTIFF=YES" -co "COMPRESS=PACKBITS"
-z 2 -s 111120 DEM.tif HILLSHADE.tif
gdaldem slope -of GTiff -co "BIGTIFF=YES" -co "COMPRESS=PACKBITS" -s
111120 DEM.tif TMP.tif
gdaldem color-relief -of GTiff -co "BIGTIFF=YES" -co
"COMPRESS=PACKBITS" TMP.tif color_slope.txt SLOPES.tif
for a in {-18..17}
do
for y1 in {-60..71}
do
x1=$(echo $a*10 | bc)
x2=$(echo $x1+10.02 | bc)
y2=$(echo $y1+1.02 | bc)
#~
echo $x1 $y2 $x2 $y1
name=E$x1"N"$y1.tif
echo $name
gdalwarp -dstnodata 255 -of GTiff -te $x1 $y1 $x2 $y2 HILLSHADE.tif
hillshade.tif
gdalwarp -dstnodata 255 -of GTiff -te $x1 $y1 $x2 $y2 SLOPES.tif
slopes.tif
composite -compose Multiply hillshade.tif slopes.tif composed.tif
mogrify -level 28%x70% composed.tif
./gdalcopyproj.py hillshade.tif composed.tif
mv composed.tif hs-tiles/$name
rm hillshade.tif slopes.tif
done
done
# sort out paletted (empty) files
for x in hs-tiles/*.tif
do
if [ $(gdalinfo $x | grep -c "ColorInterp=Palette") -ne 0 ]
then
mv $x tmp/
fi
done
gdalbuildvrt -hidenodata -vrtnodata "255" done.vrt hs-tiles/*.tif
gdal_translate -of GTiff -co "BIGTIFF=YES" -co "TILED=YES" done.vrt
wgs84.tif
gdalwarp -of GTiff -co "BIGTIFF=YES" -co "COMPRESS=PACKBITS" -co
"TILED=YES" \
-t_srs "+init=epsg:3857 +over" \
-r cubic -order 3 -multi wgs84.tif 3857.tif
_______________________________________________
dev-fr mailing list
[email protected]
http://lists.openstreetmap.org/listinfo/dev-fr
