[gdal-dev] Re: Hillshade + Topographic Map??
jluis at ualg.pt
Tue Feb 21 10:11:44 EST 2012
On 21-02-2012 14:11, Rutger wrote:
> TJMartin wrote
>> Is there a way to combine the hillshade and topo map together - ie add a
>> transparency or opacity to the hillshade.
> A common way is to convert your (i assume) RGB topo map to HSV color space
> and replace the intensity band (v) with your hillshade. This is for example
> nicely shows at Franks blog:
> I find replacing the entire intensity band always a bit harsh since for flat
> area's you would like to stick to your original color/data.
> This can be done by normalizing the hillshade band with the zenith angle
> used to create is. Following the hsv_merge.py script on Franks blog, the
> original replacement of:
> /hsv_adjusted = numpy.asarray( [hsv, hsv, v] )/
> Woud become something like:
> /v_norm = v * hsv / ( 255 * math.cos(zenith * (math.pi / 180.0)))
> hsv_adjusted = numpy.asarray( [hsv, hsv, numpy.where( v_norm> 255,
> 255, v_norm )] )/
> You have to pay a bit attention to avoid saturating the normalized band.
> This can be done by choosing a zenith angle wich gives a flat surface a
> value halfway your range of 0 and 255.
GMT does that for ages. But it doesn't replace the intensity with
shading. Instead the shading is used to change the saturation of HSV
making the facets facing illumination light become brighter and the
other become darker.
See for example
With GMT5 one can compute the shading of a grid and apply it to an
image. Here goes a small example of the GMT's test set
Even simpler is to use Mirone. We can do the above just with a couple of
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 28153 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120221/9adb79b1/Capture-0001.jpe
More information about the gdal-dev