7参数转换整体流程(不同椭球之间):
平面直角坐标系-->大地坐标系-->空间直角坐标系 转换(To) 空间直角坐标系-->大地坐标系-->平面直角坐标系
举个栗子,比如从BJ1954平面直角坐标系 转到XA1980平面直角坐标系那么需要5步:
①BJ1954平面直角坐标系 至 BJ1954大地坐标系
②BJ1954大地坐标系 至BJ1954空间直角坐标系
③BJ1954空间直角坐标系 至XA1980空间直角坐标系
④XA1980空间直角坐标系 至XA1980大地坐标系
⑤XA1980大地坐标系 至XA1980平面直角坐标系
如果从BJ1954空间直角坐标系转到XA1980平面直角坐标系只需一步:
①BJ1954空间直角坐标系 至 XA1980空间直角坐标系
1)空间直角坐标系(XYZ) 转换(To) 空间直角坐标系(XYZ):(布尔沙模型,此步重要,将一个椭球基准转换到另一个椭球基准)
其中△X,△Y,△Z是坐标平移量,R(ω)是旋转矩阵,(1+m)是比例因子
2)大地坐标系(BLH) 转换(To) 空间直角坐标系(XYZ):(相同基准)
3)空间直角坐标系(XYZ) 转换(To) 大地坐标系(BLH):(相同基准)
e2为椭球第一偏心率(a2-b2)/a2=ParaE1,e'2为椭球第二偏心率(a2-b2)/b2=ParaE2
N=double v = a / (Math.Pow(1 - Math.Pow(ParaE1, 1) * Math.Pow(Math.Sin(B), 2), 0.5));
通用的XYZ To BLH三个等价公式:(计算结果几乎无差)
4)大地坐标系(BLH) 转换(To) 平面直角坐标系(投影坐标系,xyh):(高斯投影)
计算公式见孔祥元等主编武汉大学出版社2002年出版的《控制测量学》的第72页书的的括号有问题, ( 和 [ 应该交换
5)求解七参数:(最小二乘法)
参见:
http://blog.sina.com.cn/s/blog_4c8bb86b0100bxq6.html
http://www.360doc.com/content/12/1130/11/8463843_251136066.shtml
大致流程为:
a)
b)如上5-7公式符合最小二乘法A*X=B的形式,依照最小二乘法求解:
四参数转换整体流程:(同一椭球下不同坐标系)
大地坐标系-->平面直角坐标系-->平面直角坐标系-->大地坐标系(注意与七参数的区别)
整体流程: