java 地心座標系(ECEF)和WGS-84座標系(WGS84)互轉的實現

這篇文章主要介紹了java 地心座標系(ECEF)和WGS-84座標系(WGS84)互轉的實現,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友們下面隨着小編來一起學習學習吧

本文介紹了java 地心座標系(ECEF)和WGS-84座標系(WGS84)互轉的實現,分享給大家,具體如下:

 public static String WGS84toECEF(double latitude, double longitude, double height)
  {
    double X;
    double Y;
    double Z;
    double a = 6378137;
    double b = 6356752.314245;
    double E = (a * a - b * b) / (a * a);
    double COSLAT = Math.cos(latitude * Math.PI / 180);
    double SINLAT = Math.sin(latitude * Math.PI / 180);
    double COSLONG = Math.cos(longitude * Math.PI / 180);
    double SINLONG = Math.sin(longitude * Math.PI / 180);
    double N = a / (Math.sqrt(1 - E * SINLAT * SINLAT));
    double NH = N + height;
    X = NH * COSLAT * COSLONG;
    Y = NH * COSLAT * SINLONG;
    Z = (b * b * N / (a * a) + height) * SINLAT;
    return X + "," + Y + "," + Z;
  }
 
 
  public static String ECEFtoWGS84(double x, double y, double z)
 
  {
    double a, b, c, d;
    double Longitude;//經度
    double Latitude;//緯度
    double Altitude;//海拔高度
    double p, q;
    double N;
    a = 6378137.0;
    b = 6356752.31424518;
    c = Math.sqrt(((a * a) - (b * b)) / (a * a));
    d = Math.sqrt(((a * a) - (b * b)) / (b * b));
    p = Math.sqrt((x * x) + (y * y));
    q = Math.atan2((z * a), (p * b));
    Longitude = Math.atan2(y, x);
    Latitude = Math.atan2((z + (d * d) * b * Math.pow(Math.sin(q), 3)), (p - (c * c) * a * Math.pow(Math.cos(q), 3)));
    N = a / Math.sqrt(1 - ((c * c) * Math.pow(Math.sin(Latitude), 2)));
    Altitude = (p / Math.cos(Latitude)) - N;
    Longitude = Longitude * 180.0 / Math.PI;
    Latitude = Latitude * 180.0 / Math.PI;
    return Longitude + "," + Latitude + "," + Altitude;
  }

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持神馬文庫。

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章