CN114397900B - 一种无人机航拍图片中心点经纬度误差优化方法 - Google Patents
一种无人机航拍图片中心点经纬度误差优化方法 Download PDFInfo
- Publication number
- CN114397900B CN114397900B CN202111434436.8A CN202111434436A CN114397900B CN 114397900 B CN114397900 B CN 114397900B CN 202111434436 A CN202111434436 A CN 202111434436A CN 114397900 B CN114397900 B CN 114397900B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- center point
- latitude
- longitude
- offset
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000005457 optimization Methods 0.000 title claims abstract description 23
- 230000009466 transformation Effects 0.000 claims abstract description 28
- 238000012821 model calculation Methods 0.000 claims 2
- 238000007689 inspection Methods 0.000 abstract description 16
- 238000004364 calculation method Methods 0.000 abstract description 10
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000010276 construction Methods 0.000 abstract description 3
- 238000003745 diagnosis Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0833—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using limited authority control
Landscapes
- Engineering & Computer Science (AREA)
- Computer Security & Cryptography (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种无人机航拍图片中心点经纬度误差优化方法,其包括建立飞行器姿态模型,计算出通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再构建机体坐标系和大地坐标系之间的关系;建立相机姿态模型,计算出相机垂直地面拍照中心点通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再构建机体坐标系和大地坐标系之间的关系;CGCS2000坐标系动态获取地球半径,通过CGCS2000坐标系获取当前地球半径;以及,优化图片中心点经纬度误差,通过飞行器姿态模型计算和相机姿态模型计算获取无人机航拍图片中心点经纬度;本发明巡检效率高、智能无人机机动性高、全自动飞行诊断、提升了电站巡检频次和解决了电站建设类型不同。
Description
技术领域
本发明涉及无人机航拍定点技术领域,特别是一种无人机航拍图片中心点经纬度误差优化方法。
背景技术
目前我国光伏产业规模持续扩大,行业发展总体趋好。伴随着我国光伏行业的进一步发展,预计未来几年,光伏行业市场容量将呈现出逐年增长态势。据预测,到2022年我国光伏累计装机容量将达141GW。随着并网发电的光伏电站规模越来越大,后期运营和维护业务将成为未来电站板块中增长最快的业务。
目前在光伏场站使用的无人机巡检热斑检测技术,通过红外相机拍摄的图片进行热斑检测获得热斑在图片中的位置,再根据无人机提供的经纬度信息作为无人机云台相机拍摄图片的中心点经纬度,最后获取图片中热斑实际对应的经纬度信息。以大疆M300 RTK飞行器搭载Zenmuse H20T云台为例构成一个巡检无人机,无人机高空飞行时需要不断地进行修正并产生不同程度的抖动。图1所示,O1点为飞行器4号机臂RTK天线的经纬度信息和高度信息,O2点为飞行器四号机臂在抖动时垂直到地面对应的真实位置经纬度,O3点为云台接口中心位置作为相机垂直向下拍照图片中心点位置经纬度,O4点是云台相机在抖动时拍摄的图片中心点对应的经纬度位置坐标。由设备自身误差和无人机高空定点飞行产生的抖动,对航拍图片中心点的经纬度产生一定误差。为提高最终光伏板热斑定位精度,因此无人机任意航向正确获取图片中心点经纬度位置信息O4点至关重要。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于现有技术中存在的问题,提出了本发明。
因此,本发明所要解决的技术问题是现有技术存在无人机设备本身和空中抖动,导致无人机航拍图片中心点的经纬度存在误差,造成采集结果不准确。
为解决上述技术问题,本发明提供如下技术方案:一种无人机航拍图片中心点经纬度误差优化方法,其包括建立飞行器姿态模型,计算出通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再构建机体坐标系和大地坐标系之间的关系;建立相机姿态模型,计算出相机垂直地面拍照中心点通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再构建机体坐标系和大地坐标系之间的关系;CGCS2000坐标系动态获取地球半径,通过CGCS2000坐标系获取当前地球半径;以及,优化图片中心点经纬度误差,通过飞行器姿态模型计算和相机姿态模型计算获取无人机航拍图片中心点经纬度。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:飞行器姿态模型建立:通过飞行器姿态角和RTK天线当前经纬度,结合RTK天线到云台接口中心的补偿量构建四元数,计算出通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再通过构建机体坐标系和大地坐标系之间的关系,计算出通过偏航角变换后的补偿量,及RTK天线到云台中心在大地坐标系上正北方向的偏移量和正东方向的偏移量。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:相机姿态模型建立:通过相机姿态角,结合RTK天线到地面的高度构建四元数,计算出相机垂直地面拍照中心点通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再通过构建机体坐标系和大地坐标系之间的关系,计算出通过偏航角变换后的补偿量,及相机拍照中心点位置变换后在大地坐标系上正北方向的偏移量和正东方向的偏移量。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:CGCS2000坐标系动态获取地球半径:通过CGCS2000坐标系,获取长半轴,短半轴,扁率,第一偏心率,第二偏心率参数,根据当前纬度信息动态获取当前地球半径。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:优化图片中心点经纬度误差:通过飞行器姿态模型计算的偏移量和相机姿态模型计算的偏移量,获取最终RTK天线位置到相机拍照中心点在地面上对应位置的偏移量,再结合动态获取的地球半径和当前RTK天线位置的经纬度,最终获取无人机航拍图片中心点经纬度。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:飞行器姿态模型建立包括,计算出从Z点到O点抖动之后的偏移量,构建四元数旋转向量V1并优化设备抖动误差;四元数旋转向量V1,经过构建的四元数q13旋转后,获得旋转后的四元数向量V1 /;计算飞行器经过偏航角变换后,将机体坐标系XY轴的偏移量映射到大地坐标系XgYg轴上的偏移量。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:通过将偏航角λ1取值范围[-180°,180°]分为四个区间[-180°,-90°)、[-90°,0°)、[0°,90°)、[90°,180°]四种情况,分别讨论大地坐标系上的偏移量ΔX1和ΔY1。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:相机姿态模型建立包括,构建四元数旋转向量V3=[0,0,0,-(h+Δz1)],通过俯仰角α3和侧滚角β3,构建两个四元数q31和q32,由四元数q31和q32获得新的四元数q33,四元数旋转向量V3,经过q33旋转后,获得旋转后的四元数向量V3 /,计算相机垂直地面中心点位置经过俯仰角α3和侧滚角β3变换后,将机体坐标系XY轴的偏移量映射到大地坐标系XgYg轴上的偏移量。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:CGCS2000坐标系动态获取地球半径包括,获取W-第一基本纬度函数、V-第二基本纬度函数、M-子午圈曲率半径和N-卯酉圈曲率半径;获取当前纬度地球半径R(平均曲率半径),纬圈半径r。
作为本发明所述无人机航拍图片中心点经纬度误差优化方法的一种优选方案,其中:优化图片中心点经纬度误差包括,获取Xg轴偏移量ΔX,Yg轴偏移量ΔY,获取拍摄图片中心点的经纬度坐标。
本发明的有益效果:
(1)巡检效率高,大幅缩减光伏电站巡检所需人数及巡检时间,节省人力运维成本,
(2)智能无人机机动性高,在空中飞行不受地面障碍物限制,实现自由移动。针对光伏电站幅员辽阔,地形起伏等特点,运用无人机巡检省时又省力;
(3)全自动飞行诊断,结合历史数据分析,可对光伏电站进行全面评估,还可对电站故障的出现进行有效预测;
(4)提升了电站巡检频次,有效提高电站巡检效率与精确性;
(5)解决了电站建设类型不同,组件难于巡检的难题,及以往人工巡检可能带来的人员安全问题。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为第一、二个实施例中的无人机图片中心点经纬度误差说明示意图。
图2为第一、二个实施例中的飞行器补偿示意图。
图3为第一、二个实施例中的飞行器设备偏航角[-90°,0°)偏移量计算示意图。
图4为第一、二个实施例中的相机垂直地面中心点位置俯仰和侧滚变换后偏移量计算示意图。
图5为第一、二个实施例中的相机中心点偏航角[-90°,0°)偏移量计算示意图。
图6为第一、二个实施例中的CGCS2000参考椭球示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
实施例1
参照图1~6,为本发明第一个实施例,该实施例提供了一种无人机航拍图片中心点经纬度误差优化方法,其包括建立飞行器姿态模型、建立相机姿态模型、CGCS2000坐标系动态获取地球半径和优化图片中心点经纬度误差四大步骤。
通过无人机航拍采集提供的经纬度信息,作为航拍图片中心点经纬度计算故障点经纬度的基础信息,对其进行使用和计算,可避免无人机设备产生的误差和相机抖动产生的误差,影响其计算结果精准度,进而最终获取优化后的航拍图片中心点经纬度信息。
无人机在飞行的过程中,需要建立地面坐标系0gXgYgZg与机体坐标系0XYZ来确定姿态角的信息,通过建立联系后,进而获取图片中热斑实际对应的经纬度信息,完成信息采集结果。
进一步的,无人机在飞行的过程中,俯仰角α为飞行器绕机体坐标系X轴旋转的角度,定位为水平为0,向下俯冲为正,后仰为负;侧滚角β为飞行器轴绕机体坐标系Y轴旋转的角度,定义为水平为0,左倾为负,右倾为正;偏航角λ为飞行器绕地面坐标系Zg轴旋转的角度,定义为正北为0,东半部0~180,西半部0~-180。飞行器与云台在飞行过程中,得到各自的姿态角。α1,β1,λ1为飞行器姿态角,α2,β2,λ2为云台姿态角。
飞行器姿态模型建立:通过飞行器姿态角和RTK天线当前经纬度,结合RTK天线到云台接口中心的补偿量构建四元数,计算出通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量。再通过构建机体坐标系和大地坐标系之间的关系,计算出通过偏航角变换后的补偿量,及RTK天线到云台中心在大地坐标系上正北方向的偏移量和正东方向的偏移量。
相机姿态模型建立:通过相机姿态角,结合RTK天线到地面的高度构建四元数,计算出相机垂直地面拍照中心点通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量。再通过构建机体坐标系和大地坐标系之间的关系,计算出通过偏航角变换后的补偿量,及相机拍照中心点位置变换后在大地坐标系上正北方向的偏移量和正东方向的偏移量。
CGCS2000坐标系动态获取地球半径:地球为一个椭球形,地球半径也在实时发生变化。通过CGCS2000坐标系,获取长半轴,短半轴,扁率,第一偏心率,第二偏心率参数,根据当前纬度信息动态获取当前地球半径。
优化图片中心点经纬度误差:通过飞行器姿态模型计算的偏移量和相机姿态模型计算的偏移量,获取最终RTK天线位置到相机拍照中心点在地面上对应位置的偏移量,再结合动态获取的地球半径和当前RTK天线位置的经纬度,最终获取无人机航拍图片中心点经纬度。
实施例2
参照图1~6,为本发明第二个实施例,该实施例基于上一个实施例:
飞行器姿态模型建立包括
步骤1.1:
由于飞行器四号机臂与云台接口中心之间存在补偿,如图2所示,以云台接口中心点为原点建立机体坐标系,Z点为飞行器四号机臂获取的经纬度信息,RTK天线在各轴上的补偿大小为a,b,c。通过Z点经纬度计算得出云台中心O点的经纬度信息。计算从Z点到O点抖动之后的偏移量,构建四元数旋转向量V1=[0,-a,b,-c]。
步骤1.2:
计算飞行器经过俯仰角与侧滚角变换后在机体坐标系上的偏移量,通过俯仰角α1和侧滚角β1,构建两个四元数q11和q12。旋转向量V1,先旋转q11再旋转q12,与先经过q12旋转再经过q11旋转,两次获得旋转后的向量不同。直接相乘存在顺序问题,因此将q11和q12切成n份,每次只乘q11和q12的1/n一共乘n次,不管q11在前还是q12在前,当n趋近于无穷大时,只有两端/>和/>不同,且对结果影响很小:
其中式中:q11、q12:属于四元数,代表旋转角度
单位四元数q有另外一种表示形式,其中/>是一半的旋转角,u是单位长度的旋转轴,同时也可以使用指数表示/>通过四元数相关性质,可得新的四元数q13:
q13代替q11和q12综合影响下的旋转,优化设备抖动误差。
步骤1.3:
四元数旋转向量V1,经过q13旋转后,获得旋转后的四元数向量V1 /。
V1/=q13V1q13*=q13V1q13 -1
由四元数向量V1 /=m1+a1i+b1j+c1z,可以获得飞行器四号机臂到云台接口中心的补偿在抖动后机体坐标系X轴的偏移量Δx1=a1,Y轴上的偏移量Δy1=b1,Z轴上的偏移量Δz1=c1。
步骤1.4:
计算飞行器经过偏航角变换后,将机体坐标系XY轴的偏移量映射到大地坐标系XgYg轴上的偏移量。通过将偏航角λ1取值范围[-180°,180°]分为四个区间[-180°,-90°)、[-90°,0°)、[0°,90°)、[90°,180°]四种情况,分别讨论大地坐标系上的偏移量ΔX1和ΔY1。如图3所示,以飞行器偏航角λ1∈[-90°,0°)为例,θ1=-λ1,可以获得设备抖动误差在大地坐标系上Xg轴偏移量ΔX1和Yg轴偏移量ΔY1。
ΔY1=-Δy1cos(θ1)-Δx1sin(θ1)
ΔX1=Δy1sin(θ1)-Δx1cos(θ1)
其中式中,ΔX1:代表在大地坐标系上的Xg轴偏移量;ΔY1:代表大地坐标系上Xg轴偏移量。
相机姿态模型建立包括:
步骤2.1:
通过大量数据验证和模式不同的选择,根据飞行器姿态角α1,β1,λ1和云台姿态角α2,β2,λ2,首先估算当前相机拍照时对应的真实姿态信息俯仰角α3,侧滚角β3,偏航角λ3。相机的偏移量与无人机当前高度H和设备在Z轴上的偏移量Δz1有关,构建四元数旋转向量V3=[0,0,0,-(h+Δz1)]。
步骤2.2:
计算相机中心点垂直地面的位置经过俯仰角与侧滚角变换后在机体坐标系上的偏移量。通过俯仰角α3和侧滚角β3,构建两个四元数q31和q32。由四元数q31和q32获得新的四元数q33。四元数旋转向量V3,经过q33旋转后,获得旋转后的四元数向量V3 /。
步骤2.3:
根据俯仰角α3和侧滚角β3综合变换,在机体坐标系上存在四种情况,(α3>=0,β3>=0),相机中心点落在第三象限;(α3>=0,β3<0),相机中心点落在第四象限;(α3<0,β3>=0),相机中心点落在第二象限;(α3<0,β3<0),相机中心点落在第一象限。以(α3<0,β3<0)为例,如图4所示,向量HO等于四元数向量V3,向量HA等于四元数向量V3 /。向量HA没有与XY轴平面相交,因此计算机体坐标系X轴与Y轴的偏移量分为以下几步:(1)通过向量HO与向量HA计算出(2)计算OB长度;(3)计算HB长度;(4)由向量HA与HB长度,计算向量HB;(5)向量OB=向量HB-向量HO;(6)由向量OB获得机体坐标系X轴的偏移量Δx3=OC,Y轴上的偏移量Δy3=OD。
步骤2.4:
计算相机垂直地面中心点位置经过俯仰角α3和侧滚角β3变换后,将机体坐标系XY轴的偏移量映射到大地坐标系XgYg轴上的偏移量。通过将相机偏航角λ3取值范围[-180°,180°]分为四个区间[-180°,-90°)、[-90°,0°)、[0°,90°)、[90°,180°]四种情况,分别讨论大地坐标系上的偏移量ΔX3和ΔY3。如图5所示,以飞行器偏航角λ3∈[-90°,0°)为例,θ3=-λ3,可以获得设备抖动误差在大地坐标系上Xg轴偏移量ΔX3和Yg轴偏移量ΔY3。此处与飞行器设备偏航角[-90°,0°)偏移量计算不同,由机体坐标系中心点计算变换后的相机中心点位置,因此在坐标系上偏移量的符号相反。
ΔY3=Δy3cos(θ3)+Δx3sin(θ3)
ΔX3=-Δy3sin(θ3)+Δx3cos(θ3)
其中式中,ΔX3:代表在大地坐标系上的Xg轴偏移量;ΔY3:代表大地坐标系上Xg轴偏移量。
CGCS2000坐标系动态获取地球半径包括:
CGCS2000坐标系,是我国当前最新的国家大地坐标系。如图6所示,将地球视为一个椭圆。通过CGCS2000坐标系,可以获取相关参数:长半轴a,短半轴b,扁率f,第一偏心率e1,第二偏心率e2。根据无人机飞行器4号机臂获取的纬度Lat,计算当前纬度的地球半径。
步骤3.1:
获取W-第一基本纬度函数,V-第二基本纬度函数:
步骤3.2:
获取M-子午圈曲率半径,N-卯酉圈曲率半径:
M=a*(1-e1 2)/(W*W*W)
N=M*V*V
步骤3.3:
获取当前纬度地球半径R(平均曲率半径),纬圈半径r:
r=Ncos(Lat)
优化图片中心点经纬度误差包括:
无人机获取得纬度Lat,经度Long。第一步获取设备误差在Xg轴偏移量ΔX1和Yg轴偏移量ΔY1,第二部获取相机抖动拍摄后的中心点在Xg轴偏移量ΔX3和Yg轴偏移量ΔY3。
步骤4.1:
获取Xg轴偏移量ΔX,Yg轴偏移量ΔY:
ΔY=ΔY1+ΔY3
ΔX=ΔX1+ΔX3
步骤4.2:
获取拍摄图片中心点的经纬度坐标(图1中O4点的经纬度):
Lat_new=Lat+ΔY/R*180/PI
Long_new=Long+ΔX/r*180/PI
其中式中,Long_new:代表经度坐标;
Lat_new:代表纬度坐标:
1.使用四元数算法的使用,根据四元数算法的性质,可以获取俯仰角和侧滚角综合影响变换后的状态,优化先俯仰角变换再侧滚角变换或者先侧滚角变换再俯仰角变换带来的误差;
2.对机体坐标系和大地坐标系进行转换,将姿态角分为两个部分,俯仰角和侧滚角为一组在机体坐标系上计算,偏航角为一组在机体坐标系上结合大地坐标系,来最终确定对应的偏移量,减少姿态角组合计算中带来的误差;
3.通过CGCS2000坐标系动态获取地球半径,减少地球半径计算经纬度偏移量时的误差。
本发明的有益效果:
国内光伏电站一般建设在大型山坡、戈壁、平原、沼泽、水域、厂房顶部、居民屋顶等,不同的光伏电站规模、形式、分布等均具有较大的差别。导致在后期的运维巡检过程中,存在很多不便。特别是光伏发电项目规模越大,电站巡检工作的复杂性越强,在对此种类型电站项目进行巡检工作时,一方面,常规的人工巡检方式不仅需要耗费大量的人工及时间成本,导致电站设备巡检不及时,影响电厂的经济效益;另一方面巡检人员工作也存在一定的危险性。
本发明方法通过无人机航拍采集提供的经纬度信息,作为航拍图片中心点经纬度计算故障点经纬度的基础信息,对其进行使用和计算,解决了无人机设备本身和空中抖动产生的误差和相机抖动产生的误差,最终获取了优化后的航拍图片中心点经纬度信息。
重要的是,应注意,在多个不同示例性实施方案中示出的本申请的构造和布置仅是例示性的。尽管在此公开内容中仅详细描述了几个实施方案,但参阅此公开内容的人员应容易理解,在实质上不偏离该申请中所描述的主题的新颖教导和优点的前提下,许多改型是可能的(例如,各种元件的尺寸、尺度、结构、形状和比例、以及参数值(例如,温度、压力等)、安装布置、材料的使用、颜色、定向的变化等)。例如,示出为整体成形的元件可以由多个部分或元件构成,元件的位置可被倒置或以其它方式改变,并且分立元件的性质或数目或位置可被更改或改变。因此,所有这样的改型旨在被包含在本发明的范围内。可以根据替代的实施方案改变或重新排序任何过程或方法步骤的次序或顺序。在权利要求中,任何“装置加功能”的条款都旨在覆盖在本文中所描述的执行所述功能的结构,且不仅是结构等同而且还是等同结构。在不背离本发明的范围的前提下,可以在示例性实施方案的设计、运行状况和布置中做出其他替换、改型、改变和省略。因此,本发明不限制于特定的实施方案,而是扩展至仍落在所附的权利要求书的范围内的多种改型。
此外,为了提供示例性实施方案的简练描述,可以不描述实际实施方案的所有特征(即,与当前考虑的执行本发明的最佳模式不相关的那些特征,或于实现本发明不相关的那些特征)。
应理解的是,在任何实际实施方式的开发过程中,如在任何工程或设计项目中,可做出大量的具体实施方式决定。这样的开发努力可能是复杂的且耗时的,但对于那些得益于此公开内容的普通技术人员来说,不需要过多实验,所述开发努力将是一个设计、制造和生产的常规工作。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (9)
1.一种无人机航拍图片中心点经纬度误差优化方法,其特征在于:包括,
建立飞行器姿态模型,计算出通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再构建机体坐标系和大地坐标系之间的关系;
建立相机姿态模型,计算出相机垂直地面拍照中心点通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再构建机体坐标系和大地坐标系之间的关系;
CGCS2000坐标系动态获取地球半径,通过CGCS2000坐标系获取当前地球半径;以及,
优化图片中心点经纬度误差,通过飞行器姿态模型计算和相机姿态模型计算获取无人机航拍图片中心点经纬度;
优化图片中心点经纬度误差:
通过飞行器姿态模型计算的偏移量和相机姿态模型计算的偏移量,获取最终RTK天线位置到相机拍照中心点在地面上对应位置的偏移量,再结合动态获取的地球半径和当前RTK天线位置的经纬度,最终获取无人机航拍图片中心点经纬度。
2.如权利要求1所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:飞行器姿态模型建立:
通过飞行器姿态角和RTK天线当前经纬度,结合RTK天线到云台接口中心的补偿量构建四元数,计算出通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再通过构建机体坐标系和大地坐标系之间的关系,计算出通过偏航角变换后的补偿量,及RTK天线到云台中心在大地坐标系上正北方向的偏移量和正东方向的偏移量。
3.如权利要求1所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:相机姿态模型建立:
通过相机姿态角,结合RTK天线到地面的高度构建四元数,计算出相机垂直地面拍照中心点通过俯仰角和侧滚角综合变换后在机体坐标系上的偏移量,再通过构建机体坐标系和大地坐标系之间的关系,计算出通过偏航角变换后的补偿量,及相机拍照中心点位置变换后在大地坐标系上正北方向的偏移量和正东方向的偏移量。
4.如权利要求1所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:CGCS2000坐标系动态获取地球半径:
通过CGCS2000坐标系,获取长半轴,短半轴,扁率,第一偏心率,第二偏心率参数,根据当前纬度信息动态获取当前地球半径。
5.如权利要求2所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:飞行器姿态模型建立包括,
计算出从Z点到O点抖动之后的偏移量,构建四元数旋转向量V1并优化设备抖动误差;四元数旋转向量V1,经过构建的四元数q13旋转后,获得旋转后的四元数向量V1 /;计算飞行器经过偏航角变换后,将机体坐标系XY轴的偏移量映射到大地坐标系XgYg轴上的偏移量。
6.如权利要求5所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:通过将偏航角λ1取值范围[-180°,180°]分为四个区间[-180°,-90°)、[-90°,0°)、[0°,90°)、[90°,180°]四种情况,分别讨论大地坐标系上的偏移量ΔX1和ΔY1。
7.如权利要求3所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:相机姿态模型建立包括,
构建四元数旋转向量V3=[0,0,0,-(h+Δz1)],通过俯仰角α3和侧滚角β3,构建两个四元数q31和q32,由四元数q31和q32获得新的四元数q33,四元数旋转向量V3,经过q33旋转后,获得旋转后的四元数向量V3 /,计算相机垂直地面中心点位置经过俯仰角α3和侧滚角β3变换后,将机体坐标系XY轴的偏移量映射到大地坐标系XgYg轴上的偏移量。
8.如权利要求4所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:CGCS2000坐标系动态获取地球半径包括,
获取W-第一基本纬度函数、V-第二基本纬度函数、M-子午圈曲率半径和N-卯酉圈曲率半径;获取当前纬度地球半径R(平均曲率半径),纬圈半径r。
9.如权利要求1所述的无人机航拍图片中心点经纬度误差优化方法,其特征在于:优化图片中心点经纬度误差包括,
获取Xg轴偏移量ΔX,Yg轴偏移量ΔY,获取拍摄图片中心点的经纬度坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111434436.8A CN114397900B (zh) | 2021-11-29 | 2021-11-29 | 一种无人机航拍图片中心点经纬度误差优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111434436.8A CN114397900B (zh) | 2021-11-29 | 2021-11-29 | 一种无人机航拍图片中心点经纬度误差优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114397900A CN114397900A (zh) | 2022-04-26 |
CN114397900B true CN114397900B (zh) | 2024-02-09 |
Family
ID=81225610
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111434436.8A Active CN114397900B (zh) | 2021-11-29 | 2021-11-29 | 一种无人机航拍图片中心点经纬度误差优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114397900B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115334247B (zh) * | 2022-10-11 | 2023-01-10 | 齐鲁空天信息研究院 | 摄像模块标校方法、视觉定位方法、装置及电子设备 |
CN115793698A (zh) * | 2023-02-07 | 2023-03-14 | 北京四维远见信息技术有限公司 | 自动姿态控制系统和方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101609149A (zh) * | 2009-07-13 | 2009-12-23 | 北京航空航天大学 | 一种提高机载激光雷达定姿精度的方法 |
CN103558619A (zh) * | 2013-11-06 | 2014-02-05 | 中测新图(北京)遥感技术有限责任公司 | 获取航拍像片的外方位元素的方法 |
CN106444844A (zh) * | 2016-08-31 | 2017-02-22 | 南京航空航天大学 | 一种用于多旋翼飞行器的图像程控采集方法 |
CN106500669A (zh) * | 2016-09-22 | 2017-03-15 | 浙江工业大学 | 一种基于四旋翼imu参数的航拍图像矫正方法 |
CN107247458A (zh) * | 2017-05-24 | 2017-10-13 | 中国电子科技集团公司第二十八研究所 | 无人机视频图像目标定位系统、定位方法及云台控制方法 |
CN111288890A (zh) * | 2020-02-13 | 2020-06-16 | 福建农林大学 | 基于双目摄影测量技术的道路标志尺寸及高度自动测量方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6354425B2 (ja) * | 2014-07-30 | 2018-07-11 | 株式会社デンソー | 車載カメラの取り付け姿勢検出方法およびその装置 |
-
2021
- 2021-11-29 CN CN202111434436.8A patent/CN114397900B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101609149A (zh) * | 2009-07-13 | 2009-12-23 | 北京航空航天大学 | 一种提高机载激光雷达定姿精度的方法 |
CN103558619A (zh) * | 2013-11-06 | 2014-02-05 | 中测新图(北京)遥感技术有限责任公司 | 获取航拍像片的外方位元素的方法 |
CN106444844A (zh) * | 2016-08-31 | 2017-02-22 | 南京航空航天大学 | 一种用于多旋翼飞行器的图像程控采集方法 |
CN106500669A (zh) * | 2016-09-22 | 2017-03-15 | 浙江工业大学 | 一种基于四旋翼imu参数的航拍图像矫正方法 |
CN107247458A (zh) * | 2017-05-24 | 2017-10-13 | 中国电子科技集团公司第二十八研究所 | 无人机视频图像目标定位系统、定位方法及云台控制方法 |
CN111288890A (zh) * | 2020-02-13 | 2020-06-16 | 福建农林大学 | 基于双目摄影测量技术的道路标志尺寸及高度自动测量方法 |
Non-Patent Citations (3)
Title |
---|
Research on the image rotation of the polar coordinates infrared inertial stable platform with four reflected mirrors;Qiu Zhaowen;《 Infrared and Laser Engineering》;全文 * |
基于四元数法的运动平台红外图像旋转角推导;吕波;《红外技术》;第353-358页 * |
无人机航摄系统成图精度研究;李飞达;《基础科学》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114397900A (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114397900B (zh) | 一种无人机航拍图片中心点经纬度误差优化方法 | |
CN110580717B (zh) | 一种针对电力杆塔的无人机自主巡检航线生成方法 | |
CN110554704A (zh) | 一种基于无人机的风机叶片自主巡检方法 | |
CN105698762B (zh) | 一种单机航迹上基于不同时刻观测点的目标快速定位方法 | |
CN110908401B (zh) | 一种针对未知杆塔结构的无人机自主巡检方法 | |
CN107830846B (zh) | 一种利用无人机和卷积神经网络测量通信塔天线角度方法 | |
CN110244766B (zh) | 一种光伏电站无人机巡检航线的规划方法及系统 | |
CN111652964A (zh) | 基于数字孪生的电力巡检无人机辅助定位方法和系统 | |
CN110706273B (zh) | 一种基于无人机的实时塌方区域面积的测量方法 | |
CN109272574A (zh) | 基于投影变换的线阵旋转扫描相机成像模型构建方法和标定方法 | |
CN103673995A (zh) | 一种线阵推扫式相机在轨光学畸变参数标定方法 | |
CN112710311B (zh) | 一种地形自适应无人机三维实景重建航摄点自动规划方法 | |
CN112947526B (zh) | 一种无人机自主降落方法和系统 | |
CN112102403B (zh) | 用于输电塔场景下的自主巡检无人机的高精度定位方法及其系统 | |
CN111123962A (zh) | 一种用于电力杆塔巡检的旋翼无人机重定位拍照方法 | |
CN108876828A (zh) | 一种无人机图像批处理三维重建方法 | |
CN115511956A (zh) | 一种无人机成像定位方法 | |
CN112101118A (zh) | 一种远景与近景图像结合的光伏组件定位与匹配方法 | |
CN113739767B (zh) | 针对国产面阵摆扫成像系统获取的影像生产正射影像的方法 | |
CN114564031A (zh) | 一种基于旋翼无人机实现风机巡检拍摄的路径规划方法 | |
CN117308889B (zh) | 联合拼接ccd所有芯片单元的高精度天体测量方法 | |
CN112905955B (zh) | 一种风机正反面自主巡检规划方法 | |
CN108921896A (zh) | 一种融合点线特征的下视视觉罗盘 | |
CN113295171A (zh) | 一种基于单目视觉的旋转刚体航天器姿态估计方法 | |
CN115270522B (zh) | 一种基于wgs84坐标仿真跟踪目标设备的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |