JavaScript 实现空间直角坐标系转换为大地坐标系（XYZ→BLH）

最新推荐文章于 2025-05-01 16:32:51 发布

ReBeX

最新推荐文章于 2025-05-01 16:32:51 发布

阅读量2.6k

点赞数 2

CC 4.0 BY-SA版权

文章标签： javascript

本文链接：https://blog.csdn.net/ReBeX/article/details/120857104

本文介绍如何使用JavaScript在大地测量学中，基于WGS84坐标系，从大地空间直角坐标（XYZ）转换为大地坐标（BLH）。提供了两种公式比较和迭代计算BLH坐标的详细步骤和示例代码。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一、前言

参考文章：使用Matlab对大地测量学中的大地坐标和大地空间直角坐标进行相互转换

相关文章：JavaScript 实现大地坐标系转换为空间直角坐标系（BLH→XYZ）

在相同的坐标基准下，将大地空间直角坐标系转换为大地坐标系的公式为：
在这里插入图片描述

但是不知道为什么，我用这条公式来迭代，怎么都做不对，于是我就去找课件，找到了专业课老师提供的公式：

在这里插入图片描述

二、实现代码


var mainDataXYZ = [
    [-2309213.896411623, 5396468.218629163, 2486913.5763168903],
    [-2310498.3804574762, 5395130.358054594, 2488611.242573863],
    [-2311782.0960432235, 5393791.747969109, 2490308.7006031894],
    [-2313065.032973063, 5392452.392053386, 2492005.952299308],
    [-2314347.181574534, 5391112.29519658, 2493702.9965011827]
]
//----------------------------------------------------定义坐标系椭球参数对象
let wgs84 = {
    a: 6378137,                     //长半径,单位m
    b: 6356752.3142451795,          //短半径
    f: 1/298.257223563,             //扁率
}
var square = num => num*num         //快速计算平方
let d2r = Math.PI / 180;
let r2d = 180 / Math.PI;            //角度 = 弧度 * 180 / Math.PI;

//----------------------------------------------------XYZ转换为BLH

var mainDataXYZ = [
    [-2309213.896411623, 5396468.218629163, 2486913.5763168903],
    [-2310498.3804574762, 5395130.358054594, 2488611.242573863],
    [-2311782.0960432235, 5393791.747969109, 2490308.7006031894],
    [-2313065.032973063, 5392452.392053386, 2492005.952299308],
    [-2314347.181574534, 5391112.29519658, 2493702.9965011827]
]
//----------------------------------------------------定义坐标系椭球参数对象
let wgs84 = {
    a: 6378137,                     //长半径,单位m
    b: 6356752.3142451795,          //短半径
    f: 1/298.257223563,             //扁率
}
var square = num => num*num         //快速计算平方
let d2r = Math.PI / 180;
let r2d = 180 / Math.PI;            //角度 = 弧度 * 180 / Math.PI;

//----------------------------------------------------XYZ转换为BLH
function XYZtoBLH (RE,data){
    let X = data.X;
    let Y = data.Y;
    let Z = data.Z;
    let x2y2 = square(X)+square(Y)
    let e2 = (square(RE.a)-square(RE.b))/(square(RE.a));
    //使用迭代算法来计算B的值，先求出B的初值B0

    let B0 = Math.atan(Z/Math.sqrt(x2y2)+1);
    let B1 = Math.atan(Z/Math.sqrt(x2y2));
    while (Math.abs(B1 - B0)>= 0.000000000001){
        B0 = B1;
        var N = RE.a/Math.sqrt(1-e2*square(Math.sin(B0)));
        var H = (Z/Math.sin(B0))-N*(1-e2);
        B1 = Math.atan((Z+N*e2*Math.sin(B0))/Math.sqrt(x2y2));
    }
    let L = Math.atan2(Y,X);
    let B = B1;
    
    let result = {L:L*r2d,B:B*r2d,H:H};
    return result;
}
// 读取处理数据并输出，输出结果为角度制

var outputResult2 =[];
for (var i=0;i<mainDataXYZ.length;i++){
    let XYZ = {
        X: mainDataXYZ[i][0],
        Y: mainDataXYZ[i][1],
        Z: mainDataXYZ[i][2],
    };
    let BLH = XYZtoBLH(wgs84,XYZ);
    let opBLH = [BLH.B,BLH.L,BLH.H];
    outputResult2.push(opBLH);
}
for (let i=0;i<outputResult2.length;i++){
    console.log(outputResult2[i])
}

得到数据：
在这里插入图片描述