Action: When run, this country.makefile [download][1] the raster GIS DEM sources, process them (unzip, crop, shaded relief, resize, colorize), to output 6 elegant shaded reliefs. Can work for any country or area, given you provide WNES (West/North/East/South) geocoordinates, in decimal degrees.
Output:
For similar data and similar px dimensions, file sizes are .tif:5.0MB, .png:1.6MB, .jpg:239KB. Also, whenever possible and relevant, we use the lighter format.
Direct command:
make -f country.makefile ITEM=India WEST=67.0 NORTH=37.5 EAST=99.0 SOUTH=05.0 WIDTH=1280Parameters:
Derivated from man gdal and man convert.
Z= (zFactor, integer >0, default 5): vertical exaggeration used to pre-multiply the elevationsAZ= (azimuth, integer [0-359], default: 315): azimuth of the light, in degrees. 0 if it comes from the top of the raster, 90 from the east, … The default value, 315, should rarely be changed as it is the value generally used to generate shaded maps.ALT= (altitude of light, in degrees, integer [0-90], default 60): altitude of the light, in degrees. 90 if the light comes from above the DEM, 0 if it is raking light.FUZZ= (fuzzy selection, integer [0-100], default 7): colors within this distance are considered equalSHADOW= (integer [0-100]) opacity of the hillshade upon color-relief.WIDTH=… (in px, default: 1280) : width of the bitmaps (tif, png). The EIGHT is conserved from WNES ratio and the WIDTH.Note: if the input GIS raster is in feet, then s scale should be edited. See man gdal.
[1]: the download task will automatically download the ETOPO1 (342MB) archive file, which may take from few minutes to a dozen of hours, depending on your web line.
Suggestions? Drop me a suggestion/command line @hugo_lz
<!DOCTYPE html>
<meta charset='utf-8'>
<style>
/* Style here... */
path#IND {
fill: none;
stroke: #000;
}
image.bg {
opacity: 0.2;
}
#legend, #image { border: .2em solid #C0C0C0; background: #EFEFEF; margin: 1em 8em 1em 8em; padding:0.5em; box-shadow: .1em .2em .2em .1em grey ; }
a { color: #0b5b0c;}
.code { color: white; background: #333; border: 2px solid #C0C0C0; padding: 0.4em;font-weight: 900; }
.code:first-letter { font-size:1.5em;color: #0b5b0c; }
</style>
<div id="legend" align="center" width="60%">
<big><p>Six <b>shaded reliefs</b> (aka <b>hill-shades</b>, <b>color-relief</b>) generated by single command.</p></big>
<small>
<p class="code">$make -f country.makefile ITEM=India WEST=67.0 NORTH=37.5 EAST=99.0 SOUTH=05.0 WIDTH=1280</p>
<p>Short log of past 2 weeks <a href="//twitter.com/wikiatlas">Wikimaps Atlas</a> shaded relief trials-and-errors explorations.</p>
<p>The last one named "color+hillshade_multiply" is believed to be the most elegant</p>
<p>All these layers are generated from <i>A-to-Z</i> via command lines, from the makefile.</p>
<p>Some dependencies required. Follow <a href="//twitter.com/wikiatlas">@WikiAtlas</a>/<a href="//twitter.com/hugo_lz">Hugo_lz</a>.</p></small>
</div>
<div id="image" align="center" width="60%">
<img src="./wikiatlas_relief_1280.gif" href="./wikiatlas_relief_1280.gif" width="60%"></img>
</div>65535 255 255 255
8000 224 222 216
7000 202 195 184
6000 186 174 154
5000 172 154 124
4000 170 135 83
3500 185 152 90
3000 195 167 107
2500 202 185 130
2000 211 202 157
1500 222 214 163
1000 232 225 182
600 239 235 192
400 225 228 181
200 209 215 171
100 189 204 150
50 168 198 143
20 148 191 139
1 172 208 165
0 216 242 254
-20 216 242 254
-50 198 236 255
-100 185 227 255
-250 172 219 251
-500 161 210 247
-1000 150 201 240
-2000 141 193 234
-4000 132 185 227
-7000 121 178 222
-15000 113 171 216#---- RUN
# make -f country.makefile ITEM=India WEST=67.0 NORTH=37.5 EAST=99.0 SOUTH=05.0
#---- DEFAULT VALUES (customizable):
WIDTH=1280
FUZZ=7
AZ=315
Z=5
SHADOW=50
#---- MAKEFILE
#---- End here
all: progressive_transparency merge_relief-color clean
#----PROCESSING RASTER DATA
merge_relief-color: GIS_color_relief progressive_transparency
composite -dissolve $(SHADOW) shadedrelief.trans.png color_relief-wp.jpg -alpha Set color_relief-hillshade-wp-0$(SHADOW).jpg #note: perfect
composite -dissolve 100 shadedrelief.trans.png color_relief-wp.jpg -alpha Set color_relief-hillshade-wp-100.jpg #note: perfect
convert color_relief-wp.jpg shadedrelief.trans.png -compose Multiply -composite color_relief-hillshade-wp-multiply.jpg #note: perfect
GIS_color_relief: resize
# $man gdaldem : color-relief
# Color tiff depending on color_relief.txt file. Format: elevation R G B.
# Elevation as floating point value, `nv` keyword, or percentage.
# 0% being the minimum value found in the raster, 100% the maximum value.
gdaldem color-relief crop.tiff color_relief-wikimaps.txt color_relief-wp.tiff #GIS file
#resizing
convert color_relief-wp.tiff -resize $(WIDTH) color_relief-wp.jpg #tiff:5.0MB, png:1.6MB, jpg:239KB
#--- Trans shades
progressive_transparency: grey_wiping
convert shadedrelief.grey_no.png -alpha copy -channel alpha -negate +channel shadedrelief.trans.png
grey_wiping: resize
convert shadedrelief.sized.tiff -fuzz $(FUZZ)% -transparent "#DDDDDD" shadedrelief.grey_no.png
#---- Crop, GIS Hillshade, Resize
resize: shade
#How to resize propotionally and keeping data / valid GIS tiff ?
convert shadedrelief.tiff -resize $(WIDTH) shadedrelief.sized.tiff
#gdal_translate crop.tiff -outsize $(WIDTH) 1200 resized.gdal.tif
shade: crop
gdaldem hillshade crop.tiff shadedrelief.tiff -z $(Z) -s 111120 -az $(AZ) -alt 60 -compute_edges
crop: unzip
gdal_translate -projwin $(WEST) $(NORTH) $(EAST) $(SOUTH) ETOPO1_Ice_g_geotiff.tif crop.tiff #ETOPO1_Ice_g_geotiff.tif
# ulx uly lrx lry (geodegrees) // W N E S #todo: add name parameter
#---- DOWNLOADS
unzip: download
unzip -n ../data/ETOPO1/ETOPO1.zip '*.tif'
touch ETOPO1_Ice_g_geotiff.tif
download: clean
mkdir -p ../data/ ../data/ETOPO1
curl -o ../data/ETOPO1/ETOPO1.zip -C - 'http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/ice_surface/grid_registered/georeferenced_tiff/ETOPO1_Ice_g_geotiff.zip'
echo "[fake] download by shadedrelief: done!"
.PHONY: clean
clean:
rm -f *.tiff
rm -f *.jpg
rm -f *[a-km-z0-9].png #keeps the thumbnail.png
# rm -f *.tiff
# rm -f *.dbf
# rm -f *.prj
# rm -f *.shp
# rm -f *.shx
# rm -f *.html
# rm -f *.txt
# rm -f *.json