计算三个经纬度坐标的夹角

import math

def geo2xyz(lat, lng, r=6400):
    '''
    将地理经纬度转换成笛卡尔坐标系
    :param lat: 纬度
    :param lng: 经度
    :param r: 地球半径
    :return: 返回笛卡尔坐标系
    '''
    thera = (math.pi * lat) / 180
    fie = (math.pi * lng) / 180
    x = r * math.cos(thera) * math.cos(fie)
    y = r * math.cos(thera) * math.sin(fie)
    z