学了几天mapnik

因为工作需要，要把路况信息做成切片，因为之前那种方式（直接绘制路况）效果不大好，绘制速度太慢，也太过于消耗流量了（一次就把整个市的路况信息都返回来了）。四月十三号开始接到这个任务，并着手查找关于切片的资料。
　　最开始，主任的想法是通过算法实现，很显然要用到裁剪算法。但这个方式我并不赞同，且不说裁剪算法用于此次多么复杂，而且还要考虑到效率问题，polyline排序本来就非常的复杂而且来回遍历也是很费事的。主任同意我使用一些现成的切片库，这样可以节省很多精力。

• 首先想到的arcgis engine二次开发，但不久就放弃了。网上找到一些切片demo，好像是必须要安装arcgis server，那么就真的麻烦了。而且还要安装arcgis runtime，以及授权许可，简直麻烦到家了。
• 然后想到第三方开源软件geoserver，这个可以实现切片，通过二次开发的话应该也是没问题的。但这方面的资料好少啊，简直无从下手。
• 使用gdal工具gdal2tiles.py，这个看起来不错，可以实现切片功能。需要配置gdal环境，不知怎么得就放弃了这个。
• 是看到了mapnik，才放弃了gdal2tiles.py，也许是后者资料太少吧。mapnik也是需要安装的，可以选择三种语言来开发——python、nodejs以及c++。考虑的效率问题，我最开始打算使用c++，要求版本好像挺高的，需要c++11；之后就想用nodejs，很早就像用nodejs实战一下，何况对于ajax请求以及json解析也是存在优势的，最后也不得不放弃了，因为官方放出的资料太少了，连一个demo都不给我。不得不选用python了，也不错，头一次使用python写程序，还有些小激动呢！
  　　虽然是第一次使用python，但开发起来效率还是比较高的，分为以下步骤：

1. 读取本地测试数据json，并解析
   　　中间出现过字符编码问题，纠结过一段时间。
2. 生成shapefile文件
   　　同样利用第三方库pyshp.
3. 利用网上开源代码generate_tiles.py，实现多级别切片
   　　不幸呀！出现了重叠部分，纠结了整整一天。有可能是地图投影是出现了问题，当我设置m.aspect_fix_mode属性是没有出现效果
4. 利用NSIS封装打包
   　　由于第三步出现了故障，打包推迟。
   　　老杨给出了比较好的解决方案，直接通过行列号反算经纬度，实现bbox范围的确定，解决了切面重叠问题。老杨很聪明，这一方面也不得不服。那么问题来了，路况信息5分钟刷新一次，切片也要在五分钟内完成吧。然而事实上很残酷，毕竟是到20级呀，一套下来4700秒，接近80分钟。显然不合理，20级地图占用了太多的时间，理论上是更高一级是低一级耗时的四倍，也就是说仅20级就用了一个小时。假使服务器cpu是32个核，都占用起来应该会缩短到五分钟以内，但想想机房里cpu疯狂般的不停高负荷运转也是够了。也许使用geojson，借助前端地图渲染引擎来支持或许是个好办法。