🌐 amap 根据经纬度快速查询全国行政区域信息、省市区信息
行政区域查询，根据经纬度快速地查找特定的省市区、行政区域信息。 https://wuwenrufeng.github.io/amap/
Features
- 免费，免费，免费
- 打破官方请求量限制，一次生产，终生使用
- 数据源可靠，数据源来自高德地图开放的HTTP接口
- 手把手教你自建全国行政区域信息库
- 手把手教你实现行政区域信息查询算法
- 快速地查找特定的行政区域信息
- APi接口调用
项目由来
官方服务有请求限制，又不想升级plus
配置环境
运行环境：需要运行在 python 3.7.4+
1. 安装依赖
git clone https://github.com/wuwenrufeng/amap.git
pip install -r requirements.txt快速开始
1. 下载数据库文件cn_area.sql,并导入到mysql中
链接：https://pan.baidu.com/s/1zJZTccPjrY4RD_TgAUh84A 提取码：cmuv
2. 数据库访问配置
文件地址：amap/utils/settings.py
# mysql
MYSQL_PARAM = {
'provider': 'mysql',
'host': 'localhost',
'port': 3306,
'user': 'admin',
'passwd': '123456',
'db': 'cn_area'
}3. 开启api查询服务
python amapSearchApi.py4. 开始查询
URL: http://127.0.0.1/search?point=117.195907,39.118327
使用手册
第一步：高德地图-初始化出各省级(有边界信息)、市级和区级的初始信息
python amapInit.py第二步：高德地图-补充市、区的边界信息 (ps: 高德不提供有关街道的边界顶点信息)
python amapPopulate.py第三步：计算出province, city, district表中的最大最小边界经纬度，并填充到数据库中
python amapMaxMin.py第四步：开启服务，查询给定坐标点所属的行政区域信息
python amapSearchApi.py点在区域内的判定算法
使用外包矩形+闭合路径包含算法
一级过滤：外包矩形算法
def in_box(point, maxP, minP):
"""
判断点是否在多边形的外包矩形中
point: 待判定的坐标点
maxP: 最大坐标点
minP: 最小坐标点
"""
return (minP.x <= point.x <= maxP.x) and (minP.y <= point.y <= maxP.y)二级过滤：闭合路径包含算法
def in_polyline(point, polyline):
"""
根据经纬度，判断是否在区域内
point: 待判定的坐标点
polyline: 边界坐标
"""
lats = []
lngs = []
for pol in polyline:
lats.append(pol.x)
lngs.append(pol.y)
xc = np.array(lats)
yc = np.array(lngs)
xycrop=np.vstack((xc,yc)).T
pth=Path(xycrop,closed=False)
mask=pth.contains_points([[point.x,point.y]])
if mask:
return True
return False测试结果
测试环境：个人PC
1.查询耗时
2.查询准确率
PS: 准确率测试是使用区级center进行查询的，所得准确率应该是100%，丢失的准确率是因为官方将center标记错误了！
比如抚顺县的实际地理区域是不包含抚顺县的，但官方的center却标到了抚顺县。
后续更新。。。
GPS坐标、mapbar坐标、baidu坐标转换成高德坐标
逆地址地理信息解析
交流 & 讨论
如果您对amap有任何疑问，欢迎添加作者微信(微信号：wuwenrufeng)来交流讨论。或者，您可以扫下方二维码给作者打赏去升级其他地图api套件或买一杯咖啡哈
