求问已经(x,y,z)空间直角坐标,怎么得到大地坐标,就是B,H根据公式如何求...
发布网友
发布时间:2024-05-11 18:53
我来回答
共1个回答
热心网友
时间:2024-06-17 16:21
这是大地测量学的内容。经度L你会求吧?这个简单,直接根据反三角函数公式就可求出。关键是如何求纬度B和大地高H。可以采用迭代法。
先求初始值B0:B0=arctan[Z/sqrt(X^2+Y^2)](直接舍去了后面含N和H的接近于1的因子,因而可以得到近似值),其中sqrt代表开根号;
然后,利用B0,求初始值N0,H0:根据卯酉圈曲率半径计算公式N=a/W=a/sqrt(1-e^2*sinB*sinB)得,N0=a/sqrt(1-e^2*sinB0*sinB0),其中,a代表椭球长半轴,e代表离心率,对于一个给定的参考椭球,a与e是已知的。则H0=sqrt(X^2+Y^2)/cosB0-N0;
再把N0和H0的值代入完整的B的计算公式:得B1=arctan[Z/sqrt(X^2+Y^2)*(1-1/(e^2*N0/(N0+H0)))].再利用B1算出N1、H1……如此循环,直到B(n)与B(n-1)差值的绝对值满足限差要求即可。
以上过程,可以利用循环语句编程求解。(数学符号键盘输入不便,望谅解,也希望你看懂)