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

本文主要是介绍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}});

显示效果如下
前台展示