PostGIS：使用 Lats/Lngs 进行空间连接求和

Question

我有一个名为 pets 的表格，其中包含 lats 和 lng 坐标和变量，还有一个通过上传名为 city 的形状文件创建的表格，其中包含 left, right, top, bottom, lat, lng, 和 geom. 列

    lat    lng   cats  dogs
   -99.4   42.1   2     0
   -98.1   44.3   1     1
   -99.7   43.6   0     3 
   -99.8   42.0   2     1

我想知道如何构建一个查询来汇总每个边界内的猫和狗的数量。通常我会在 QGIS 中进行空间连接，但我的数据集太大，在成功连接两者之前会崩溃。

我对 QGIS 和 Postgres 非常熟悉，但对 PostGIS 很陌生。我已经构建了一个查询，但未能获得所需的答案：

select sum(pets.dogs)
from pets, city
WHERE ST_Within(pets.dogs, city.geom);

非常感谢为新手编写的 PostGIS 上的任何指针或资源。谢谢！

Answer 1

pets 表需要一个代表 lat/lng 的 GeometryColumn。像这样添加它：

psql> SELECT AddGeometryColumn ('my_schema','pets','geom',4326,'POINT',2);
psql> UPDATE pets SET geom=ST_SetSRID(ST_MakePoint(lng, lat), 4326);

对于大型数据集，建议在列中添加spatial index：

psql> CREATE INDEX geom_idx ON pets USING GIST (geom);

这应该允许快速查询，例如：

psql> SELECT sum(pets.dogs)
      FROM pets, city 
      WHERE ST_Within(pets.geom, city.geom);

这将对任何城市几何体中的所有狗进行一次计数。按城市名称（或 ID）分组以按城市获取计数：

psql> SELECT city.id, city.name, sum(pets.dogs)
      FROM pets, city 
      WHERE ST_Within(pets.geom, city.geom)
      GROUP BY city.id;