如何将单个经纬度点转为PostGIS 30米范围多边形（适配谷歌地图）

我刚好之前踩过一模一样的坑！核心原因是你没注意到坐标系的单位差异——直接用WGS84（EPSG:4326）的经纬度做缓冲，单位是「度」而不是「米」，这就难怪范围会超大了（1度大概相当于111公里，想想你缓冲0.00027度才是30米，而且不同纬度还不准）。下面给你两个靠谱的解决方案，都是PostGIS+Ruby的实操步骤：

PostGIS里的ST_Buffer函数是基于当前坐标系的单位计算的：

• 当你用经纬度（EPSG:4326，地理坐标系）时，缓冲的单位是度，不是米
• 要实现米级缓冲，必须先把点转换到投影坐标系（单位是米），缓冲后再转换回经纬度坐标系

谷歌地图默认使用Web墨卡托（EPSG:3857）作为底图坐标系，这个方案简单直接，适合大多数场景：

对应的PostGIS SQL语句 #

SELECT ST_AsGeoJSON(
  ST_Transform(
    ST_Buffer(
      ST_Transform(
        ST_SetSRID(ST_MakePoint(:longitude, :latitude), 4326),
        3857 -- 转换为Web墨卡托（单位：米）
      ),
      30 -- 这里的30就是你要的米级范围
    ),
    4326 -- 转换回经纬度坐标系，方便谷歌地图解析
  )
) AS polygon_geojson;

语句拆解 #

1. ST_MakePoint(:longitude, :latitude)：用传入的经纬度创建点
2. ST_SetSRID(..., 4326)：标记这个点的坐标系是WGS84经纬度
3. ST_Transform(..., 3857)：把点转换为Web墨卡托投影（单位变为米）
4. ST_Buffer(..., 30)：围绕点生成30米半径的缓冲多边形
5. ST_Transform(..., 4326)：把多边形转换回经纬度坐标系
6. ST_AsGeoJSON(...)：输出GeoJSON格式，方便谷歌地图直接渲染

Web墨卡托在高纬度地区会有明显的拉伸变形，如果你的业务覆盖全球，推荐用UTM坐标系（每个UTM带对应6度经度范围，精度更高）：

对应的PostGIS SQL语句 #

SELECT ST_AsGeoJSON(
  ST_Transform(
    ST_Buffer(
      ST_Transform(
        ST_SetSRID(ST_MakePoint(:longitude, :latitude), 4326),
        -- 动态计算当前点对应的UTM带SRID（北半球用326xx，南半球用327xx）
        ST_UTMZone(ST_SetSRID(ST_MakePoint(:longitude, :latitude), 4326)) + CASE WHEN :latitude >= 0 THEN 32600 ELSE 32700 END
      ),
      30
    ),
    4326
  )
) AS polygon_geojson;

关键说明 #

ST_UTMZone函数会返回当前点所在的UTM带编号（比如中国大部分地区是50带），再根据纬度加上对应的基数（北半球32600，南半球32700），就能得到正确的UTM坐标系SRID。

# 假设你有传入的经纬度参数
target_lon = 116.397228
target_lat = 39.9075

# 执行SQL查询
result = ActiveRecord::Base.connection.execute(
  "SELECT ST_AsGeoJSON(ST_Transform(ST_Buffer(ST_Transform(ST_SetSRID(ST_MakePoint($1, $2), 4326), 3857), 30), 4326)) AS polygon",
  target_lon, target_lat
).first

# 获取生成的GeoJSON字符串，直接传给谷歌地图渲染
polygon_geojson = result['polygon']

1. 确保你的PostGIS版本支持ST_UTMZone函数（PostGIS 2.2及以上版本支持）
2. 谷歌地图的google.maps.Polygon组件可以直接解析GeoJSON的坐标数组，不需要额外转换
3. 如果是批量处理点，建议把坐标系转换的逻辑封装成PostGIS函数，提高效率

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