通过多个经纬度计算出该经纬度范围的中心点

业务上有一个需求，在地图上画出小区的面积范围，在范围的中心点画上名字。因此封装了一个方法，传入多组经纬度，根据多组经纬度算出多组经纬度组成的面积中心点。

该方法传入一个二维数组，数组的子项都是有经纬度组成的一个数组；实现方法如下。

/**
 * 根据坐标范围，得出中心点坐标
 * [[lon,lat],[lon,lat],[lon,lat]]
 */
export const getRangeCenterPoint = (points) => {
    let point_num = points.length; // 坐标点个数
    let X = 0, Y = 0, Z = 0;

    for (let i = 0; i < point_num; i++) {
        let point = points[i];

        if (!point || point.length < 2) {
            continue; // 如果点无效或格式不正确，跳过
        }

        let lon = point[0] * Math.PI / 180; // 经度转换为弧度
        let lat = point[1] * Math.PI / 180; // 纬度转换为弧度

        let x = Math.cos(lat) * Math.cos(lon);
        let y = Math.cos(lat) * Math.sin(lon);
        let z = Math.sin(lat);

        X += x;
        Y += y;
        Z += z;
    }

    // 计算平均 X, Y, Z
    X = X / point_num;
    Y = Y / point_num;
    Z = Z / point_num;

    // 通过 X, Y, Z 反算经纬度
    let tmp_lon = Math.atan2(Y, X);
    let tmp_lat = Math.atan2(Z, Math.sqrt(X * X + Y * Y));

    // 将弧度转换为角度并返回
    return [tmp_lon * 180 / Math.PI, tmp_lat * 180 / Math.PI];
}