如何在SQL中对地理位置坐标进行网格分组统计_利用经纬度换算函数

轻萱同学_7684

轻萱同学_7684

2026-04-27

786人浏览

原创

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

如何在sql中对地理位置坐标进行网格分组统计_利用经纬度换算函数

为什么直接用 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 指的是“度”,毫无意义。关键在于先转投影、再网格化、最后转回便于阅读。

拍我AI
拍我AI

一款AI视频创作工具,主要用于AI视频生成平台PixVerse的国内版本,适合需要提升相关任务效率的用户。

下载
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

网格精度和投影选择不是技术炫技,而是直接影响统计结论是否可信。尤其当业务涉及“热力图”“服务覆盖半径”“竞品网点密度”时,差一个投影可能让决策南辕北辙。

相关专题

更多
数据分析工具有哪些
数据分析工具有哪些

数据分析工具有Excel、SQL、Python、R、Tableau、Power BI、SAS、SPSS和MATLAB等。详细介绍:1、Excel,具有强大的计算和数据处理功能;2、SQL,可以进行数据查询、过滤、排序、聚合等操作;3、Python,拥有丰富的数据分析库;4、R,拥有丰富的统计分析库和图形库;5、Tableau,提供了直观易用的用户界面等等。

2023.10.12

4463

8

SQL中distinct的用法
SQL中distinct的用法

SQL中distinct的语法是“SELECT DISTINCT column1, column2,...,FROM table_name;”。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

2023.10.27

931

4

SQL中months_between使用方法
SQL中months_between使用方法

在SQL中,MONTHS_BETWEEN 是一个常见的函数,用于计算两个日期之间的月份差。想了解更多SQL的相关内容,可以阅读本专题下面的文章。

2024.02.23

1169

5

SQL出现5120错误解决方法
SQL出现5120错误解决方法

SQL Server错误5120是由于没有足够的权限来访问或操作指定的数据库或文件引起的。想了解更多sql错误的相关内容,可以阅读本专题下面的文章。

2024.03.06

6421

10

sql procedure语法错误解决方法
sql procedure语法错误解决方法

sql procedure语法错误解决办法:1、仔细检查错误消息;2、检查语法规则;3、检查括号和引号;4、检查变量和参数;5、检查关键字和函数;6、逐步调试;7、参考文档和示例。想了解更多语法错误的相关内容,可以阅读本专题下面的文章。

2024.03.06

3203

4

oracle数据库运行sql方法
oracle数据库运行sql方法

运行sql步骤包括:打开sql plus工具并连接到数据库。在提示符下输入sql语句。按enter键运行该语句。查看结果,错误消息或退出sql plus。想了解更多oracle数据库的相关内容,可以阅读本专题下面的文章。

2024.04.07

6380

11

sql中where的含义
sql中where的含义

sql中where子句用于从表中过滤数据,它基于指定条件选择特定的行。想了解更多where的相关内容,可以阅读本专题下面的文章。

2024.04.29

8721

6

sql中删除表的语句是什么
sql中删除表的语句是什么

sql中用于删除表的语句是drop table。语法为drop table table_name;该语句将永久删除指定表的表和数据。想了解更多sql的相关内容,可以阅读本专题下面的文章。

2024.04.29

1190

5

sql中删除一列的命令是什么
sql中删除一列的命令是什么

在sql中,使用alter table语句可以删除一列,语法为:alter table table_name drop column column_name。想了解更多sql的相关内容,可以阅读本专题下面的文章。

2024.04.29

1032

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程