ST_HexagonGrid — 返回一组六边形和单元格索引,它们完全覆盖几何参数的边界。
setof record ST_HexagonGrid(
float8 size, geometry bounds)
;
从平面六边形平铺的概念开始。(不是地球的六边形平铺,这不是H3平铺方案。)对于给定的平面 SRS 和给定的边长,从 SRS 的原点开始,平面有一个唯一的六边形平铺,平铺(SRS,大小)。此函数回答了以下问题:给定平铺(SRS,大小)中的哪些六边形与给定边界重叠。
输出六边形的 SRS 由边界几何提供。
将六边形的边长加倍或三倍会生成一个新的父平铺,它与原始平铺相吻合。不幸的是,无法生成子瓦片完美拟合的父六边形平铺。
可用性:3.1.0
要对六边形平铺进行点汇总,请使用点的范围作为边界生成六边形网格,然后进行空间连接到该网格。
SELECT COUNT(*), hexes.geom FROM ST_HexagonGrid( 10000, ST_SetSRID(ST_EstimatedExtent('pointtable', 'geom'), 3857) ) AS hexes INNER JOIN pointtable AS pts ON ST_Intersects(pts.geom, hexes.geom) GROUP BY hexes.geom;
如果我们为每个多边形边界生成一组六边形,并过滤掉那些与它们的六边形不相交的六边形,我们最终会得到每个多边形的平铺。
平铺状态导致每个州的六边形覆盖,以及多个六边形在州之间的边界处重叠。
当引用 FROM 列表中的先前表时,集返回函数隐含 LATERAL 关键字。因此,CROSS JOIN LATERAL、CROSS JOIN 或仅使用逗号是此示例的等效构造。 |
SELECT admin1.gid, hex.geom FROM admin1 CROSS JOIN ST_HexagonGrid(100000, admin1.geom) AS hex WHERE adm0_a3 = 'USA' AND ST_Intersects(admin1.geom, hex.geom)