055-【保姆级教程】手把手教你如何抓取高德全部POI-基于ArcGIS 渔网分割-六虎

这是坚持技术写作计划（含翻译）的第55篇，定个小目标999，每周最少2篇。

工作原因需要获取本地全部poi，直接调用高德api会有数量限制（理论上超过1000条appear就容易丢失部分数据，而且这种丢失很隐蔽，难以校验）

最终经过一番周折，基于网上资料进行了整理

注册高德账apple watch号，获取高德key
下载安装ArcMap/ArcGIS
将抓取区域等距程序员的黑科技系统分割并导出多边形经纬度
基于程序员切割后的经纬度重新轮高德api

前置准备

注册高德账号，并创建应用，获得key ，这个看http协议官网文档或者网上自己搜就行，不多枸杞说了

获取行政区划点位，lbs.amap.com/api/webserv…

复制出polyline,用任何一个现代化编辑器(vs程序员是做什么的 code,notepad++等)将 ;替换成换行符,在第一行插入 x,y保存成 jn.csv

根据 ArcMap 0 （ArcGIS10.2安装(完善版–能解决常见问题)）下载并安装ArcMaapproachp 10.2,并破解，但别汉化（后边点集转线的时候会报错）

获取分割后的坐标点位

导入点集

File->Add Data->Add XY Data…

建议如果不会用arcmap切换文件夹的，直接将文件放入 文档ArcGI程序员工资一般多少S 下,自动带入 x,y字段，点OK即可

出现行政区划点图像了

点集转线

ArcToolbox -> Data Management Tools -> Features -&apple tvgt;Points To Line

使用渔网功能等分切割行政区域

cell size widthapp安装下载/height这俩别选，因为高德导入的是火星坐标，不属于标准gps坐标，ArcGis 不识别

所以只能自行计算，假设我需要将济南切割成3kapple idm X 3http://192.168.1.1登录km的小格httpwatch子，参考经纬度1度等于多少米Go，经度1度=85.39km，纬度1度 = 大约111km。那就拿出计算器算呗，(Right-Left)*85.39/3=49.78826191,取整用http://www.baidu.com50就行，同理 (Topappleid-Bottom)*111/3=程序员怎么学57.060475,取整用57就行，我截图用的58，懒得换了。

注意别选 Create Label Points(optional) 选项,因为用不到

将行政区划填充颜色

ArcToolbox -> Data Management Tools -> Features -> Featapple idurapproache To Polygonhttp://192.168.1.1登录

取渔网和行政区划交集（擦除行工作细胞政区划外的渔网网格部分）

将httpwatch其导出后边要用

获取网格点位的经appreciate纬度

参考大佬的文章 Python3爬取高德地图POI ，下载作者开发的渔网对角坐标获取工具
链接：链接：pan.baidu.com/s/11nwWHDf7…
提取码：rdu程序员怎么学v

获application取对角坐标

抓取高德POI

98%都是抄自 Python3爬取高德地图POappreciateI 代码，只是稍微处理了下从dict获取字段防止抛异常部分，以及加了省字段

import requests
import json
from papple tvymongo import MongoClient
import time
client = MongoClient('localhost',27017)
db = client.POI_Jinanshi
collection = db.table_1
polygon_listapp安装下载 = list()
with opappearen("jn2.txt",程序员工资一般多少 'r', encoding='UTF-8') as txt_file:
for each_line in txt_file:
if each_line != "" and each_line != "n":
fields = each_line.split("n")
polygon = fields[0]
polygon_list.append(papple watcholygon)
def getjson(polygon, page):
headers = {
'User-Agent': 'Mozilla/5.0 (Windows NT 10.0; WOW64) AppleWebKit/537.36 (KHTML, like程序员那么可爱 G程序员需要什么学历ecko) Chrome/70.0.3538.67 SafaAppleri/537.36'
}
pa = {
'key': 'xxxxxxxxxxxxxxxxxxxxxxhttp代理', #从控制台申请
'polygon': polygon,
'app storetypes':'970000|990000',    #不要过多
'httpclientcity':'0531',
'offset': 20,
'page': page,
'extensions': 'allgoogleplay',
'output': 'JSON'
}
r = requests.get('https://restapi.amap.com/v3/place/polygon?', params=pa, headers=headers)
decodejson = json.loads(r.text)
return decodejson
for each_po狗狗币lygon in polygon_list:
not_last_Gopage = True
page宫外孕 = 1
while not_la程序员是做什么的st_page:
decodejson = getjson(eahttpwatchch_polygon, page)
print(decodejson)
count = decodejson.get('HTTPcount',0)
print(each_polygon, page)
if decodejson['pois']:
for eachone in decodejson['pois']:
data={
'name':eachone.get('name',None), #POI名称
'typeshttp协议':eachone.get('type',None), #POI所属类别
'address':eachone.get('address',None), #POI地址
'location':eappearachone.get('lapple watchocation',Non程序员工资一般多少钱一个月e), #POI坐标
'city':eachone.get('cityname',None), #城市
'co程序员的黑科技系统unty':eachone.get('adname',None), #区县
'province':eachone.get('pna工作总结me',None), #程序员工资一般多少钱一个月 省份
'count'apple pay:count,http协议 # 条数
'polygon':each_polygon # 多https和http的区别边形经纬度，便于后边再次抓取
}
collection.apple watchinsert_one(datHTTPa)
time.sleep(0.2)
page += 1
else:
not_last_page = Fals程序员的黑科技系统e

招聘小广告

山东济南的小伙伴欢迎投简历啊加入我们 , 一apple store起搞app id注册事情。
长期招聘，Java程序员，大数据工程师，运维工程师，前APP端工程师。

参考资料

我的博客
我的
ArcMap 0 （ArcGIS10.2安装(完善版–能解决常见问题)）
Python3爬appleid取高德地图POI
pythapp id注册on获取高德papple id密码重置oi
经纬度1度等于多少米

055-【保姆级教程】手把手教你如何抓取高德全部POI-基于ArcGIS 渔网分割

前置准备

获取分割后的坐标点位

导入点集