一种基于无人机的实时塌方区域面积的测量方法与流程

本发明涉及测量领域，尤其涉及一种基于无人机的实时塌方区域面积的测量方法。

背景技术：

应急测绘是为各类突发公共事件提供地理信息和现代测绘技术的基础支撑，是国家突发事件应急体系的重要组成部分，是指挥决策和抢险救灾的保障和依据，其中及时获取灾害发生地的实时测绘信息对于开展突发事件应急救援具有重要意义，目前常用的方法是使用无人机到地质灾害现场拍摄正射/三维影像，待飞机飞回指挥中心后，利用高性能计算机对现场进行正射成像/三维重建，在正射影像/三维模型上进行灾害区域的面积测量，但是存在的问题是需要飞机飞回指挥中心后，再提取相关资料进行计算，实时性差，往往实现某个区域的正射影像图需要几个小时甚至几天，大大影响灾害评估的时效性，从而影响救援速度。

技术实现要素：

本发明的目的就在于为了解决上述问题设计了一种基于无人机的实时塌方区域面积的测量方法。

本发明通过以下技术方案来实现上述目的：

一种基于无人机的实时塌方区域面积的测量方法，包括以下步骤：

s1、无人机搭载光电吊舱采集数据，并通过地面站进行图像数据和光电吊舱参数的保存；

s2、读取地面站保存的图像数据，并解析出图像数据中保存的相机内参数和采集数据时的外参数；

s3、将读取的图像数据进行图像坐标系转换为世界坐标系；

s4、根据转换后的世界坐标系参数得出待测塌方区域面积的相关情况。

本发明的有益效果在于：无人机在地质灾害发生的第一时间到达现场，利用高清图传实时回传现场高清图片和飞机姿态、相机姿态信息，实时测量塌方区域面积，实现灾害现场实时测绘的方法，提高灾害评估的时效性，提高救援速度。

附图说明

图1是本发明一种基于无人机的实时塌方区域面积的测量方法的流程图；

图2是本发明一种基于无人机的实时塌方区域面积的测量方法的坐标示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。

因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。

在本发明的描述中，需要理解的是，术语“上”、“下”、“内”、“外”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系，或者是该发明产品使用时惯常摆放的方位或位置关系，或者是本领域技术人员惯常理解的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

此外，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

在本发明的描述中，还需要说明的是，除非另有明确的规定和限定，“设置”、“连接”等术语应做广义理解，例如，“连接”可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接连接，也可以通过中间媒介间接连接，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。

下面结合附图，对本发明的具体实施方式进行详细说明。

一种基于无人机的实时塌方区域面积的测量方法，包括以下步骤：

s1、无人机搭载光电吊舱采集数据，并通过地面站进行图像数据和光电吊舱参数的保存；

s2、读取地面站保存的图像数据，并解析出图像数据中保存的相机内参数和采集数据时的外参数，内参数包括摄像机的焦距f、成像平面上的单位距离的像素sx和sy、成像分辨率(w,h)；外参数包括无人机飞行的滚转角r、方位角a、俯仰角p、高度(height)；

s3、将读取的图像数据进行图像坐标系转换为世界坐标系；

s31、根据采集到的图像数据选取建立图像坐标系为uov，并根据成像分辨率的(w,h)的求得图像的中心点坐标为(u0，v0)，以摄像机的光心点为坐标原点建立摄像机坐标系为oc-xcyczc，摄像机的光轴作为yc轴，xc和zc分别与成像的u和v轴平行，世界坐标系为ow-xwywzw；

s32、假定摄像机的姿态向量b＝(a,p,r)^t，位置向量pc＝(xwc,ywc,zwc)^t,点ui＝(ui,1,vi)^t,根据pi＝ycrt^-1(ui-u0)+pc求得ui所在的实际坐标点pi＝(xwi,ywi,zwi)^t，其中u0＝(u0,0,v0)^t，t为摄像机坐标系到成像平面的投影矩阵，r为摄像机坐标系与世界坐标系转换的旋转矩阵，zwc为无人机飞行的高度，ui在光轴yc上的投影为其中rij为r矩阵中第i行第j列元素，(ui,vi)为图像上待求点坐标；

s33、根据实际坐标点通过世界坐标系中平面点坐标，求得pi和pj两点间距离，

s4、根据转换后的世界坐标系参数得出待测塌方区域面积的相关情况，通过多边形面积计算公式，结合在世界坐标系中每个点的坐标信息，最终得出待测塌方区域的面积m：

本发明的技术方案不限于上述具体实施例的限制，凡是根据本发明的技术方案做出的技术变形，均落入本发明的保护范围之内。