前段时间由于项目需要,自己开发了通过查阅资料,编写了经纬高转化成XYZ坐标系的程序,并测试通过!有需要的可以看一下(转载请注明出处)
#define PI 3.14159265358979//*---------------------------------------------------------------函数: LoLaHtoDxSpaceXYZ(double L,double B,double Height,double &X,double &Y,double &Z)参数: L,B,Heith -- 地心坐标系下的经纬度坐标 单位:度X,Y,Z-- 地心空间坐标系下的坐标 单位:米功能: WGS84经纬度转换为地心空间直角坐标----------------------------------------------------------------*/{double a = 6378137.0; double e2 = 0.00669438002290;//转换为弧度 L = L * PI/180;B = B * PI/180;double fac1 = 1- e2*sin(B)*sin(B);double N = a/sqrt(fac1); //卯酉圈曲率半径double Daita_h = 0;//高程异常,默认为0double h = Daita_h + Height;X = (N+h)*cos(B)*cos(L);Y = (N+h)*cos(B)*sin(L);Z = ( N*(1-e2)+ h ) * sin(B);}