具体实施方式

为使本发明实施方式的目的、技术方案和优点更加清楚，下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整地描述。

本发明提供了一种基于高精度POS构架航线的辅助空三加密联合平差方法，具体包括如下步骤：

外业操作：进行多次分辨率精度递增的航摄操作；外业操作航摄航线设计具体为：根据地形均匀设置基线，以垂直基线的方向进行构架航线布设设计，其中，平均跨度不大于六条基线就补充一条构架航线。具体包括如下子步骤：

采用有人机在第一绝对航高获取测区第一轮数码影像；

采用无人机在第二绝对航高获取测区第二轮数码影像。

其中，第一绝对航高的高度大于第三绝对航高。第一轮数码影像的地面分辨率优于5cm，所述第二轮数码影像的地面分辨率优于3cm。

其中，有人机搭载焦距为79.8mm的UCD航摄仪，所述无人机采用搭载焦距为35mm的神思P1航摄仪的大疆经纬M300或搭载焦距为8.8mm的FC6310R航摄仪的大疆精灵4RTK。

内业操作：将经过预处理后的POS数据依序进行影像点匹配、检查、自检校和再平差后输出。具体包括如下子步骤：

1)对POS数据进行预处理；具体为：对POS数据的平面坐标与高程值进行转换和精化处理，其中，所述平面坐标采用高斯-克吕格直角平面投影，3度分带，中央经线为117度；所述高程值采用测区起算控制点所求的七参数进行转换，或采用似大地水准精化面改算。

2)进行数据导入与参数配置；

3)进行影像同名点的自动匹配；具体为：获取第一轮数码影像和第二轮数码影像，对第一轮数码影像和第二轮数码影像中的影像同名点进行匹配连接。

4)经过人工检查并剔除飞点与粗差点；

5)对航摄相机进行自检校得出相对检校前准确的相机参数；

6)将相机参数更新导入，再平差优化所有像片的对齐方式：

7)重复执行步骤4)-6)直至满足预设条件后输出平差结果。其中，预设条件具体为：平差精度满足1/2像素的要求。

实施例：

测试测区位于福建省沿海地区，地形平坦，以城镇居民区和新建新区为主，面积约40km²。需测绘1：500地形图。测区交通便利，植被较多。采用有人机搭载UCD航摄仪获取测区优于5cm分辨率的数码影像，虽然有带POS数据，但坐标定位精度较差，不能满足空三加密。因此采用大疆精灵4RTK与大疆经纬M300分别进行补充构架航线，航摄情况如表1：

表1航摄情况表

带高精度POS的构架航线是为减少像控点的布设，因此构架航线设计一般垂直于基本航向，也可根据测区形状灵活设计，也可参照像控布设的规范要求来设计，构架航线航摄时一定要采用带GPS-RTK功能的无人机平台搭载航摄仪进行，确保航摄过程中所获取的每一张相片都能参与空三加密平差使用。目的就是完全替代像片控制点，并且通过增加构架航线密度和高精度POS辅助空三加密，对于平差结果精度和区域网之间的连续性具有较好的效果。

本测区根据像控布设的规范要求来设计构架航线，平均跨度6条基线就补充一条构架航线，基本满区域网布控的原则。构架航线示意图如图1所示。

内业数据处理首先对POS数据的平面坐标与高程值进行转换和精化处理。平面坐标直接采用高斯－克吕格直角平面投影，3度分带，中央经线为117度，高程值采用测区起算控制点所求的七参数进行转换，这部分的POS数据处理工作至关重要。具体包括如下步骤，如图2所示：1)对POS数据进行预处理；2)进行数据导入与参数配置；3)进行影像同名点的自动匹配；4)经过人工检查并剔除飞点与粗差点；5)对航摄相机进行自检校得出相对检校前准确的相机参数；6)将相机参数更新导入，再平差优化所有像片的对齐方式：7)重复执行步骤4)-6)直至满足预设条件后输出平差结果。本测区空三加密平差采用photoscan软件进行，此软件影像同名点匹配精度高，匹配点连接效果好，准确度高，能满足此方法作业的需求。空三作业中每步聚都非常关键，每个参数配置都不能有误，不同时期的影像存在色差、纹理、地物要素等不同，匹配误差相应的增大，因此更需要高精度的POS数据作为辅助，才能提高匹配点的可靠性和提升匹配的速度。同时需要对POS数据的权重进行设置，每张构架航线的像主点坐标值的权重参数如图3所示。影像匹配点首先需要经过人工检查并剔除飞点与粗差点，同时对航摄相机进行自检校得出相对准确的相机参数，再将相机参数更新导入，再平差优化所有像片对齐方式，反复前述工作流程直到平差精度满足1/2像素的要求。构架航线和基本航线相机参数更新配置如图4所示。

本测区中，构架航线和基本航线间航高不一致，航摄时间不同，影像色差大，但是通过高精度的POS数据辅助空三匹配后，影像匹配连接点较多，匹配精度高，满足空三加密精度指标要求。构架航线和基本航线匹配连接点效果图如图5所示，构架航线和基本航线影像匹配效果图如图6所示。

精度结果分析：本测区空三平差过程需要对相机参数进行自检校后重新配置相机的参数。因此我们首先需查看相机参数畸变误差。本测区的相机畸变误差分析图如图7所示。

根据本测区的实际情况，在构架航线之间均匀布设了23处检查区，每个检查区采用外业碎部测量的方式采集特征点坐标，通过平面坐标转换和高程精化导入空三做为检查点。检查点分布图如图8：航测法成图内业空三加密点平差精度要求如表2所示：

表2内业加密点精度要求

本测区空三加密经过多次优化平差，最终检查点精度值见图9，可看出检查点精度满足表2要求。空三成果输出后，进行立体测图，外业巡图，编辑建库形成本项目最终DLG成果。本测区DLG成果通过野外精度验证，分别统计出平面绝对精度中误差±11.42cm，平面相对精度中误差±8.07cm，高程点精度中误差±11.66cm。统计表如图10所示。

以上所述为本发明的实施例，对于本领域的普通技术人员而言，根据本发明的教导，在不脱离本发明的原理和精神的情况下凡依本发明申请专利范围所做的均等变化、修改、替换和变型，皆应属本发明的涵盖范围。