一、前言

參考文章：使用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, //長半徑,單位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, 2492005.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轉換為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的初值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]) }

得到數據：

總結

以上是生活随笔為你收集整理的JavaScript 实现空间直角坐标系转换为大地坐标系（XYZ→BLH）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。