本文介绍了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;
}
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持萬仟网。
希望与广大网友互动??
点此进行留言吧!