PostGIS边界框查询返回奇怪的结果

时间:2014-11-04 19:39:40

标签: postgresql geospatial postgis

我尝试了以下SQL命令:

CREATE TABLE places(
    lat_lng geography(Point,4326),
    place_name varchar(50)
);

CREATE INDEX places_lat_lng_idx ON places USING gist(lat_lng);

INSERT INTO places values ('POINT(-126.4 45.32)', 'Food Bar1');
INSERT INTO places values ('POINT(-126.4 47.32)', 'Food Bar2');
INSERT INTO places values ('POINT(-125.4 47.42)', 'Food Bar3');

SELECT place_name, ST_AsText(lat_lng) as point FROM places WHERE places.lat_lng && ST_MakeEnvelope(-130.0, 44.0, -100.0, 46.7, 4326);

结果是:

 place_name |        point        
------------+---------------------
 Food Bar1  | POINT(-126.4 45.32)
 Food Bar2  | POINT(-126.4 47.32)
 Food Bar3  | POINT(-125.4 47.42)

这对我来说不合适,因为ymax是46.7,但“Food Bar2”和“Food Bar3”的ymax值分别为47.32和47.42。问题在哪里?

1 个答案:

答案 0 :(得分:11)

Here is your geography envelope

url_mapper

在这里,您的查询点。使用ST_Segmentize:

将信封展平为笛卡尔空间
SELECT ST_Segmentize(ST_MakeEnvelope(-130.0, 44.0, -100.0, 46.7, 4326)::geography,50000);

point_in_bbox

所以你是正确的,这些点应该在地理范围内,但你使用了&&边界框操作符,它忽略了几何形状。地理信封的边界框如下所示:

point_in_bbox_w_env

显示边界框中的所有点。

使用以下内容修复查询:

SELECT place_name, ST_AsText(lat_lng) as point
FROM places
WHERE ST_Intersects(ST_MakeEnvelope(-130.0, 44.0, -100.0, 46.7, 4326), places.lat_lng)