因为工作需要,要把路况信息做成切片,因为之前那种方式(直接绘制路况)效果不大好,绘制速度太慢,也太过于消耗流量了(一次就把整个市的路况信息都返回来了)。四月十三号开始接到这个任务,并着手查找关于切片的资料。
最开始,主任的想法是通过算法实现,很显然要用到裁剪算法。但这个方式我并不赞同,且不说裁剪算法用于此次多么复杂,而且还要考虑到效率问题,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,但开发起来效率还是比较高的,分为以下步骤:
- 读取本地测试数据json,并解析
中间出现过字符编码问题,纠结过一段时间。 - 生成shapefile文件
同样利用第三方库pyshp. - 利用网上开源代码generate_tiles.py,实现多级别切片
不幸呀!出现了重叠部分,纠结了整整一天。有可能是地图投影是出现了问题,当我设置m.aspect_fix_mode
属性是没有出现效果 - 利用NSIS封装打包
由于第三步出现了故障,打包推迟。
老杨给出了比较好的解决方案,直接通过行列号反算经纬度,实现bbox范围的确定,解决了切面重叠问题。老杨很聪明,这一方面也不得不服。那么问题来了,路况信息5分钟刷新一次,切片也要在五分钟内完成吧。然而事实上很残酷,毕竟是到20级呀,一套下来4700秒,接近80分钟。显然不合理,20级地图占用了太多的时间,理论上是更高一级是低一级耗时的四倍,也就是说仅20级就用了一个小时。假使服务器cpu是32个核,都占用起来应该会缩短到五分钟以内,但想想机房里cpu疯狂般的不停高负荷运转也是够了。也许使用geojson,借助前端地图渲染引擎来支持或许是个好办法。