經(jīng)緯度和球體三維坐標(biāo)換算
在球面坐標(biāo)系上，坐標(biāo)是用經(jīng)緯度標(biāo)示，緯線和經(jīng)線一樣是為度量方便而假設(shè)出來的輔助線，定義為地球表面某點(diǎn)隨地球自轉(zhuǎn)所形成的軌跡。
任何一根緯線都是圓形而且兩兩平行。緯線的長(zhǎng)度是赤道的周長(zhǎng)乘以緯線的緯度的余弦，所以赤道最長(zhǎng)，離赤道越遠(yuǎn)的緯線，周長(zhǎng)越短，到了兩極就縮為0。從赤道向北和向南，各分90°，稱為北緯和南緯，分別用“N”和“S”表示。
在實(shí)際的計(jì)算中，經(jīng)緯度一般是不能拿來直接計(jì)算的，都需要進(jìn)行一定的轉(zhuǎn)換，比方說計(jì)算距離等，假設(shè)把地球是一個(gè)正圓的，放在一個(gè)很大三維坐標(biāo)系中，原心為地球的球心，坐標(biāo)換算成三維的（笛卡爾）坐標(biāo)數(shù)字，就方便計(jì)算了。
C#算法代碼
下邊分享一下C#版的，經(jīng)緯度和立體三維坐標(biāo)的相互轉(zhuǎn)換，都是一個(gè)數(shù)學(xué)換算，簡(jiǎn)單的進(jìn)行注釋：
//設(shè)置地球半徑
?private? double WGS_84_RADIUS_EQUATOR = 6378137.0;?
//經(jīng)緯度轉(zhuǎn)立體坐標(biāo)
? ? ?public CVector WGS84ToWorld(CVector _wgscood)
? ? ? ? {
? ? ? ? ? ? if (_wgscood.VZ == 0 || Double.IsNaN(_wgscood.VZ))
? ? ? ? ? ? {
? ? ? ? ? ? ? ? _wgscood.VZ = 1;
? ? ? ? ? ? }
? ? ? ? ? ? _wgscood.VZ += WGS_84_RADIUS_EQUATOR;
? ? ? ? ? ? double tempx = 0;
? ? ? ? ? ? double tempy = 0;
? ? ? ? ? ? double tempz = 0;
? ? ? ? ? ? tempy = _wgscood.VZ * Math.Sin(_wgscood.VY / 180 * Math.PI);
? ? ? ? ? ? tempx = _wgscood.VZ * Math.Cos(_wgscood.VY / 180 * Math.PI) * Math.Sin(_wgscood.VX / 180 * Math.PI);
? ? ? ? ? ? tempz = _wgscood.VZ * Math.Cos(_wgscood.VY / 180 * Math.PI) * Math.Cos(_wgscood.VX / 180 * Math.PI);
? ? ? ? ? ? CVector rexyz = new CVector { VX = tempx, VY = tempy, VZ = tempz };
? ? ? ? ? ? return rexyz;
? ? ? ? }
? ? ?//立體坐標(biāo)轉(zhuǎn)回經(jīng)緯度
? ? ? ? public CVector WorldToWGS84(CVector _worldcood)
? ? ? ? {
? ? ? ? ? ? double tempx = 0;
? ? ? ? ? ? double tempy = 0;
? ? ? ? ? ? double tempz = 0;
? ? ? ? ? ? if (_worldcood.VX >= 0)
? ? ? ? ? ? {
? ? ? ? ? ? ? ? tempx = Math.Acos(_worldcood.VZ / Math.Sqrt(_worldcood.VX * _worldcood.VX + _worldcood.VZ * _worldcood.VZ)) * 180 / Math.PI;
? ? ? ? ? ? ? ? tempx = tempx < 0 ? tempx + 180 : tempx;
? ? ? ? ? ? }
? ? ? ? ? ? else
? ? ? ? ? ? {
? ? ? ? ? ? ? ? tempx = Math.Acos(_worldcood.VZ / Math.Sqrt(_worldcood.VX * _worldcood.VX + _worldcood.VZ * _worldcood.VZ)) * 180 / Math.PI;
? ? ? ? ? ? ? ? tempx = tempx < 0 ? -(tempx + 180) : -tempx;
? ? ? ? ? ? }
? ? ? ? ? ? tempy = Math.Asin(_worldcood.VY / Math.Sqrt(_worldcood.VX * _worldcood.VX + _worldcood.VY * _worldcood.VY + _worldcood.VZ * _worldcood.VZ)) * 180 / Math.PI;
? ? ? ? ? ? tempy = tempy > 90 ? 90 - tempy : tempy;
? ? ? ? ? ? tempz = (Math.Sqrt(_worldcood.VX * _worldcood.VX + _worldcood.VY * _worldcood.VY + _worldcood.VZ * _worldcood.VZ)) - WGS_84_RADIUS_EQUATOR;
? ? ? ? ? ? CVector relonlat = new CVector { VX = tempx, VY = tempy, VZ = tempz };
? ? ? ? ? ? return relonlat;
? ? ? ? }
https://mp.weixin.qq.com/s?__biz=MzU2ODYzNzc4OQ==&mid=2247485621&idx=1&sn=579bfb0f484e66894343197d4cac8e59&chksm=fc8ba887cbfc219186479a71c82c6a6409c7df2a423252ac090ca24c841e1fb147d4ec9cd02f&token=1732935870&lang=zh_CN#rd
