python 经纬度计算相关

最新推荐文章于 2023-01-03 13:48:56 发布

原创最新推荐文章于 2023-01-03 13:48:56 发布 · 888 阅读

2 ·

CC 4.0 BY-SA版权

文章标签：

#python

python 专栏收录该内容

23 篇文章

订阅专栏

本文介绍了如何使用高德API获取地址的经纬度信息，并提供了两种基于经纬度计算两点间距离的方法。此外，还探讨了经纬度划分区域的技术，用于大数据量下的区域划分，提高数据处理效率。

1，高德接口返回地址经纬度

import requests

def getcode(site):
    parameters = { 'address' : site, 'key': '******高德平台申请******'}
    base = 'http://restapi.amap.com/v3/geocode/geo'
    response = requests.get(base, parameters)
    info_site = response.json()
    lng = info_site['geocodes'][0]['location'].split(',')[0]
    lat = info_site['geocodes'][0]['location'].split(',')[1]
    return [lng,lat]

2，利用经纬度计算距离

from math import *

def Distance1(Lat_A,Lng_A,Lat_B,Lng_B): #第一种计算方法
    ra=6378.140 #赤道半径
    rb=6356.755 #极半径 （km）
    flatten=(ra-rb)/ra  #地球偏率
    rad_lat_A=radians(Lat_A)
    rad_lng_A=radians(Lng_A)
    rad_lat_B=radians(Lat_B)
    rad_lng_B=radians(Lng_B)
    pA=atan(rb/ra*tan(rad_lat_A))
    pB=atan(rb/ra*tan(rad_lat_B))
    xx=acos(sin(pA)*sin(pB)+cos(pA)*cos(pB)*cos(rad_lng_A-rad_lng_B))
    c1=(sin(xx)-xx)*(sin(pA)+sin(pB))**2/cos(xx/2)**2
    c2=(sin(xx)+xx)*(sin(pA)-sin(pB))**2/sin(xx/2)**2
    dr=flatten/8*(c1-c2)
    distance=ra*(xx+dr)
    return distance
def Distance2(lat1,lng1,lat2,lng2):# 第二种计算方法
    radlat1=radians(lat1)  
    radlat2=radians(lat2)  
    a=radlat1-radlat2  
    b=radians(lng1)-radians(lng2)  
    s=2*asin(sqrt(pow(sin(a/2),2)+cos(radlat1)*cos(radlat2)*pow(sin(b/2),2)))  
    earth_radius=6378.137  
    s=s*earth_radius  
    if s<0:  
        return -s  
    else:  
        return s

#验证
Lat_A=32.060255; Lng_A=118.796877 # 南京
Lat_B=39.904211; Lng_B=116.407395 # 北京
r_id = 'nanjing'; t_id = 'beijing'
distance=Distance(t_id,r_id,lat_a,lng_a,lat_b,lng_b)
Lat_A = [32.060255,31.060255];Lng_A=[118.796877]
Lat_B = [39.060255,39.060255];Lng_B=[116.796877]
print('(Lat_A, Lng_A)=({0:.6f},{1:.6f})'.format(Lat_A,Lng_A))
print('(Lat_B, Lng_B)=({0:.6f},{1:.6f})'.format(Lat_B,Lng_B))
print('Distance1={0:.3f} km'.format(distance))
print('Distance2={0:.3f} km'.format(Distance2(Lat_A,Lng_A,Lat_B,Lng_B)))

3，经纬度划分区域（数据量多的时候可以用来判断数据分成那几块区域）

_base32 = '0123456789bcdefghjkmnpqrstuvwxyz'
#10进制和32进制转换，32进制去掉了ailo
_decode_map = {}
_encode_map = {}
for i in range(len(_base32)):
    _decode_map[_base32[i]] = i
    _encode_map[i]=_base32[i]
del i



def encode(lat,lon,precision=7):
    '''
    编码位数对应区域大小（精度）
    precision =  5 ±2.4km
    precision =  6 ±0.61km
    precision =  7 ±0.076km
    precision =  8 ±0.019km
    '''
    lat_range, lon_range = [-90.0, 90.0], [-180.0, 180.0]
    geohash=[]
    code=[]
    while len(geohash)<precision:
#         print(code,lat_range,lon_range,geohash)
        lat_mid=sum(lat_range)/2
        lon_mid=sum(lon_range)/2
        #经度
        if lon<lon_mid:
            code.append(0)
            lon_range[1]=lon_mid
        else:
            code.append(1)
            lon_range[0]=lon_mid
        #纬度
        if lat<lat_mid:
            code.append(0)
            lat_range[1]=lat_mid
        else:
            code.append(1)
            lat_range[0]=lat_mid
        ##encode
        if len(code)>=5:
            geohash.append(_encode_map[int(''.join(map(str,code[:5])),2)])
            code=code[5:]
    return ''.join(geohash)