我的数据库结构是普通Postgres:
CREATE TABLE IF NOT EXISTS locations (
name text NOT NULL,
lat double precision NOT NULL,
lng double precision NOT NULL,
);
CREATE INDEX ON locations(lat,lng);
当我要计算边界框中左下角和右上角的所有位置时,请使用以下查询:
SELECT * FROM locations
WHERE lat >= min_lat AND
WHERE lat <= max_lat AND
WHERE lng >= min_lng AND
WHERE lng <= max_lng;
现在,我想生成一个给定点的边界框,并在locations查询中使用边界框结果。我正在使用以下PostGIS查询来生成边界框:
SELECT
ST_Extent(
ST_Envelope(
ST_Rotate(
ST_Buffer(
ST_GeomFromText('POINT (-87.6297982 41.8781136)',4326)::GEOGRAPHY,160934)::GEOMETRY,0)));
结果:BOX(-89.568160053866 40.4285062983089,-85.6903925527536 43.3273499289221)
但是,我不确定一次调用如何将PostGIS查询边界框中的结果用于原始lat / lng Postgres查询中。关于如何将两者合并的任何想法?最好这样以保留索引。
答案 0 :(得分:3)
如果要获取bbox坐标作为单独的值,则可以查看ST_XMax
,ST_YMax
,ST_XMin
,ST_YMin
。嵌入您的查询的以下CTE应该给您一个提示:
WITH j (geom) AS (
SELECT
ST_Extent(ST_Envelope(
ST_Rotate(ST_Buffer(
ST_GeomFromText('POINT(-87.6297982 41.8781136)',4326)::GEOGRAPHY,160934)::GEOMETRY,0)))
)
SELECT
ST_XMax(geom),ST_YMax(geom),
ST_XMin(geom),ST_YMin(geom)
FROM j
st_xmax | st_ymax | st_xmin | st_ymin
-------------------+-----------------+-------------------+------------------
-85.6903925527536 | 43.327349928921 | -89.5681600538661 | 40.4285062983098
附带说明:将几何值存储为数字可能看起来很简单,但这绝不是更好的选择-特别是在处理多边形时!因此,我真的建议您将这些值存储为geometry
或geography
,乍一看似乎很复杂,但从长远来看肯定会有所收获。
此答案可能有助于您了解涉及多边形的距离/包含性查询:Getting all Buildings in range of 5 miles from specified coordinates