原文地址:http://blog.thematicmapping.org/2012/07/terrain-mapping-with-mapnik.html
作者:作者:Bj?rnSandvik(Google Plus:https://plus.google.com/118196887774002693676)
在之前的三篇博文中,我们使用DEM数据分别创建了三个GroTiff的数据,分别是:山体阴影,坡度阴影和彩色地形图,如图1所示。在这篇文章中,我们将把这三个数据合并为一个数据,即Jotunheimen的地形图。
图1 之前三篇博文中生成的三个数据
上面图 1中的三个数据下载地址如下:
- 山体阴影:http://thematicmapping.org/playground/terrain/jotunheimen_hillshade.png
- 坡度阴影:http://thematicmapping.org/playground/terrain/jotunheimen_slopeshade.png
- 彩色地形图:http://thematicmapping.org/playground/terrain/jotunheimen_color_relief.png
Mapnik是一个用于地图渲染和绘制强大的开源工具包。接下来我们将使用RasterSymbolizer来对GeoTiff数据进行合并。样式和图层(影像数据)可以使用xml(jotunheimen_terrain.xml)文件来进行配置,文件内容如下:
在这个xml文件的第一行,定义了一个地图投影(UTM 32N)。Mapnil使用PROJ.4库,你可以在网站http://spatialreference.org/ref/epsg/32632/proj4/找到这个地图投影对应的PROJ4的格式。Mapnik不能修改图像的投影,所以在使用Mapnik之间把数据都转为所需的投影。
接下来,我定义了此地图的两个风格的RasterSymbolizer。山体阴影来和彩色地形图文件使用乘法混合模式进行混合。山体阴影的每个像素乘以对应的彩色地形图的像素值。同时还设置了山体阴影的不透明度为0.6,减少其的贡献量。接下来是地图图层的样式的定义。Mapnik将按照这个顺序进行渲染,首先从最顶层的图层开始。
你可以使用一个简单的python脚本(jotunheimen_terrain.py)来渲染地图图像(你也可以使用像TileCache和MapProxy工具)。
#!/usr/bin/env python
import mapnik
map = mapnik.Map(3134, 3134)
mapnik.load_map(map, 'jotunheimen_relief.xml')
map.zoom_all()
mapnik.render_to_file(map,'jotunheimen_relief.png')
上面的脚本创建一个PNG图像,使用方法是:
python jotunheimen_terrain.py
得到的结果如下图所示:
接下来,我们在上面的xml中添加坡度阴影数据(jotunheimen_terrain2.xml),改后的xml内容如下:
上面我调整了坡度阴影和山体阴影的透明度来生成新的图像。处理的结果如下图所示:
上面两个图像那个好是完全是由个人的喜好决定的。第一个图中,坡面比较亮,第二个图中,陡峭的山比较明显,而且坡面的细节信息更丰富。
在这两个图中,我很想念Jotunheimen地区的两个重要的地类——湖泊和冰川。在下一个博文中将讨论怎么添加土地覆盖数据。
………………………………………………………分割线……………………………………………………………
关于Bjørn Sandvik系列的博文,暂时翻译这四篇,剩下的不是我关心的,我就不翻译了,有兴趣的童鞋可以去作者的博客看看。有空写个程序,把上面的流程整理一下。