PostGIS如何判断geography类型地理点是否位于多边形范围内

查询返回错误结果或无结果，核心是两个基础错误：

1. 点坐标构造顺序错误
   空间函数遵循WKT标准，ST_POINT的参数顺序为经度在前，纬度在后，写入数据时把latitude放在第一个参数、longitude放在第二个参数，所有点的坐标完全错位。
2. 空间类型与坐标系不匹配
   存储字段是geography类型，构造多边形时用ST_GEOMFROMTEXT默认生成无坐标系声明的geometry类型，两边类型、SRID不统一，空间运算无法得到正确结果。

1. 修正存量geo_point数据 #

重新执行更新语句，调整经纬度顺序，显式指定WGS84坐标系（SRID=4326，对应通用GPS经纬度）：

UPDATE addresses SET geo_point = ST_POINT(longitude, latitude, 4326);

2. 执行正确的范围查询 #

优先使用ST_Covers替代ST_Contains，避免多边形边界上的点被漏判，同时统一使用geography类型做球面运算，精度更高：

SELECT * FROM addresses a 
WHERE ST_Covers(
    ST_GeogFromText('SRID=4326;POLYGON((
        -104.95768433837907 39.71956203404201,
        -104.95768433837907 39.75883237170433,
        -105.02291566162134 39.75883237170433,
        -105.02291566162134 39.71956203404201,
        -104.95768433837907 39.71956203404201
    ))'),
    a.geo_point
);

• 如果表数据量超过1万条，给geo_point字段创建GIST空间索引，查询性能会提升10~100倍：

CREATE INDEX idx_addresses_geopoint ON addresses USING GIST(geo_point);

• 如果习惯用ST_Contains，需要把两边的空间类型、SRID完全对齐，写法如下：

SELECT * FROM addresses a 
WHERE ST_Contains(
    ST_GeomFromText('POLYGON((
        -104.95768433837907 39.71956203404201,
        -104.95768433837907 39.75883237170433,
        -105.02291566162134 39.75883237170433,
        -105.02291566162134 39.71956203404201,
        -104.95768433837907 39.71956203404201
    ))', 4326),
    a.geo_point::geometry
);

内容的提问来源于stack exchange，提问作者engineer715154