java 地心坐标系(ECEF)和WGS-84坐标系(WGS84)互转的实现

2022-10-16,

本文介绍了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;
  }

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持。

《java 地心坐标系(ECEF)和WGS-84坐标系(WGS84)互转的实现.doc》

下载本文的Word格式文档,以方便收藏与打印。