地理网格模型

简介

地理网格模型是一种统一、简单的地理空间划分和定位参考系统，依据统一规则，将地面区域按照一定经纬度或地面距离进行连续分割，并将空间不确定性控制在一定范围内，形成规则多边形，每个多边形均称为格网单元，从而构成分级、分层次的多级格网体系，实现地面空间离散化，并赋予统一编码。

Ganos GeomGrid是对象关系型数据库PostgreSQL兼容版本（PolarDB PostgreSQL版）的一个时空引擎扩展，用于表达和处理地理网格对象。此模型主要提供了地理网格以及地理网格数组类型的空间操作函数，空间关系判断函数，以及与Ganos中其他数据类型的转换计算函数

GeomGrid在Ganos中指代地理网格数据类型。地理网格是一种再现地球表面的多边形网格单元集合，可以用于表示地物在地理空间中的位置信息，融合其他各类时空数据。每个网格单元都会进行编码，网格与编码是一一对应的。三维地理网格不只考虑经纬度，还把高度维纳入剖分和编码范围。

GeomGrid模块目前涵盖GeoSOT和H3两种地理网格。GeoSOT是中国提出的一套地球空间剖分理论，并在此基础上发展出的一种离散化、多尺度区域位置标识体系。H3是Uber研发的一种覆盖全球表面的二维地理网格，采用的基本网格是正六边形。

功能概述

GeomGrid支持地理网格的输入输出，空间操作和空间关系判断：

• 支持Text、Bytea、Geometry、Meshgeom、Sfmesh、Vomesh等类型与GeomGrid类型的互相转换。

• 支持网格层级获取，网格坐标获取，父子节点计算，网格合并和网格剖分等空间查询操作。

• 支持geomgrid与geometry、meshgeom、sfmesh、vomesh以及gemogrid类型之间的空间关系判断操作。

更多内容请参见GeomGrid SQL参考。

主要业务场景

Ganos GeomGrid可支持多种业务场景，最常见的包括以下几个：

• 网格查询

  给定一个网格或网格数组，查询是否与数据库中存储的网格相交。支持退化网格，即可以通过子网格查出父网格，反之亦然。例如在停放共享单车时，停放区的网格码存在数据库中，而用户的位置可以转换为网格码到库中搜索，找出最近的停放区。网格查询也常用于关联不同的图层，把不同图层的地物转化为网格存储，然后根据某一网格查出所有相关图层的地物。

• 网格聚合

  物流、快送等行业中经常以网格为单位统计时空要素。例如船或车的轨迹数据可经过网格聚合生成热力图，先把轨迹转化为网格表达，然后统计每个网格中的轨迹数量作为热力进行绘制。再例如以街道为单位收集的人口数据可分配至网格进行可视化表达。

  

• 网格寻路

  无人机在复杂环境下可以借助三维网格进行航线规划，通过网格设置空间障碍物的穿越开销，借助路径规划算法可得到网格路径。

  

概述

基于GeoSOT和H3地理网格的理论框架，在Ganos GeomGrid模块中，两种网格分别被实现成GeomGrid数据类型和H3Grid数据类型。

编码

GeoSOT

GeoSOT将地球通过简单投影变换到平面，将180°×360° 地球表面空间扩展为512°×512°，并将该空间作为第0级剖分面片，且面片中心与赤道和本初子午线的交点重合。在此基础上，对第0级剖分面片进行递归四叉剖分，直到32级。这样，整个地球表面经纬度空间在经线方向和纬线方向可严格的整型二分，由此将整个地球分割为大到全球、小到厘米级的整度、整分和整秒的层次网格。

GeoSOT二维网格最长的编码位为32位四进制数值编码。第1-9位是度级网格编码，第10-15位是分级网格编码，第16-21位是秒级网格编码，第22-32位是秒以下网格编码，编码长度即为网格层级。

经度和纬度均采用A°B′C.D″的表达形式，将度由十进制数转换成8bit的二进制数A2；将分由十进制数转换成6bit的二进制数B2；将秒由十进制数转换成6bit的二进制数C2；将秒的小数部分由十进制数转换成11bit的二进制数D2；将二进制数A2、B2、C2 以及D2顺次拼接为31bit的二进制数。将纬度前置，经度后置，采用莫顿交叉的方式生成62位的混合代码。在混合代码前加上G0、G1、G2或G3即可得到网格代码。G0对应东北半球，G1对应西北半球，G2对应东南半球，G3对应西南半球。

GeoSOT三维网格在二维网格的基础上增加了高度域。为保证高度域的划分与经纬度一致，对应地表1º 网格，地面向上延伸至高空256个网格，向下延伸至地心256个网格。而后对应不同层级进行二分处理得到垂直方向的层数，并以此进行二进制编码，与经纬度编码莫顿交叉形成三维网格编码。

H3

H3采用的投影方式是以二十面体近似地球球体，如下图所示，这个二十面体的每一个面都是球面三角形，有12个顶点，称为球形二十面体。

如下图所示，H3在每个三角面上都划分相同排列方式的六边形，进而将全球划分成122个基本单元，然后对每个基本单元再进行六边形剖分。

H3从122个基本单元开始，对每个单元递归地细分为7个子单元，如下图所示，以编号为20的基本单元为例进行分解。

H3的编码值最多占63个比特位，可以用一个长整型表示，通常序列化为十六进制字符串。H3编码的位布局提供了一个相对紧凑的结构，用于存储有关给定索引的信息并定义其地理位置，如下图所示。

索引

索引通过将数据组织到搜索树中来加速搜索，避免大数据量时全局“顺序扫描”带来的延迟。GeomGrid模块适用的索引主要包括以下几种：

简介

快速入门文档帮助用户快速理解 Ganos GeomGrid 引擎的基本用法，包括扩展创建、建表、插入数据、计算网格码、创建索引、查询等内容。

语法说明

• 创建扩展。

  CREATE EXTENSION Ganos_GeomGrid CASCADE;

  说明

  建议将扩展安装在public模式下，避免权限问题。

  CREATE extension ganos_geometry WITH schema public cascade;

• GeoSOT网格操作。

  • 创建具有网格码的表。

     -- 可以为几何
     CREATE TABLE t_grid(id integer,
                         geom geometry,  -- 几何对象
                         grid1 geomgrid[], -- 网格码，精度1
                         grid2 geomgrid[], -- 网格码，精度2
                         grid3 geomgrid[] -- 网格码，精度3
                        );

  • 插入数据。

     --插入点数据
     INSERT INTO t_grid(id, geom)
     VALUES (1, ST_GeomFromText('POINT(116.31522216796875 39.910277777777778)', 4490)),
            (2, ST_GeomFromText('POINT(116.31522217796875 39.910277776777778)', 4490)),
            (3, ST_GeomFromText('POINT(116.31522217797875 39.910277776787778)', 4490)),
            (4, ST_GeomFromText('POINT(116.31522227796875 39.910277776775778)', 4490));
    
     -- 插入面数据
     INSERT INTO t_grid(id, geom) VALUES(5, 'SRID=4490;POLYGON((-0.08077 -0.02814, 0.0482 -0.03, 0.07426 0.03724, -0.08077 -0.02814))'::geometry);
    
     -- 插入三维数据
     INSERT INTO t_grid(id, geom) VALUES(6, 'SRID=4490;CIRCULARSTRING Z (-63.597471 44.8071 20,-63.597 44.807 0,-63.5974 44.807 40)'::geometry);

  • 计算网格码。

     -- 创建不同精度的网格码
     UPDATE t_grid
     SET grid1 = ST_AsGrid(geom, 10),
         grid2 = ST_AsGrid(geom, 15),
         grid3 = ST_AsGrid(geom, 26);
    
     -- 生成退化网格码
     UPDATE t_grid SET grid1 = st_asgrid(geom, 18, true) WHERE id = 5;
    
     -- 生成三维网格码
     UPDATE t_grid SET grid1 = st_as3dgrid(geom, 25) WHERE id=6;

  • 对网格码创建索引。

     -- 针对不同精度的网格码创建GIN索引
     CREATE INDEX idx_grid_gin1
     ON t_grid
     USING GIN(grid1);
    
     CREATE INDEX idx_grid_gin2
     ON t_grid
     USING GIN(grid2);
    
     CREATE INDEX idx_grid_gin3
     ON t_grid
     USING GIN(grid3);

  • 查询。

     -- 完全在某个网格中
     SELECT id
     FROM t_grid
     WHERE grid2 = ARRAY[ST_GridFromText('G001310322230230')];
    
     -- 和某个网格相交
     SELECT id
     FROM t_grid
     WHERE grid3 @> ARRAY[ST_GridFromText('G00131032223023031031033223')];
    
     -- 和某些网格相交
     SELECT id
     FROM t_grid
     WHERE grid3 && ARRAY[ST_GridFromText('G00131032223023031031211001'),
                              ST_GridFromText('G00131032223023031031211111')];
    
     -- 和某个几何对象在网格上相交
     SELECT id
     FROM t_grid
     WHERE grid3 &&
     ST_AsGrid(
       ST_GeomFromText('LINESTRING(116.31522216796875 39.910277777777778, 116.31522217797875 39.910277776787778)', 4490), 26);

• H3网格操作。

  • 建表。

     CREATE TABLE h3_grid(
      id integer,
      geom geometry,  -- 几何对象
      h3 h3grid[]     -- H3对象类型
     );

  • 插入数据。

     INSERT INTO h3_grid VALUES (1, 'POINT(102.5 25.7)'::geometry);
     INSERT INTO h3_grid VALUES (2, 'POLYGON((-0.08077 -0.02814, 0.0482 -0.03, 0.07426 0.03724, -0.08077 -0.02814)'::geometry);

  • 计算H3网格码。

     -- 普通网格码
     UPDATE h3_grid SET h3 = ST_AsH3Grid(geom, 7);
    
     -- 退化网格码
     UPDATE h3_grid SET h3 = ST_AsH3Grid(geom, 7, true);

  • 建立索引。

    CREATE INDEX h3_grid_btree ON h3_grid(h3);

  • 查询。

     -- 显示H3编码
     SELECT st_astext(h3[1]) FROM h3_grid;
    
     -- 搜索
     SELECT * FROM h3_grid WHERE h3 > ARRAY[ST_H3FromText('884a126689fffff')];

• 删除扩展（可选）。

  DROP EXTENSION Ganos_GeomGrid CASCADE;

SQL参考

详细SQL手册请参见GeomGrid SQL参考。