python实现经纬度转笛卡尔坐标（czml文件制作）

在学习mars3d的时候，想用czml制作动画效果，但我只有时间点，经纬度（wgs84），海拔的表格数据，而编写czml文件需要使用笛卡尔坐标系。然后我参考Cesium中笛卡尔坐标系到底是什么鬼_coslatitude-CSDN博客对Cartesian3.fromRadians的解释用python写了经纬度转为笛卡尔坐标系的代码

Cartesian3.fromRadian如下：

Cartesian3.fromRadians = function(longitude, latitude, height, ellipsoid, result) {
    //>>includeStart('debug', pragmas.debug);
    Check.typeOf.number('longitude', longitude);
    Check.typeOf.number('latitude', latitude);
    //>>includeEnd('debug');

    height = defaultValue(height, 0.0);
    var radiiSquared = defined(ellipsoid) ? ellipsoid.radiiSquared : wgs84RadiiSquared;

    var cosLatitude = Math.cos(latitude);
    scratchN.x = cosLatitude * Math.cos(longitude);
    scratchN.y = cosLatitude * Math.sin(longitude);
    scratchN.z = Math.sin(latitude);
    scratchN = Cartesian3.normalize(scratchN, scratchN);

    Cartesian3.multiplyComponents(radiiSquared, scratchN, scratchK);
    var gamma = Math.sqrt(Cartesian3.dot(scratchN, scratchK));
    scratchK = Cartesian3.divideByScalar(scratchK, gamma, scratchK);
    scratchN = Cartesian3.multiplyByScalar(scratchN, height, scratchN);

    if (!defined(result)) {
        result = new Cartesian3();
    }
    return Cartesian3.add(scratchK, scratchN, result);
};

 python代码：

import math



def degrees_to_radians(degrees):
    return math.radians(degrees)


def cartesian3_from_radians(longitude, latitude, height=0.0):
    """
    根据经度、纬度和高度以及椭球体参数计算三维笛卡尔坐标。

    参数:
    longitude (float): 经度
    latitude (float): 纬度
    height (float): 高度，默认为0.0

    返回:
    list: 计算出的三维笛卡尔坐标
    """
    scratchN = [0, 0, 0]

    wgs84RadiiSquared = [6378137.0 ** 2, 6378137.0 ** 2, 6356752.3142451793 ** 2]  # WGS84椭球体半径平方
    ellipsoid = wgs84RadiiSquared

    latitude = degrees_to_radians(latitude)
    longitude = degrees_to_radians(longitude)

    cos_latitude = math.cos(latitude)
    scratchN[0] = cos_latitude * math.cos(longitude)
    scratchN[1] = cos_latitude * math.sin(longitude)
    scratchN[2] = math.sin(latitude)

    # 归一化向量
    norm = sum([n ** 2 for n in scratchN])
    scratchN = [n / norm for n in scratchN]

    # 计算乘积
    scratchK = [a * b for a, b in zip(ellipsoid, scratchN)]

    # 计算gamma
    gamma = math.sqrt(sum([a * b for a, b in zip(scratchN, scratchK)]))
    scratchK = [k / gamma for k in scratchK]

    # 高度乘以归一化向量
    scratchN = [n * height for n in scratchN]


    # 结果相加
    result = [k + n for k, n in zip(scratchK, scratchN)]

    return result

运行结果

和cesium的计算结果一样。
cesium将经纬度转为世界坐标（笛卡尔坐标）_cesium 经纬度转 cartesian3-CSDN博客