CN115164944B - 一种多波束测深异常的角度和传播时间相关改正方法 - Google Patents
一种多波束测深异常的角度和传播时间相关改正方法 Download PDFInfo
- Publication number
- CN115164944B CN115164944B CN202211099138.2A CN202211099138A CN115164944B CN 115164944 B CN115164944 B CN 115164944B CN 202211099138 A CN202211099138 A CN 202211099138A CN 115164944 B CN115164944 B CN 115164944B
- Authority
- CN
- China
- Prior art keywords
- sounding
- depth
- deviation
- angle
- propagation time
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 43
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000010410 layer Substances 0.000 claims description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 18
- 238000010276 construction Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000002344 surface layer Substances 0.000 claims description 3
- 238000012876 topography Methods 0.000 description 6
- 230000002159 abnormal effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 230000003313 weakening effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52004—Means for monitoring or calibrating
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Manufacturing & Machinery (AREA)
- Computer Networks & Wireless Communication (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明涉及水文测量领域,尤其涉及一种多波束测深异常的角度和传播时间相关改正方法。
背景技术
水下地形是最基础地理要素,在河床演变、水域治理、水资源利用等邻域具有重要的应用价值。水下地形常借助单波束测深系统或多波束测深系统通过水下地形测量方式来获得,相较单波束测深系统,多波束测深系统一次可以获得100~500个测深点,可实现条带式全覆盖扫测,显著提高了测深效率,以成为当前水下地形获取的主要手段。
多波束测深系统是一个由多元传感器组成的综合系统,姿态传感器安装偏差、换能器安装偏差、罗经校准偏差、声速误差等均会给最终的测深结果带来显著的系统性误差,这种误差随着波束入射角和深度的增大而增大,在测深数据的断面地形中呈现边缘波束测深异常,导致相邻测线公共覆盖区测深数据不一致,进而引起地形的不合理突变。针对该问题,现有的解决方案主要有两种方法,即强制压制法和趋势面拟合法。强制压制法是根据相邻测线公共覆盖区测深数据,寻求公共位置的测深数据,取分别属于两个测线的测点深度的平均值作为该位置最终的深度,通过这种强制压制的方式,实现相邻测线公共位置测深数据的一致。强制压制法实施简单,但缺少原理,虽然实现了相邻测线公共覆盖区测深数据的一致,但不能正确的反映真实的地形变化。趋势面拟合法利用多波束的测深特点,即中央波束测深数据受诸因素影响较小,边缘波束受影响较大,利用相邻测线的中央波束测深数据,构建地形趋势面,以该趋势面为参考,内插边缘波束测量的地形趋势,将相邻测线的边缘波束压制到该参考面上,进而实现诸因素对边缘波束造成的测深异常改正。趋势面法较强制压制法而言,顾及了诸因素对多波束测深影响的机理,但仅适用于地形变化相对平缓的水下地形,而对于地形变化比较复杂的水下地形,则基于相邻测线中央波束构建的地形趋势则不能真实的反映边缘波束扫测位置的水下地形趋势,由此实施的边缘波束测深异常压制则会造成虚假地形。
发明内容
本发明的目的在于针对上述现有技术的不足,提供了一种多波束测深异常的角度和传播时间相关改正方法,用于获取高精度的多波束测深数据和水下地形,从机理上解决多波束测深带来系统性误差的影响。
为实现上述目的,本发明采用了如下技术方案:
本发明提供了一种多波束测深异常的角度和传播时间相关改正方法,包括以下步骤:
S1、多波束测深偏差计算:多波束测深结果与单波束测深结果在同名点位置的深度偏差;
S2、多波束测深偏差序列构建:根据所述同名点位置上的深度偏差、多波速侧深点对应的波速角和传播时间,构建所述深度偏差与波束角和传播时间的数据序列;
S5、多波束测深偏差修正:利用所述S4中的多波束测深数据改正模型,根据各个波束的传播时间和波束入射角,计算改正量,修正原始测深数据。
进一步,所述S1中,多波束测深偏差计算的模型为:
进一步,所述S2中,数据序列为:
进一步,所述S3中,构建自变量表达式包括以下步骤:
S301、多波束测深计算模型为:
S302、根据Snell法则,下列关系成立:
S305、简化公式(4)和公式(5):
进一步,所述S4中,多波束的测深误差模型用多项式关系模型表达为:
根据所述S2中的数据序列,
公式(10)的矩阵形式为:
其中,
进一步,根据构建公式(9)的角度和深度相关的多项式关系模型,对所有的多波束测深数据,根据波束入射角和测量深度,代入公式(9)中,计算改正量,利用,对边缘波束,即波束角大于45度的多波束测深点进行改正,进而消除多波束测深偏差的影响:
本发明的有益效果为:通过研究发现了多波束边缘波束测深异常与两个因素相关,即波束入射角和传播时间,以相同位置的单波束测深系统测深结果为参考,寻求出多波束扫测断面地形在不同入射角和传播时间下的偏差量,构建与波束入射角度、传播时间相关的测深改正模型,实现多波束边缘波束测深异常的削弱。该方法从机理上彻底解决多波束边缘波束测深异常的削弱难题,无论边缘波束地形变化平缓还是复杂,均能实现测深异常的削弱和真实地形的获取。
附图说明
图1 为本发明一种多波束测深异常的角度和传播时间相关改正方法的流程图;
图2为同名点的测深偏差示意图;
图3为测深偏差与入射角的示意图;
图4为测深偏差与传播时间的关系图;
图6为处理前的断面图;
图7为处理后的断面图;
图8为原始测深数据误差分布图;
图9为处理后的误差分布图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,下面结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
一种多波束测深异常的角度和传播时间相关改正方法,包括以下步骤:
S1、多波束测深偏差计算:多波束测深结果与单波束测深结果在同名点位置的深度偏差;
S2、多波束测深偏差序列构建:根据所述同名点位置上的深度偏差、多波速侧深点对应的波速角和传播时间,构建所述深度偏差与波束角和传播时间的数据序列;
S5、多波束测深偏差修正:利用所述S4中的多波束测深数据改正模型,根据各个波束的传播时间和波束入射角,计算改正量,修正原始测深数据。
所述S1中,多波束测深偏差计算的模型为:
所述S2中,数据序列为:
所述S3中,构建自变量表达式包括以下步骤:
S301、多波束测深计算模型为:
S302、根据Snell法则,下列关系成立:
S305、简化公式(4)和公式(5):
所述S4中,多波束的测深误差模型用多项式关系模型表达为:
根据所述S2中的数据序列,
公式(10)的矩阵形式为:
其中,
根据构建公式(9)的角度和深度相关的多项式关系模型,对所有的多波束测深数据,根据波束入射角和测量深度,代入公式(9)中,计算改正量,利用,对边缘波束,即波束角大于45度的多波束测深点进行改正,进而消除多波束测深偏差的影响:
实施例一
将该技术方案应用于南海某水域,该水域水深范围104m~110m。对该水域施测了多波束测深,为消除多波束边缘波束异常测深的影响,在该区域的中央部分施测了单波束水下地形测量;
S1,根据S1中式(1)给出的方法,寻找单波束与多波束相同位置的测深点(同名的),计算多波束测深结果与单波束测深结果在同名点位置的深度偏差如图2所示。
S2、根据S2,分别构建多波束测深偏差与波束入射角序列、测深偏差与波束传播时间数据序列。两个序列的分布图如图3和图4所示。
利用最小二乘原理,解算方程组,即式(11),得到多波束测深修正模型的多项式模型系数,建立多波束测深数据改正模型。
解算得到的模型系数如表1所示。
表1 回归模型的参数
S5、根据S5所述方法和式(12),利用S4中建立的多波束测深数据改正模型,根据多波束各个波束的传播时间和波束入射角,计算改正量,修正原始多波束测深数据,得到修正后的测深数据。图6为原始多波束测深数据形成的地形断面,图7为修正后测深数据形成的地形断面。
图5直观的给出了修正前、后地形断面,结果表明利用建立的多波束测深异常修正模型对原始多波束测深数据修正,实现了相邻多波束测线(L1、L2、L3)测深点对地形的一致性表达,也表明S1~S5给出的多波束测深异常修正方法的正确性。
图5仅用于直观视觉定性分析,为定性分析,对全区域改正前、后的同名点测深不符值进行统计,结果如图8、图9和表2所示。
表2 多波束测深数据修正前后精度统计结果
图8、图9和表2表明,在整个区域,利用本专利给出的S1~S5方法和过程,实现了多波束异常测深结果的修正,在不同多波束测量线的同名测深点上,经过异常测深修正,测深精度提供了2倍多,进一步表明了本专利方法的正确性和实用性。
以上所述实施例仅表达了本发明的实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求。
Claims (2)
1.一种多波束测深异常的角度和传播时间相关改正方法,其特征在于,包括以下步骤:
S1、多波束测深偏差计算:多波束测深结果与单波束测深结果在同名点位置的深度偏差;
S2、多波束测深偏差序列构建:根据所述同名点位置上的深度偏差、多波速侧深点对应的波速角和传播时间,构建所述深度偏差与波束角和传播时间的数据序列;
S5、多波束测深偏差修正:利用所述S4中的多波束测深数据改正模型,根据各个波束的传播时间和波束入射角,计算改正量,修正原始测深数据;
所述S1中,多波束测深偏差计算的模型为:
所述S2中,数据序列为:
所述S3中,构建自变量表达式包括以下步骤:
S301、多波束测深计算模型为:
S302、根据Snell法则,下列关系成立:
S305、简化公式(4)和公式(5):
所述S4中,多波束的测深误差模型用多项式关系模型表达为:
根据所述S2中的数据序列,
公式(10)的矩阵形式为:
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211099138.2A CN115164944B (zh) | 2022-09-09 | 2022-09-09 | 一种多波束测深异常的角度和传播时间相关改正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211099138.2A CN115164944B (zh) | 2022-09-09 | 2022-09-09 | 一种多波束测深异常的角度和传播时间相关改正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115164944A CN115164944A (zh) | 2022-10-11 |
CN115164944B true CN115164944B (zh) | 2022-12-02 |
Family
ID=83482311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211099138.2A Active CN115164944B (zh) | 2022-09-09 | 2022-09-09 | 一种多波束测深异常的角度和传播时间相关改正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115164944B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115930916B (zh) * | 2022-12-06 | 2023-06-06 | 长江水利委员会水文局 | 一种船体涌动频谱相关性的多波束条纹测深影响消除方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269811A (zh) * | 2010-06-04 | 2011-12-07 | 国家海洋局第一海洋研究所 | 基于垂直波束水深的边缘波束水深数据声速改正方法 |
CN104613906A (zh) * | 2015-02-06 | 2015-05-13 | 长江水利委员会水文局 | 基于声线跟踪的库区深水水深测量方法 |
CN105258684A (zh) * | 2015-10-30 | 2016-01-20 | 山东科技大学 | 基于激光点云为约束的多波束低掠射角波束归位方法 |
CN111079080A (zh) * | 2019-12-17 | 2020-04-28 | 中国人民解放军海军大连舰艇学院 | 基于卡尔曼滤波的多波束数据声速整体误差改正方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103344954B (zh) * | 2013-07-08 | 2014-07-02 | 国家海洋局第二海洋研究所 | 一种基于多来源水深数据融合的海底地形地貌构建方法 |
US10008364B2 (en) * | 2015-02-27 | 2018-06-26 | Kla-Tencor Corporation | Alignment of multi-beam patterning tool |
CN110146895B (zh) * | 2019-05-16 | 2021-04-20 | 浙江大学 | 基于倒置式多波束回声仪的声速剖面反演方法 |
-
2022
- 2022-09-09 CN CN202211099138.2A patent/CN115164944B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102269811A (zh) * | 2010-06-04 | 2011-12-07 | 国家海洋局第一海洋研究所 | 基于垂直波束水深的边缘波束水深数据声速改正方法 |
CN104613906A (zh) * | 2015-02-06 | 2015-05-13 | 长江水利委员会水文局 | 基于声线跟踪的库区深水水深测量方法 |
CN105258684A (zh) * | 2015-10-30 | 2016-01-20 | 山东科技大学 | 基于激光点云为约束的多波束低掠射角波束归位方法 |
CN111079080A (zh) * | 2019-12-17 | 2020-04-28 | 中国人民解放军海军大连舰艇学院 | 基于卡尔曼滤波的多波束数据声速整体误差改正方法 |
Non-Patent Citations (3)
Title |
---|
Multi-beam Data Automatic Filtering Technology;Yu Yan 等;《2022 3rd International Conference on Geology, Mapping and Remote Sensing (ICGMRS)》;20220424;第1-4页 * |
多波束与单波束测深数据的融合处理技术;黄谟涛等;《测绘学报》;20011130;第30卷(第04期);第299-303页 * |
多波束与单波束测深数据融合处理方法;罗君;《海洋测绘》;20181231;第38卷(第006期);第21-24页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115164944A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zabotin et al. | NeXtYZ: Three-dimensional electron density inversion for dynasonde ionograms | |
CN109344865B (zh) | 一种多数据源的数据融合方法 | |
CN115164944B (zh) | 一种多波束测深异常的角度和传播时间相关改正方法 | |
CN110031877B (zh) | 一种基于grnn模型的区域nwp对流层延迟改正方法 | |
CN106872963A (zh) | 一种多组多线激光雷达的自动标定算法 | |
CN110081864B (zh) | 一种顾及水深值的水深测量综合延迟改正方法 | |
US20200225112A1 (en) | Four-dimensional imaging method for structural damage based on time-invariant characteristic signal | |
CN104459696A (zh) | 一种基于平地相位的sar干涉基线精确估计方法 | |
CN111008355B (zh) | 一种基于信任传播的气象地面要素插值方法 | |
CN101216555A (zh) | Rpc模型参数提取方法和几何纠正方法 | |
CN107679361A (zh) | 一种测风数据处理方法 | |
CN112381337B (zh) | 一种多源气象数据融合处理方法、系统、终端及介质 | |
CN106709460A (zh) | 一种高频底座力天平的动力校准方法 | |
CN115687983A (zh) | 一种桥梁健康状态监测方法、系统及电子设备 | |
CN109084708A (zh) | 颗粒二维轮廓表面整体粗糙度计算方法 | |
CN109359360B (zh) | 一种基于局部特征的结构应力处理方法 | |
CN115690106A (zh) | 一种基于计算机视觉的深埋锚具封闭检测方法 | |
Zhao et al. | From 10 m to 11000 m, automatic processing multi-beam bathymetric data based on PGO method | |
CN107831516A (zh) | 融合gnss和地面监测网获取大坝实时高精度位移的方法 | |
CN111505619A (zh) | 一种时空分布不规则不均匀卫星高度计测高数据网格化处理方法 | |
CN109269453B (zh) | Pbf增材制造高频超声在线检测单次检测层数确定方法 | |
Regier et al. | Observations of the power and directional spectrum of ocean surface waves | |
CN113723541A (zh) | 一种基于混合智能算法的边坡位移预测方法 | |
CN113111529A (zh) | 融合数值模式和卫星微波云反演数据的红外亮温模拟方法 | |
CN110348040B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |