关于超图模型坐标系转换参数精度评定方法探索

关于超图模型坐标系转换参数精度评定方法探索

以超图帮助文件例子为例：

超图判断精度有效坐标对方法：
计算得到重合点坐标残差（Dx、Dy、DH)及点位中误差，程序判定坐标残差大于3倍点位中误差的重合点为不可用点，即在“可用”一栏显示为“否”；反之满足精度要求。
为提高转换参数的精度，用户可将不可用的重合点剔除或设置为不采用，重新计算转换参数。

相关误差判断公式：
$重合点残差：$
$$V=重合点转换坐标值-重合点已知坐标值$$
$点位中误差计算公式：$
$$M_p=\sqrt{M_X^2+M_Y^2+M_Z^2}$$
$$M_X:空间直角坐标X残差中误差，M_X=\pm \sqrt{\frac{[VV{]}_X}{n-1}}$$
$$M_Y:空间直角坐标Y残差中误差，M_Y=\pm \sqrt{\frac{[VV{]}_Y}{n-1}}$$
$$M_Z:空间直角坐标残差中误差，M_Z=\pm \sqrt{\frac{[VV{]}_Z}{n-1}}$$
$平面点位中误差：$
$$M_p=\sqrt{M_X^2+M_Y^2}$$
$$M_X:空间直角坐标X残差中误差，M_X=\pm \sqrt{\frac{[VV{]}_X}{n-1}}$$
$$M_Y:空间直角坐标Y残差中误差，M_Y=\pm \sqrt{\frac{[VV{]}_Y}{n-1}}$$

以上述例子来说：
12组点对，第一组点对R551的中误差计算就是
$$M_p=\sqrt{M_X^2+M_Y^2+M_Z^2}$$
$$M_X=\pm \sqrt{\frac{[VV{]}_X}{n-1}}=\pm \sqrt{\frac{D_X\ast D_X}{12-1}}$$
$$M_Y=\pm \sqrt{\frac{[VV{]}_Y}{n-1}}=\pm \sqrt{\frac{D_Y\ast D_Y}{12-1}}$$
$$M_Z=\pm \sqrt{\frac{[VV{]}_Z}{n-1}}=\pm \sqrt{\frac{D_Z\ast D_Z}{12-1}}$$

这里换算理解：以单点残差 除以 总数量 来计算点位中误差

经过测试发现：针对规则程序判定坐标残差大于3倍点位中误差的重合点为不可用点，超图实现只计算X，Y残差分量与点位中误差的关系，不考虑Z分量残差

针对模型转换的理解：
布尔莎模型针对的是地心笛卡尔积空间坐标系XYZ，所以需要从投影坐标系转换为大地坐标系，再转换为地心笛卡尔积空间坐标系进行换算