• 周一. 10 月 7th, 2024

5G编程聚合网

5G时代下一个聚合的编程学习网

热门标签

百万级要素优化:基于PG与PostGIS的实时矢量瓦片服务

admin

11 月 28, 2021

背景介绍

对于某些应用场景,数据量很容易达到百万级要素,数据每天都在改变。而且还要考虑在1-6级这种小级别时候的数据展示,那么这个时候仅仅是PG与PostGIS简单的矢量瓦片服务比较困难达到客户想要的效果,这个时候就要考虑做一些优化了。

优化方法

百万级要素优化对于不同的业务来说,有相应不同的解决办法,没有最好,只有更好,这个时候要做考量,下面介绍几个常用的会考虑使用的方法:

1、依据业务,对要素进行抽稀过滤

根据业务层面,在不同场景浏览不同的要素,虽然总的要素在几百万,但是这个时候过滤一部分后剩余要素可以进行实时矢量瓦片展示。

2、几种抽稀过滤的算法

除第一种方法外,其余几种方法包含一些个人的想法与思路,其它小伙伴有其他好的思路,或者对我的想法思路有疑问的话,欢迎在评论区留言。

1)道格拉斯-普客算法(DP)

该算法个人理解为是简化的算法,对于线或者面的边界,可以使用该算法对节点进行抽稀,达到对线或者面的简化。但是为什么在点要素抽稀过滤中也介绍该算法呢?

答案:在轨迹数据中,每一次存储的都是点位信息(点),所以对于该种数据可以使用该算法进行过滤,但是对于一个点即代表一种信息的数据来说,该算法不太适用。

2)基于网格的抽稀过滤算法

该算法结果为均匀分布。

  1. 根据数据范围以及阈值将范围划分网格。
  2. 遍历网格,利用网格与数据求交(使用PostgreSQL的gist索引,加快速度),计算交到的要素。
  3. 将交到要素随机选取一个,其余要素舍弃。
  4. 遍历完网格后,即可完成要素过滤。

本人写了一个简单的实现,供各位参考。

--drop FUNCTION point_simplify_grid;
CREATE
OR REPLACE FUNCTION point_simplify_grid (
    PAR_minx DOUBLE PRECISION,
    PAR_miny DOUBLE PRECISION,
    PAR_maxx DOUBLE PRECISION,
    PAR_maxy DOUBLE PRECISION,
    PAR_threshold DOUBLE PRECISION
) RETURNS geometry AS $$
DECLARE REC_points RECORD ; ARR_result_points geometry []; GEO_point geometry ; grid_key VARCHAR (10) ; GEOM_result geometry ; grid_num_x INTEGER ; grid_num_y INTEGER ; grid_json jsonb ;
BEGIN
    grid_num_x := CEIL (
        (PAR_maxx - PAR_minx) / PAR_threshold
    ) ; grid_num_y := CEIL (
        (PAR_maxy - PAR_miny) / PAR_threshold
    ) ; grid_json := '{}' :: jsonb ; ARR_result_points := '{}' :: geometry []; FOR REC_points IN SELECT
        *
    FROM
        tycd_view
    WHERE
        geom && st_makeenvelope (
            PAR_minx,
            PAR_miny,
            PAR_maxx,
            PAR_maxy,
            4326
        ) LOOP grid_key := CEIL (
            (
                st_x (REC_points.geom) - PAR_minx
            ) / PAR_threshold
        ) || ',' || CEIL (
            (
                st_y (REC_points.geom) - PAR_miny
            ) / PAR_threshold
        ) ;
    IF NOT (grid_json ? grid_key) THEN
        ARR_result_points := array_append(
            ARR_result_points,
            st_setsrid (REC_points.geom, 4326)
        ) ; grid_json := grid_json || ('{"' || grid_key || '":true}') :: jsonb ;
    END
    IF ;
    END loop ; GEOM_result := st_union (ARR_result_points) ; RETURN GEOM_result ;
    END ; $$ LANGUAGE plpgsql;

--select st_asgeojson(point_simplify_grid(73,3,135,53,0.2))

github链接:https://github.com/MrSmallLiu/point_dilution

3)基于距离的抽稀过滤算法

该算法结果为均匀分布。

  1. 设置距离阈值。
  2. 选取初始点,计算其它点到该点的距离,阈值范围内的点,舍弃。
  3. 从上次结果中,选取另一个点作为基准点,再次执行(2)的步骤,只需计算保留下来的点。
  4. 全部点计算完成后,即可完成要素抽稀过滤。

4)基于随机数的抽稀过滤算法

该算法结果可以保持数据疏密。

对于大数据量时的随机数,可以达到伪随机性,所以可以利用该算法进行随机性的抽稀。

  1. 设置经验阈值,例如r=0.5,同时因为该算法效率可以实时进行,可以考虑基于每一次查询矢量切片时,根据级别设置阈值。
  2. 查询的sql语句中或者其它方式后面添加where条件,where random()>0.5时保留要素。
  3. 基于以上就可以完成50%的抽稀率。

使用的局限:

  1. 对于小数据量可能效果达不到。

优势:

  1. 效率很高,可以不用提前处理数据。
  2. 可以将random()做到数据字段中,使用btree加速。
  3. 设置字段中后,可以保证例如小级后的结果例如3级,在后续级别中保留例如4-18级。

5)基于网格与随机数抽稀过滤算法

该算法结果可以保持数据疏密。

  1. 根据范围与网格大小阈值,划分数据网格。
  2. 遍历网格,利用网格与数据求交(使用PostgreSQL的gist索引,加快速度),计算交到的要素。
  3. 在网格内根据随机数,将网格内数据随机丢弃一些。
  4. 网格遍历完成后,即完成数据抽稀过滤。

该算法弥补了单纯基于数据数的一些局限性,不会导致数据存在缺块(对于大数据量,可能性微乎其微)。

6)基于距离并保持疏密程度的抽稀过滤算法

该算法结果可以保持数据疏密。

  1. 设置距离阈值,设置丢弃比例。
  2. 选取初始点,遍历其它要素,计算其它要素到该要素的距离,小于距离阈值的归为一片。
  3. 将该片数据,基于丢弃比例,舍弃一部分要素。
  4. 继续下一个点为基准点,循环(2)(3)的步骤。
  5. 全部要素完成后,即可完成数据抽稀过滤。

以上即为介绍的几种抽稀过滤算法,可以在合适的场景使用合适的方式。其中虽然有些算法存在局限性(基于随机数的算法),但是对于大数据量来说,可能性很低,并且效率极高,可以做到实时,以下效果也是采用该方法。

效果展示

示例:

  • 数据量: 本示例数据量为300万点要素;
  • 抽稀过滤方法:基于随机数抽稀过滤,将random字段添加到字段中,同时利用该方法结合经验值,实时抽稀过滤展示(例如:<4级:0.2,4-6级:0.4,7-9:0.7,>10不抽稀过滤);
  • 数据展示方法: 实时矢量切片,结合根据级别抽

发表回复