CN104835186B - 一种历史地震图波形自动跟踪方法 - Google Patents
一种历史地震图波形自动跟踪方法 Download PDFInfo
- Publication number
- CN104835186B CN104835186B CN201510214008.2A CN201510214008A CN104835186B CN 104835186 B CN104835186 B CN 104835186B CN 201510214008 A CN201510214008 A CN 201510214008A CN 104835186 B CN104835186 B CN 104835186B
- Authority
- CN
- China
- Prior art keywords
- mrow
- seismogram
- path
- point
- curve
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
一种历史地震图波形自动跟踪方法,属于油气地球物理勘探技术领域。其特征包括:(1)利用微分几何曲线提取算法从地震图中提取出地震波曲线点,以及地震波曲线的局部方向;(2)依据已经提取出来的曲线局部方向将提取出来的曲线点,连成多条曲线段;(3)最后依据地震波曲线的自身特性,如地震波曲线在其拐点处呈现出的对称特性以及地震波曲线围绕其中轴线做往返运动的震动特性,以回溯的思路对地震波曲线段进行跟踪。本发明的效果和益处是提高了地震波波形自动跟踪的精度,无需人工干预,能够处理有交叉线和断点的地震图。
Description
技术领域
本发明属于油气地球物理勘探技术领域,涉及到对历史地震图的数字化,特别涉及到数字化过程中波形自动跟踪的方法。
背景技术
在数字化地震观测系统建成之前的百年间,人们通过使用模拟笔的方法对地震图进行记录,留下了许多宝贵的无可复制的地震资料。这些历史地震图一般使用纸张、胶带等进行记录,具有明显的保存时间限制,为了方便使用现代地震学方法对其进行分析同时也为了延长这些资料的保存时间,我们需要对其进行数字化。
近年来,越来越多的学者专家投入到历史地震图的数字化工作中。其中如何高效并且自动的跟踪地震波成为历史地震图数字化的关键步骤。目前已有的跟踪方法主要分为两大类。第一类跟踪算方法属于半自动范畴,这些方法需要大量的人机交互过程,过程很繁琐,比较耗时,如2003年Bromirki等设计的SeisDig系统,以及2005年Pintore等设计的Teseo2系统。第二类跟踪方法可以实现对地震图波形的全自动跟踪,然而这些方法一般依赖于地震图图像的局部特征,对于交叉线和断点较多的地震图,他们不能给出很好的结果,如Church等人在2013年提出的基于形态学的波形自动跟踪算法,Mao等人在2014年提出的基于逐点跟踪策略的波形自动跟踪算法。
发明内容
本发明的目的是提供了一种历史地震图波形自动跟踪方法,该方法将利用地震波曲线的自身特性,如地震波曲线在其拐点处呈现出的对称特性以及地震波曲线围绕其中轴线做往返运动的震动特性,采用回溯的思路完成对波形的自动跟踪。该方法提高了地震波波形自动跟踪的精度,无需人工干预,能够处理有交叉线和断点的地震图。
本发明的技术方案是:
一、技术方案的原理
(1)技术方案的基本概念。
历史地震图:以滚筒的形式由模拟笔记录到纸张、胶带上的原始地震资料。
目标地震图:历史地震图经过扫描仪得到的光栅图像,是本发明的处理对象。
跟踪:从目标地震图中,提取并且识别出不同的地震波曲线。
地震波曲线:目标地震图中,地震波波形曲线往往具有一定的线宽,本方案中,地震波曲线指的是地震波波形曲线的中心线,并且是本方案的跟踪对象。
曲线局部走向:即地震波曲线的切线方向。
拐点:即地震波曲线走向发生变化的点。
地震波中轴线:地震波曲线由模拟笔记录完成,模拟笔在记录过程中沿其平衡位置上下摆动,本方案称这个平衡位置为地震波中轴线。
(2)技术方案的具体原理。
本发明方法由三部分组成,首先依据微分几何曲线提取算法从目标地震图中提取出地震波曲线点以及地震波曲线的局部走向;然后结合地震波的局部走向将提取出来的地震波曲线点连接成多条曲线段;最后依据地震波曲线的自身特性,如地震波曲线在其拐点处呈现出的对称特性以及地震波曲线围绕其中轴线做往返运动的震动特性,以回溯的思路对地震波曲线段进行跟踪。下面将对上述每一步的技术原理进行详细介绍,并在介绍地震波跟踪算法之前,将地震波的一些性质与其计算方法也给与相应的介绍。
1)地震波曲线点以及曲线局部走向的提取
对于目标地震图I,沿地震波波形曲线的切线方向做灰度剖面,该灰度剖面的极值点即为地震波曲线点。由极值点的性质可知,对于地震波曲线点来说,它沿切线方向的一阶方向导数为零二阶方向导数的绝对值较大。本发明取I的二阶方向导数最大的方向为地震波波形曲线的切线方向,记上述切线方向为其中为该方向与X轴的夹角。I沿任意方向(cosθ,sinθ)的二阶方向导数为:Iθθ=Ixxcos2θ+2Ixycosθsinθ+Iyysin2θ。令可以求得I的二阶方向导数最大的方向,该方向即为法线方向,如公式(1)所示:
其中Ixx,Ixy,Iyx,Iyy为I的各阶偏导。一旦地震波波形曲线的法线方向被计算出来,与该方向垂直的方向即为地震波曲线的局部走向。
由于目标地震图是离散的,其最小计量单位为像素,在本发明中认为像素为正方形,对于某一像素点(x,y),它的定义范围为直接计算上述灰度剖面的极值点,仅能得到像素级别的精度,为了得到亚像素级别的精度,我们需要对每个像素点进行插值。图像处理中泰勒展开式是获取亚像素的一种重要手段。对像素点(x,y)沿方向做泰勒展开:
其中为估计的亚像素点,t为亚像素点与展开点(x,y)的距离,由于在一个像素内部做展开故t<0.5。这样亚像素级别的极值点即地震波曲线点可以通过公式(3)、(4)求得。
其中γ1为用户设置的阈值。
2)地震波曲线点的连接策略
记已提取的地震波曲线点的集合为C,从中选取二阶导数较大的点作为连线的初始种子点。对于每一个初始种子点,本发明将按照图3所示的策略进行连线。从初始种子点开始,沿当前地震波曲线方向搜索初始种子点的三个邻域像素内是否存在地震波曲线点,例如初始种子点为(x,y),当前地震波曲线方向为θ且θ∈[-22.5°,22.5°],本发明将搜索(x+1,y-1),(x+1,y),(x+1,y+1)这三个邻域像素内是否存在地震波曲线点,如果不存在,本次连线结束,重新选择新的初始种子点重复连线过程;如果存在,标记这些曲线点为候选种子点,按公式(5)计算这些候选种子点与初始种子点的代价函数值,删除代价大于阈值的候选种子点,记该阈值为γ2;删除后如果一个候选种子点都没剩下,声明本次连线结束,重新选择新的初始种子点重复连线过程;删除操作后如果还有多个候选种子点,取代价最小的候选种子点作为最优连线点,而后更新初始种子点为已经找到的最优连线点;重复上述策略。本过程是一个迭代过程,当搜索不到候选种子点时,该迭代过程结束。
式中Cost表示代价函数,α1表示地震波曲线在初始种子点的方向,α2表示地震波曲线在候选种子点的方向。
3)地震波曲线的自身性质以及计算方法
在地震波曲线点提取完成后,本方案对提取出来的地震波曲线点进行了初步连接,将其连成一段段小的曲线段。为了完成地震图波形的跟踪,需要设计算法对这些小的曲线段进行跟踪,即那些属于同一条地震波曲线的曲线段被划分为一类。由于地震波曲线的自身性质在跟踪算法中极为重要,本发明在介绍跟踪算法之前,首先介绍地震波曲线的一些重要的性质。本发明通过对大量地震图的研究得出地震波曲线的以下三点重要性质。
性质1:对于一笔画下来的相邻的两个地震波曲线段来说,这两条曲线段曲线走向相差很小。
曲线段走向可按如下策略计算:对于任一条曲线线段L,对其首尾两端M个曲线点的方向进行加权平均计算其首尾两端的曲线走向,如公式(6)所示,取L的最左端为其首部,最右端为其尾部。
其中β1,β2为曲线首部和尾部走向,d(L(k))表示种子线段L在第k个点的局部方向,N表示种子线段L包含的曲线点数。
性质2:地震波曲线由模拟笔的摆动绘制而成,存在大量的拐点。对于这些拐点来说,位于拐点处两端的地震波曲线段具有一定的对称性和走向相反的特性。
关于曲线走向是否相反的判断,通过公式(6)计算曲线的走向并作出判断。关于两条曲线对称性的判断,本发明采用如下策略:记这两条曲线为L1,L2。将这两条曲线分别等分为N份,计算切分后的每一小段曲线的斜率。而后按公式(7)判断两条曲线的对称性。公式(7)中K(i)为曲线L1切分后的第i条曲线段的斜率,k(i)为曲线L2切分后的第i条曲线段的斜率,γ3是用来判断对称性的阈值,由人工设定。
性质3:每条地震波曲线都有一条属于自己的中轴线,地震波曲线围绕其做往返运动,并且不同的地震波曲线具有不同的中轴线。
本发明采用IRLS迭代加权最小二乘法对地震波曲线的中轴线进行估计。记当前地震波曲线为Γ,关于Γ的中轴线具体估计步骤如下:
Step1:对Γ的纵坐标进行平均求和,给出中轴线的初始估计,记为Linital。
Step2:对Γ的纵坐标进行加权求和,Γ中每个点的权值按公式(8)计算,记该次估计的中心线为Lnext。
Step3:判断|Lnext-Linital|<1是否成立,如果成立转向Step5,否则转入Step4。
Step4:更新Linital为Lnext,转入Step2.
Step5:中轴线估计结束。
公式(8)中,ωi为Γ中第i个点纵坐标的权值,yi为Γ中第i个点纵坐标,σ是一个很小的值,主要是为了防止除数为零的情况发生,取0.001。
3)基于回溯的跟踪方法
回溯算法也叫做试探法,它是一种系统的搜索问题解的方法。回溯的基本思想为:从一条路线出发开始跟踪,能进则进,不能进则退回来,换一条路再试试。本发明将回溯的思路引入到地震图波形的跟踪上。
将整幅地震图看作是一个有向图G;记初步连接地震波曲线点过程中得到的地震波曲线段为G的有向边。则地震图波形的跟踪问题,可以看做是从有向图G中按照地震波波形的特点寻找最优路径的问题。对于有向图中最优路径的选取可以采用回溯的算法给予实现,其基本思路为:从G中的一条有向边出发往前走,碰到交叉点,首先判断该点是否为拐点。如果不是拐点选择与当前有向边走向相差最小的一条有向边作为最合适的路径并往前走;如果是拐点执行回溯的思路,在拐点处对跟踪路径进行分叉,沿拐点处的每条有向边都分出一条新的路径,记这些路径为Ψ,按照IRLS迭代加权最小二乘法对每条路径的曲线的中轴线进行估计并计算这些估计的中轴线与当前跟踪的地震波中轴线的差值,取差值最小的路径最为最合适的路径并往前走。整个跟踪过程是个迭代过程,每次碰到交叉点,就按照上述规则进行路径选取,当在G中搜索不到新的有向边时,迭代过程结束。
上述跟踪过程中需要用到当前跟踪的地震波中轴线。我们利用已经跟踪到的曲线段,通过IRLS迭代加权最小二乘法对该中轴线进行估计。
二、技术方案的步骤
S1获取目标地震图的灰度图片记为I,I的像素灰度级为0到255;
S2将I与高斯模糊核进行卷积,用于压制I中噪声,通过sobel算子计算I的各阶偏导数Ix,Iy,Ixx,Ixy,Iyy。
S3依据公式(1)计算地震波波形曲线的切线方向,取与该方向垂直的方向为地震波曲线的局部走向。设置公式(4)中的阈值,依据公式(3)、(4)提取出地震波曲线点。
S4记步骤S2中提取的曲线点的集合为Φ,对于Φ中每一点沿公式(1)给出的方向计算其二阶方向导数,对这些二阶方向导数进行归一化处理,取二阶方向导数大于0.2的曲线点作为连线的起始种子点。对于每一个起始种子点,按图3所示的流程连接成曲线段,流程图中的阈值设定为0.12*π。
S5将整幅地震图看作是一个有向图G;记步骤S3中得到的所有地震波曲线段为G的有向边,按照公式(6)计算G中所有有向边的首尾走向。从G中选取用于跟踪地震波曲线的初始跟踪路径,初始跟踪路径按照如下方案选取:将有向图G中所有的有向边按照它们的起始点的横坐标大小从小到大进行排序,取排在前10%的边作为初始跟踪路径。
S5.1从步骤S5中取出第一条初始跟踪路径。
S5.2在初始跟踪路径尾部节点处搜索有向边作为候选的跟踪路径,该搜索区域为半径为2.5倍的地震波波形曲线宽度的圆形区域,如果一条候选的跟踪路径都没找到跳到步骤S5.8;如果找到了候选路径,进入步骤S5.3。
S5.3利用本次跟踪过程已经跟踪过的路径,结合IRLS迭代加权最小二乘法对当前跟踪的地震波曲线的中轴线进行估计。
S5.4按照公式(7)计算所有候选跟踪路径与初始路径之间的相似性,计算过程中N取4。如果这些候选跟踪路径中存在某条路径和初始路径相似并且这两条路径走向相反,则认为本次跟踪过程中碰到了拐点,并转入步骤S5.6,否则转入步骤S5.5。
S5.5从步骤S5.2中所有候选跟踪路径中,选取首部走向和初始跟踪路径尾部走向差值最小的路径,作为后续跟踪路径,并转入步骤S5.7。
S5.6执行回溯策略,在拐点处进行路径分叉,拐点处的每条候选跟踪路径都作为后续跟踪路径分出多条新的路径,记这些路径为Ψ,采用IRLS迭代加权最小二乘法对每条路径的中轴线进行估计并计算这些中轴线与当前跟踪的地震波中轴线的差值,取差值最小的路径作为本次跟踪过程中的后续跟踪路径,舍弃其他路径。
S5.7更新初始跟踪线段为找到的后续跟踪路径,转入步骤S5.2.
S5.8本次跟踪过程结束,保存本次跟踪过程中的路径,存入计算机。
S5.9按顺序取步骤S5中未执行跟踪过程的起始种子线段,重复步骤S5.2-5.8。
本发明的和效果和益处是,本发明将地震图的局部特征和地震波波形的宏观特性相结合并借鉴回溯的思想,提高了地震波波形自动跟踪的精度,无需大量的人工干预,可以较好的处理交叉线断点较多的地震图。
附图说明
附图1是经过扫描仪得到的目标地震示意图。
附图2是历史地震图波形自动跟踪算法流程图。
附图3是地震波曲线点连接成线的流程图。
具体实施方式
以下结合技术方案和附图详细叙述本发明的具体实施方式。
实施例:
本例以冰岛地震图的一部分为实验目标,该地震图由大连地震台记录于2000年7月21号。实验平台CPU为I3-2120,运行环境matlab 2012a。
步骤1、获取实例地震图的灰度图片记为I,I的像素灰度级为0到255
步骤2、将I与σ=1.4的高斯模糊核与进行卷积,用于压制I中噪声,采用sobel算子求取卷积后的图像的各阶偏导数,Ix,Iy,Ixx,Ixy,Iyy。
步骤3、依据公式(1)计算地震波波形曲线的切线方向,取与该方向垂直的方向为地震波曲线的局部走向。设置公式(4)中的阈值γ1=1.2,依据公式(3)、(4)提取出地震波曲线点,对每一个曲线点沿公式(1)给出的方向向计算其二阶方向导数,从这些曲线点中删除二阶方向导数小于0的点,余下的点作为最终提取出来的地震波曲线点。
步骤4、对步骤3中地震波曲线点的二阶方向导数进行归一化,取二阶方向导数大于0.2的曲线点作为连线的起始种子点。对于每一个起始种子点,按图3所示的流程连接成曲线段,连线过程的γ3设定为0.12π。
步骤5、将整幅地震图看作是一个有向图G;记步骤S3中得到的地震波曲线段为G的有向边,按照公式(6)计算G中所有有向边的首尾走向,公式(6)中M取8。从G中选取用于跟踪地震波曲线的初始跟踪路径,初始跟踪路径按照如下方案选取:将有向图G中所有的有向边按照它们的起始点的横坐标大小从小到大进行排序,取排在前10%的有向边作为初始跟踪路径。
步骤5.1、从步骤5中的取出第一条初始路径。
步骤5.2、在初始路径尾部节点处搜索有向边作为候选的跟踪路径,该搜索区域为半径为15个像素的圆形区域,如果一条候选的跟踪路径都没找到跳到步骤5.8;如果找到了候选跟踪路径,进入步骤5.3。
步骤5.3利用本次跟踪过程已经跟踪过的路径,结合IRLS迭代加权最小二乘法对当前跟踪的地震波曲线的中轴线进行估计。
步骤5.4按照公式(7)计算所有候选跟踪路径与初始路径之间的相似性,计算过程中N取4,γ2取0.1。如果这些候选跟踪路径中存在某条路径和初始路径的相似并且这两条路径走向相反,则认为本次跟踪过程中遇到了拐点,转入步骤5.6,否则转入步骤5.5。
步骤5.5从步骤5.2中所有候选跟踪路径中,选取首部走向和初始跟踪路径尾部走向差值最小的路径,作为后续跟踪路径,并转入步骤5.7。
步骤5.6执行回溯策略,在拐点处进行路径分叉,拐点处的每条候选跟踪路径都作为后续跟踪路径分出多条新的路径,记这些路径为Ψ,按照IRLS迭代加权最小二乘法对每条路径的中轴线进行估计并计算这些估计的中轴线与当前跟踪的地震波中轴线的差值,取差值最小的路径作为本次跟踪过程中的后续跟踪路径,舍弃其他路径。
步骤5.7更新初始跟踪线段为找到的后续跟踪路径,转入步骤5.2.
步骤5.8声明本次跟踪过程结束,保存本次跟踪过程中的路径,存入计算机。
步骤5.9按顺序取步骤5中未执行跟踪过程的跟踪路径,重复步骤5.2-5.9。
Claims (1)
1.一种历史地震图波形自动跟踪方法,利用地震波曲线的自身特性并借鉴回溯的思想,完成对历史地震图波形的自动跟踪,该方法由三部分组成,首先依据微分几何曲线提取算法从目标地震图中提取出地震波曲线点以及地震波曲线的局部走向;然后结合地震波曲线的局部走向将提取出来的地震波曲线点连接成多条曲线段;最后依据地震波曲线的自身特性,地震波曲线在其拐点处呈现出的对称特性以及地震波曲线围绕其中轴线做往返运动的震动特性,以回溯的思路对地震波曲线段进行跟踪,其特征在于如下步骤:
S1将整幅地震图看作是一个有向图G;本方法前两部分得到的所有地震波曲线段为G的有向边,按照公式(1)计算G中所有有向边的首尾走向;从G中选取用于跟踪地震波曲线的初始跟踪路径,初始跟踪路径按照如下方案选取:将有向图G中所有的有向边按照它们的起始点的横坐标大小从小到大进行排序,取排在前10%的有向边作为初始跟踪路径;
<mrow>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>M</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>N</mi>
<mo>-</mo>
<mi>M</mi>
</mrow>
<mi>N</mi>
</munderover>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(1)中β1,β2为有向边的首部和尾部走向,L为有向边,d(L(k))表示有向边L在第k个点的切线方向,N表示L中包含的曲线点数,M表示采用多少个点做加权平均;
S1.1从步骤S1中取出第一条初始跟踪路径;
S1.2在初始跟踪路径尾部节点处搜索有向边作为候选跟踪路径,搜索区域为半径为2.5倍的地震波波形曲线宽度的圆形区域,如果一条候选跟踪路径都没找到跳到步骤S1.6;如果找到了候选跟踪路径进入步骤S1.3;
S1.3利用本次跟踪过程已经跟踪过的路径,采用IRLS迭代加权最小二乘法对当前跟踪的地震波曲线的中轴线进行估计;
S1.4按照公式(2)计算所有候选跟踪路径与初始跟踪路径之间的相似性,计算过程中N取4,如果这些候选跟踪路径中存在某条路径和初始跟踪路径相似并且这两条路径走向相反,则认为本次跟踪过程中遇到了拐点,并转入步骤S1.6,否则转入步骤S1.5;
<mrow>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>|</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo><</mo>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
公式(2)中,N表示将这两条待比较对称性的曲线等分为多少份,K(i)为其中一条曲线切分后的第i条曲线段的斜率,k(i)为另外一条曲线切分后的第i条曲线段的斜率,γ2是用来判断对称性的阈值,由人工设定;
S1.5从步骤S1.2中所有候选跟踪路径中,选取首部走向和初始跟踪路径尾部走向差值最小的路径,作为后续跟踪路径,并转入步骤S1.7;
S1.6执行回溯策略,在拐点处进行路径分叉,将拐点处的每条候选跟踪路径都作为后续跟踪路径分出多条新的路径,记这些路径为Ψ,按照IRLS迭代加权最小二乘法对每条路径的中轴线进行估计并计算这些中轴线与当前跟踪的地震波中轴线的差值,取差值最小的路径作为本次跟踪过程中的后续跟踪路径,舍弃其他路径;
S1.7更新初始跟踪线段为找到的后续跟踪路径,转入步骤S1.2;
S1.8声明本次跟踪过程结束,保存本次跟踪过程中的路径,存入计算机;
S1.9取步骤S1中未执行跟踪过程中的起始种子线段,重复步骤S1.2-1.8,所述的起始种子线段为步骤S1中排在前10%的有向边。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510214008.2A CN104835186B (zh) | 2015-04-30 | 2015-04-30 | 一种历史地震图波形自动跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510214008.2A CN104835186B (zh) | 2015-04-30 | 2015-04-30 | 一种历史地震图波形自动跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104835186A CN104835186A (zh) | 2015-08-12 |
CN104835186B true CN104835186B (zh) | 2017-11-10 |
Family
ID=53813051
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510214008.2A Expired - Fee Related CN104835186B (zh) | 2015-04-30 | 2015-04-30 | 一种历史地震图波形自动跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104835186B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111247874B (zh) * | 2018-04-27 | 2022-03-18 | 上海趋视信息科技有限公司 | 照明控制系统和方法 |
CN116520049B (zh) * | 2023-04-13 | 2024-05-17 | 广州致远仪器有限公司 | 波形输出设置方法、装置、设备以及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551467A (zh) * | 2009-04-23 | 2009-10-07 | 中国石油化工股份有限公司胜利油田分公司物探研究院 | 一种基于边界检测的初至波自动拾取方法 |
EP2594964A2 (en) * | 2011-11-21 | 2013-05-22 | CGGVeritas Services S.A. | Device and method for computing depth velocity variations |
CN103592681A (zh) * | 2013-09-16 | 2014-02-19 | 电子科技大学 | 一种基于信号分类的地震图像层位追踪方法 |
-
2015
- 2015-04-30 CN CN201510214008.2A patent/CN104835186B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551467A (zh) * | 2009-04-23 | 2009-10-07 | 中国石油化工股份有限公司胜利油田分公司物探研究院 | 一种基于边界检测的初至波自动拾取方法 |
EP2594964A2 (en) * | 2011-11-21 | 2013-05-22 | CGGVeritas Services S.A. | Device and method for computing depth velocity variations |
CN103592681A (zh) * | 2013-09-16 | 2014-02-19 | 电子科技大学 | 一种基于信号分类的地震图像层位追踪方法 |
Non-Patent Citations (4)
Title |
---|
Automatic Misalignment Correction of Seismograms Using Low-Rank Matrix Recovery;Lifu Wang et al;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20130430;第10卷(第2期);352-356 * |
On the digitizing of historical seismograms;D Grabrovec et al;《GEOFIZIKA》;19941231;第11卷;27-31 * |
历史地震图波形曲线提取方法研究;孟涛 等;《中国地球科学联合学术年会 2014》;20141231;865 * |
地震波的计算机自动提取—智能跟踪与细化;肖锦玉 等;《地球物理学报》;19920430;第35卷(第2期);262-265 * |
Also Published As
Publication number | Publication date |
---|---|
CN104835186A (zh) | 2015-08-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Perceiving spectral variation: Unsupervised spectrum motion feature learning for hyperspectral image classification | |
Zhang et al. | Boundary-constrained multi-scale segmentation method for remote sensing images | |
CN101364259B (zh) | 多层次知识驱动的全色遥感影像的道路变化信息提取方法 | |
CN103400134B (zh) | 非接触式手指静脉样本感兴趣区域提取方法 | |
Pan et al. | A generic FCN-based approach for the road-network extraction from VHR remote sensing images–using openstreetmap as benchmarks | |
CN103034982A (zh) | 一种基于变焦视频序列的图像超分辨率重建方法 | |
CN104867137A (zh) | 一种基于改进ransac算法的图像配准方法 | |
CN107704837A (zh) | 道路网拓扑和几何信息的提取方法 | |
CN104835186B (zh) | 一种历史地震图波形自动跟踪方法 | |
CN106600965B (zh) | 基于尖锐度的交通流早晚高峰时段自动识别方法 | |
CN105005995A (zh) | 一种计算三维点云模型骨骼的方法 | |
CN107179543A (zh) | 对叠前数据进行规则化的方法和装置 | |
CN102136133B (zh) | 一种图像处理方法及图像处理装置 | |
Zhang et al. | Enhanced interpreter-aided salt boundary extraction using shape deformation | |
Kahraman et al. | Road detection from high satellite images using neural networks | |
US11636649B2 (en) | Geospatial modeling system providing 3D geospatial model update based upon predictively registered image and related methods | |
CN102867174B (zh) | 一种人脸特征定位方法及装置 | |
Seo et al. | Ortho-image analysis for producing lane-level highway maps | |
Lei et al. | Pedestrian lane detection for assistive navigation of vision-impaired people: Survey and experimental evaluation | |
CN109781126A (zh) | 一种基于矢量法的全局静态路径规划方法 | |
CN113633375A (zh) | 一种虚拟支气管镜的构建方法 | |
Zhang et al. | Road topology extraction from satellite imagery by joint learning of nodes and their connectivity | |
CN106796725A (zh) | 一种血管脊线追踪方法及装置 | |
US20220215619A1 (en) | Geospatial modeling system providing 3d geospatial model update based upon iterative predictive image registration and related methods | |
CN115619900A (zh) | 基于距离地图和概率路图的点云地图拓扑结构提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171110 |