CN112465983A - 一种高程地表模型点线状障碍物修正方法和装置 - Google Patents
一种高程地表模型点线状障碍物修正方法和装置 Download PDFInfo
- Publication number
- CN112465983A CN112465983A CN202011197881.2A CN202011197881A CN112465983A CN 112465983 A CN112465983 A CN 112465983A CN 202011197881 A CN202011197881 A CN 202011197881A CN 112465983 A CN112465983 A CN 112465983A
- Authority
- CN
- China
- Prior art keywords
- elevation
- obstacle
- point
- work area
- model
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 74
- 230000001788 irregular Effects 0.000 claims abstract description 14
- 238000004590 computer program Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 4
- 230000004927 fusion Effects 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims 1
- 238000003379 elimination reaction Methods 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract description 9
- 230000002159 abnormal effect Effects 0.000 abstract description 8
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 18
- 230000014509 gene expression Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Road Repair (AREA)
Abstract
本发明提供一种高程地表模型点线状障碍物修正方法和装置，该方法包括：获取地表高程数据，建立地表高程模型；将采集观测系统和工区道路加载到地表高程模型；将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；融合工区道路进入网格；针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；对剔除点进行薄板样条插值，建立不规则三角网格。针对道路两旁的线状障碍物能够快速定位，精准剔除点线状异常值，运用薄板样条插值方法将障碍物处高程平滑的填充，实现点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，能够真实反映地表三维模型，有利于采集施工。
Description
技术领域
本发明涉及地球物理勘探数据处理技术领域，尤其涉及一种高程地表模型点线状障碍物修正方法和装置。
背景技术
本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
在石油地震勘探领域中利用无人机获取的地表高精度高程数据经过底层图像处理等一系列处理最终实现三维地表建模，来指导观测系统的布设。
但是在复杂地表情况下，在工区道路一定距离内一些小的地表障碍物如点状的电线杆、树木，线状的广告牌等，会对三维建模产生影响。
现在没有好的办法和技术来迅速的消除这些尖锐的异常值，从而影响了三维地表建模的精度，给观测系统布设带来困难。
因此，如何提供一种新的方案，其能够解决上述技术问题是本领域亟待解决的技术难题。
发明内容
本发明实施例提供一种高程地表模型点线状障碍物修正方法，实现了点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，有利于采集施工，该方法包括：
获取地表高程数据，建立地表高程模型；
将采集观测系统和工区道路加载到地表高程模型；
将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
融合工区道路进入网格；
针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
对剔除点进行薄板样条插值，建立不规则三角网格。
本发明实施例还提供一种高程地表模型点线状障碍物修正装置，包括：
地表高程模型建立模块，用于获取地表高程数据，建立地表高程模型；
加载模块，用于将采集观测系统和工区道路加载到地表高程模型；
网格化模块，用于将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
融合模块，用于融合工区道路进入网格；
点线状障碍判断剔除模块，用于针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
薄板样条插值模块，用于对剔除点进行薄板样条插值，建立不规则三角网格。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述一种高程地表模型点线状障碍物修正方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述一种高程地表模型点线状障碍物修正方法的计算机程序。
本发明实施例提供的一种高程地表模型点线状障碍物修正方法和装置，首先获取地表高程数据，建立地表高程模型；然后将采集观测系统和工区道路加载到地表高程模型；接着将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；下一步融合工区道路进入网格；再下一步针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；最后对剔除点进行薄板样条插值，建立不规则三角网格。本发明实施例在勘探地震勘探高程数据中勘探工区道路时，针对道路两旁一定距离内的电线杆、树木、广告牌等点线状障碍物能够快速定位，精准剔除点线状异常值，运用薄板样条插值方法将障碍物处高程平滑的填充，实现点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，插值后高程数据精确，能够快速清晰表达地表高程，这样能够真实反映地表三维模型，有利于采集施工的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例一种高程地表模型点线状障碍物修正方法示意图。
图2为本发明实施例一种高程地表模型点线状障碍物修正方法的观测系统与工区道路的平面图。
图3为本发明实施例一种高程地表模型点线状障碍物修正方法的网格化后的平面图。
图4为本发明实施例一种高程地表模型点线状障碍物修正方法的莱茵达准则示意图。
图5为本发明实施例一种高程地表模型点线状障碍物修正方法的薄板样条曲面图。
图6为本发明实施例一种高程地表模型点线状障碍物修正方法的平滑后的不规则三角网格图。
图7为运行本发明实施的一种高程地表模型点线状障碍物修正方法的计算机装置示意图。
图8为本发明实施例一种高程地表模型点线状障碍物修正装置示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1为本发明实施例一种高程地表模型点线状障碍物修正方法示意图，如图1所示，本发明实施例提供一种高程地表模型点线状障碍物修正方法，实现了点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，有利于采集施工，该方法包括：
步骤101：获取地表高程数据，建立地表高程模型；
步骤102：将采集观测系统和工区道路加载到地表高程模型；
步骤103：将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
步骤104：融合工区L道路进入网格；
步骤105：针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
步骤106：对剔除点进行薄板样条插值，建立不规则三角网格。
本发明实施例提供的一种高程地表模型点线状障碍物修正方法，首先获取地表高程数据，建立地表高程模型；然后将采集观测系统和工区道路加载到地表高程模型；接着将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；下一步融合工区道路进入网格；再下一步针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；最后对剔除点进行薄板样条插值，建立不规则三角网格。本发明实施例在勘探地震勘探高程数据中勘探工区道路时，针对道路两旁一定距离内的电线杆、树木、广告牌等点线状障碍物能够快速定位，精准剔除点线状异常值，运用薄板样条插值方法将障碍物处高程平滑的填充，实现点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，插值后高程数据精确，能够快速清晰表达地表高程，这样能够真实反映地表三维模型，有利于采集施工的目的。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法，可以包括：
获取地表高程数据DEM(Digital Elevation Model)，建立地表高程模型；
将采集观测系统S(x，y)和工区道路L加载到地表高程模型；
将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
融合工区道路进入网格；
针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达(Pauta)准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
对剔除点进行薄板样条插值，建立不规则三角网格。
图2为本发明实施例一种高程地表模型点线状障碍物修正方法的观测系统与工区道路的平面图，横轴为东坐标，单位m，纵轴为北坐标，单位m；如图2所示，具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，前述的将采集观测系统和工区道路加载到地表高程模型，包括：
将采集观测系统S(x，y)以散点的形式加载到地表高程模型；
将工区道路L以线段形式加载到地表高程模型。
图3为本发明实施例一种高程地表模型点线状障碍物修正方法的网格化后的平面图，横轴为东坐标，单位m，纵轴为北坐标，单位m；如图3所示，具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，将加载后的地表高程模型网格化，结合采集观测系统散点构建网格，主要包括：结合采集观测系统S散点构建网格，以实现将高程数据网格化。
图4为本发明实施例一种高程地表模型点线状障碍物修正方法的莱茵达准则示意图，如图4所示，具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，按照如下过程判定是否有点线状障碍：
从融合后的网格中获取高程数据列；
根据高程数据列，计算算数平均值；
根据算数平均值，计算剩余误差；
根据贝塞尔法，结合算数平均值和剩余误差，计算均方根偏差；
根据剩余误差和均方根偏差，判断高程数据中是否有点线状障碍。
实施例中，针对融合后网格对工区道路周围设定距离范围内的高程点，例如道路两旁一定距离内的电线杆、树木、广告牌等点线状障碍物，采用莱茵达准则能够快速定位和判定，具体过程包括：
从融合后的网格中获取高程数据列X1、X2、X2、…、Xn；根据高程数据列，计算算数平均值
根据算数平均值
计算剩余误差Vi；根据贝塞尔法，结合算数平均值
和剩余误差Vi，计算均方根偏差σX；根据剩余误差Vi和均方根偏差σX，判断高程数据中是否有点线状障碍。
本实施例在勘探地震勘探高程数据中勘探工区道路时，针对道路两旁一定距离内的电线杆、树木、广告牌等点线状障碍物采用莱茵达准则能够快速定位，精准剔除点线状异常值。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，按照如下方式，计算算数平均值：
其中，
为算数平均值；高程数据列为X1、X2、X2、…、Xn。
前述提到的计算算数平均值的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，按照如下方式，计算剩余误差：
其中，Vi为剩余误差；
为算数平均值。
前述提到的计算剩余误差的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，按照如下方式，计算均方根偏差：
其中，σX为均方根偏差；
为剩余误差。
前述提到的计算均方根偏差的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，按照如下方式，判断高程数据中是否有点线状障碍：
其中，σX为均方根偏差；
为剩余误差。
前述提到的判断高程数据中是否有点线状障碍的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
图5为本发明实施例一种高程地表模型点线状障碍物修正方法的薄板样条曲面图，如图5所示，具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正方法时，在一个实施例中，按照如下方式，对剔除点进行薄板样条插值：
其中，z0为待插值点的高程值；x0和y0为待插值点的横坐标和纵横坐标；a、b、c为标量；n为插值点周围控制点数量，其坐标为(xi,yi,zi)，i＝1、2、…、n；A∈R1*n，*为点乘；dj为当前待插值点与已知点集中第j个点之间的距离。
实施例中，前述的对剔除点进行薄板样条插值的表达式，也称作薄板样条插值公式。前述提到的对剔除点进行薄板样条插值的表达式为举例说明，本领域技术人员可以理解，在实施时还可以根据需要对上述公式进行一定形式的变形和添加其它的参数或数据，或者提供其它的具体公式，这些变化例均应落入本发明的保护范围。
运用薄板样条插值方法将障碍物处高程平滑的填充，实现点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，插值后高程数据精确，能够快速清晰表达地表高程，这样能够真实反映地表三维模型，有利于采集施工的目的。
在前述的对剔除点进行薄板样条插值的公式⑥中，共有N+3个系数，系数有下面的公式⑦确定：
薄板样条插值公式的求解步骤包括：
1、求解已知点与插值点的距离dij(i,j＝0、1、...、n)；
2、根据薄板样条插值公式计算公式中A的系数值Rij，即dij 2logdij。
3、根据公式1联立矩阵方程求解插值系数A如下：
4、用最小二乘法求解上述矩阵方程A＝R-1Z。得到系数A。
5、求插值点的值
插值效果如图5所示。
在对剔除点进行薄板样条插值后，地表高程数据已修正完好，接下来，图6为本发明实施例一种高程地表模型点线状障碍物修正方法的平滑后的不规则三角网格图，如图6所示，根据修正好的地表高程数据建立不规则三角网格(TIN网格，TriangulatedIrregular Network)。
本发明实施例主要目的是针对高精度DEM高程数据；结合地震采集观测系统与工区道路特征，通过采用莱茵达准则，快速定位需要剔除的点线状障碍物，进行障碍物高程剔除；利用薄板样条插值算法对此障碍物处的高程近行填充，实现点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，有利于采集施工。
本发明实施例采用了基于地表高程数据DEM，根据加载的观测系统后进行网格化，融入工区道路，然后通过莱茵达算法来对工区道路附近一定距离内的点线状障碍物进行快速定位与剔除，最后经过薄板样条插值算法插出异常值所在点的高程数据。
本发明实施例可以在地震勘探高程数据中勘探工区道路一定范围内快速定位点线状障碍物。能够精准剔除点线状异常值剔。插值后高程数据精确。能够快速清晰表达地表高程。
本发明实施例读取某一地表高程数据，将采集观测系统以散点形式加载到地表高程模型中，并结合观测系统以一定间距对地表高程进行网格化，将工区道路加载到地表高程模型中，并融进网格中，运用莱茵达准则针对道路两旁一定距离内的电线杆，树木，广告牌等点线状障碍物能够快速定位，并将其高程值剔除，运用薄板样条插值方法将障碍物处高程平滑的填充，这样能够真实反映地表三维模型，有利于采集施工的目的。
图7为运行本发明实施的一种高程地表模型点线状障碍物修正方法的计算机装置示意图，如图7所示，本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述一种高程地表模型点线状障碍物修正方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行实现上述一种高程地表模型点线状障碍物修正方法的计算机程序。
本发明实施例中还提供了一种高程地表模型点线状障碍物修正装置，如下面的实施例所述。由于该装置解决问题的原理与一种高程地表模型点线状障碍物修正方法相似，因此该装置的实施可以参见一种高程地表模型点线状障碍物修正方法的实施，重复之处不再赘述。
图8为本发明实施例一种高程地表模型点线状障碍物修正装置示意图，如图8所示，本发明实施例还提供一种高程地表模型点线状障碍物修正装置，具体实施时可以包括：
地表高程模型建立模块801，用于获取地表高程数据，建立地表高程模型；
加载模块802，用于将采集观测系统和工区道路加载到地表高程模型；
网格化模块803，用于将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
融合模块804，用于融合工区道路进入网格；
点线状障碍判断剔除模块805，用于针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
薄板样条插值模块806，用于对剔除点进行薄板样条插值，建立不规则三角网格。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的加载模块，具体用于：
将采集观测系统以散点的形式加载到地表高程模型；
将工区道路以线段形式加载到地表高程模型。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的点线状障碍判断剔除模块，具体用于：
从融合后的网格中获取高程数据列；
根据高程数据列，计算算数平均值；
根据算数平均值，计算剩余误差；
根据贝塞尔法，结合算数平均值和剩余误差，计算均方根偏差；
根据剩余误差和均方根偏差，判断高程数据中是否有点线状障碍。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的点线状障碍判断剔除模块，还用于按照如下方式，计算算数平均值：
其中，
为算数平均值；高程数据列为X1、X2、X2、…、Xn。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的点线状障碍判断剔除模块，还用于按照如下方式，计算剩余误差：
其中，Vi为剩余误差；
为算数平均值。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的点线状障碍判断剔除模块，还用于按照如下方式，计算均方根偏差：
其中，σX为均方根偏差；
为剩余误差。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的点线状障碍判断剔除模块，还用于按照如下方式，判断高程数据中是否有点线状障碍：
其中，σX为均方根偏差；
为剩余误差。
具体实施本发明实施例提供的一种高程地表模型点线状障碍物修正装置时，在一个实施例中，前述的薄板样条插值模块，具体用于按照如下方式，对剔除点进行薄板样条插值：
其中，z0为待插值点的高程值；x0和y0为待插值点的横坐标和纵横坐标；a、b、c为标量；n为插值点周围控制点数量，其坐标为(xi,yi,zi)，i＝1、2、…、n；A∈R1*n，*为点乘；dj为当前待插值点与已知点集中第j个点之间的距离。
综上，本发明实施例提供的一种高程地表模型点线状障碍物修正方法和装置，首先获取地表高程数据，建立地表高程模型；然后将采集观测系统和工区道路加载到地表高程模型；接着将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；下一步融合工区道路进入网格；再下一步针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；最后对剔除点进行薄板样条插值，建立不规则三角网格。本发明实施例在勘探地震勘探高程数据中勘探工区道路时，针对道路两旁一定距离内的电线杆、树木、广告牌等点线状障碍物能够快速定位，精准剔除点线状异常值，运用薄板样条插值方法将障碍物处高程平滑的填充，实现点线状障碍物的修正，达到观测系统及道路沿线的线路平滑，插值后高程数据精确，能够快速清晰表达地表高程，这样能够真实反映地表三维模型，有利于采集施工的目的。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (18)
1.一种高程地表模型点线状障碍物修正方法，其特征在于，包括：
获取地表高程数据，建立地表高程模型；
将采集观测系统和工区道路加载到地表高程模型；
将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
融合工区道路进入网格；
针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
对剔除点进行薄板样条插值，建立不规则三角网格。
2.如权利要求1所述的方法，其特征在于，将采集观测系统和工区道路加载到地表高程模型，包括：
将采集观测系统以散点的形式加载到地表高程模型；
将工区道路以线段形式加载到地表高程模型。
3.如权利要求1所述的方法，其特征在于，按照如下过程判定是否有点线状障碍：
从融合后的网格中获取高程数据列；
根据高程数据列，计算算数平均值；
根据算数平均值，计算剩余误差；
根据贝塞尔法，结合算数平均值和剩余误差，计算均方根偏差；
根据剩余误差和均方根偏差，判断高程数据中是否有点线状障碍。
4.如权利要求3所述的方法，其特征在于，按照如下方式，计算算数平均值：
其中，
为算数平均值；高程数据列为X1、X2、X2、…、Xn。
5.如权利要求3所述的方法，其特征在于，按照如下方式，计算剩余误差：
其中，Vi为剩余误差；
为算数平均值。
6.如权利要求3所述的方法，其特征在于，按照如下方式，计算均方根偏差：
其中，σX为均方根偏差；
为剩余误差。
7.如权利要求3所述的方法，其特征在于，按照如下方式，判断高程数据中是否有点线状障碍：
其中，σX为均方根偏差；
为剩余误差。
8.如权利要求1所述的方法，其特征在于，按照如下方式，对剔除点进行薄板样条插值：
其中，z0为待插值点的高程值；x0和y0为待插值点的横坐标和纵横坐标；a、b、c为标量；n为插值点周围控制点数量，其坐标为(xi,yi,zi)，i＝1、2、…、n；A∈R1*n，*为点乘；dj为当前待插值点与已知点集中第j个点之间的距离。
9.一种高程地表模型点线状障碍物修正装置，其特征在于，包括：
地表高程模型建立模块，用于获取地表高程数据，建立地表高程模型；
加载模块，用于将采集观测系统和工区道路加载到地表高程模型；
网格化模块，用于将加载后的地表高程模型网格化，结合采集观测系统散点构建网格；
融合模块，用于融合工区道路进入网格；
点线状障碍判断剔除模块，用于针对融合后网格对工区道路周围设定距离范围内的高程点采用莱茵达准则进行判定是否有点线状障碍，对判定出的点线状障碍进行剔除；
薄板样条插值模块，用于对剔除点进行薄板样条插值，建立不规则三角网格。
10.如权利要求9所述的装置，其特征在于，加载模块，具体用于：
将采集观测系统以散点的形式加载到地表高程模型；
将工区道路以线段形式加载到地表高程模型。
11.如权利要求9所述的装置，其特征在于，点线状障碍判断剔除模块，具体用于：
从融合后的网格中获取高程数据列；
根据高程数据列，计算算数平均值；
根据算数平均值，计算剩余误差；
根据贝塞尔法，结合算数平均值和剩余误差，计算均方根偏差；
根据剩余误差和均方根偏差，判断高程数据中是否有点线状障碍。
12.如权利要求11所述的装置，其特征在于，点线状障碍判断剔除模块，还用于按照如下方式，计算算数平均值：
其中，
为算数平均值；高程数据列为X1、X2、X2、…、Xn。
13.如权利要求11所述的装置，其特征在于，点线状障碍判断剔除模块，还用于按照如下方式，计算剩余误差：
其中，Vi为剩余误差；
为算数平均值。
14.如权利要求11所述的装置，其特征在于，点线状障碍判断剔除模块，还用于按照如下方式，计算均方根偏差：
其中，σX为均方根偏差；
为剩余误差。
15.如权利要求11所述的装置，其特征在于，点线状障碍判断剔除模块，还用于按照如下方式，判断高程数据中是否有点线状障碍：
其中，σX为均方根偏差；
为剩余误差。
16.如权利要求9所述的装置，其特征在于，薄板样条插值模块，具体用于按照如下方式，对剔除点进行薄板样条插值：
其中，z0为待插值点的高程值；x0和y0为待插值点的横坐标和纵横坐标；a、b、c为标量；n为插值点周围控制点数量，其坐标为(xi,yi,zi)，i＝1、2、…、n；
*为点乘；dj为当前待插值点与已知点集中第j个点之间的距离。
17.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至8任一项所述方法。
18.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行实现权利要求1至8任一项所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197881.2A CN112465983A (zh) | 2020-10-30 | 2020-10-30 | 一种高程地表模型点线状障碍物修正方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011197881.2A CN112465983A (zh) | 2020-10-30 | 2020-10-30 | 一种高程地表模型点线状障碍物修正方法和装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112465983A true CN112465983A (zh) | 2021-03-09 |
Family
ID=74834347
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011197881.2A Pending CN112465983A (zh) | 2020-10-30 | 2020-10-30 | 一种高程地表模型点线状障碍物修正方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112465983A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112883339A (zh) * | 2021-03-11 | 2021-06-01 | 北京市地震局 | 一种地震有感范围确定方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101008394B1 (ko) * | 2010-10-12 | 2011-01-14 | (주)동광지엔티 | 수치표면모형의 경계선을 이용한 수치표고모형의 건물영역 평탄화 장치 |
| CN103399990A (zh) * | 2013-07-18 | 2013-11-20 | 北京工业大学 | 城市排水仿真模拟系统的道路精细离散网格构建方法 |
| KR101548647B1 (ko) * | 2014-12-19 | 2015-08-31 | 한국건설기술연구원 | 3차원 지형정보 가시화 프로세서 및 그 동작방법 |
| CN106017472A (zh) * | 2016-05-17 | 2016-10-12 | 成都通甲优博科技有限责任公司 | 全局路线规划方法、全局路线规划系统及无人机 |
| CN106228521A (zh) * | 2016-07-25 | 2016-12-14 | 哈尔滨工业大学 | 一种基于薄板样条插值的障碍物特征提取方法 |
-
2020
- 2020-10-30 CN CN202011197881.2A patent/CN112465983A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101008394B1 (ko) * | 2010-10-12 | 2011-01-14 | (주)동광지엔티 | 수치표면모형의 경계선을 이용한 수치표고모형의 건물영역 평탄화 장치 |
| CN103399990A (zh) * | 2013-07-18 | 2013-11-20 | 北京工业大学 | 城市排水仿真模拟系统的道路精细离散网格构建方法 |
| KR101548647B1 (ko) * | 2014-12-19 | 2015-08-31 | 한국건설기술연구원 | 3차원 지형정보 가시화 프로세서 및 그 동작방법 |
| CN106017472A (zh) * | 2016-05-17 | 2016-10-12 | 成都通甲优博科技有限责任公司 | 全局路线规划方法、全局路线规划系统及无人机 |
| CN106228521A (zh) * | 2016-07-25 | 2016-12-14 | 哈尔滨工业大学 | 一种基于薄板样条插值的障碍物特征提取方法 |
Non-Patent Citations (1)
| Title |
|---|
| 石翠萍;: "利用等高线制作格网DEM方法的研究", 黑龙江工程学院学报, no. 04, 25 August 2014 (2014-08-25) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112883339A (zh) * | 2021-03-11 | 2021-06-01 | 北京市地震局 | 一种地震有感范围确定方法及系统 |
| CN112883339B (zh) * | 2021-03-11 | 2023-11-28 | 北京市地震局 | 一种地震有感范围确定方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109033538B (zh) | 一种基于实测结构面参数的裂隙岩体渗透张量的计算方法 | |
| CN112381937B (zh) | 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法 | |
| CN102506824B (zh) | 一种城市低空无人机系统生成数字正射影像图的方法 | |
| CN104536043B (zh) | 一种深度域整体速度模型融合方法及装置 | |
| CN107527382B (zh) | 数据处理方法以及装置 | |
| CN109917454B (zh) | 基于双基准面的真地表叠前深度偏移成像方法及装置 | |
| CA2876632C (en) | System and method for migration velocity modeling | |
| CN104181593B (zh) | 一种三维无射线追踪回折波层析成像方法及装置 | |
| CN102944896A (zh) | 表层调查数据的模型法静校正方法 | |
| CN101984366A (zh) | 基于共孔径面的起伏地表叠前时间偏移方法及设备 | |
| WO2012001388A2 (en) | Gravity survey data processing | |
| CN104865597A (zh) | 一种深度域层速度初始模型的建模方法 | |
| CN112465983A (zh) | 一种高程地表模型点线状障碍物修正方法和装置 | |
| CN105353406B (zh) | 一种生成角道集的方法和装置 | |
| CN111239822B (zh) | 一种消除井震闭合差的断层下速度建模方法及处理终端 | |
| CN111257969B (zh) | 一种断层控制下高精度速度建模方法及处理终端 | |
| CN112505780B (zh) | 地层深度数据的校正方法及装置 | |
| CN112363221B (zh) | Walkaway VSP测线的布设方法及装置 | |
| CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
| NO20181090A1 (en) | Visualizing attributes of multiple fault surfaces in real time | |
| CN105549075B (zh) | 求取浅层速度场的方法及装置 | |
| CN114966852B (zh) | 双界面控制浅深层速度融合与静校正量计算方法及装置 | |
| JP7122901B2 (ja) | ブレイクライン抽出プログラム、ブレイクライン抽出システム | |
| CN112379420B (zh) | 高精度弯曲测线叠前时间域成像方法及装置 | |
| CN111198400B (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 |