Chapter 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

请注意，为了回答这个问题，我们必须计算每个多边形的面积。如果我们经常这样做，那么向表中添加一个可为性能编制索引的Area列将是有意义的。通过将结果按降序排列，并使用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
...

这个查询需要一段时间，因为表中的每条道路都被汇总到最终结果中(示例表中大约有25万条道路)。对于较小的数据集(数百条记录中的数千条记录)，响应速度可能非常快。

这是一个“叠加”的例子，它接受两个表，并输出一个由空间剪裁或剪切结果组成的新表。与上面演示的“空间连接”不同，该查询创建了新的几何图形。覆盖类似于增强型空间连接，对于更精确的分析工作非常有用：

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