# 代码格式需要自行优化# 通过读取GeoJSON格式的上海边界数据# 获取上海边界外接矩形的最小经度、最大经度、最小纬度、最大纬度import geopandas as gp# 结果变量resultMinLon = 99999.0# 最小经度resultMaxLon = 0.0# 最大经度resultMinLat = 99999.0# 最小纬度resultMaxLat = 0.0# 最大纬度# 读取上海边界数据loadGeoData = gp.read_file("D:\\上海边界数据.geojson")for i in range(0, len(loadGeoData)):# 读取几何数据    loadGeometry = loadGeoData.loc[i, "geometry"]# 上海边界的几何数据是由多个多边形组成，如崇明岛、长兴岛、上海市区等# 所以需要遍历所有多边形# 多边形的类型是Polygon或者MultiPolygonfor j in range(0, len(loadGeometry)):        loadPolygon = loadGeometry[j]# 读取多边形的边界# 边界的类型是LineString(单线)或MultiLineString(多线)if loadPolygon.boundary.geom_type == 'LineString':# 读取边界中的每个坐标点for z in range(0, len(loadPolygon.boundary.coords)):                loadLon = loadPolygon.boundary.coords[z][0]  # 经度                loadLat = loadPolygon.boundary.coords[z][1]  # 纬度# 进行大小判断，并替换结果if loadLon <= resultMinLon:                    resultMinLon = loadLonif loadLon >= resultMaxLon:                    resultMaxLon = loadLonif loadLat <= resultMinLat:                    resultMinLat = loadLatif loadLat >= resultMaxLat:                    resultMaxLat = loadLatelif loadPolygon.boundary.geom_type == 'MultiLineString':# 如果边界是MultiLineString，则需要遍历其中的每条LineStringfor z in range(0, len(loadPolygon.boundary[0].coords)):                loadLon = loadPolygon.boundary[0].coords[z][0]                loadLat = loadPolygon.boundary[0].coords[z][1]if loadLon <= resultMinLon:                    resultMinLon = loadLonif loadLon >= resultMaxLon:                    resultMaxLon = loadLonif loadLat <= resultMinLat:                    resultMinLat = loadLatif loadLat >= resultMaxLat:                    resultMaxLat = loadLat# 结果输出print("上海最小经度：", resultMinLon)print("上海最大经度：", resultMaxLon)print("上海最小纬度：", resultMinLat)print("上海最大纬度：", resultMaxLat)

步骤2输出结果

步骤3：以四叉树索引的方式不断四分矩形，直到所有矩形符合POI数量限制

# 代码格式需要自行优化import requests# 上海边界经纬度值minLon = 120.859465  # 最小经度maxLon = 122.122756  # 最大经度minLat = 30.677191  # 最小纬度maxLat = 31.868458  # 最大纬度# 存放符合POI数量限制多边形的列表resultPolygonList = []# 用来存放不符合POI数量限制多边形的列表currentPolygonList = [[minLon, maxLon, minLat, maxLat]]# 以四叉树的方式不断四分外接矩形for currentLevel in range(1, 9999999):# 用来存放不符合POI数量限制多边形的临时列表tempPolygonList = []for z in range(0, len(currentPolygonList)):loadMinLon = currentPolygonList[z][0]loadMaxLon = currentPolygonList[z][1]loadMinLat = currentPolygonList[z][2]loadMaxLat = currentPolygonList[z][3]# 判断当前矩形是否符合POI数量限制# 如果符合则放入resultPolygonList# 如果不符合则进行四分，并将四分结果放入tempPolygonListifSatisfy = judgeIfSatisfy_GaoDe([loadMinLon, loadMaxLon, loadMinLat, loadMaxLat])if ifSatisfy == True:resultPolygonList.append([loadMinLon, loadMaxLon, loadMinLat, loadMaxLat])else:# 进行四叉树tempQuadtree = executeQuadtree(loadMinLon, loadMaxLon, loadMinLat, loadMaxLat)# 左下角tempPolygonList.append(tempQuadtree[0])# 右下角tempPolygonList.append(tempQuadtree[1])# 右上角tempPolygonList.append(tempQuadtree[2])# 左上角tempPolygonList.append(tempQuadtree[3])currentPolygonList = tempPolygonList.copy()# 如果没有多边形进入下一层级，则跳出print("第", currentLevel, "层", "......", "符合要求的矩形：", len(resultPolygonList), "......", "剩余矩形：",len(currentPolygonList))if len(currentPolygonList) == 0:break

步骤3结果

点击文末的阅读原文，可在电子地图中可查看步骤3生成的矩形

步骤4：遍历所有矩形，获取上海瑞幸POI

# 代码格式需要自行优化# 遍历所有矩形，获取上海瑞幸POI数据savePoiList = []  # 存放POI数据的结果列表for i inrange(0, len(resultPolygonList)):# 读取每个矩形    loadPolygon = resultPolygonList[i]# 获取矩形中的POI数据，并添加到结果列表中    savePoiList.extend(getPoiFromPolygon(loadPolygon, gaodeKey))

下面代码为获取单个矩形中POI数据的方法

# 代码格式需要自行优化# 从单个矩形获取瑞幸POI的方法# 其中key为高德地图API密钥，需要自行申请defgetPoiFromPolygon(inputPolygon, key):# 存放POI数据结果的列表    resultList = []# 由于POI数量限制是100，每页返回POI数量是20# 因此，最大的页数为5for currentPage inrange(1, 6):# 组成Url        currentUrl = "https://restapi.amap.com/v3/place/polygon?polygon=" + str(inputPolygon[0]) + "," + str(            inputPolygon[2]) + "|" + str(inputPolygon[1]) + "," + str(            inputPolygon[2]) + "|" + str(inputPolygon[1]) + "," + str(inputPolygon[3]) + "|" + str(            inputPolygon[0]) + "," + str(inputPolygon[3]) + "|" + str(inputPolygon[0]) + "," + str(            inputPolygon[2]) + "&offset=20&page=" + str(            currentPage) + "&keywords=瑞幸&types=050000&output=json&key=" + key# 发送Get请求，并接收返回内容        response = requests.get(currentUrl, stream=True, verify=False, timeout=60)        returnData = response.json()        returnPoiList = returnData['pois']iflen(returnPoiList) > 0:for i inrange(0, len(returnPoiList)):                saveName = returnPoiList[i]['name']                saveType = returnPoiList[i]['type']                saveAddress = returnPoiList[i]['address']                saveLocation = returnPoiList[i]['location']                saveProvince = returnPoiList[i]['pname']                saveCity = returnPoiList[i]['cityname']                saveArea = returnPoiList[i]['adname']                resultList.append(                    [saveName, saveType, saveType, saveAddress, saveLocation, saveProvince, saveCity, saveArea])print(saveName, saveType, saveType, saveAddress, saveLocation, saveProvince, saveCity, saveArea)else:# 如果当前页POI数量为0，则返回已获取的POI数据并跳出return resultListreturn resultList

上海瑞幸POI数据

获取方式

第一步：点击下方链接或扫描二维码，关注地学大数据公众号

第二步：在公众号对话框中输入资料获取，即可获取

END

本文来自网友投稿或网络内容，如有侵犯您的权益请联系我们删除，联系邮箱：wyl860211@qq.com 。

【数据技术】利用Python获取高德地图POI数据——以上海瑞幸咖啡门店为例

最新文章

热门文章

随机文章

【数据技术】利用Python获取高德地图POI数据——以上海瑞幸咖啡门店为例

Python实现发票信息识别

使用PyCharm+Python程序控制舵机转动

最新文章

热门文章

随机文章