摘 要：针对多镜头测绘级无人机价格昂贵、技术门槛较高的问题，本文提出利用消费级无人机获取的倾斜影像制作大比例尺地形图的方法。该方法运用消费级单镜头无人机，在ＲＴＫ （ ｒｅａｌ － ｔｉｍｅ ｋｉｎｅｍａｔｉｃ ）辅助的前提下，采用低空飞行、井字形航线设计、高重叠率的作业方案，配以高精度像控点、立体测图方法改进，从而实现大比例尺数字测图。以莆田市农村地籍与房屋调查工程为实践，对该方法的精度进行了分析，结果表明，该方法精度完全能够满足１∶５００数字地形图测绘的需要，且作业效率高。
关键词：消费级单镜头无人机；倾斜摄影测量；大比例尺测图；精度分析随着数字测图的人工成本不断攀升和小型旋翼无人机的快速发展，基于倾斜三维模型的立体测图方法受到各个测绘单位的欢迎。采用单镜头相机进行高精度的无人机倾斜摄影测量，将在小范围大比例尺地形图修补测、城市三维建模 、超高层建筑竣工测量、“多测合一”二、三维校核等领域发挥重要作用。大比例尺测图方面，黎其添提出了一种利用倾斜摄影测量技术快速测绘大比例尺地形图的方法，平面点位中误差为０．１５ｍ ；吴献文等利用精灵４Ｐｒｏ无人机获取的倾斜影像制作实景三维模型测绘大比例尺地形图，相对航高１００ｍ ，地面分辨率２．９９ｃｍ ，最终检查结果的平面中误差为±０.０６５４ｍ ，达到了１∶５００数字地形图精度，但无法满足地籍测量精度；段宗恩等采用多旋翼无人机搭载单镜头相机，对小范围、低成本的实景三维建模进行研究，构建贵阳市花溪区某试验区实景三维模型，通过单点精度检验，符合大比例尺航空摄影测量内业规范要求。本文基于莆田市农村地籍与房屋调查工程，探讨了使用消费级、单镜头、低成本无人机进行倾斜摄影测量制作大比例尺地形图的相关技术，并通过改进相关作业方法，提高数字测图工作效率，保证测图精度，降低无人机倾斜摄影测量的应用门槛。研究区位于福建省第三大岛南日岛，全岛面积５２ｋｍ２ ，海岸线总长６６．４ｋｍ 。岛上全年４级风以上时间达３００ｄ以上，为消费级无人机的航摄增加了较大困难；但岛内房屋错落有致，植被稀疏，为倾斜摄影测量提供了有利条件。
大疆精灵４ＲＴＫ 版无人机具有 ＲＴＫ 导航定位系统、相机微秒级同步、 ＡＰＰ航线规划等功能，具备ｃｍ 级导航定位系统和高性能成像系统，配备千寻位置公司集成的千寻知寸网络 ＲＴＫ 模块，有利于弱化姿态控制单元的重要性，减少荷载，提高定位精度。在 ＲＴＫ 的辅助下，采用井字形航线设计、低空飞行、高重叠率的作业方案获取影像数据。为了得到与多镜头云台一样的航测成果，单相机倾斜摄影需要更多的航带数和照片数 。航线规划涉及地面分辨率、航高、航速、相机倾角、重叠度等。对于单镜头无人机，常用的飞行方案为多架次“之”字形方案，即在航向规划软件中依架次设置不同的航线角度、相机倾角，来模拟多镜头的航飞效果。也可采用“井”字形航线，即固定相机倾角不变，采用“井”字形反复。飞行的多角度飞行方案，既保证同一地物能被３个以上的不同方向航拍到，又保持相机倾角的稳定。本文采用如图１所示的“井”字形航线。由于大疆精灵４ＲＴＫ的单块电池续航时间低于１５ｍｉｎ ，为提高外业采集效率，一般需配备１６块电池或６块电池加便携式充电站。在满负荷１６架次飞行的情况下，单天可完成１ｋｍ２ 的倾斜航摄。研究区面积为０．３７ｋｍ２ ，实际飞行８个架次，加上起降、换电池的时间损耗，在１８０ｍｉｎ内完成了数据采集工作，共获取照片２８３６张。
图 １ 航线设计图
考虑到现场和设备情况，采用低空飞行（设置默认航高为低空８０ｍ ），在没有高压电塔、风力发电塔或房屋较低矮、密集的区域，将航高降低到 ７０ｍ 。顾及地表高低起伏，相对航高一般５０多米，局部相对航高只有３０多米。航向重叠度８０％ ，旁向重叠度８０％ ，云台俯仰角度－６０°，地面分辨率优于２．１９ｃｍ 。相对航高低于文献［ ５ ］，而地面分辨率优于文献［５ ］。
像控点布设是空中三角测量的基础，其测量精度、布设密度和分布直接影响成果精度。本文按实地间隔１００ｍ 均匀布设像控点，像控点标志按“ Ｌ ”型布设，在Ｌ 型标志的外角处打入水泥钉，并以红漆绘出点号，便于后期检测精度，且兼作图根控制点。空中三角测量的精度最弱点位于测区边缘，因此要加强测区边角的控制，以改善区域四周的精度 ［８ ］ 。以研究区内的东岱村为例，东岱村房屋密集区面积０．３７ｋｍ２ ，三面临海。测区内共布设了２１个像控点，如图２所示。每个像控点独立观测３次，在 ＲＴＫ固定解收敛至 ｍｍ 级后开始观测，每次观测历元数３０个，采样间隔２ｓ 。观测成果平面坐标分量 较 差 均 小 于 ３ｃｍ ，垂 直 坐 标 分 量 较 差 小 于５ｃｍ ，取中数作为该点测量的最终成果。
图２ 像控点布置图
空中三角测量利用外业采集的影像、像控点成果、 ＰＯＳ信息，通过建模获得所有影像的外方位元素，再通过多视影像密集匹配，自动获取高密度三维点云信息 ［８ ］ 。空中三角测量步骤主要分为多视影像密集匹配及自由网空三、像控点判刺、绝对网空三解算等。多视影像密集匹配及自由网空三步骤基于同名特征点信息进行自动匹配解算，基本无需人工干预。个别情况下，空三结果可能会出现丢片、分层等异常情况，可能是由于气象条件恶劣导致像片抖动剧烈或过曝光、 ＲＴＫ 信号失锁、ＰＯＳ数据异常、植被和水系区域弱纹理等 ［９ ］ 。
本文空中三角测量采用 ＣｏｎｔｅｘｔＣａｐｔｕｒｅ软件进行处理。研究区域共判刺２１个像控点，密集匹配５７５ ９６５个连接点，连接点如图３所示。
图 ３ 空三加密网连接点分布图
为比较不同像控点密度对空三精度的影响，将２１个像控点均匀抽稀到８个，得到两种像控点密度下的空三控制点精度（表１ ），可以看出像控点减少，空三精度下降。
三维自动建模基于高密度点云构建不规则三角网，对不规则三角网进行自动纹理映射，生成基于真实影像纹理的三维实景模型。生成的模型包含了高密度点云、不规则三角网、真实影像纹理３种信息。生成的倾斜三维模型成果的精度直接决定了最终成图的精度 ［９ ］ 。因此，需要对模型的内符合精度进行检验。如图 ４ 所示，对同一个房角点在不同高度取点对比得到 Δ Ｘ ＝０．０２ｍ ， Δ Ｙ ＝０．０１ｍ 。类似的，在测区内随机选择３０栋混房的角点，分别在同一墙体角点的不同高度取两点，得到如图５所示的点位误差分布图，根据 （ Δ Ｓ ×Δ Ｓ ）／（ ｎ －１槡）计算得到内符合精度为±０．０３４ｍ ，其中， Δ Ｓ 表示测量误差。
图 ４ 同一房角点平面坐标差异分析
图５ 三维模型内符合精度图
基于倾斜三维模型进行的立体测图也称“裸眼测图”。采用人机交互的方式直接从倾斜三维模型上获取特征点坐标，更直观、更快捷。绘图者不需要长时间的专业培训，即可快速上手，因此受到了各测绘单位的欢迎。目前，常用的立体测图的软件有武汉天际航公司的 ＤＰ Ｍｏｄｅｌｅｒ软件系统、北京清华山维的ＥＰＳ软件系统等。为进一步节省倾斜摄影测量的成本，降低立体测图的门槛，提高立体测图效率，且充分考虑不同区域房屋构造不同的特性，开发了基于 ＡｕｔｏＣＡＤ 的立体测图插件 ［１０ ］ 。该插件利用 Ｗｉｎｄｏｗｓ ＡＰＩ 函数，捕获并锁定“ Ａｃｕｔｅ３Ｄ窗口选点”和“ＣＡＳＳ窗口绘图”两个窗口，通过两个窗口间的信息交互来进行立体测图。即实时将 Ａｃｕｔｅ３Ｄ 测量窗口中显示的点位坐标信息同步展绘到 ＣＡＳＳ中。利用该插件，考虑到房屋自身的直角、对称、镜像等特点，可以充分利用 ＡｕｔｏＣＡＤ的构造线、垂线、延长线等工具进行辅助绘图。针对模型上房角点变形的问题，分别在房角点以及房角点所在的墙面取点，并兼顾房屋平行于相邻的房屋、加固坎、围墙、马路等特点，绘制构造线进行印证，交会出最理想的房角点位置。开发的插件与ＣＡＳＳ原有的各工具按钮、快捷命令完全不冲突，使用方便。图６为自主研发的立体测图插件的界面。
图 ６ 自主研发的立体测图插件
为了验证消费级无人机获取倾斜影像制作大比例尺地形图的精度，将试验区域内全站仪实测检查点坐标与立体测图成果同名点平面坐标进行比较，差值见表２ 。由表２可知点位精度为０．０４２ｍ。
为了更加系统、全面地检验精度，在试验区选取了８个村的立体测图成果进行打点检验，如表３所示。由表 ３ 可知，采用本文方案及改进措施，能达到农村地籍与房屋调查工程的精度要求 ［１１ ］ 。
本文分析了基于大疆精灵４ＲＴＫ 版无人机的单镜头倾斜摄影测量方法，提出了低空飞行、井字形飞行、较高的重叠度、改进的立体测图方法等方案。以莆田市南日岛８个村为试验区域，结果证明，该方案可以满足地籍图测绘的精度要求。本文方法将测图的主要工作从原来的全站仪外业打点转变为室内立体测图，并且整套装备价格低、工作效率高、易于上手，大大缩短了测图人员的培养周期。但该方法超低空飞行受测区房屋高低影响较大，航片数量巨大，像控点密度较高，建模效率不高。下一步将研究如何提高建模效率、减少像控点数量以及全自动的三维重建工作等问题。