一、前言

参考文章：使用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, 249.952299308],[-2314347.181574534, 5391112.29519658, 2493702.9965011827]]//----------------------------------------------------定义坐标系椭球参数对象let wgs84 = {a: 6378137, //长半径,单位mb: 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转换为BLHvar mainDataXYZ = [[-2309213.896411623, 5396468.218629163, 2486913.5763168903],[-2310498.3804574762, 5395130.358054594, 2488611.242573863],[-2311782.0960432235, 5393791.747969109, 2490308.7006031894],[-2313065.032973063, 5392452.392053386, 249.952299308],[-2314347.181574534, 5391112.29519658, 2493702.9965011827]]//----------------------------------------------------定义坐标系椭球参数对象let wgs84 = {a: 6378137, //长半径,单位mb: 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转换为BLHfunction 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的初值B0let 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])}

得到数据：