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

Jus*_*der 3 postgresql postgis geospatial

我尝试了以下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);
Run Code Online (Sandbox Code Playgroud)

结果是:

 place_name |        point        
------------+---------------------
 Food Bar1  | POINT(-126.4 45.32)
 Food Bar2  | POINT(-126.4 47.32)
 Food Bar3  | POINT(-125.4 47.42)
Run Code Online (Sandbox Code Playgroud)

这对我来说不合适,因为ymax是46.7,但"Food Bar2"和"Food Bar3"的ymax值分别为47.32和47.42.问题出在哪儿?

Mik*_*e T 11

这是您的地理信封:

url_mapper

在这里,您的查询点.使用ST_Segmentize将信封展平为笛卡尔空间:

SELECT ST_Segmentize(
   ST_MakeEnvelope(-130.0, 44.0,
                   -100.0, 46.7, 4326)::geography,50000);
Run Code Online (Sandbox Code Playgroud)

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)
Run Code Online (Sandbox Code Playgroud)