贾延青 温立文

(北京星天地信息科技有限公司 北京 102299)

随着遥感技术和计算机技术的不断更新、无人机技术的迅猛发展，无人机摄影测量技术也应运而生。无人机摄影测量系统相比较于传统的航空摄影测量系统，其独特的优势在于：（1）平台构建容易，维护方便，成本在大幅度降低；（2）无人机摄影测量系统体积小、重量轻、操作方便，飞机在起飞降落时不需要专门的起降机场，且可用多种方式起飞和降落，机动灵活[1]；（3）受天气条件和地面地貌状况影响较小，作业方式相对比较安全，人工在地面控制飞行，使无人机可进入危险地带开展作业，从而降低工作人员在作业过程中的风险。正是无人机摄影测量系统的这些优点，使其成为补充传统航天摄影测量的一种重要手段，逐步得到测绘、地理信息等领域的广泛关注，并被普遍应用在天气的监测和预报、国土资源环境的调查、城市管理、海事信息化的建设与动态管理、灾害预报、灾害监测与评估、国家海洋权益保障、农业监测、电力维修、水利勘察等各个领域[2]。

1 无人机影像获取流程和技术路线

该文采用南方测绘“天行”八旋翼无人机搭载Sony RX1RM2 数码相机，在北京市某镇进行了生产作业，获取无人机低空航摄遥感影像并完成1∶500 大比例尺地形图的成图。该次作业选用的数据处理软件是Smart 3D 三维实景建模软件和清华山维EPS 裸眼3D测图软件。无人机影像获取流程具体见图1。

图1 无人机影像获取流程

1.1 测区概况

该次的作业区位于北京市某镇；测区属于温暖温带半湿润半干旱季风气候，夏季高温多雨，冬季寒冷干燥，春、秋短促；全年无霜期180～200天；该次作业任务是为了该地区前期的旅游开发项目做一个前期的勘察工作，项目面积约为3 km2左右。

起飞前，需要对测区已有资料进行收集。该次作业测区的已有资料包括从相关部门收集到该地区已有3个控制点，可以作为布设像控点的起算点，该次作业控制点坐标系采用1980 年西安坐标系，1985 高程基准。提前与甲方做好沟通并提供测区范围线，并且告知地区的难度重点，以及最大的海拔高差等有关信息。

1.2 航线设计

航线设计确保无人机能够按照预定轨道进行飞行作业，并保证飞行过程中的各种参数满足相关轨迹[3]。航线设计是否合理关乎最终成果，要严格按照相关规范进行设计，并做好复查工作，确保飞行过程的安全和飞行数据的可靠，该节主要是按照相关要求，计算出航线设计的所需参数[4]。

1.3 外业像控点测量

为保证该次测量作业的精确，该次野外像控点的布设严格遵守野外像片控制点的布设原则，共布设了76 个野外像控点，该次坐标系采用1980 年西安坐标系，1985高程基准，中央子午线为114°。其中3个已知控制点作为起算点，使用GPS-RTK 结合温州CORS 账号进行平滑采集，计算结果平面、高程精度均满足规范要求。

2 内业数据处理

得到两个测区的影像数据后，采用Smart 3D 软件进行内业数据处理，该软件是一款集数据导入、空三加密、三维建模于一体的自动化航空影像处理软件，该软件自动化操作程度高，能够生成应用于EPS裸眼3D测图软件的OSGB 瓦片格式的数据[5]。通过该软件生成的三维模型精度高、效果好，可以满足大比例尺成图所需要的模型的要求，Smart 3D 软件处理数据的流程具体见图2。

图2 Smart 3D 的作业流程

2.1 航测原始数据的导入

首先将该次航测获取到的记录飞行姿态的POS数据、影像数据和像控点数据在数据导入界面导入到Smart 3D中，完成航测数据的录入工作。

2.2 空中三角测量

设置好像片的相关属性后，选择提交空中三角测量，可以设置空三的名称、定位、参考方式以及其他设置等。如果对一次空三加密的计算结果不满意，Smart3D可以进行多次提交空三加密，直到符合精度要求[6]。该次试验共进行了两次空三加密。

2.3 构建三维模型

空三加密完成之后，开始进行三维模型的构建，Smart 3D 的三维模型构建高度自动化，该次构建需要生成的是OSGB 格式瓦片模型，为了加快生成模型的速度和精度，需要对整个测区的空间框架进行构建，选取生成的模型范围。根据电脑性能和测区的实际情况，将整个测区共分为49 个瓦片。随后提交项目生产。Smart 3D可以生成三维网格、三维点云、正射影像和数字表面模型，根据不同生产项目的需要可以生产不同类型的数据。

2.4 数据成果

经过上述步骤，就完成了内业处理的整个步骤，再经过对后续成果整理和修饰，就可得到测区的数字三维模型。具体见图3。

图3 测区数字三维模型

3 精度分析

该文精度分析主要是利用GNSS-RTK系统对测区内明显的地物点实测其平面坐标和高程，然后在数字线划图上选取该点的图上坐标和高程，最后利用中误差计算公式计算其平面精度和高程精度。检查点主要选取一些比较明显、在图片上比较容易辨别的点，比如道路上斑马线或者交通标志的垂直拐点、测区内平房的房屋拐点、硬化水泥路十字路口的交点等。

3.1 平面精度分析

利用GNSS-RTK 系统对测区内20 个明显的地物点实测其平面坐标和高程，然后在数字线划图上选取该点的图上坐标和高程，最后利用中误差计算公式计算其平面精度和高程精度，可以得到测区检查点的平面精度统计结果，具体见表1。

表1 测区检查点的平面精度统计（单位：m）

从表1数据可以得出测区平面精度统计结果：X方向中误差为0.065 m，最大误差为-0.108 m，Y方向中误差为0.089 1 m，最大误差为0.154 m，平面坐标中误差为0.112 m，最大误差为0.198 m，均满足规范要求。

3.2 高程精度分析

同样可以算出测区的高程误差和精度，具体数据见表2。

表2 测区检查点的高程精度统计（单位：m）

从表2 数据可知，测区的高程中误差为0.274 m，最大误差为0.438 m，该测区地形起伏不大，大多属于平地和丘陵地，满足1∶500地形图对高程精度要求。

4 结语

系统阐述了无人机摄影测量系统在大比例成图中的完整流程，通过实例精度分析，可以得到结论：在控制好成图过程误差传播，严格按照规范要求获取影像数据，数据处理过程严格把关，无人机航测大比例尺地形图可达到1∶500 比例尺的成图精度，尤其是在平原地区平面精度和高程精度都很高，可以在该地区应用于生产实践。