在Postgis 2.5中查询点是否在面内

Question

我在地图中有一个定义为GeoJSON的多边形。需要知道某个点是否在该多边形内。

多边形:

ST_AsGeoJSON( 
 'POLYGON ((
    -69.62607888946572 -22.389720278609026, 
    -69.62632833490409 -22.39009971947581,
    -69.6252903200153 -22.3904196394022,
    -69.62501673469583 -22.390032759398007,
    -69.62553708324471 -22.389888919122058,
    -69.62607888946572 -22.389720278609026))')

点：

ST_Point(-69.625137, -22.389777)

需要知道问Postgis这个查询的正确方法。

我们测试了不同的替代方案，其中一些返回正确答案，但不知道它们有何不同以及哪种是“正确的方法”（解释ST_Intersects、ST_Contains和ST_DWithin之间的差异将不胜感激）。我是否需要在查询中“设置SRID”点和多边形？

1：

select
ST_SetSRID(ST_GeomFromGeoJSON(ST_AsGeoJSON( 
 'POLYGON ((
    -69.62607888946572 -22.389720278609026, 
    -69.62632833490409 -22.39009971947581,
    -69.6252903200153 -22.3904196394022,
    -69.62501673469583 -22.390032759398007,
    -69.62553708324471 -22.389888919122058,
    -69.62607888946572 -22.389720278609026))')), 4326) &&
ST_SetSRID(ST_Point(-69.625137, -22.389777), 4326)

2：ST_Intersects，返回false，CORRECT

select 
ST_Intersects(
ST_SetSRID(
 'POLYGON ((
    -69.62607888946572 -22.389720278609026, 
    -69.62632833490409 -22.39009971947581,
    -69.6252903200153 -22.3904196394022,
    -69.62501673469583 -22.390032759398007,
    -69.62553708324471 -22.389888919122058,
    -69.62607888946572 -22.389720278609026))':: geography, 4326),
 ST_SetSRID(ST_Point(-69.625137, -22.389777):: geography, 4326))

3：ST_Contains，返回错误，正确

select 
ST_Contains(
ST_SetSRID(ST_GeomFromGeoJSON(ST_AsGeoJSON( 
 'POLYGON ((
    -69.62607888946572 -22.389720278609026, 
    -69.62632833490409 -22.39009971947581,
    -69.6252903200153 -22.3904196394022,
    -69.62501673469583 -22.390032759398007,
    -69.62553708324471 -22.389888919122058,
    -69.62607888946572 -22.389720278609026))':: geography)), 4326),
ST_SetSRID(ST_Point(-69.625137, -22.389777), 4326))

4：ST_DWithin，返回错误，正确

select 
ST_DWithin( 
 'POLYGON ((
    -69.62607888946572 -22.389720278609026, 
    -69.62632833490409 -22.39009971947581,
    -69.6252903200153 -22.3904196394022,    
    -69.62501673469583 -22.390032759398007,
    -69.62553708324471 -22.389888919122058,
    -69.62607888946572 -22.389720278609026))':: geography,
 ST_Point(-69.625137, -22.389777):: geography,0)

如果您想复制/粘贴geojson.io：

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "Point",
        "coordinates": [
          -69.625137,
          -22.389777
        ]
      }
    },
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "Polygon",
        "coordinates": [
          [
            [
              -69.62607888946572, 
              -22.389720278609026
            ],
            [
              -69.62632833490409,
              -22.39009971947581
            ],
            [
              -69.6252903200153,
              -22.3904196394022
            ],
            [
              -69.62501673469583,
              -22.390032759398007
            ],
            [
              -69.62553708324471,
              -22.389888919122058
            ],
            [
              -69.62607888946572,
              -22.389720278609026
            ]
          ]
        ]
      }
    }
  ]
}

Answer 1

<代码>

st_interacts不言自明。

对于面操作中的点，st_contains类似（尽管面边上的点不被视为包含在面中）。线与面或多边形与面之间非常不同，因为第一个要素必须位于另一个要素内部（而对于交点，如果 geom 1 的一部分位于 geom2 中，则相交为真）。

st_dwithin检查 2 个几何之间的最小距离是否小于或等于指定的距离。值为 0 时，它与 st_intersects 相同。