2000与WGS84大地坐标系转换为空间直角坐标系

 

图1中，P点为测量点，P0为其在地球椭球上的投影点。由图1测量点P点大地坐标P(B,L,Hn) 与空间直角坐标P(X,Y,Z)的几何关系，通过简单矢量运算和三角函数运算即可得出大地坐标(B,L,H)转换到空间直角坐标(X,Y,Z)的公式。

相同的基准，大地坐标系转换为空间直角坐标系，其公式为

X=(N+H)·cosBcosL

Y=(N+H)·cosBsinL

Z=[N(1-e2)+H]·sinB

其中，N为卯酉圈曲率半径；e为椭球第一偏心率；e2为 “e的平方”。

注：卯酉圈上一点的曲率半径(N)为：N=a(1-e^2·sin^2 B)^(-1/2)，

       e2=(a^2-b^2)/2

       式中：a—椭球长半径；e—椭球第一偏心率；B—大地纬度。

我国主要采用的，2000国家大地坐标系与WGS-84坐标系的地球椭球参数如表1所示。

JS代码如下

// WGS84坐标系转空间直角坐标
function TransRect(lon,lat){
	 var mm = 6378137;
     var b = 6356752.3142;
     var f = 1 / 298.257223563;
     var b = mm * (1 - f);
     var eeee = Math.sqrt(mm * mm - b * b);
     var e = eeee / mm;
     var N = mm / Math.sqrt(1 - e * e * Math.sin(lat * Math.PI / 180) * Math.sin(lat * Math.PI / 180));
     var WGS84_X = (N) * Math.cos(lat * Math.PI / 180) * Math.cos(lon * Math.PI / 180);
     var WGS84_Y = (N) * Math.cos(lat * Math.PI / 180) * Math.sin(lon * Math.PI / 180);
     return [WGS84_X, WGS84_Y]
}

参考文献：

[1]仝巧珍,赵计环.谈大地坐标系和空间直角坐标系转换的研究[J].山西建筑,2011,37(08):194-195.