计算每个多边形中的点数

Question

我有两张桌子，

points：有两列 latitude(float) 和 longitude(float)。

polygon：它有两列 polygon_id(int) 和 geom(多边形的形状)

我想打印出每个 polygon_id 中包含的点数。

我可以弄清楚如何获取包含一个点但无法进一步发展的多边形 ID。

SELECT polygon.id
FROM polygon
WHERE st_contains(geom, ST_GeomFromText('POINT(120.32155 -33.843322)', 4326)) = true

我想使用点表中的纬度和经度，但是，

ST_GeomFromText('POINT(points.latitude points.longitude)

没用。

预期结果是具有 polygon_id 以及该 polygon_id 包含的每一行的点。因为有 100 个多边形，所以我希望有 100 行。

Answer 1

你缺少的是ST_MakePoint。

考虑您的表格结构和以下示例，它们由一个多边形和两个重叠点组成..

CREATE TABLE polygon (id int, geom geometry);
INSERT INTO polygon VALUES (1,'POLYGON((120.20 -33.75,120.43 -33.75,120.43 -33.88,120.20 -33.88,120.20 -33.75))'::geometry);

CREATE TABLE points (lon numeric, lat numeric);
INSERT INTO points VALUES (120.32,-33.84);
INSERT INTO points VALUES (120.39, -33.85);

.. 执行以下查询：

SELECT polygon.id,COUNT(points.lon)
FROM polygon,points
WHERE ST_Contains(geom, ST_MakePoint(lon,lat))
GROUP BY polygon.id

 id | count 
----+-------
  1 |     2
(1 row)

我建议您不要将纬度和经度存储在不同的列中。将点存储为geometry 或geography 始终是PostgreSQL/PostGIS 中的最佳选择。