第 5 章 空间查询

您可以使用非常简单的 SQL 语句来回答这个问题

SELECT sum(ST_Length(geom))/1000 AS km_roads FROM bc_roads;

km_roads
------------------
70842.1243039643

此查询将属性条件（在市政区域名称上）与空间计算（多边形面积）结合起来。

SELECT
  ST_Area(geom)/10000 AS hectares
FROM bc_municipality
WHERE name = 'PRINCE GEORGE';

hectares
------------------
32657.9103824927

此查询使用空间度量作为排序值。解决此问题的方法有很多，但最有效的方法如下

SELECT
  name,
  ST_Area(geom)/10000 AS hectares
FROM bc_municipality
ORDER BY hectares DESC
LIMIT 1;

name           | hectares
---------------+-----------------
TUMBLER RIDGE  | 155020.02556131

请注意，为了回答此查询，我们必须计算每个多边形的面积。如果我们经常这样做，那么将面积列添加到表中以用于索引性能将是有意义的。通过以降序排列结果，然后使用 PostgreSQL 的“LIMIT”命令，我们可以轻松地选择最大值，而无需使用像 MAX() 这样的聚合函数。

这是一个“空间连接”的例子，它使用空间交互（“包含”）作为连接条件（而不是通常基于共同键的关联方法）将来自两个表的数据连接在一起。

SELECT
  m.name,
  sum(ST_Length(r.geom))/1000 as roads_km
FROM bc_roads AS r
JOIN bc_municipality AS m
  ON ST_Contains(m.geom, r.geom)
GROUP BY m.name
ORDER BY roads_km;

name                        | roads_km
----------------------------+------------------
SURREY                      | 1539.47553551242
VANCOUVER                   | 1450.33093486576
LANGLEY DISTRICT            | 833.793392535662
BURNABY                     | 773.769091404338
PRINCE GEORGE               | 694.37554369147
...

此查询需要一段时间，因为表中的每条道路都会汇总到最终结果中（示例表中约有 250,000 条道路）。对于较小的数据集（几千条记录，几百条记录），响应速度非常快。

这是一个“叠加”的例子，它接收两个表并输出一个新的表，该表包含空间裁剪或切割的结果。与上面演示的“空间连接”不同，此查询创建新的几何图形。叠加就像一个增强的空间连接，对于更精确的分析工作很有用。

CREATE TABLE pg_roads as
SELECT
  ST_Intersection(r.geom, m.geom) AS intersection_geom,
  ST_Length(r.geom) AS rd_orig_length,
  r.*
FROM bc_roads AS r
JOIN bc_municipality AS m
  ON ST_Intersects(r.geom, m.geom)
WHERE
  m.name = 'PRINCE GEORGE';

SELECT
  sum(ST_Length(r.geom))/1000 AS kilometers
FROM bc_roads r
JOIN bc_municipality m
  ON ST_Intersects(m.geom, r.geom
WHERE
  r.name = 'Douglas St'
  AND m.name = 'VICTORIA';

kilometers
------------------
4.89151904172838

SELECT gid, name, ST_Area(geom) AS area
FROM bc_municipality
WHERE ST_NRings(geom) > 1
ORDER BY area DESC LIMIT 1;

gid  | name         | area
-----+--------------+------------------
12   | SPALLUMCHEEN | 257374619.430216