根据经纬度坐标获得省市区县行政区划城市名称，自建数据库 java python php c# .net 均适用

在LBS应用中，根据坐标来解析获得对应是哪个城市是一个很常见的功能，比如App里面通过手机定位自动选择城市；本文介绍的是通过自己建的数据库，利用SQL空间查询来进行坐标解析得到对应的省市区，绝大部分支持空间数据类型(Spatial)的数据库均支持，包括但不限于：MySQL、SQL Server、Oracle、PostgreSQL等；开发语言不限，只要能进行数据库查询就都支持，用Java、Python、PHP、C#均能很简单的实现。

省市区边界数据在线预览、下载：https://xiangyuecn.gitee.io/areacity-jsspider-statsgov/
GitHub地址：https://github.com/xiangyuecn/AreaCity-JsSpider-StatsGov

通过坐标来获取地址：百度地图提供的接口叫 “地址逆解析”，高德地图提供的接口叫 “地理逆编码”，它们开放平台均提供了前后端接口，发个http请求就可以拿到数据，相关的接口调用请自行阅读平台开发文档，使用起来很方便。

当然这篇文章不会去介绍怎么搞街道门牌号这些完整地址，也不讲怎么样调用人家的接口，只介绍坐标对应的省市区名称的获取，自建数据库写SQL进行空间查询，完全自己实现，比调开放平台接口相对会复杂些。

由于全国区县每年都会有比较多的变更，所以省市区边界数据需要经常去同步维护，好在这上面这个开源库有在持续的长期维护，新数据发布后更新相对容易很多。由于开源库更新维护数据比较及时，所以只要开源库没有被关闭，本文介绍的提取方法就一直适用； 比那些上传到下载平台就万年不更新的数据优秀很多。

查询效果展示：
数据库查询效果

从这里随便拿的一个坐标：
坐标位置

直观的效果如上图，在百度地图里面随便点选一下（或App定位获得的坐标）得到一个坐标，然后到数据库里面利用空间查询SQL就能查询到坐标所在的城市。

步骤一、下载省市区边界数据

到开源库里面下载最新的 ok_geo.csv.7z 文件（13M大小，解压130M+），点此下载。下载好后解压得到 ok_geo.csv，此文件内包含了最新全国所有的省市区县坐标边界矢量数据。

注：这个文件只包含三级（省市区）数据不含第四级（乡镇街道），如需乡镇坐标边界数据，可以请点此下载 ok_geo4_*.csv 文件（90MB+压缩包导出后300M+）。

步骤二、解析CSV文件导入数据库

下载好的文件 ok_geo.csv 是纯文本文件，可以自行编写脚本进行解析，然后导进数据库中，自行解析处理比较复杂，请参考开源库内文档；开源库内提供了一个格式转换工具，此工具支持将CSV数据导入数据库，因此我们直接在下载数据时顺带把工具下载好，点此下载。

此转换工具除了支持将 ok_geo.csv 导入数据库外，还支持导出：sql、shp、geojson，和坐标系转换；还可执行自定义 JavaScript 脚本，扩展出丰富功能；软件是Windows版，如果需要在MacOs中用，可以用虚拟机。

导入操作

转换工具执行导入数据库操作：

点击 选择ok_geo.csv文件 按钮，选择解压出来的CSV文件；
数据库设置中选择要导入的数据库类型，这里选的是MySQL，再填写数据库连接，包括：端口、数据库名称、账号密码；
点击导入数据库按钮，等待一会，大约3分钟左右，所有数据就都导入到了数据库按今天日期新建的表里面。

注：csv文件内的边界数据默认是高德地图GCJ-02火星坐标系，如果需要其他坐标系，比如百度的BD-09或GPS的WGS-84，可以通过高级脚本中的坐标系转换插件进行转换，选择好对应的插件后，点击应用就可以了，在导入数据库时会自动进行坐标系转换。

注：这个工具限制每次操作只可导出一个城市和它的下一级数据，导出少量数据还是很轻松的，所以我们可以多操作几次，将需要的城市数据全部导入数据库；比如要深圳的所有区县数据：第一遍导入全国所有的省，第二遍在城市名前缀中填写广东省（结尾带一个空格）导入广东所有的市，第三遍在城市名前缀中填写广东省深圳市（结尾带一个空格）导入深圳所有的区。如果在密钥输入框中填写了密钥，此工具就没有这些限制了，开源库里面会不定期发放密钥搞福利，点击一次操作即可导出全国所有的省市区三级数据。

表结构和空间字段（MySQL版，其他数据库类似）：

CREATE TABLE Areacity_Geo_20220216 (
  id int NOT NULL, --城市id
  pid int NOT NULL, --上级城市id
  deep int NOT NULL, --层级：0省、1市、2区
  name varchar(250) NOT NULL, --城市名称：`深圳市`
  ext_path varchar(255) NOT NULL, --省市区三级完整名称：广东省 深圳市 罗湖区
  geo geometry NOT NULL, --城市中心坐标，空间数据格式
                --，ST_AsText转成WKT文本后：`POINT EMPTY`、`POINT (123.456 34.567)`
  polygon geometry NOT NULL --城市边界范围图形，空间数据格式
                --，ST_AsText转成WKT文本后：`POLYGON EMPTY`、`POLYGON ((123.456 34.567,...))`、`MULTIPOLYGON (((123.456 34.567,...)),...)`
)

对空间字段的查询，需要用`ST_AsText()`方法才能查询出字符串文本(WKT: Well Known Text)，否则查询出来的是二进制数据
-- MySQL版：
SELECT id, name, ST_AsText(geo) AS geo, ST_AsText(polygon) AS polygon FROM 表名
-- SQL Server版：
SELECT id, name, geo.STAsText() AS geo, polygon.STAsText() AS polygon FROM 表名

步骤三、在程序中根据坐标解析获得城市

省市区边界导入到了数据库后，我们就可以在在 Java、Python、PHP、C# 等程序中对数据库进行查询，通过SQL的空间计算函数ST_Intersects来查询一个坐标在哪些边界范围内，就能得到对应的省市区信息了。

空间查询SQL语句

比如要查询坐标`lng:113.929976 lat:22.529497`是在哪个城市
-- MySQL版：
SELECT id,deep,name FROM 表名 WHERE ST_Intersects(polygon, ST_GeomFromText('POINT(113.929976 22.529497)',0))=1
-- SQL Server版：
SELECT id,deep,name FROM 表名 WHERE polygon.STIntersects(geometry::STGeomFromText('POINT(113.929976 22.529497)',0))=1