CN108256049B - 一种不同区块多波束数据接边处理方法 - Google Patents
一种不同区块多波束数据接边处理方法 Download PDFInfo
- Publication number
- CN108256049B CN108256049B CN201810031754.1A CN201810031754A CN108256049B CN 108256049 B CN108256049 B CN 108256049B CN 201810031754 A CN201810031754 A CN 201810031754A CN 108256049 B CN108256049 B CN 108256049B
- Authority
- CN
- China
- Prior art keywords
- block
- water
- depth
- edge fit
- data
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control For Baths (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Bathtub Accessories (AREA)
Abstract
一种不同区块多波束数据接边处理方法,工作流程和步骤为:(1)将区块A和区块B的数据分别进行网格化处理;(2)在步骤(1)重叠区中间分别提取区块A和区块B同一位置剖面的水深值;(3)以步骤(2)提取的剖面为界,删除重叠在对方区域的数据;(4)计算上述提取的两个剖面的同点水深差值h;(5)结合区域水深观测值Z0,利用趋势面计算处理,得出区块水深校正值Z;通过本发明对不同区块数据进行接边处理,可提高多波束测深数据的精度,能够很好地解决多波束数据在接边时地形趋势不协调、等深线走向异常的问题。
Description
技术领域
本发明涉及一种不同区块多波束数据接边处理方法。
背景技术
近年来,我国海洋环境、生态和资源都发生了重大变化,海洋开发活动不断增加,以及海洋调查、海洋工程建设、海洋科学研究的发展,多波束测深系统的应用备受关注,在海洋基础测绘中更是广泛应用。多波束测深系统已经成为目前海底地形地貌测量的主要设备,其具有全覆盖、高精度、高效率、高分辨率的特点,且能够全面真实准确的反映海底地形地貌。
在大规模的海底地形地貌调查项目中,需要将调查区域进行分区块,再由不同的单位分区块完成。最终的成果集成时,因为不同区块的数据由不同单位的不同设备测量提供,在测量精度和分辨率方面存在差异,以及系统误差的存在,导致在接边时会出现同一个位置的水深不一致,等深线出现拐角的问题,从而要对不同区块的多波束测深数据进行接边处理。而目前处理这类问题常用拟合的方法来融合多区块的数据,该方法主要包括有移动平均法、反距离加权法和最小曲率法等,其主要是依靠数学加权的方式对不同区块重叠区数据进行压制平滑,处理后的地形仍然会表现出趋势不协调,等深线走向不自然,人为斧凿修改的痕迹明显。
发明内容
本发明的目的在于提供一种不同区块多波束数据接边处理方法,以解决现有的多波束数据接边时仍然存在的地形趋势不协调、等深线走向异常的问题。
为了解决上述的技术问题,本发明采用的技术方案如下:
为方便表达本发明接边处理的方法,以下都以两个区块(即区块A和区块B)的接边来说明,且这两个区块数据都符合接边标准。
所述一种不同区块多波束数据接边处理方法,工作流程和步骤为:
(1)将区块A和区块B的数据分别进行网格化处理;
(2)在步骤(1)重叠区中间分别提取区块A和区块B同一位置剖面的水深值;
(3)以步骤(2)提取的剖面为界,删除重叠在对方区域的数据;
(4)计算上述提取的两个剖面的同点水深差值h;
(5)结合区域水深观测值Z0,利用趋势面计算处理,得出区块水深校正值Z;
所述步骤(4)计算两个剖面的同点水深差值h的具体方法如下:
定义区块A、B中提取的同位置的剖面水深值分别为hA、hB,则同点水深差值h=hA-hB;
所述步骤(5)利用趋势面计算处理,得出区块水深校正值Z的具体方法如下:
将所述的接边的两个区块沿着与所述步骤(2)提取的剖面垂直的方向,进行条带化,使得每个区块划分成若干细小的条带,然后对每个条带中的水深值进行校正计算,条带划分的越细小,精度相对越高;所述的水深校正值Z表示为:Z=Z0+f(x,h/2);
其中x是校正点到提取的剖面的距离;h/2表示相接边的两个细小条带在剖面处均分的水深差值;f(x,h/2)表示沿地形趋势面变化,校正点按距离远近吸收的水深差值。
本发明的有益效果是:它提供了一种不同区块多波束数据接边处理方法,通过该方法对不同区块数据进行接边处理,可提高多波束测深数据的精度,能够很好地解决多波束数据在接边时地形趋势不协调、等深线走向异常的问题。
附图说明
图1为本发明方法实施的流程图。
图2为多波束测深数据分布图。
图3为多波束数据网格化后的地形及等深线图。
图4为在重叠区分别提取同一位置的剖面示意图。
图5为删除处理后的多波束测深数据分布图。
图6为删除处理后的多波束测深数据网格化后的地形及等深线图。
图7为提取的两条剖面水深差值图。
图8为趋势面接边示意图。
图9为经本发明接边方法处理后的地形及等深线图。
具体实施方式
为了使本发明的目的、技术方案以及优点显示的更加清楚明白,下面结合附图1-9对本发明做进一步的详细说明。
如图1所示,一种不同区块多波束数据接边处理方法,本发明的方法流程主要包括以下几个步骤:(1)将区块A和区块B的数据分别进行网格化处理;(2)在重叠区中间分别提取区块A和区块B同一位置剖面的水深值;(3)以上述提取的剖面为界,删除重叠在对方区域的数据;(4)计算上述提取的两个剖面的同点水深差值h;(5)结合区域水深观测值Z0,利用趋势面计算处理,得出区块水深校正值Z。
以长江口海域两个区块A和B的多波束测深数据为例,对不同区块多波束数据接边处理方法进行实证分析,图2为两个不同区块A和B的多波束测深数据分布图,这些数据都是经过预处理后符合接边标准的,从图2中可以很明显的看出A、B两个区块相邻分布且有一定的重叠区。
在图3中,对应本发明所述工作流程中的步骤(1),是对两个区块A和B多波束数据利用传统拟合方法进行网格化处理形成的地形图,图3中可以清晰的看出传统方法融合的两个区块重叠处地形趋势不协调、不光滑,甚至出现地形突变的现象。
在图4中,对应本发明所述工作流程中的步骤(2),选取两个区块A和B的重叠区中间位置作一条剖面,分别在这一剖面上提取两个区块接边区同一位置的水深值。
在图5和6中,对应本发明所述工作流程中的步骤(3),是以步骤(2)中提取的剖面为界,删除重叠在对方区域的数据后形成的数据分布图以及地形图,可以更加明显的看出两个区块A和B的拼接处水深不一致,两区域接边位置等深线错位。
在图7中,对应本发明所述工作流程中的步骤(4),是根据步骤(3)中提取两个区块A和B接边区同一位置的水深值,相减计算出两个区块接边位置的水深相差值h。
在图8中,对应本发明所述工作流程中的步骤(5),是结合区域水深观测值Z0,利用趋势面计算处理,得出区块水深校正值Z。根据本发明提出的方法,将所述的接边的两个区块A和B沿着与所述工作流程中第(2)步提取的剖面垂直的方向,对两个区块数据进行条带化划分,使得每个区块划分成若干个细小的条带,从而每个细小条带与剖面垂直相交位置所要吸收的误差值即为该位置两个区块接边位置的水深相差值的一半(h/2),然后对每个条带中的水深值进行校正计算,条带划分的越细小,精度相对越高;所述的水深校正值Z表示为:Z=Z0+f(x,h/2),其中x是校正点到提取的剖面的距离;h/2表示相接边的两个细小条带在剖面处均分的水深差值;f(x,h/2)表示沿地形趋势面变化,校正点按距离远近吸收的水深差值。
在图9中,是经过本发明所述的接边处理方法处理后形成的地形图,与图6相比地形趋势更加平滑,等深线走向更加合理。
综上所述,本发明是适用于不同区块多波束数据接边处理的,通过实例验证,解决了传统的利用拟合的方法融合多波束数据后,依然存在的地形趋势不协调,等深线走向异常的问题,提高了多波束测深数据的精度。
Claims (1)
1.一种不同区块多波束数据接边处理方法,其特征是:工作流程和步骤为,
(1)将区块A和区块B的数据分别进行网格化处理;
(2)在步骤(1)网格化重叠区中间分别提取区块A和区块B同一位置剖面的水深值;
(3)以步骤(2)提取的剖面为界,删除重叠在对方区域的数据;
(4)计算上述提取的两个剖面的同点水深差值h;
(5)结合区域水深观测值Z0,利用趋势面计算处理,得出区块水深校正值Z;
所述步骤(4)计算两个剖面的同点水深差值h的具体方法如下:
定义区块A、B中提取的同位置的剖面水深值分别为hA、hB,则同点水深差值h=hA-hB;
所述步骤(5)利用趋势面计算处理,得出区块水深校正值Z的具体方法如下:
将所述的接边的两个区块沿着与所述步骤(2)提取的剖面垂直的方向,进行条带化,使得每个区块划分成若干细小的条带,然后对每个条带中的水深值进行校正计算,条带划分的越细小,精度相对越高;所述的水深校正值Z表示为:Z=Z0+f(x,h/2);
其中x是校正点到提取的剖面的距离;h/2表示相接边的两个细小条带在剖面处均分的水深差值;f(x,h/2)表示沿地形趋势面变化,校正点按距离远近吸收的水深差值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810031754.1A CN108256049B (zh) | 2018-01-12 | 2018-01-12 | 一种不同区块多波束数据接边处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810031754.1A CN108256049B (zh) | 2018-01-12 | 2018-01-12 | 一种不同区块多波束数据接边处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108256049A CN108256049A (zh) | 2018-07-06 |
CN108256049B true CN108256049B (zh) | 2019-10-11 |
Family
ID=62726398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810031754.1A Active CN108256049B (zh) | 2018-01-12 | 2018-01-12 | 一种不同区块多波束数据接边处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108256049B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100898617B1 (ko) * | 2008-11-20 | 2009-05-27 | 주식회사 범아엔지니어링 | LiDAR 측량값 및 MBES 측량값의 불규칙삼각망(TIN) 보간 성과 검증을 통한 육상부와 수심부가 함께 공존하는 지역의 수치표고모델을 구축하는 방법 |
CN102927974A (zh) * | 2012-10-31 | 2013-02-13 | 山东科技大学 | 多波束测深系统检测方法 |
CN104949692A (zh) * | 2015-06-25 | 2015-09-30 | 江苏中海达海洋信息技术有限公司 | 多波束测深系统安装校正值的计算方法 |
CN106060493A (zh) * | 2016-07-07 | 2016-10-26 | 广东技术师范学院 | 多源投影无缝边缘拼接方法及系统 |
-
2018
- 2018-01-12 CN CN201810031754.1A patent/CN108256049B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100898617B1 (ko) * | 2008-11-20 | 2009-05-27 | 주식회사 범아엔지니어링 | LiDAR 측량값 및 MBES 측량값의 불규칙삼각망(TIN) 보간 성과 검증을 통한 육상부와 수심부가 함께 공존하는 지역의 수치표고모델을 구축하는 방법 |
CN102927974A (zh) * | 2012-10-31 | 2013-02-13 | 山东科技大学 | 多波束测深系统检测方法 |
CN104949692A (zh) * | 2015-06-25 | 2015-09-30 | 江苏中海达海洋信息技术有限公司 | 多波束测深系统安装校正值的计算方法 |
CN106060493A (zh) * | 2016-07-07 | 2016-10-26 | 广东技术师范学院 | 多源投影无缝边缘拼接方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108256049A (zh) | 2018-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Diegel et al. | Cenozoic structural evolution and tectono-stratigraphic framework of the northern Gulf Coast continental margin | |
Bird et al. | Gulf of Mexico tectonic history: Hotspot tracks, crustal boundaries, and early salt distribution | |
Kauahikaua et al. | Deep magmatic structures of Hawaiian volcanoes, imaged by three-dimensional gravity models | |
Moore et al. | The evolution and ages of Makgadikgadi palaeo-lakes: consilient evidence from Kalahari drainage evolution south-central Africa | |
McNeill et al. | Tectonics of the Neogene Cascadia forearc basin: Investigations of a deformed late Miocene unconformity | |
Gao et al. | Fine crustal structure beneath the junction of the southwest Tian Shan and Tarim Basin, NW China | |
CN103472483B (zh) | 基于真地表或浮动基准面的速度建模方法 | |
CN110838175B (zh) | 注气开发油藏的地质模型建立方法 | |
CN109919070B (zh) | 一种剖面形态自适应拟合的海岸线遥感推算方法 | |
CN107101617B (zh) | 一种基于遥感监测的冲沟地形下采场覆岩活动控制方法 | |
CN110501743B (zh) | 一种储层构型约束的砂岩厚度等值线图编制方法 | |
CN106097448A (zh) | 一种多特征约束下的盐腔围岩地质三维建模方法 | |
CN106569260B (zh) | 针对高陡构造地区二维地震变速空校成图方法 | |
CN111681316B (zh) | 一种高精度河道地形插值方法 | |
CN103425801A (zh) | 一种面向数字地球的数据融合方法和装置 | |
CN108896004A (zh) | 一种裂缝面粗糙度各向异性表征方法 | |
CN103941286A (zh) | 一种基于平面导航的快速三维层位解释方法 | |
CN102419436A (zh) | 基于总传播误差滤波器的多波束数据处理方法 | |
CN112394404A (zh) | 一种渐进式储层精细表征方法 | |
CN103941288A (zh) | 一种基于平面导航的三维目标地质体解释方法 | |
Hynes et al. | Lithoprobe line 55: integration of out-of-plane seismic results with surface structure, metamorphism, and geochronology, and the tectonic evolution of the eastern Grenville Province | |
CN111257870A (zh) | 一种利用InSAR监测数据的采煤沉陷积水区水下地形反演方法 | |
CN102454400A (zh) | 碳酸盐岩缝洞型储层识别方法 | |
CN107908901B (zh) | 一种基坑施工中承压水降水对周边环境影响的计算方法 | |
CN106199704A (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: No. 178 University Road, Siming District, Xiamen City, Fujian Province, 361000 Applicant after: Third Institute of Oceanography, Ministry of Natural Resources Address before: No. 178 University Road, Siming District, Xiamen City, Fujian Province, 361000 Applicant before: Oceanography Inst. No.3, State Bureau of Oceanography |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |