CN108507576B - 一种中高山区航空重力测量三维航迹规划方法 - Google Patents
一种中高山区航空重力测量三维航迹规划方法 Download PDFInfo
- Publication number
- CN108507576B CN108507576B CN201810256318.4A CN201810256318A CN108507576B CN 108507576 B CN108507576 B CN 108507576B CN 201810256318 A CN201810256318 A CN 201810256318A CN 108507576 B CN108507576 B CN 108507576B
- Authority
- CN
- China
- Prior art keywords
- flight path
- cutting line
- line
- planning method
- path planning
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/16—Measuring gravitational fields or waves; Gravimetric prospecting or detecting specially adapted for use on moving platforms, e.g. ship, aircraft
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Automation & Control Theory (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明涉及一种中高山区航空重力测量三维航迹规划方法，根据航空重力仪对飞行起伏坡度和起伏周期的要求，以测区地形数据为基础，求取满足一定要求的极大值和求取测线与切割线交叉点处均值相结合的方式获得控制点，然后对控制点进行线性插值并整体抬升一定高度，得到三维航迹曲线。本发明有效降低了平均飞行高度，保证了测线和切割线交叉点位置高度在理论上重合，从而确保了交叉点上的重力场差值尽可能小。同时，还可将计算得到的输出控制点提供给导航设备，从而在航空重力缓起伏飞行中突破了传统目视导航方法，实现了基于导航设备的精确导航，确保飞行高度和起伏坡度满足重力测量需求，减少了人为因素造成的飞机俯仰角过大而导致的虚假重力异常问题。
Description
技术领域
本发明属于航空重力勘探技术领域，涉及一种中高山区航空重力测量三维航迹规划方法。
背景技术
航空重力测量是将航空重力仪等设备装载在飞行器上，在测量地区上空按照预先设定的航迹对地球重力场进行测量的地球物理方法。航空重力测量航迹规划通常以平飞为主，在中高山区(高差大于500m、海拔高度大于1000m的复杂山地)采取平飞航空重力测量，导致飞行高度增大，使得重力仪数据中采集到的有效异常信息大大减弱，降低了地质勘探的精准度。周锡华等在《起伏飞行在航空重力测量的应用研究》(物探与化探，2015，39：98-104)中从理论上探讨了航空重力起伏飞行的可行性和有效性，然而并未提出起伏飞行的三维航迹规划方法。当前，航空重力起伏飞行仅能依靠目视进行，由于缺少精确的三维航迹规划，一是容易造成飞机俯仰角过大，导致出现虚假重力异常问题，严重干扰了测量结果的准确度；二是无法保证测线和切割线交点位置的飞行高度重合，不利于全区重力场调平。
发明内容
本发明的目的在于克服现有技术的缺陷，提供一种计算步骤简单、运算效率高，可满足中高山区航空重力测量需求的三维航迹规划方法。
为实现上述目的，本发明采用了如下技术方案：
一种中高山区航空重力测量三维航迹规划方法，包括以下步骤：(1)求取测线和切割线上对应地形的最高峰、极大值；(2)以最高峰为参考点，删除测线和切割线上坡度大于3°的极大值；(3)删除剩余极大值中间距小于窗口距离的极大值；剩余的极大值作为测线和切割线的一般控制点；(4)对一般控制点进行线性插值，形成测线和切割线的初始航迹曲线；(5)求取测线和切割线交叉点位置高差，当高差均小于阈值时，所述一般控制点即为输出控制点，转步骤(7)；(6)当高差超出阈值时，求取交叉点位置航迹均值作为测线和切割线的特殊控制点，然后删除测线和切割线上相邻控制点坡度大于3°的一般控制点；对测线和切割线的特殊控制点和未删除的一般控制点进行线性插值，更新测线和切割线航迹，再次求取测线和切割线交叉点位置高差，直至高差均小于阈值；特殊控制点和未删除的一般控制点即为输出控制点；(7)对输出控制点进行插值，形成三维航迹曲线。
进一步，所述步骤(1)之前还包括以下步骤：布置测网。
进一步，步骤(3)中所述删除剩余极大值中距离小于窗口距离的极大值的方法为：以最高峰为起点，分别向两端以固定窗口距离进行滑动，仅保留窗口距离内最大的极大值。
进一步，所述步骤(3)中窗口距离的宽度为2km。
进一步，所述步骤(5)中高差阈值为0.01米。
进一步，所述步骤(7)中还包括：对形成的三维航迹曲线进行整体抬升，使航迹曲线高出最高峰100米。
本发明一种中高山区航空重力测量三维航迹规划方法，根据航空重力仪对飞行起伏坡度和起伏周期的要求，以测区地形数据为基础，求取满足一定要求的极大值和求取测线与切割线交叉点处均值相结合的方式获得控制点，然后对控制点进行线性插值并整体抬升一定高度，得到三维航迹曲线。本发明有效降低了平均飞行高度，且兼顾测线和切割线上的地形差异，保证了测线和切割线交叉点位置高度在理论上重合，从而确保了交叉点上的重力场差值尽可能小，利于全区重力场调平。同时，本发明还可将计算得到的输出控制点提供给导航设备，从而在航空重力缓起伏飞行中突破了传统目视导航方法，实现了基于导航设备的精确导航，不仅保证了人身安全，还可以确保飞行高度和起伏坡度满足重力测量需求，减少了人为因素造成的飞机俯仰角过大而导致的虚假重力异常问题。
附图说明
图1是实施例1中一种中高山区航空重力测量三维航迹规划方法的整体流程示意图；
图2是实施例2中根据实施例1所述方法规划的一条初始航迹曲线。
具体实施方式
以下结合附图1和2，进一步说明本发明一种中高山区航空重力测量三维航迹规划方法的具体实施方式。本发明一种中高山区航空重力测量三维航迹规划方法不限于以下实施例的描述。
实施例1：
本实施例给出中一种高山区航空重力测量三维航迹规划方法，如图1所示，包括以下步骤：
1、收集测区地形数据，根据高山区航空重力测量需求，布置测网(包括测线和切割线)，获取测线和切割线的水平投影坐标。
2、求取每条测线和切割线上对应地形的最高峰、极大值。所述最高峰仅有一个，所述极大值有多个，所述最高峰、极大值均为三维坐标值。
3、以最高峰为参考点，删除测线和切割线上坡度大于3°的极大值。具体的，以最高峰和一个极大值的三维坐标建立直线，若该直线的斜率大于3°，则删除该极大值。该步骤的目的是，确保最终规划的航迹中，其飞行坡度始终不大于3°，以确保飞行过程中的起伏坡度始终满足重力测量需求。
4、删除剩余极大值中间距小于窗口距离的极大值；剩余的极大值作为测线和切割线的一般控制点。具体方法为：以最高峰为起点，以2km宽度的窗口距离，分别沿测线或切割线向两端以固定窗口距离进行滑动，仅保留窗口距离内最大的极大值。该步骤的目的是，尽量忽略或减少飞行中的连续起伏变化，确保飞行过程中的起伏周期始终满足重力测量需求。
5、对一般控制点进行线性插值，形成测线和切割线的初始航迹曲线。为获得更好地飞行效果，本步骤也可采用其他差值拟合方法。
6、求取测线和切割线交叉点位置高差，当高差均小于阈值(可设定为0.01米)时，所述一般控制点即为输出控制点，转步骤8。当高差均小于阈值时，说明测线和切割线交叉点位置基本重合，满足将来对测线和切割线采集的测量数据的重力场差值要求，在进行全区重力场调平时可以将误差控制在允许范围内；反之，则说明测线和切割线交叉点之间的高差过大，不能满足测量要求，还需要对一般控制点进行进一步调整。
7、当高差超出阈值时，求取交叉点位置航迹均值作为测线和切割线的特殊控制点，然后删除测线和切割线上相邻控制点坡度大于3°的一般控制点；对测线和切割线的特殊控制点和未删除的一般控制点进行线性插值，更新测线和切割线航迹，再次求取测线和切割线交叉点位置高差，直至高差均小于阈值；特殊控制点和未删除的一般控制点即为输出控制点。本步骤的目的在于，将测线和切割线交叉点位置作为特殊控制点，并作为最终输出控制点的重要参照，对一般控制点进行调整，最终达到交叉点位置高程一致、满足测量需求的目标。
8、对输出控制点进行插值，形成三维航迹曲线。对形成的三维航迹曲线进行整体抬升，例如使航迹曲线高出最高峰100米，以确保飞行安全。
实施例2：
本实施例给出根据实施例1所述方法具体规划出的一条初始航迹曲线。
如图2所示，横坐标为测线或切割线延伸方向的投影距离，纵坐标为高度。第一实线101为地形轮廓线，倒三角形标记1-73为该地形轮廓线上的极大值，其中63为最高峰。根据实施例1所述方法，首先删除坡度大于3°的极大值，剩余的极大值用圆形图案进行标记；然后删除剩余极大值中间距小于窗口距离的极大值，剩下的极大值用*形图案进行标记，即为一般控制点；最后，对一般控制点进行线性插值，得到的第二实线102即为根据一般控制点规划出的初始航迹曲线。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。
Claims (6)
1.一种中高山区航空重力测量三维航迹规划方法，其特征在于：包括以下步骤：
(1)求取测线和切割线上对应地形的最高峰、极大值；
(2)以最高峰为参考点，删除测线和切割线上坡度大于3°的极大值；
(3)删除剩余极大值中间距小于窗口距离的极大值；剩余的极大值作为测线和切割线的一般控制点；
(4)对一般控制点进行线性插值，形成测线和切割线的初始航迹曲线；
(5)求取测线和切割线交叉点位置高差，当高差均小于阈值时，所述一般控制点即为输出控制点，转步骤(7)；
(6)当高差超出阈值时，求取交叉点位置航迹均值作为测线和切割线的特殊控制点，然后删除测线和切割线上相邻控制点坡度大于3°的一般控制点；对测线和切割线的特殊控制点和未删除的一般控制点进行线性插值，更新测线和切割线航迹，再次求取测线和切割线交叉点位置高差，直至高差均小于阈值；特殊控制点和未删除的一般控制点即为输出控制点；
(7)对输出控制点进行线性插值并整体抬升一定高度，得到三维航迹曲线。
2.根据权利要求1所述的中高山区航空重力测量三维航迹规划方法，其特征在于：所述步骤(1)之前还包括以下步骤：
布置测网。
3.根据权利要求2所述的中高山区航空重力测量三维航迹规划方法，其特征在于：步骤(3)中所述删除剩余极大值中距离小于窗口距离的极大值的方法为：以最高峰为起点，分别向两端以固定窗口距离进行滑动，仅保留窗口距离内最大的极大值。
4.根据权利要求3所述的中高山区航空重力测量三维航迹规划方法，其特征在于：所述步骤(3)中窗口距离的宽度为2km。
5.根据权利要求4所述的中高山区航空重力测量三维航迹规划方法，其特征在于：所述步骤(5)中高差阈值为0.01米。
6.根据权利要求5所述的中高山区航空重力测量三维航迹规划方法，其特征在于：所述步骤(7)中还包括：对形成的三维航迹曲线进行整体抬升，使航迹曲线高出最高峰100米。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810256318.4A CN108507576B (zh) | 2018-03-27 | 2018-03-27 | 一种中高山区航空重力测量三维航迹规划方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810256318.4A CN108507576B (zh) | 2018-03-27 | 2018-03-27 | 一种中高山区航空重力测量三维航迹规划方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108507576A CN108507576A (zh) | 2018-09-07 |
| CN108507576B true CN108507576B (zh) | 2020-06-30 |
Family
ID=63378677
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810256318.4A Expired - Fee Related CN108507576B (zh) | 2018-03-27 | 2018-03-27 | 一种中高山区航空重力测量三维航迹规划方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108507576B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109471191B (zh) * | 2018-11-13 | 2021-07-06 | 中国地震局地震研究所 | 一种航空重力测量的方法及设备 |
| CN113311505B (zh) * | 2021-04-13 | 2021-10-22 | 中国自然资源航空物探遥感中心 | 航空重力仪器飞行高度质量评估方法及装置 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5402340A (en) * | 1994-05-02 | 1995-03-28 | Textron, Inc. | Gradiometer based terrain estimation |
| CN101231167A (zh) * | 2008-02-20 | 2008-07-30 | 刘雁春 | 海洋测线网系统误差的检测与调整方法 |
| CN101692315B (zh) * | 2009-09-25 | 2011-08-10 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
| CN104122597B (zh) * | 2013-08-22 | 2016-09-14 | 中国科学院遥感与数字地球研究所 | 一种无人机航磁探测系统及方法 |
| CN105068131A (zh) * | 2015-08-03 | 2015-11-18 | 中国科学院电子学研究所 | 一种航磁数据调平方法 |
| CN105571588A (zh) * | 2016-03-10 | 2016-05-11 | 赛度科技（北京）有限责任公司 | 一种无人机三维空中航路地图构建及其航路显示方法 |
| CN106226828B (zh) * | 2016-07-18 | 2017-11-24 | 吉林大学 | 时间域航空电磁数据约束多项式拟合调平方法 |
| CN106970426B (zh) * | 2017-03-31 | 2018-02-13 | 吉林大学 | 一种基于测线差分与主成分分析的航空电磁数据调平方法 |
-
2018
- 2018-03-27 CN CN201810256318.4A patent/CN108507576B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN108507576A (zh) | 2018-09-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111927552B (zh) | 基于天空地勘察技术的复杂山区长大深埋隧道勘察方法 | |
| CN106093963B (zh) | 一种提高铁路车载激光雷达扫描数据精度的方法 | |
| CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
| CN106087621A (zh) | 一种基于车载激光雷达技术的铁路既有线复测方法 | |
| CN103020966B (zh) | 一种基于建筑物轮廓约束的航空与地面LiDAR数据自动配准方法 | |
| CN102494667B (zh) | 一种表征地面沉降的方法 | |
| CN102445186A (zh) | 一种激光雷达扫描生成公路设计地表信息方法 | |
| CN104807440B (zh) | 一种用于机场的全自动高频率地表沉降监测方法 | |
| CN113189615A (zh) | 一种使用垂直起降固定翼无人机对输电线路进行巡检的方法 | |
| CN109508508A (zh) | 一种露天矿山治理勘查设计方法 | |
| CN110057367A (zh) | 一种无人机的航线规划的方法及装置 | |
| CN108507576B (zh) | 一种中高山区航空重力测量三维航迹规划方法 | |
| CN107132578A (zh) | 一种微地震地面监测速度模型校正算法 | |
| CN105067004A (zh) | 一种基于地形的路径搜索方法 | |
| CN111854692A (zh) | 一种在道路测设中的无人机影像匹配点云的测量方法 | |
| CN111045091A (zh) | 一种玄武岩覆盖下古河道的识别定位方法 | |
| CN113742437B (zh) | 地图更新方法、装置、电子设备和存储介质 | |
| CN114897777A (zh) | 顾及空间关系的接触网支撑设施激光点云全自动提取方法 | |
| CN117607893A (zh) | 基于无人机非接触测量的铁路既有线复测方法 | |
| CN108596397A (zh) | 一种基于水文模型多机航线规划方法 | |
| CN109242247A (zh) | 一种煤矿区地面塌陷破坏程度评价方法 | |
| CN107886573A (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
| CN116794651A (zh) | 基于航测探地雷达与影像数据的滑坡特征解译方法及系统 | |
| CN103543744A (zh) | 一种无人驾驶飞艇航空磁测飞行航线布置的方法 | |
| CN111897015A (zh) | 基于航空电磁法的深大断裂三维探测方法 |
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 | ||
| CB02 | Change of applicant information |
Address after: 100083 No. 31, Haidian District, Beijing, Xueyuan Road Applicant after: CHINA AERO GEOPHYSICAL SURVEY AND REMOTE SENSING CENTER FOR NATURAL RESOURCES Address before: 100083 No. 31, Haidian District, Beijing, Xueyuan Road Applicant before: CHINA AERO GEOPHYSICAL SURVEY & REMOTE SENSING CENTER FOR LAND AND RESOURCES |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200630 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |