wgc84 笛卡尔_从经纬度到笛卡尔坐标的转换

料青山看我应如是

转换理论GPS(WGS84)到笛卡尔坐标 https：/en.wikipara.org/wiki/地理_坐标_转换#从_测地线_到_ECEF_坐标以下是我使用的内容：GPS(WGS 84)和笛卡尔坐标中的经度是相同的。所需纬度由WGS 84椭球参数转换为半长轴6378137米，平坦度倒数为298.257223563。我附上了一个VB代码我写道：Imports System.Math

'Input GPSLatitude is WGS84 Latitude,h is altitude above the WGS 84 ellipsoid

Public Function GetSphericalLatitude(ByVal GPSLatitude As Double, ByVal h As Double) As Double

        Dim A As Double = 6378137 'semi-major axis 

        Dim f As Double = 1 / 298.257223563  '1/f Reciprocal of flattening

        Dim e2 As Double = f * (2 - f)

        Dim Rc As Double = A / (Sqrt(1 - e2 * (Sin(GPSLatitude * PI / 180) ^ 2)))

        Dim p As Double = (Rc + h) * Cos(GPSLatitude * PI / 180)

        Dim z As Double = (Rc * (1 - e2) + h) * Sin(GPSLatitude * PI / 180)