ST_intersects地理学

Question

我在想办法，我在postgres 9.6上使用postgis 2.3

地理上的ST_intersects返回 false，尽管两个地理相交。为了调试，我使用几何尝试了相同的查询，两者相交。

有人能帮我解释一下吗？

示例：

我正在运行此查询：

select ST_intersects(ST_GeogFromText('POLYGON((12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890))'), ST_GeogFromText('POLYGON((-170 -47,174 -47,174 77,-170 77,-170 -47))'))  from table limit 1;

我得到：

 st_intersects 
---------------
 f
(1 row)

当我跑的时候：

 select ST_intersects(ST_GeomFromText('POLYGON((12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890))'), ST_GeomFromText('POLYGON((-170 -47,174 -47,174 77,-170 77,-170 -47))')) from  table limit 1;

我得到：

 st_intersects
 ---------------
  t
 (1 row)

Answer 1

ST_Intersects -- 如果几何/地理在 2D 空间上相交，则返回 True。如果从多边形（在具有 5 个顶点的情况下）转换为 LINESTRING，则查询应返回相同的结果。测试：

SELECT 
 ST_Intersects(
 ST_GeomFromText('SRID=4326;LINESTRING(12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890)'), 
 ST_GeomFromText('SRID=4326;LINESTRING(-170 -47,174 -47,174 77,-170 77,-170 -47)'));


SELECT 
 ST_Intersects(
 ST_GeogFromText('SRID=4326; LINESTRING(12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890)'), 
 ST_GeogFromText('SRID=4326; LINESTRING(-170 -47,174 -47,174 77,-170 77,-170 -47)'));

更新:

实际上，对于< code>ST_Intersects，使用< code>Polygon就可以了。我意识到< code>Linestring不是您想要的——返回false，因为一个完全在另一个内部。问题是这两个多边形在球面坐标系中相距甚远。检查:

SELECT 
 ST_Distance(
 ST_GeomFromText('SRID=4326;POLYGON((12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890))'), 
 ST_GeomFromText('SRID=4326;POLYGON((-170 -47,174 -47,174 77,-170 77,-170 -47))'));

==

SELECT 
 ST_Distance(
 ST_GeogFromText('SRID=4326;POLYGON((12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890))'), 
 ST_GeogFromText('SRID=4326;POLYGON((-170 -47,174 -47,174 77,-170 77,-170 -47))'));

==

当然，您可以将< code>geography转换为< code>geometry类型，如下所示，以获得所需的输出:

SELECT 
 ST_Intersects(
 ST_GeogFromText('SRID=4326;POLYGON((12.856 41.890,12.856 41.935,12.811 41.935,12.811 41.890,12.856 41.890))')::geometry, 
 ST_GeogFromText('SRID=4326;POLYGON((-170 -47,174 -47,174 77,-170 77,-170 -47))')::geometry);