现代测量数据处理理论 - 粗差检测Python代码

Python版本：Python3.13

实验要求

今有一模拟的立体像对（f为100.5mm，摄影比例尺为1:25000，像点坐标量测中误差为±2.8μm，数据文件见附件），在y坐标中人为加入了50个4σ_0~100σ_0的粗差，请用本课程所学知识将这些粗差逐一找出

数据介绍

• （1）地面控制点坐标文件 PRECKI.DAT：包含6个控制点，坐标系为大地测量坐标系下的三维坐标，为右手系。
• （2）相机参数文件 PRECAI.DAT：包含成像时的摄影机主距、像主点坐标 等内方位元素，像主点坐标已归零化；
• （3）像点量测坐标文件 PREPHIy10.DAT：由测量点点号、横坐标、纵坐标组成。像点量测坐标文件以“-99”为标记界限，划分出205021、205022、205023三幅影像

核心思想

采用基于验后方差估计原理的选权迭代法结合立体像对相对定向理论实现粗差检测，其基本思想是在最小二乘的基础上，对定权的方法加以改进，通过自适应调整观测值权值，削弱粗差对参数估计的影响，最终识别出超出允许范围的粗差

更新时间

最后更新时间：2025年12月1日