如何使用PostGIS查询指定经纬度2英里半径内的相交记录

你之前的查询无法得到正确结果，源于三个常见的PostGIS使用误区：

• 误用PostgreSQL原生CIRCLE类型：这不是PostGIS支持的空间对象，不携带SRID信息，和空间字段做运算时不会做坐标参考匹配，结果无参考价值。
• 经纬度坐标系单位错误：SRID 4269（NAD 83）是地理坐标系，坐标单位为度，不是米/英里这类长度单位，直接传入数值作为半径生成范围，得到的不是实际地面距离对应的圆形区域。
• 坐标顺序错误：PostGIS的ST_Point函数遵循经度在前、纬度在后的传参规则，你最开始的查询把纬度放在第一位，生成的点坐标位置完全错误。

优先使用ST_DWithin函数实现范围查询，相比先生成缓冲圆再做相交判断的写法，ST_DWithin不需要生成临时几何对象，可以直接命中空间索引，查询效率高1~2个数量级，且计算精度可靠。
由于SRID 4269是地理坐标系，直接转成GEOGRAPHY类型做大地测量计算即可，不需要手动匹配区域投影，写法最简便：

SELECT * FROM us_census_zcta 
WHERE ST_DWithin(
  geom::geography,
  ST_SetSRID(ST_Point(-73.935242, 40.730610), 4269)::geography,
  2 * 1609.344 -- 2英里换算为米，GEOGRAPHY类型计算距离的默认单位为米
);

注：1英里对应1609.344米，2英里对应半径为3218.688米。如果你的表还没给geom字段建空间索引，先执行CREATE INDEX idx_us_census_zcta_geom ON us_census_zcta USING GIST(geom);，查询速度会有质的提升。

如果你确实需要显式构造圆形几何对象做相交判断，必须先把坐标转换为对应区域的平面投影坐标系（单位为长度单位），生成缓冲后再转回原表SRID做匹配，示例（纽约区域适用NAD83 / UTM 18N投影，EPSG编码为26918）：

SELECT * FROM us_census_zcta 
WHERE ST_Intersects(
  geom,
  ST_Transform(
    ST_Buffer(
      ST_Transform(ST_SetSRID(ST_Point(-73.935242, 40.730610), 4269), 26918),
      2 * 1609.344
    ),
    4269
  )
);

• 所有PostGIS空间运算要求参与计算的几何对象SRID完全一致，否则会返回空结果或直接报错。
• 禁止在单位为度的地理坐标系下直接传入长度数值做缓冲、距离计算：1度经度对应的实际地面距离随纬度变化从0到111公里不等，直接计算得到的范围不是正圆，结果完全失真。
• 你之前调整后能返回1条记录的点查询，本质是判断中心点落在哪个ZCTA范围内，没有带2英里半径约束，不符合查询需求。

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