简单的经纬度点聚类
import logging logging.basicConfig(level=logging.INFO, format='%(asctime)s-%(filename)s[line:%(lineno)d]-%(levelname)s:%(message)s', datefmt='%Y-%m-%d %H:%M:%S') from geopy.distance import geodesic def stop2label(pnts, near=300): # 基于300m的半径对位置进行合并 # pnts [(lng,lat), (), *] res = [] pfun = lambda p1, p2: round(geodesic(p1[::-1], p2[::-1]).m, 1) count = 1 while pnts.__len__(): ipnt = pnts[0] flags = [pfun(ipnt, _) < near for _ in pnts] for idx, flag in enumerate(flags): if flag: res.append((pnts[idx], count)) logging.info((flags, len(pnts))) pnts = [pnts[_] for _, flag in enumerate(flags) if not flag] count += 1 return res
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步