//空间直角坐标系(xyz)转换大地坐标系(经纬度)
double a = 6378137.000;//长半轴
double f = 1 / 298.3;//扁率
double b = a - a * f;//短半轴
double ec = 1 - Math.pow(b, 2) / Math.pow(a, 2);//第一偏心率平方
double ecc= ec / (1-ec);//第二偏心率平方
double X =-2098813.2676;
double Y =4804967.1055;
double Z =3619537.1465;
double L = Math.atan(Y / X);
if (X < 0)
L += Math.PI;
double r = Math.sqrt(X * X + Y * Y);
double B1 = Math.atan(Z / r);
double B2;
while (true)
{
double W1 = Math.sqrt(1 - ec * (Math.sin(B1) * Math.sin(B1)));
double N1 = a / W1;
B2 = Math.atan((Z + N1 * ec * Math.sin(B1)) / r);
if (Math.abs(B2 - B1) <= 0.0000000001)
break;
B1 = B2;
}
double B = B2;
double W = Math.sqrt(1 - ec * (Math.sin(B2) * Math.sin(B2)));
double N = a / W;
double H = r / Math.cos(B2) - N;
System.out.println(L+","+B+","+H);
//大地坐标系(经纬度)转换空间直角坐标系(xyz)
double L=113.608437;
double B=34.802459;
double Height=0;
double a = 6378137.0;
double e2 = 0.00669438002290; //转换为弧度
L = L * PI/180;//经度
B = B * PI/180;//纬度
double fac1 = 1- e2*Math.sin(B)*Math.sin(B);
double N = a/Math.sqrt(fac1); //卯酉圈曲率半径
double Daita_h = 0; //高程异常,默认为0
double h = Daita_h + Height;
double X = (N+h)*Math.cos(B)*Math.cos(L);
double Y = (N+h)*Math.cos(B)*Math.sin(L);
double Z = ( N*(1-e2)+ h ) * Math.sin(B);
System.out.println(X+","+Y+","+Z);
该方法转换的经纬度有问题,因为我需要的是高度差,所以我只关心了高的标准值。