根据经纬度计算范围_亿级经纬度距离计算88.73秒，秒杀VBA！

最新推荐文章于 2023-03-20 13:07:26 发布

weixin_39640085

最新推荐文章于 2023-03-20 13:07:26 发布

阅读量395

点赞数

文章标签：根据经纬度计算范围

本文链接：https://blog.csdn.net/weixin_39640085/article/details/111706817

版权

本文介绍了一种高效的方法，能够在88.73秒内完成亿级经纬度范围的计算，对比传统的VBA计算方式，速度有了显著提升。

摘要由CSDN通过智能技术生成

计算经纬度的代码网上一搜一大把，通常是单点距离的计算，无法实现批量计算，本文将利用pandas实现亿级经纬度距离代码的实现。最短距离计算建议参考下文，mapinfo能够很好的实现。
MAPINFO 最小站间距统计本文将实现两张表的任意点之间100、200、300、500、800、1000米范围内的距离计算。首先导入需要使用的包

1import pandas as pd2import numpy as np3from math import radians, cos, sin, asin, sqrt, ceil4import math5import time

经纬度计算自定义函数

1def geodistance(lng1,lat1,lng2,lat2):2    lng1, lat1, lng2, lat2 = map(radians, [float(lng1), float(lat1), float(lng2), float(lat2)]) 3    # 经纬度转换成弧度4    dlon=lng2-lng15    dlat=lat2-lat16    a=sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**27    distance=2*asin(sqrt(a))*6371*1000 # 地球平均半径，6371km8    distance=round(distance,0)9    return distance

实现不同范围内的距离计算，例如100、200、300、500、800、1000，适合做成一张参数表。由于地球是球形，不同纬度下，同一经度差值对