PostGIS计算矢量切片(二)--按值渲染

方案背景

今年三月份写了一篇postgis计算矢量切片，参考了网上资料给出了一份很粗糙的相关方案(文章写的更粗糙)。当时的方案中只能针对gis形状进行渲染，而不能用属性渲染.针对这个情况，本文进行相对应的修改。

前期准备

软件是用的是Qgis和Postgis(Postgis版本为2.42 ,Postgresql版本为10.0)。qgis使用Ramdom points in extent 用来生成测试数据，数据范围参数使用
120.6327590942379970,120.8625335693359943,31.2309341430663991,31.4467678070068004（xmin,xmax,ymin,ymax）,数据量设置为十万。生成好数据，再建好索引,添加字段v,根据奇偶性给字段赋值。

后台梳理

先看了ST_AsMVT和ST_AsMVTGeom帮助文档，一头雾水.使用谷歌大法，发现了这个资料参考资料

只要通过上图的方式就能查询到我想要的带属性值的矢量切片，代码如下。只要先输入需要查找的范围，用ST_AsMVTGeom查询出范围内矢量并转换成屏幕坐标，最后用ST_AsMVT压缩数据成mvt格式，和之前文章有所不同的是，查询过程中加入了属性的查询，压缩之后矢量切片就会带上我们想要的属性值（w.v就是我们想要的属性值）

SELECT ST_AsMVT(tile,'points',4096,'geom') tile  FROM (SELECT w.v,ST_AsMVTGeom(w.the_geom,Box2D(ST_MakeEnvelope(119.531250,30.751278,120.937500,31.952162, 4326)),4096, 0, true)	 AS geom FROM pnt w) AS  tile where  tile.geom is not null

做完这些我们只要用后台语言写一个服务,将前台请求的切片坐标z/x/y这些转换成经纬度得到最大最小经纬度，放到前面提到的sql中去查询就可以。

前台展示

前台展示还是选用mapbox,添加自定义矢量切片数据源，修改下按值渲染颜色

map.addLayer({
               "id": "custom-go-vector-tile-layer",
                "type": "circle",
                "source": "custom-go-vector-tile-source",
                "source-layer": "points",
                paint: {
                    'circle-radius': {
                        stops: [
                            [8, 0.1], [11, 0.5],[15, 3],[20, 60]
                        ]
                    },
                    'circle-color': {
                        property: 'v',
                        stops: [
                            [0, '#990055'],
                            [1, '#2a55b9']
                        ]
                    },
                    'circle-opacity': 1
                }
            });

显示效果如下
前台展示

性能对比

一百万数据，页面加载结束5s;十万数据，页面加载结束1.86s.同样数据，百万数据，geoserver的(openlayer)页面加载结束20.20s;十万数据时，geoserver的的(openlayer)页面加载结束4.32s

附代码地址：

https://github.com/tpolong/postgisvectortile