CN102722640B - 一种顾及邻近波形信息的机载激光波形数据分解算法 - Google Patents
一种顾及邻近波形信息的机载激光波形数据分解算法 Download PDFInfo
- Publication number
- CN102722640B CN102722640B CN201210162472.8A CN201210162472A CN102722640B CN 102722640 B CN102722640 B CN 102722640B CN 201210162472 A CN201210162472 A CN 201210162472A CN 102722640 B CN102722640 B CN 102722640B
- Authority
- CN
- China
- Prior art keywords
- waveform
- pulse crest
- crest
- pulse
- candidate
- 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
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种顾及邻近波形信息的机载激光波形数据分解方法,包括A、波形数据的预处理;B、波形数据的脉冲波峰探测;C、基于广义高斯模型的显著脉冲波峰提取;D、顾及邻近波形的候选脉冲波峰提取;E、基于区域生长策略的待分析波形列表生成;F、波形分解结果输出。本发明采用高、低双阈值,可以克服单阈值难以合适选取的缺陷,有效考虑波形数据中的微弱脉冲信号;波形分解过程中,通过区域生长的方式,不断利用邻近波形显著脉冲波峰所提供的可靠信息,对待分析波形的候选脉冲波峰进行分析,可以实现这类微弱脉冲波峰信息的稳健提取,获取更加完整的点云信息。
Description
技术领域
本发明涉及测绘科学与技术领域,涉及一种顾及邻近波形信息的机载激光波形数据分解算法,该方法能对波形数据中微弱的脉冲信号进行有效检测,以获取更加完整的点云,非常适用于植被覆盖严重的地区。
背景技术
机载LiDAR(Light Detection And Ranging),又称机载激光扫描技术,作为一种三维空间信息的实时获取手段,可以快速、精确的获取大范围区域的地表信息。根据回波记录方式的不同,机载激光雷达系统可以分为离散激光雷达系统和全波形激光雷达系统。前者记录有限个离散的回波信号,而后者以很小的采样间隔对来自目标的激光反射信号进行采样记录,形成一个随时间变化的回波信号。相对于离散激光雷达而言,全波形雷达可以提供更多的细节信息,具有更强的地物区分能力和更优的植被穿透特性等。同时,也给数据处理和信息提取提出了更高的要求。
如何从波形采样数据中提取出高质量的点云是数据处理中的一个关键问题。Hofton等提出机载激光雷达波形数据高斯分解算法,Wagner等从理论角度阐述了机载全波形激光雷达的成像机理,并使用高斯模型提取波形数据峰值点。这些方法大多数情况下可以获取较好结果,但是由于各种干扰因素,回波信息并不是高斯函数的精确表达。Chauve等提出利用广义高斯模型来拟合波形数据,并通过迭代的方式实现波形数据的分解。马洪超等采用Expectation-Maximization算法(EM)进行机载雷达波形数据高斯混合模型参数的最大似然估计,取得较好结果但计算量较大。通过分析发现,目前算法中主要存在如下几点缺陷:
第一,噪声阈值难以确定。目前,为排除信号噪声的影响,大多算法采用经验取值的方法,并指定为较高的数值。
第二,微弱信号难以可靠检测。当利用较高阈值对波形数据进行预处理后,一些微弱脉冲信号由于信号强度不够而无法得到考虑,从而不能进行检测。这种现象在植被茂密地区尤为严重,由于植被遮挡造成地面点反射脉冲信号比较微弱,给地面信息的提取造成较大挑战。
发明内容
本发明的目的是在于针对现有技术存在的上述问题,提供一种顾及邻近波形信息的机载激光波形数据分解方法,利用高、低双阈值数据预处理,同时通过顾及邻近波形显著脉冲波峰信号所提供的信息,对比较微弱的脉冲信息进行分析,融入区域生长策略,最终实现波形数据中脉冲波峰信号的完整提取,可以得到更加精确、完整的点云数据。
本发明的上述目的通过以下技术方案实现:
一种顾及邻近波形信息的机载激光波形数据分解方法,包括以下步骤:
步骤1、波形数据的预处理:利用一维高斯模板对原始波形数据进行平滑处理;利用一阶微分算子计算经过平滑处理后的波形数据的一阶导数;计算每一个波形数据对应的近似物方平面坐标,并根据近似物方平面坐标建立k-d树索引;
步骤2、波形数据的脉冲波峰探测:通过计算步骤1中得到的一阶导数零交叉点,确定波形的脉冲波峰中心的初始位置,采用高、低双阈值对波峰强度进行分析,确定出显著脉冲波峰和候选脉冲波峰,并剔除伪脉冲波峰;
步骤3、基于广义高斯模型的显著脉冲波峰提取:对当前所确定的显著脉冲波峰,采用广义高斯模型和非线性最小二乘Levenburg-Marquardt算法进行波形数据的拟合,提取显著脉冲波峰的参数,并统计波形数据拟合的单位权中误差;
步骤4、顾及邻近波形的候选脉冲波峰提取:以含有候选脉冲波峰的波形为中心,通过k-d树索引搜索设定半径范围内的邻近波形,并计算邻近波形的显著脉冲波峰提供的高程值与设置为中心的波形的候选脉冲波峰的高程值的高程差,若高程差小于给定阈值,则认为该候选脉冲波峰与邻域环境相容,将该候选脉冲波峰重新认定为显著脉冲波峰,采用广义高斯模型和非线性最小二乘Levenburg-Marquardt算法进行重新认定的显著脉冲波峰的参数的提取,并将重新认定为显著脉冲波峰的对应的波形加入到已更新波形列表中;
步骤5、基于区域生长策略的待分析波形列表生成:以已更新波形列表中的波形为中心,通过k-d树搜索设定半径范围内的邻近波形,如果邻近波形中含有候选脉冲波峰,将该邻近波形加入到待分析波形列表中,当待分析波形列表不为空时,返回步骤4;当待分析波形列表为空时,则进入步骤6;
步骤6、波形分解结果输出:利用所获取的显著脉冲波峰的参数,计算其对应激光点的三维空间坐标,输出点云。
如上所述的步骤1中一维高斯模板,其连续情况下的数字形式如下:
其中,σ为高斯模板的标准偏差,π为圆周率。
如上所述的σ是基于以下公式:
其中,FWHM为激光器波形的半宽参数。
如上所述的步骤1中建立k-d树索引包括以下步骤:
步骤4.1、计算波形对应的近似物方平面坐标,第k个波形近似物方平面坐标的计算,采用如下公式:
其中,E0,N0为第k个波形第一个采样点所对应的平面坐标,dE,dN为第k个波形的平面坐标微分,WFOFFSET为第k个波形的回波数据第一个采样点到发射波形第一次采样点的偏移量;
步骤4.2、利用波形对应的近似物方平面坐标建立k-d树索引。
如上所述的步骤2中采用高、低双阈值对波峰强度进行分析,确定出显著脉冲波峰和候选脉冲波峰,并剔除伪脉冲波峰包括以下步骤:
设定高阈值为HighThred,低阈值为LowThred,脉冲波峰位置i处信号强度为WVal[i];
当WVal[i]>HighThred时,该脉冲波峰为显著脉冲波峰;
当LowThred≤WVal[i]≤HighThred时,该脉冲波峰为候选脉冲波峰,将候选脉冲波峰所在的波形作为待分析波形加入待分析波形列表中;
当WVal[i]<LowThred时,该脉冲波峰为伪脉冲波峰,剔除伪脉冲波峰;
若波形中候选脉冲波峰个数不为零,则将该波形加入待分析波形列表中。
如上所述的步骤4中高程值是基于以下公式:
H=H0+dH×(WFOFFSET+t)
其中,H为高程值,H0为显著脉冲波峰/候选脉冲波峰所在波形第一个采样点所对应的高程坐标,dH为波形数据单元的高程坐标微分,WFOFFSET为回波数据第一个采样点到发射波形第一次采样点的偏移量,t为显著脉冲波峰/候选脉冲波峰对应的位置。
如上所述的步骤6中点云的计算基于以下公式:
E=E0+dE×(WFOFFSET+t)
N=N0+dN×(WFOFFSET+t)
H=H0+dH×(WFOFFSET+t)
其中,E、N、H为波形点云三维坐标,E0、N0、H0为波形第一次采样点所对应的平面坐标,dE、dN、dH为波形的平面坐标微分,WFOFFSET为回波数据第一次采样点到发射波形第一次采样点的偏差量,t为显著脉冲波峰对应的位置。
本发明与现有技术相比,具有以下优点:
1、采用高、低双阈值,可以克服单阈值难以合适选取的缺陷,有效考虑波形数据中的微弱脉冲信号;
2、波形分解过程中,通过区域生长的方式,不断利用邻近波形显著脉冲波峰所提供的可靠信息,对待分析波形的候选脉冲波峰进行分析,可以实现这类微弱脉冲波峰信息的稳健提取,获取更加完整的点云信息。
3、利用该算法提取的点云相比于相关商业软件生成的点云,提取的有效地面点云信息可达到30%以上,成效显著。
附图说明
图1是本发明的流程示意图。
具体实施方式
以下结合附图对本发明的技术方案作进一步详细描述。
实施例1
如图1,一种顾及邻近波形信息的机载激光波形数据分解方法,包括以下步骤:
步骤1、波形数据的预处理:利用一维高斯模板对原始波形数据进行平滑处理;利用一阶微分算子计算经过平滑处理后的波形数据的一阶导数;计算每一个波形数据对应的近似物方平面坐标,并根据近似物方平面坐标建立k-d树索引;
基于高斯模板的波形数据平滑处理:利用数字图像处理中常用的一维高斯模板,其连续情况下的数字形式如下:
其中,σ为高斯模板的标准偏差,π为圆周率。σ与信号滤波平滑的程度密切相关,σ越大,高斯滤波器的频带就越宽,平滑的程度就越强;反之,平滑的程度就越弱。为了在平滑噪声的同时尽量保留信号特征,本实施例中σ的选取根据激光器波形的半宽参数(full width at half maximum,FWHM)确定,如果激光器波形的半宽参数厂商没有提供,可以利用平地处所反射激光波形的半宽参数来代替。当得到半宽参数FWHM后,依据公式(2)可以计算得到高斯模板的标准偏差σ:
高斯滤波模板窗口的半径取为3σ,按照一定间隔,对连续的一维高斯函数进行离散化和归一化处理,即可获取所需的高斯模板。之后,利用高斯模板对原始波形数据进行卷积运算,即可实现波形数据中噪声的有效抑制。
计算波形数据的一阶导数值:假设高斯滤波后的波形数据为WVal[i](i=0,…,m),其中m为波形数据的采样数,对于位置i处波形数据对应的一阶导数Grad[i],采用如下公式计算:
Grad[i]=WVal[i+1]-WVal[i] (3)
建立波形数据的k-d树索引:首先计算波形对应的近似物方平面坐标,第k个波形近似物方平面坐标()的计算,采用如下公式:
式中,(E0,N0)为第k个波形第一个采样点所对应的平面坐标,(dE,dN)为第k个波形的平面坐标微分,WFOFFSET为第k个波形的回波数据第一个采样点到发射波形第一次采样点的偏移量,这些参数均可以从后缀名为lgc的文件中读取。
然后利用波形对应的平面坐标建立k-d树索引,k-d树是由Bentley于1975年提出并逐渐发展成为一种比较好的多维空间索引,特别适合于空间点状目标。它通过超平面将一个空间递归划分为两个子空间来实现搜索,其实现算法比较稳定和成熟,本专利中采用了马里兰大学David M.Mount和香港理工大学Sunil Arya提供的最近邻算法软件库(Approximate Nearest Neighbor,ANN)实现,该软件包提供了最临近查询和范围查询两种方式,其中最临近查询方式为查询给定点周围若干最临近点,范围查询方式为查询给定点给定半径大小范围内的所有点。至此,可以利用k-d树索引和波形对应的近似平面坐标完成邻域波形的查询。
步骤2、波形数据的脉冲波峰探测:通过计算步骤1中得到的一阶导数零交叉点,确定波形的脉冲波峰中心的初始位置,采用高、低双阈值对波峰强度进行分析,确定出显著脉冲波峰和候选脉冲波峰,并剔除伪脉冲波峰;主要包括以下步骤:
设定高阈值为HighThred,低阈值为LowThred,脉冲波峰位置i处信号强度为WVal[i];
当WVal[i]>HighThred时,该脉冲波峰为显著脉冲波峰;
当LowThred≤WVal[i]≤HighThred时,该脉冲波峰为候选脉冲波峰,将候选脉冲波峰所在的波形作为待分析波形加入待分析波形列表中;
当WVal[i]<LowThred时,该脉冲波峰为伪脉冲波峰,剔除伪脉冲波峰;
若波形中候选脉冲波峰个数不为零,则将该波形加入待分析波形列表中。遍历每一个波形数据,并进行上述处理。
步骤3、基于广义高斯模型的显著脉冲波峰提取:对当前所确定的显著脉冲波峰,采用广义高斯模型和非线性最小二乘Levenburg-Marquardt算法进行波形数据的拟合,提取显著脉冲波峰的参数,并统计波形数据拟合的单位权中误差;
当波形的显著脉冲波峰个数不为零时,利用当前确定的显著脉冲波峰的近似参数,通过广义高斯模型和非线性最小二乘算法Levenberg-Marquardt对波形进行分解。
波形分解的目的就是要将波形分解成多个脉冲信号的集合,通过记录各个脉冲信号的特征来描述不同地物的位置和反射强度等信息。本实施例中将每个脉冲信号视为广义高斯函数,则波形的形式如下:
式(5)中,y为波形值,n代表脉冲信号个数,取值为当前所确定的显著脉冲波峰个数,fi(x)为广义高斯函数,其具体表达形式如下:
式(6)中,Pi为脉冲波形的振幅,ti为脉冲波形对应波峰的位置,si为脉冲波形的标准方差,ri为形状参数。容易发现,当时,相当于标准的高斯函数;广义高斯模型通过在标准高斯函数的基础上增加了形状参数ri,达到改变波形形状的目的,从而可以更好的拟合波形的形状。
在建立波形分解的数学模型后,波形分解的过程其实就是各个脉冲波形参数的优化计算过程,具体采用Levenberg-Marquardt算法实现。该算法是数学中最常见的优化算法,是一种介于牛顿法与梯度下降法之间的一种非线性优化算法,其优点在于对过参数化问题不敏感,能够有效处理冗余参数带来的强相关问题,是迭代优化陷于局部极小值的机会大大减少。
利用Levenberg-Marquardt算法实现波形数据分解的步骤如下:
1、参数初始值的确定:对于广义高斯函数参数,在首次迭代中,ti为利用波形一阶导数的零交叉点所确定的波峰位置,Pi为波峰位置处信号强度值WVal[ti],缺省认为脉冲为标准高斯函数形状,故ri取为此外,脉冲信号半宽参数FWHM一般为3,故si取为1.5;在随后的迭代中,参数初始值取为前一次迭代平差所计算的新参数值。对于Levenberg-Marquardt算法中所使用的变量,k为迭代次数,λk为阻尼系数,v为阻尼系数调整常数,本实施例中k初始化为0,λk初始值λ0设为0.01,v取值为10。
2、计算拟合函数的雅可比矩阵Jk:
式(7)中,n为当前探测出的显著波峰个数,m为所使用的波形数据点个数,pk=[P1 t1 s1 r1 … Pn tn sn rn]T为待解算的未知数,为各脉冲信号的波形参数。
然后,计算海塞矩阵H:
并构造增量正规方程
式中,δk表示待解算的未知数pk的增量向量,εk表示利用未知数当前值所计算的波形值与其对应观测值的误差向量。利用公式(9)进行求解,可以得到δk。
3、阻尼系数调整:
3.1、利用δk更新参数向量并计算新的误差向量,如果||y-f(pk+δk)||<||εk||(||·||表示向量的L2-范数),即误差向量变小,则令pk+1=pk+δk;如果||εk||小于为给定收敛阈值,满足收敛条件,则停止迭代,输出结果,否则令λk+1=λk/v,转到步骤2中进行迭代。
3.2、如果||y-f(pk+δk)||>||εk||,则令λk+1=λk·v,返回到步骤1重新解算。
计算拟合中误差σ:当迭代满足收敛条件后,利用最终得到的误差向量εk计算拟合中误差
步骤4、顾及邻近波形的候选脉冲波峰提取:以含有候选脉冲波峰的波形为中心,通过k-d树索引搜索设定半径范围内的邻近波形,并计算邻近波形的显著脉冲波峰提供的高程值与设置为中心的波形的候选脉冲波峰的高程值的高程差,若高程差小于给定阈值,则认为该候选脉冲波峰与邻域环境相容,将该候选脉冲波峰重新认定为显著脉冲波峰,采用广义高斯模型和非线性最小二乘Levenburg-Marquardt算法进行重新认定的显著脉冲波峰的参数的提取,并将重新认定为显著脉冲波峰的对应的波形加入到已更新波形列表中;
首先,利用步骤1中波形数据预处理中所建立的k-d树索引,和待分析波形的近似平面位置,查询待分析波形一定半径大小范围内的所有波形数据,半径取值为点云平均距离的3倍。然后,利用如下公式计算邻近波形中显著脉冲波峰所对应的地面高程H:
H=H0+dH×(WFOFFSET+t) (10)
式中,H0为显著脉冲波峰所在波形第一个采样点所对应的高程坐标,dH为波形数据单元的高程坐标微分,WFOFFSET为回波数据第一个采样点到发射波形第一次采样点的偏移量,t为显著脉冲波峰对应的位置,这些参数均可以从后缀名为lgc文件获取得到。
同样,采用公式(10),计算待分解波形中候选脉冲波峰对应的高程值,并将其与邻近波形中显著脉冲波峰对应的高程值求差,如果差值小于给定高差阈值,则认为该候选脉冲与周围环境是相容的,并将其确定为显著脉冲波峰。
如果有候选脉冲波峰被认定为显著脉冲波峰,则重新对该波形进行基于广义高斯模型的波形分解,精确提取当前显著脉冲波峰的参数信息。波形分解成功后,将该波形加入到已更新波形列表中。
步骤5、基于区域生长策略的待分析波形列表生成:以已更新波形列表中的波形为中心,通过k-d树搜索设定半径范围内的邻近波形,如果邻近波形中含有候选脉冲波峰,将该邻近波形加入到待分析波形列表中,当待分析波形列表不为空时,返回步骤4;当待分析波形列表为空时,则进入步骤6;
对于已更新波形列表,由于增加了新的显著脉冲波峰,将对其所在邻域中含有候选脉冲波峰的波形数据造成影响。为此,对已更新波形列表中的波形,逐一通过前面建立的k-d树索引,搜索给定半径大小范围内的邻近波形,半径取值为点云平均距离的3倍。如果邻近波形中存在候选脉冲波峰,将该波形加入到待分析波形列表中。
当得到待分析波形列表后,采用步骤4中的顾及邻近波形的候选脉冲波峰提取方法,逐一对待分析波形的候选脉冲波峰进行分析和探测。波形分解完成后,重新生成已更新波形列表。同样,通过建立的k-d树索引,根据新产生的已更新波形列表,生成新的待分析波形列表。可以发现,这种方式非常类似于区域生长策略,仅不断考虑有改变的波形的邻域,这样不仅减少了计算量,而且实现候选脉冲波峰的递进提取。
步骤6、波形分解结果输出:利用所获取的显著脉冲波峰的参数,计算其对应激光点的三维空间坐标,输出点云。
点云的计算基于以下公式:
E=E0+dE×(WFOFFSET+t)
N=N0+dN×(WFOFFSET+t)
H=H0+dH×(WFOFFSET+t)
其中,E、N、H为波形点云三维坐标,E0、N0、H0为波形第一次采样点所对应的平面坐标,dE、dN、dH为波形的平面坐标微分,WFOFFSET为回波数据第一次采样点到发射波形第一次采样点的偏差量,t为显著脉冲波峰对应的位置。此外,根据实际需要,可以输出激光点所对应脉冲波峰的强度、波形宽度等信息。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (6)
1.一种顾及邻近波形信息的机载激光波形数据分解方法,其特征在于:包括以下步骤:
步骤1、波形数据的预处理:利用一维高斯模板对原始波形数据进行平滑处理;利用一阶微分算子计算经过平滑处理后的波形数据的一阶导数;计算每一个波形数据对应的近似物方平面坐标,并根据近似物方平面坐标建立k-d树索引;
步骤2、波形数据的脉冲波峰探测:通过计算步骤1中得到的一阶导数零交叉点,确定波形的脉冲波峰中心的初始位置,采用高、低双阈值对波峰强度进行分析,确定出显著脉冲波峰和候选脉冲波峰,并剔除伪脉冲波峰;
步骤3、基于广义高斯模型的显著脉冲波峰提取:对当前所确定的显著脉冲波峰,采用广义高斯模型和非线性最小二乘Levenburg-Marquardt算法进行波形数据的拟合,提取显著脉冲波峰的参数,并统计波形数据拟合的单位权中误差;
步骤4、顾及邻近波形的候选脉冲波峰提取:以含有候选脉冲波峰的波形为中心,通过k-d树索引搜索设定半径范围内的邻近波形,并计算邻近波形的显著脉冲波峰提供的高程值与设置为中心的波形的候选脉冲波峰的高程值的高程差,若高程差小于给定阈值,则认为该候选脉冲波峰与邻域环境相容,将该候选脉冲波峰重新认定为显著脉冲波峰,采用广义高斯模型和非线性最小二乘Levenburg-Marquardt算法进行重新认定的显著脉冲波峰的参数的提取,并将重新认定为显著脉冲波峰的对应的波形加入到已更新波形列表中;
步骤5、基于区域生长策略的待分析波形列表生成:以已更新波形列表中的波形为中心,通过k-d树搜索设定半径范围内的邻近波形,如果邻近波形中含有候选脉冲波峰,将该邻近波形加入到待分析波形列表中,当待分析波形列表不为空时,返回步骤4;当待分析波形列表为空时,则进入步骤6;
步骤6、波形分解结果输出:利用所获取的显著脉冲波峰的参数,计算其对应激光点的三维空间坐标,输出点云;
所述的步骤2中采用高、低双阈值对波峰强度进行分析,确定出显著脉冲波峰和候选脉冲波峰,并剔除伪脉冲波峰包括以下步骤:
设定高阈值为HighThred,低阈值为LowThred,脉冲波峰位置i处信号强度为WVal[i];
当WVal[i]>HighThred时,该脉冲波峰为显著脉冲波峰;
当LowThred≤WVal[i]≤HighThred时,该脉冲波峰为候选脉冲波峰,将候选脉冲波峰所在的波形作为待分析波形加入待分析波形列表中;
当WVal[i]<LowThred时,该脉冲波峰为伪脉冲波峰,剔除伪脉冲波峰;
若波形中候选脉冲波峰个数不为零,则将候选脉冲波峰所在的波形加入待分析波形列表中。
2.根据权利要求1所述的一种顾及邻近波形信息的机载激光波形数据分解方法,其特征在于,所述的步骤1中一维高斯模板,其连续情况下的数字形式如下:
其中,σ为高斯模板的标准偏差,π为圆周率。
3.根据权利要求2所述的一种顾及邻近波形信息的机载激光波形数据分解方法,其特征在于,所述的σ是基于以下公式:
其中,FWHM为激光器波形的半宽参数。
4.根据权利要求1所述的一种顾及邻近波形信息的机载激光波形数据分解方法,其特征在于,所述的步骤1中建立k-d树索引包括以下步骤:
步骤4.1、计算波形对应的近似物方平面坐标,第k个波形近似物方平面坐标的计算,采用如下公式:
其中,E0,N0为第k个波形第一个采样点所对应的平面坐标,dE,dN为第k个波形的平面坐标微分,WFOFFSET为第k个波形的回波数据第一个采样点到发射波形第一次采样点的偏移量;
步骤4.2、利用波形对应的近似物方平面坐标建立k-d树索引。
5.根据权利要求1所述的一种顾及邻近波形信息的机载激光波形数据分解方法,其特征在于,所述的步骤4中高程值是基于以下公式:
H=H0+dH×(WFOFFSET+t)
其中,H为高程值,H0为显著脉冲波峰/候选脉冲波峰所在波形第一个采样点所对应的高程坐标,dH为波形数据单元的高程坐标微分,WFOFFSET为回波数据第一个采样点到发射波形第一次采样点的偏移量,t为显著脉冲波峰/候选脉冲波峰对应的位置。
6.根据权利要求1所述的一种顾及邻近波形信息的机载激光波形数据分解方法,其特征在于,所述的步骤6中点云的计算基于以下公式:
E=E0+dE×(WFOFFSET+t)
N=N0+dN×(WFOFFSET+t)
H=H0+dH×(WFOFFSET+t)
其中,E、N、H为波形点云三维坐标,E0、N0、H0为波形第一次采样点所对应的平面坐标,dE、dN、dH为波形的平面坐标微分,WFOFFSET为回波数据第一次采样点到发射波形第一次采样点的偏差量,t为显著脉冲波峰对应的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210162472.8A CN102722640B (zh) | 2012-05-24 | 2012-05-24 | 一种顾及邻近波形信息的机载激光波形数据分解算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210162472.8A CN102722640B (zh) | 2012-05-24 | 2012-05-24 | 一种顾及邻近波形信息的机载激光波形数据分解算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102722640A CN102722640A (zh) | 2012-10-10 |
CN102722640B true CN102722640B (zh) | 2014-10-15 |
Family
ID=46948397
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210162472.8A Active CN102722640B (zh) | 2012-05-24 | 2012-05-24 | 一种顾及邻近波形信息的机载激光波形数据分解算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102722640B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104639277B (zh) * | 2015-03-14 | 2017-07-28 | 河北工业大学 | 一种高斯光谱半宽复用通信系统的解调方法 |
CN106525000B (zh) * | 2016-10-31 | 2017-07-11 | 武汉大学 | 基于激光扫描离散点强度梯度的道路标线自动化提取方法 |
CN106990401B (zh) * | 2017-05-24 | 2019-12-10 | 武汉大学 | 基于全波形机载激光雷达数据二类高程误差修正方法 |
CN109872281A (zh) * | 2018-12-29 | 2019-06-11 | 中国科学院遥感与数字地球研究所 | 一种波形信息辅助下的渐进加密三角网LiDAR点云滤波方法 |
CN111551952B (zh) * | 2020-05-08 | 2022-10-18 | 中国地质大学(武汉) | 一种水深测量激光数据的提取方法、水深测量方法及装置 |
CN112711001B (zh) * | 2020-12-17 | 2023-08-01 | 中国科学院空天信息创新研究院 | 一种精细去噪辅助的激光雷达波形分解方法 |
CN115345208B (zh) * | 2022-10-19 | 2023-02-03 | 成都理工大学 | 基于顶帽变换的中子-伽马脉冲堆积甄别方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101196562A (zh) * | 2007-12-14 | 2008-06-11 | 武汉大学 | 一种基于改进的em算法的激光雷达波形数据分解的方法 |
US7656526B1 (en) * | 2006-07-21 | 2010-02-02 | University Corporation For Atmospheric Research | Lidar system for remote determination of calibrated, absolute aerosol backscatter coefficients |
CN101908068A (zh) * | 2010-08-03 | 2010-12-08 | 武汉大学 | 一种基于四叉树索引的海量激光扫描点云实时绘制方法 |
CN102393180A (zh) * | 2011-10-19 | 2012-03-28 | 中国林业科学研究院资源信息研究所 | 一种从LiDAR点云数据中自动提取林分上层树木参数的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8144048B2 (en) * | 2009-03-25 | 2012-03-27 | Honeywell International Inc. | Systems and methods for gaussian decomposition of weather radar data for communication |
-
2012
- 2012-05-24 CN CN201210162472.8A patent/CN102722640B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7656526B1 (en) * | 2006-07-21 | 2010-02-02 | University Corporation For Atmospheric Research | Lidar system for remote determination of calibrated, absolute aerosol backscatter coefficients |
CN101196562A (zh) * | 2007-12-14 | 2008-06-11 | 武汉大学 | 一种基于改进的em算法的激光雷达波形数据分解的方法 |
CN101908068A (zh) * | 2010-08-03 | 2010-12-08 | 武汉大学 | 一种基于四叉树索引的海量激光扫描点云实时绘制方法 |
CN102393180A (zh) * | 2011-10-19 | 2012-03-28 | 中国林业科学研究院资源信息研究所 | 一种从LiDAR点云数据中自动提取林分上层树木参数的方法 |
Non-Patent Citations (16)
Title |
---|
A marked Point Process for Modeling Lidar Waveforms;clement mallet,et al;《IEEE TRANSACTIONS ON IMAGE PROCESSING》;20101231;第19卷(第12期);3204~3221 * |
BretarF et al.MANAGING FULL WAVEFORM LIDAR DATA:A CHALLENGING TASK FOR THE FORTHCOMING YEARS.《International Archives of Photogrammetry and Remote Sensing |
ChauveA et al.PROCESSING FULL-WAVEFORM LIDAR DATA:MODELLING RAW SIGNALS.《International Archives of Photogrammetry |
clement mallet,et al.A marked Point Process for Modeling Lidar Waveforms.《IEEE TRANSACTIONS ON IMAGE PROCESSING》.2010,第19卷(第12期), |
Decomposition of laser altimeter waveforms;Hofton, M.A.,et al;《IEEE Transactions on Geoscience and Remote Sensing》;20000731;第38卷(第4期);1989~1996 * |
Gaussian decomposition and calibration of a novel small-footprint full-waveform digitising airborne laser scanner;Wolfgang Wagnera,et al;《ISPRS Journal of Photogrammetry and Remote Sensing》;20060430;第60卷(第2期);100~112 * |
Hofton M.A. |
MANAGING FULL WAVEFORM LIDAR DATA:A CHALLENGING TASK FOR THE FORTHCOMING YEARS;Bretar F,et al;《International Archives of Photogrammetry and Remote Sensing,2008》;20080514;1~5 * |
PROCESSING FULL-WAVEFORM LIDAR DATA:MODELLING RAW SIGNALS;Chauve A,et al;《International Archives of Photogrammetry, Remote Sensing and Spatial Information Sciences 2007》;20071231;102~107 * |
Wagner W,et al.WAVEFORM ANALYSIS TECHNIQUES IN AIRBORNE LASER SCANNING.《International Archives of Photogrammetry and Remote Sensing》.2007,第36卷(第3期), |
WAVEFORM ANALYSIS TECHNIQUES IN AIRBORNE LASER SCANNING;Wagner W,et al;《International Archives of Photogrammetry and Remote Sensing》;20071231;第36卷(第3期);413~418 * |
Wolfgang Wagnera,et al.Gaussian decomposition and calibration of a novel small-footprint full-waveform digitising airborne laser scanner.《ISPRS Journal of Photogrammetry and Remote Sensing》.2006,第60卷(第2期), |
吴华意等.机载激光雷达系统的应用与数据后处理技术.《测绘与空间地理信息》.2006,第29卷(第3期), |
徐光彩.机载LIDAR波形数据处理及分类研究.《中国优秀硕士学位论文全文数据库基础科学辑(月刊)》.2011,(第5期), * |
机载激光雷达系统的应用与数据后处理技术;吴华意等;《测绘与空间地理信息》;20060630;第29卷(第3期);58~63 * |
覃驭楚等.小光斑激光雷达全波形数据递进分解与相对辐射校正.《中国科学:地球科学》.2011,第41卷(第1期), * |
Also Published As
Publication number | Publication date |
---|---|
CN102722640A (zh) | 2012-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102722640B (zh) | 一种顾及邻近波形信息的机载激光波形数据分解算法 | |
Zhou et al. | Gaussian inflection point selection for LiDAR hidden echo signal decomposition | |
CN103236063B (zh) | 基于多尺度谱聚类及决策级融合的sar图像溢油检测方法 | |
CN101196562B (zh) | 一种基于改进的em算法的激光雷达波形数据分解的方法 | |
CN106500671B (zh) | 一种基于lm算法分解激光雷达波形确定海水深度的方法 | |
Waters et al. | Data assimilation of partitioned HF radar wave data into Wavewatch III | |
CN103941254A (zh) | 一种基于地质雷达的土壤物性类别识别方法和装置 | |
CN103217679A (zh) | 一种基于遗传算法的全波形激光雷达回波数据的高斯分解法 | |
CN104614718A (zh) | 基于粒子群算法的激光雷达波形数据分解的方法 | |
CN109709527B (zh) | 一种全波形激光测高回波信号中高斯分解的高斯波峰法 | |
CN105093203A (zh) | 星体次表层遥感探测雷达回波模拟及参数反演方法 | |
Rossi et al. | Coupled inverse modeling of a controlled irrigation experiment using multiple hydro-geophysical data | |
CN105676205A (zh) | 一种机载LiDAR波形数据高斯分解方法 | |
JP5719075B1 (ja) | 空洞厚探査方法 | |
CN105259537A (zh) | 基于频移迭代的多普勒谱中心频率估计方法 | |
CN105527617A (zh) | 一种基于鲁棒主元分析的探地雷达数据背景去除方法 | |
JP5629840B1 (ja) | 空洞厚探査方法 | |
KR101426226B1 (ko) | 레이더의 신호 처리 방법 | |
CN110135299B (zh) | 一种用于浅水测深的单波段蓝绿激光波形分析方法及系统 | |
CN115422981A (zh) | 面向单频机载激光测深数据的水陆分类方法、系统及应用 | |
Luo et al. | Accurate tree roots positioning and sizing over undulated ground surfaces by common offset GPR measurements | |
Li et al. | A fuzzy mean-shift approach to lidar waveform decomposition | |
CN105738881A (zh) | 基于波形分类的近海Envisat卫星测高波形重定方法 | |
Zhao et al. | A Comprehensive Horizon‐Picking Method on Subbottom Profiles by Combining Envelope, Phase Attributes, and Texture Analysis | |
Xiong et al. | Automated reconstruction of subsurface interfaces in Promethei Lingula near the Martian south pole by using SHARAD data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |