CN111273282B - 一种无人机机载雷达地面扫描覆盖区域计算方法 - Google Patents
一种无人机机载雷达地面扫描覆盖区域计算方法 Download PDFInfo
- Publication number
- CN111273282B CN111273282B CN202010105096.3A CN202010105096A CN111273282B CN 111273282 B CN111273282 B CN 111273282B CN 202010105096 A CN202010105096 A CN 202010105096A CN 111273282 B CN111273282 B CN 111273282B
- Authority
- CN
- China
- Prior art keywords
- point
- scanning
- axis
- coordinate system
- calculating
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000012800 visualization Methods 0.000 abstract description 4
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/60—Rotation of whole images or parts thereof
- G06T3/604—Rotation of whole images or parts thereof using coordinate rotation digital computer [CORDIC] devices
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种无人机机载雷达地面扫描覆盖区域计算方法,确定整个扫描覆盖区域地理范围,建立一个参空间坐标系,用于边界点的相对位置计算,计算雷达扫描区域在OXYZ坐标系Y轴方向的有效范围,计算任意边界点在OXYZ坐标系下的坐标,计算扫描区域边界线的地理范围,通过坐标转换,最终获得一组采用经纬度表示的扫描区域边界点。本发明针对一种具有圆锥状扫描范围的机载雷达,提出了其地面扫描覆盖区域地理位置的计算方法,可以实现雷达在任意扫描角度时地面扫描覆盖区域的计算,对于无人机任务规划或者飞行监控等任务中,雷达扫描区域的可视化提供了计算基础。
Description
技术领域
本发明涉及无人机应用领域,尤其是一种无人机机载雷达监控方法。
背景技术
机载雷达是一种常用的无人机目标侦察手段,使用装备在无人机上的雷达对指定地面区域进行扫描,可以发现该区域感兴趣目标并获取其准确地理位置。为了使该类侦察无人机准确高效地完成指定区域的侦察任务,需要根据侦察雷达的工作方式,预先完成相应的任务规划,并在飞行过程中实时监控侦察雷达的扫描区域。由于侦察雷达工作方式的独特性,针对本发明所涉及的侦察无人机,尚无有效的任务规划和任务监控辅助手段。
发明内容
为了克服现有技术的不足,本发明提供一种无人机机载雷达地面扫描覆盖区域计算方法。为了解决该类无人机的任务规划和监控问题,本发明可准确计算雷达扫描覆盖区域的地理范围,实现扫描覆盖区域的可视化,用于完成任务规划的有效性评估以及侦察任务的辅助监控。
本发明解决其技术问题所采用的技术方案的计算步骤如下:
步骤1:,在以无人机地面投影点为中心的参考地平面上,将雷达地面扫描覆盖区域边界线表征为该平面上的一组离散边界点,通过计算这一组离散边界点的经纬度,实现整个扫描覆盖区域地理范围的计算;
步骤2:建立一个参空间坐标系,用于边界点的相对位置计算;假设雷达执行扫描工作时,飞机飞行平稳,俯仰角和滚转角为0°,以飞机位置的地面投影点为坐标原点,飞机航向方向为Y轴,垂直于地面方向为Z轴建立空间直角坐标系OXYZ,则扫描覆盖区域位于XY平面,且关于Y轴对称;
步骤3:计算雷达扫描区域在OXYZ坐标系Y轴方向的有效范围;
根据飞机飞行高度H、雷达框架角α(雷达扫描中心线与飞机中轴线夹角)、雷达扫描视场角θ以及雷达最大扫描距离L0,计算覆盖区域在Y轴最远点p1,最近点p2以及最远有效点p0,各点的Y坐标计算方法如下:
步骤4:计算任意边界点在OXYZ坐标系下的坐标;
在步骤3计算出的Y轴有效扫描覆盖范围内任取一点d,过点d作X轴的平行线,与XY平面上的边界曲线相交于d1和d2点,d1和d2点关于Y轴对称;过d点做扫描圆锥体的横截面,e1和e2为截面圆与YZ平面的相交线的端点,直线pp1和直线pp2为扫描圆锥面与YZ平面的相交线,e1e2为过直线d1d2的圆锥截面圆与YZ平面的相交线;首先根据直线pp1、pp2以及e1e2的相交关系,根据空间几何关系计算e1,e2的坐标,计算公式如下:
然后计算截面圆圆心点e和点d的距离Led:
步骤5:计算扫描区域边界线的地理范围;
在步骤2计算所得的Y轴有效扫描覆盖范围内均取N个离散点,按照步骤4计算每一个离散点d对应边界点d1,d2的在OXYZ坐标系下的坐标;结合飞机位置经纬度和航向角计算出各边界点在高斯-克吕格投影直角坐标系下的地理坐标:
其中(xp,yp)为无人机地面投影点在高斯-克吕格投影直角坐标系下的地理坐标,φ为无人机航向角,(xd1,yd1)为边界点d1在OXYZ坐标系的XY平面坐标,(x′d1,y′d1)为边界点d1在高斯-克吕格投影直角坐标系下的地理坐标;
根据高斯-克吕格投影直角坐标系与大地坐标系的投影关系,将该高斯直角坐标转换为对应的经纬度坐标,最终获得一组采用经纬度表示的扫描区域边界点。
本发明有益效果针对一种具有圆锥状扫描范围的机载雷达,提出了其地面扫描覆盖区域地理位置的计算方法,可以实现雷达在任意扫描角度时地面扫描覆盖区域的计算。本发明对于无人机任务规划或者飞行监控等任务中,雷达扫描区域的可视化提供了计算基础。
附图说明
图1是本发明机载雷达扫描区域示意图。
图2是本发明机载雷达扫描区域解析图。
图3是本发明机载雷达扫描区域截面图。
图4是本发明本发明方法的流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明涉及的无人机机载雷达工作过程示意图如图1所示。雷达安装在无人机机头位置,其立体扫描范围为圆锥状,扫描中心线位于飞机侧平面上。无人机飞行过程中雷达按一定扫描角度对地面进行持续扫描,扫描到地面目标的有效区域可表示为扫描圆锥体与地平面相切部分,基本扫描形状如图1所示。为了在无人机任务规划,或者飞行监控过程中实现雷达扫描区域的可视化,就需要实时计算无人机飞行过程中,雷达地面扫描覆盖区域的准确地理位置。本发明针对这类问题提出了一种适于计算机实现的计算方法。
现结合附图说明本发明所述方法的具体实现过程。假设雷达工作时飞机飞行平稳,飞机在某一时刻位姿表示为(B,L,H,φ,0,0),其中B为纬度,L为经度,H为高度,φ为航向角,飞机俯仰和滚转角为0°。根据本发明,雷达扫描覆盖区域计算的具体实施过程如下:
步骤1:,在以无人机地面投影点为中心的参考地平面上,将雷达地面扫描覆盖区域边界线表征为该平面上的一组离散边界点,通过计算这一组离散边界点的经纬度,实现整个扫描覆盖区域地理范围的计算;
步骤2:建立一个参空间坐标系,用于边界点的相对位置计算;假设雷达执行扫描工作时,飞机飞行平稳,俯仰角和滚转角为0°,以飞机位置的地面投影点为坐标原点,飞机航向方向为Y轴,垂直于地面方向为Z轴建立空间直角坐标系OXYZ,则扫描覆盖区域位于XY平面,且关于Y轴对称,扫描区域基本形状如图2所示。
步骤3:计算雷达扫描区域在OXYZ坐标系Y轴方向的有效范围;
根据飞机飞行高度H、雷达框架角α(雷达扫描中心线与飞机中轴线夹角)、雷达扫描视场角θ以及雷达最大扫描距离L0,计算覆盖区域在Y轴最远点p1,最近点p2以及最远有效点p0,如图2所示,各点的Y坐标计算方法如下:
步骤4:计算任意边界点在OXYZ坐标系下的坐标;
在步骤3计算出的Y轴有效扫描覆盖范围内任取一点d,过点d作X轴的平行线,与XY平面上的边界曲线相交于d1和d2点,如图2所示,d1和d2点关于Y轴对称;过d点做扫描圆锥体的横截面如图3所示,e1和e2为截面圆与YZ平面的相交线的端点,如图2所示,直线pp1和直线pp2为扫描圆锥面与YZ平面的相交线,e1e2为过直线d1d2的圆锥截面圆与YZ平面的相交线;首先根据直线pp1、pp2以及e1e2的相交关系,根据空间几何关系计算e1,e2的坐标,计算公式如下:
然后计算截面圆圆心点e和点d的距离Led:
步骤5:计算扫描区域边界线的地理范围;
在步骤2计算所得的Y轴有效扫描覆盖范围内均取N个离散点,按照步骤4计算每一个离散点d对应边界点d1,d2的在OXYZ坐标系下的坐标;结合飞机位置经纬度和航向角计算出各边界点在高斯-克吕格投影直角坐标系下的地理坐标:
其中(xp,yp)为无人机地面投影点在高斯-克吕格投影直角坐标系下的地理坐标,φ为无人机航向角,(xd1,yd1)为边界点d1在OXYZ坐标系的XY平面坐标,(x′d1,y′d1)为边界点d1在高斯-克吕格投影直角坐标系下的地理坐标;
根据高斯-克吕格投影直角坐标系与大地坐标系的投影关系,将该高斯直角坐标转换为对应的经纬度坐标,最终获得一组采用经纬度表示的扫描区域边界点。
实施例如下:
1)以飞机投影点为原点建立空间直角坐标系OXYZ,如图2所示。飞机位于p点,坐标原点O为飞机位置p垂直于地面的投影点,Y轴与飞机的航向平行,Z轴垂直于地面向上。雷达工作时其扫描覆盖区域如图2中虚线区域所示,该区域位于坐标系OXYZ的XY平面上,且关于Y轴对称。
2)计算雷达扫描区域在OXYZ坐标系Y轴方向的有效范围。设雷达框架角(雷达扫描轴线与飞机中轴线夹角)为α(逆时针为正),扫描视场范围角θ,雷达最大扫描距离为L0,飞机当前高度为H。此时,扫描区域在Y轴的最远点表示为p1,最近点表示为p2,与雷达最大扫描距离对应的有效扫描距离点表示为p0。记直线pp1与Y轴的夹角为β1=α+θ/2,直线pp2与Y轴的夹角为β2=α-θ/2,直线pp0于与Y轴的夹角为β0=-arcsin(H/L0)。则对于雷达扫描覆盖区域的中轴线上任意一点d,其在Y轴方向有效范围为:
其中yd表示扫描区域内任意一点d的Y轴坐标,H为飞机当前高度,β1,β2,β0分别为直线pp1,pp2,pp0与Y轴的夹角。
3)计算雷达扫描覆盖区域边界线任意一点的坐标。对于扫描区域中轴线上任意一点d,其对应的边界线上的点d1,d2如图2所示。点d1,d2的Y坐标即为yd,Z坐标为0,X坐标计算方法如下:
过点d做垂直于扫描圆锥体中心线的截面如图3所示,截面为圆形,中心点记为e,则边界点d1,d2为截面与XY平面相交线的两端点。记该截面与YZ平面相交线的两个端点为e1,e2,则e1,e2分为位于直线pp1和pp2上。直线pp1、pp2和e1e2在YZ平面的直线方程为:
则e1,e2在YZ平面分别为直线pp1和直线pp2与直线e1e2的交点,根据上述直线方程可得点e1在YZ平面的坐标为:
点e2在YZ平面的坐标为:
根据e1,e2两点坐标计算得图3中截面圆的半径R:
线段ed的长度D为:
4)计算扫描区域边界线的地理范围。在步骤(2)计算的扫描区域Y轴范围内等距取N个离散点,对于任一离散点d,根据步骤(3)计算其对应的扫描区域边界点d1,d2在XY平面的坐标(xd1,yd1)和(xd2,yd2):
其中yd表示扫描区域内任意一点d的Y轴坐标。
设无人机在点p处的位姿参数为(B,L,H,φ,0,0),首先根据大地坐标系与高斯-克吕格投影直角坐标系的投影关系,计算其在高斯-克吕格投影直角坐标系下的坐标(xp,yp)(地理坐标系转换通用计算公式),然后计算Y轴上的离散点d对应的两个边界点d1,d2在高斯-克吕格投影直角坐标系下的坐标(x′d1,y′d1),(x′d2,y′d2):
再将该直角坐标转换为对应的经纬度坐标。
对于扫描区域在Y轴有效范围内等距取的N个离散点,按照上述计算方法依次计算每一个离散点d对应的两个扫描区域边界点d1,d2的经纬度坐标,最终生成由2N的边界点组成的扫描区域边界点集,根据每一个边界点的经纬度可在地图上绘制出相应的扫描覆盖区域图形。
Claims (1)
1.一种无人机机载雷达地面扫描覆盖区域计算方法,其特征在于包括下述步骤:
步骤1:在以无人机地面投影点为中心的参考地平面上,将雷达地面扫描覆盖区域边界线表征为该平面上的一组离散边界点,通过计算这一组离散边界点的经纬度,实现整个扫描覆盖区域地理范围的计算;
步骤2:建立一个参空间坐标系,用于边界点的相对位置计算;假设雷达执行扫描工作时,飞机飞行平稳,俯仰角和滚转角为0°,以飞机位置的地面投影点为坐标原点,飞机航向方向为Y轴,垂直于地面方向为Z轴建立空间直角坐标系OXYZ,则扫描覆盖区域位于XY平面,且关于Y轴对称;
步骤3:计算雷达扫描区域在OXYZ坐标系Y轴方向的有效范围;
根据飞机飞行高度H、雷达框架角α即雷达扫描中心线与飞机中轴线夹角、雷达扫描视场角θ以及雷达最大扫描距离L0,计算覆盖区域在Y轴最远点p1,最近点p2以及最远有效点p0,各点的Y坐标计算方法如下:
步骤4:计算任意边界点在OXYZ坐标系下的坐标;
在步骤3计算出的Y轴有效扫描覆盖范围内任取一点d,过点d作X轴的平行线,与XY平面上的边界曲线相交于d1和d2点,d1和d2关于Y轴对称;过d点做扫描圆锥体的横截面,e1和e2为截面圆与YZ平面的相交线的端点,直线pp1和直线pp2为扫描圆锥面与YZ平面的相交线,e1e2为过直线d1d2的圆锥截面圆与YZ平面的相交线;首先根据直线pp1、pp2以及e1e2的相交关系,根据空间几何关系计算e1,e2的坐标,计算公式如下:
然后计算截面圆圆心点e和点d的距离Led:
步骤5:计算扫描区域边界线的地理范围;
在步骤2计算所得的Y轴有效扫描覆盖范围内均取N个离散点,按照步骤4计算每一个离散点d对应边界点d1,d2的在OXYZ坐标系下的坐标;结合飞机位置经纬度和航向角计算出各边界点在高斯-克吕格投影直角坐标系下的地理坐标:
根据高斯-克吕格投影直角坐标系与大地坐标系的投影关系,将高斯直角坐标转换为对应的经纬度坐标,最终获得一组采用经纬度表示的扫描区域边界点。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010105096.3A CN111273282B (zh) | 2020-02-20 | 2020-02-20 | 一种无人机机载雷达地面扫描覆盖区域计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010105096.3A CN111273282B (zh) | 2020-02-20 | 2020-02-20 | 一种无人机机载雷达地面扫描覆盖区域计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111273282A CN111273282A (zh) | 2020-06-12 |
CN111273282B true CN111273282B (zh) | 2023-04-14 |
Family
ID=70999124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010105096.3A Active CN111273282B (zh) | 2020-02-20 | 2020-02-20 | 一种无人机机载雷达地面扫描覆盖区域计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111273282B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112206441B (zh) * | 2020-10-12 | 2022-02-08 | 江西省智能产业技术创新研究院 | 一种消防机器人调度系统的协同调度方法 |
CN112698317A (zh) * | 2020-12-04 | 2021-04-23 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种机载雷达探测区域在数字地图上的投影方法及装置 |
CN115350547B (zh) * | 2022-10-24 | 2023-01-06 | 湖南九九智能环保股份有限公司 | 一种基于激光雷达定位的喷雾除尘方法及系统 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1342677A (en) * | 1970-06-08 | 1974-01-03 | Communications Satellite Corp | Conical torus antenna |
JP2012242368A (ja) * | 2011-05-24 | 2012-12-10 | Nsk Ltd | 移動環境認識装置及び方法 |
CN106885531A (zh) * | 2017-04-20 | 2017-06-23 | 河北科技大学 | 基于二维激光雷达的车斗形容器三维扫描系统标定方法 |
CN109752703A (zh) * | 2019-03-18 | 2019-05-14 | 上海盎维信息技术有限公司 | 具有锥形扫描区域的激光雷达及扫描方法 |
CN109932701A (zh) * | 2019-04-02 | 2019-06-25 | 哈尔滨工程大学 | 一种模拟船用雷达的目标船回波2d成像方法 |
CN209878987U (zh) * | 2019-03-18 | 2019-12-31 | 上海盎维信息技术有限公司 | 具有锥形扫描区域的激光雷达 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2888944B1 (fr) * | 2005-07-20 | 2007-10-12 | Eurocopter France | Procede de detection par telemetrie d'objets filaires suspendus |
-
2020
- 2020-02-20 CN CN202010105096.3A patent/CN111273282B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1342677A (en) * | 1970-06-08 | 1974-01-03 | Communications Satellite Corp | Conical torus antenna |
JP2012242368A (ja) * | 2011-05-24 | 2012-12-10 | Nsk Ltd | 移動環境認識装置及び方法 |
CN106885531A (zh) * | 2017-04-20 | 2017-06-23 | 河北科技大学 | 基于二维激光雷达的车斗形容器三维扫描系统标定方法 |
CN109752703A (zh) * | 2019-03-18 | 2019-05-14 | 上海盎维信息技术有限公司 | 具有锥形扫描区域的激光雷达及扫描方法 |
CN209878987U (zh) * | 2019-03-18 | 2019-12-31 | 上海盎维信息技术有限公司 | 具有锥形扫描区域的激光雷达 |
CN109932701A (zh) * | 2019-04-02 | 2019-06-25 | 哈尔滨工程大学 | 一种模拟船用雷达的目标船回波2d成像方法 |
Non-Patent Citations (3)
Title |
---|
Imaging for a Forward Scanning Automotive Synthetic Aperture Radar;Gishkori, Shahzad等;《IEEE Transactions on Aerospace & Electronic Systems》;20180920;第55卷(第03期);第1420-1434页 * |
基于激光雷达三维点云数据的配准技术;王小清;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180415(第04期);第I136-2498页 * |
机载激光雷达测深系统定位模型与视准轴误差影响分析;俞家勇 等;《红外与激光工程》;20190630;第48卷(第06期);第337-345页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111273282A (zh) | 2020-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111273282B (zh) | 一种无人机机载雷达地面扫描覆盖区域计算方法 | |
CN106774392B (zh) | 一种电力线路巡检过程中飞行路径的动态规划方法 | |
US10017271B2 (en) | Methods of three dimensional (3D) airflow sensing and analysis | |
WO2021082396A1 (zh) | 基于低空空域限制条件的无人机飞行网络建模方法 | |
US20170269612A1 (en) | Flight control methods for operating close formation flight | |
CN104298248A (zh) | 旋翼无人机精确视觉定位定向方法 | |
CN112789672B (zh) | 控制和导航系统、姿态优化、映射和定位技术 | |
CN107084725A (zh) | 一种多旋翼无人机电力巡检三维航线规划方法 | |
CN109031312A (zh) | 适用于烟囱内部作业的飞行平台定位装置和定位方法 | |
CN113963276A (zh) | 一种输电线路无人机自主巡检方法及系统 | |
CN114490815B (zh) | 一种基于地形通视的无人机任务载荷侦察覆盖率计算方法 | |
CN103954953A (zh) | 一种基于数据驱动的机载激光雷达盲源误差补偿方法 | |
CN112597664B (zh) | 一种基于无人机的铁路既有线精细实景三维建模航带设计方法 | |
CN112947526B (zh) | 一种无人机自主降落方法和系统 | |
CN108627142A (zh) | 一种结合离线高程和机载光电吊舱的目标定位方法 | |
CN113093801B (zh) | 一种多无人机协同编队的航路文件制作方法 | |
CN111208497A (zh) | 一种机载激光雷达系统平差处理方法 | |
CN109444915B (zh) | 一种基于激光雷达数据的危险区域预判方法 | |
CN112987781A (zh) | 一种无人机航线生成方法及装置 | |
CN116758442A (zh) | 一种基于激光扫描技术的输电线路数据处理方法 | |
CN117288166A (zh) | 一种基于bim的无人机测绘装置及方法 | |
CN110148218B (zh) | 一种大批量机载LiDAR点云数据整体优化的方法 | |
CN114018244B (zh) | 一种基于无人机光电平台的目标跟踪航路生成方法 | |
CN115525062A (zh) | 面向返回式航天器测量任务无人机实时在线航迹规划方法 | |
CN106371096B (zh) | 机载双天线InSAR三维构像模型构建方法 |
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 |