直接四舍五入经纬度小数位会导致高纬度网格面积缩水、密度失真,因经度距离随纬度升高而减半;正确做法是先投影(如epsg:3857)再按米级步长网格化。

为什么直接用 ROUND(lat, 2) + ROUND(lng, 2) 做网格会出错
因为经纬度的物理距离不均匀:赤道上经度1°≈111km,但到纬度60°时,1°经度只剩约55km;而纬度1°在任何地方都接近111km。直接四舍五入经纬度小数位,会导致高纬度区域网格实际面积大幅缩水,统计密度严重失真。
正确做法是先将经纬度转为平面近似坐标(如 Web Mercator 或等距矩形投影),再按固定米级步长切分。常见且轻量的方案是使用 ST_GeomFromText + ST_SnapToGrid(PostGIS),或手动实现“地理网格编码”逻辑。
- PostGIS 用户优先用
ST_SnapToGrid(ST_Transform(geom, 3857), 1000):先转 Web Mercator(EPSG:3857),再按1000米网格对齐 - MySQL 8.0+ 可用
ST_SRID(ST_Transform(geom, 3857), 3857)配合ST_X/ST_Y手动取整,但需注意 MySQL 对 3857 支持有限,建议降级用等距法 - 纯 SQL(无空间扩展)场景:用公式
FLOOR(lat / 0.01)+FLOOR(lng / 0.01)仅适用于小范围(如城市内),且需明确标注“非等面积”,避免误读密度
PostGIS 中用 ST_SnapToGrid 实现 500 米精度网格聚合
ST_SnapToGrid 本身不理解地理单位,必须配合坐标系转换。直接对 WGS84(EPSG:4326)几何调用它,参数 500 指的是“度”,毫无意义。关键在于先转投影、再网格化、最后转回便于阅读。
SELECT ST_AsText(ST_Transform(ST_SnapToGrid(ST_Transform(geom, 3857), 500), 4326)) AS grid_wkt, COUNT(*) AS cnt FROM locations WHERE geom IS NOT NULL GROUP BY ST_SnapToGrid(ST_Transform(geom, 3857), 500);
- 内部
ST_Transform(geom, 3857)把经纬度转成米制坐标(x/y 单位≈米) -
ST_SnapToGrid(..., 500)将每个点吸附到最近的 500×500 米网格左下角 - 外层
ST_Transform(..., 4326)是为了把网格角点转回经纬度方便调试,生产中可只保留网格唯一 hash(如MD5(ST_AsBinary(...))) - 注意:3857 在两极严重畸变,若数据含纬度 > 85° 或
MySQL 8.0 手动模拟地理网格(无 ST_Transform 时)
MySQL 8.0 内置空间函数不支持动态投影转换,ST_Transform 仅限于少数预设 CRS 且 3857 不在其中。此时稳妥做法是用“纬度校正系数”近似等距网格:
SELECT FLOOR(latitude / 0.0045) AS lat_grid, FLOOR(longitude / (0.0045 / COS(RADIANS(latitude)))) AS lng_grid, COUNT(*) AS cnt FROM locations GROUP BY lat_grid, lng_grid;
- 0.0045° ≈ 500 米(按赤道估算),
COS(RADIANS(latitude))补偿经度方向收缩 - 该公式在 ±60° 内误差
- 不能用于跨纬度大范围聚合(如全国),否则高纬网格被过度细分,低纬则过粗
- 若表中
latitude/longitude为空或非法,务必加WHERE latitude BETWEEN -90 AND 90 AND longitude BETWEEN -180 AND 180
网格 ID 设计:别用字符串拼接,用整型哈希
常见错误是写 CONCAT(FLOOR(lat/0.01), '_', FLOOR(lng/0.01)) 生成网格 ID —— 字符串 group by 效率低,且无法做范围查询(比如“查东边3个网格内的所有点”)。更优解是构造唯一整型 key。
- PostGIS:用
ST_XMin(grid_geom)::bigint * 1000000 + ST_YMin(grid_geom)::bigint(需确保坐标已转 3857) - 通用方案:
(FLOOR(lat / 0.0045) ,左移避免冲突,适合 INT32 - 若后续要叠加地图瓦片(如 XYZ),直接复用
tile_x, tile_y, zoom三元组,比自定义网格更易集成 - 永远不要在网格字段上建 TEXT 索引——改为 GENERATED COLUMN + STORED + INDEX
网格精度和投影选择不是技术炫技,而是直接影响统计结论是否可信。尤其当业务涉及“热力图”“服务覆盖半径”“竞品网点密度”时,差一个投影可能让决策南辕北辙。











