CN105137411B - 基于地物散射的全波LiDAR波形分解方法 - Google Patents
基于地物散射的全波LiDAR波形分解方法 Download PDFInfo
- Publication number
- CN105137411B CN105137411B CN201510564627.4A CN201510564627A CN105137411B CN 105137411 B CN105137411 B CN 105137411B CN 201510564627 A CN201510564627 A CN 201510564627A CN 105137411 B CN105137411 B CN 105137411B
- Authority
- CN
- China
- Prior art keywords
- scattering
- point
- wave lidar
- effective
- scattering unit
- 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
Classifications
-
- 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/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于地物散射的全波LiDAR波形分解方法,包括:激光雷达发射信号到地物目标上反馈回来的信号叠加形成全波LiDAR波形;采用有效散射单元来描述全波LiDAR波形:将全波LiDAR波形转化为各采样点对应的具有相同散射能量的有效散射单元;采用阶梯函数对全波LiDAR波形中的有效散射单元个数建模;利用RJMCMC算法对全波LiDAR波形进行分解,得到地物目标相对高程以及散射特性。本发明针对采用有效散射单元来描述全波LiDAR波形,将全波LiDAR波形转化为各采样点对应的具有相同散射能量的有效散射单元,能够将地物散射特性与波形分解相结合,并提出采用阶梯函数对有效散射单元个数建模,打破了使用传统高斯统计分布模型来拟合原始波形的思想。
Description
技术领域
本发明所属全波LiDAR数据处理技术领域,具体是一种基于地物散射的全波LiDAR波形分解方法。
背景技术
波形分解是全波LiDAR数据研究的主要方面,是全波LiDAR数据应用的前期工作,因此波形分解成为全波LiDAR数据处理领域研究的热点及难点问题。波形分解是指用统计分布模型拟合原始波形数据的过程。
LiDAR波形分解包括两个关键问题:(1)实现波形数据的最佳分解;(2)分解结果能够反映被测目标的散射特性。目前波形分解方法中多数采用高斯分布来拟合原始波形数据,尽管高斯分解模型具有模型简单,拟合容易等优点,但在实际激光反射中,当激光光束照到较为复杂的地物目标时,波形会发生倾斜、拖尾等形变,在这种情况下,采用高斯函数拟合原始波形数据将产生较大的拟合误差。广义高斯模型、Nakagami和Burr等统计模型也常用来拟合原始波形数据。
但这些波形分解方法仅仅考虑波形拟合,并没有从地物散射机理上考虑波形的形成,不能充分的反映被测目标的散射特性。并且,目前的波形分解方法仅仅针对波形数据本身,并没有考虑研究区域内地物目标及其散射特征等特性,因此难于挖掘蕴含于波形数据中的地物信息。
发明内容
针对现有技术存在的不足,本发明提供一种基于地物散射的全波LiDAR波形分解方法。
本发明的技术方案是:
基于地物散射的全波LiDAR波形分解方法,包括:
步骤1、激光雷达发射信号到地物目标上反馈回来的信号叠加形成全波LiDAR波形,包括分布在有效高程区域内的采样点及各采样点处的散射能量;
步骤2、采用有效散射单元来描述全波LiDAR波形:将全波LiDAR波形转化为各采样点对应的具有相同散射能量的有效散射单元;
步骤3、采用阶梯函数对全波LiDAR波形中的有效散射单元个数建模:用阶梯的高度表示各个高程点对应的有效散射单元个数的均值,阶梯的位置表示地物目标的相对高程;
步骤4、利用RJMCMC(可逆跳转马尔科夫链蒙特卡罗)算法对全波LiDAR波形进行分解,即求解阶梯函数,得到地物目标相对高程以及散射特性。
所述步骤2的具体步骤如下:
步骤2-1、设各采样点处的散射能量中的散射能量最小值c=min(yi)为一个有效散射单元;
步骤2-2、计算每个采样点对应的有效散射单元个数zi=yi/c,i=1,2,...,n;即采样点处的散射能量yi除以有效散射单元的散射能量c;每个采样点对应的有效散射单元个数最多为zmax;
步骤2-3、设有效高程区域内每两个采样点间有zmax个高程点,且zi个有效散射单元随机分布在zmax个高程点上,则每个高程点对应的有效散射单元个数wj为0或1,其中,j=1,2,...,t,j为高程点索引,t为高程点个数,即:t=zmax×n。
所述步骤3的具体步骤如下:
步骤3-1、将全波LiDAR波形分解为m+1个阶梯,包括m个跳变点和m+1个阶梯高度;
步骤3-2、定义全波LiDAR波形或然率服从均值是阶梯函数的泊松分布,建立全波LiDAR波形或然率模型,即参数集Φ={s,h,m}条件下高程点对应有效散射单元个数的条件概率,其中,s为跳变点,h为阶梯高度,m为跳变点个数;
步骤3-3、根据贝叶斯定理构建已知高程点对应的有效散射单元个数的条件下参数集Φ={s,h,m}的联合后验概率密度函数。
有益效果:
本发明针对采用有效散射单元来描述全波LiDAR(Light Detection AndRanging)波形,将全波LiDAR波形转化为各采样点对应的具有相同散射能量的有效散射单元,能够将地物散射特性与波形分解相结合,并提出采用阶梯函数对有效散射单元个数建模,打破了使用传统高斯统计分布模型来拟合原始波形的思想,为全波LiDAR波形建模提供了有效的新思路。
附图说明
图1为本发明具体实施方式的基于地物散射的全波LiDAR波形分解方法流程图;
图2为本发明具体实施方式的步骤2的具体流程图;
图3为本发明具体实施方式的步骤4-1的具体流程图;
图4为本发明具体实施方式的步骤4-2的具体流程图;
图5为本发明具体实施方式的增加阶梯操作实例;
图6为本发明具体实施方式的仿真实验所采用的数据集以及实验区域对应的遥感影像,其中,图6(a1)至(d1)为原始波形数据,图6(a2)至(d2)为实验区域分别对应的遥感影像,其中不同灰度代表不同高程的地物目标;
图7为本发明具体实施方式的由原始数据得到的每个采样点对应有效散射单元个数均值以及由本发明得到的阶梯函数的叠加图,其中,图7(a)至(d)分别为图6(a1)至(d1)的分解结果;
图8为本发明具体实施方式的由数据计算得到的每个阶梯对应散射单元个数与本发明仿真实验结果得到的每个阶梯对应散射单元个数的叠加图,(a)~(d)分别为不同实验区域的叠加图;
图9为本发明具体实施方式的10000次迭代过程中后验概率取对数后的变化情况,(a)~(d)分别为不同实验区域后验概率取对数后的变化情况。
具体实施方式
下面结合附图详细说明本发明的具体实施方式。
基于地物散射的全波LiDAR波形分解方法,如图1所示,包括:
步骤1、激光雷达发射信号到地物目标上反馈回来的信号叠加形成全波LiDAR波形,包括分布在有效高程区域内的采样点及各采样点处的散射能量。实验数据采用ICESat-GLAS测高数据,每个波形有544个采样点,前151个采样点间隔增加到4ns(60cm),后393个采样点间隔为1ns(15cm),有效高程区间为[0,L],L=150m;
x={xi;i=1,2,...,n},y={yi;i=1,2,...,n},其中,i为采样点索引,xi为分布在有效高程区间[0,L]内的采样点,yi为在采样点xi处的散射能量,n=544为采样点个数。
步骤2、采用有效散射单元来描述全波LiDAR波形:将全波LiDAR波形转化为各采样点对应的具有相同散射能量的有效散射单元;
根据雷达方程原理将全波LiDAR波形用有效散射单元描述,如图2所示,具体如下:
步骤2-1、设各采样点处的散射能量中的散射能量最小值c=min(yi)为一个有效散射单元,c为一个有效散射单元的散射能量;
步骤2-2、计算每个采样点xi对应的有效散射单元个数zi=yi/c,i=1,2,...,n;即采样点处的散射能量yi除以有效散射单元的散射能量c;每个采样点对应的有效散射单元个数最多为zmax;
步骤2-3、设有效高程区域内每两个采样点间有zmax个高程点,且zi个有效散射单元随机分布在zmax个高程点上,则每个高程点对应的有效散射单元个数wj为0或1,其中,j=1,2,...,t,j为高程点索引,t为高程点个数,即:t=zmax×n。
步骤3、采用阶梯函数对全波LiDAR波形中的有效散射单元个数建模:用阶梯的高度表示各个高程点对应的有效散射单元个数的均值,阶梯的位置表示地物目标的相对高程;
步骤3-1、将全波LiDAR波形分解为m+1个阶梯,包括m个跳变点和m+1个阶梯高度;
步骤3-1-1、确定跳变点个数m。
假设跳变点个数m服从均值为λ的泊松分布,其概率密度函数为:其中,m≤mmax,mmax为预先指定的最大阶梯数;
步骤3-1-2、生成m个跳变点s={sk;k=1,2,...,m}。
假设跳变点sk(k=1,2,...,m)服从均匀分布在有效高程区间[0,L]上的2m+1个随机点的偶数有序统计分布,各跳变点的联合概率密度函数为:p(s)=p(s1,s2,...,sm)=(2m+1)!L-2m-1s1(s2-s1)...(sm-sm-1)(L-sm),即:在有效高程区间[0,L]上随机生成2m+1个随机点,把其中偶数点的位置作为跳变点s,跳变点设为0(s0)<s1<s2<...<sm<L(sm+1),其中,k为跳变点索引;
步骤3-1-3、生成m+1个阶梯高度h={hl;l=0,1,...,m}。
跳变点子区间[s0,s1)对应的阶梯高度为h0,跳变点子区间[sk,sk+1),k=1,2,...,m对应的阶梯高度为hk,假设阶梯高度hl(l=0,1,...,m)服从形状参数和尺度参数分别为α和β的Gamma分布且相互独立,各阶梯高度的联合概率密度函数为:
步骤3-2、建立全波LiDAR波形或然率模型,即参数集Φ={s,h,m}条件下高程点对应有效散射单元个数的条件概率;
定义或然率p(wj|Φ)(j=1,2,...,t)为服从均值为g(·)的泊松分布,g(·)为定义在区间[0,L]上的阶梯函数,在[0,L]区间内有m+1个阶梯,包括m个跳变点s以及m+1个阶梯高度h,其中,参数集Φ={s,h,m},s={s1,s2,...,sm},h={h0,h1,...,hm}。或然率可表示为:
步骤3-3、根据贝叶斯定理构建已知高程点对应的有效散射单元个数w的条件下参数集Φ={s,h,m}的联合后验概率密度函数。
根据Bayes定理,已知高程点对应的散射单元个数w条件下参数集Φ={s,h,m}的后验概率密度函数为,
步骤4、利用RJMCMC(可逆跳转马尔科夫链蒙特卡罗)算法对全波LiDAR波形进行分解,即求解阶梯函数,得到地物目标的相对高程以及散射特性。
步骤4-1,更新参数集Φ={s,h,m},包括更新跳变点s和更新阶梯高度h,如图3所示,具体步骤如下:
步骤4-1-1,更新跳变点s。
从当前m个跳变点中随机抽取一个跳变点sk作为候选跳变点sk *,k∈{1,2,...,m},该操作只改变该跳变点的位置,而保持其它参数不变,候选跳变点sk *服从区间[sk-1,sk+1]上的均匀分布,hl和hl+1对应的候选跳变点区间分别为[sk-1,sk *]和[sk *,sk+1];
步骤4-1-2,计算更新跳变点s的接受率其中,在[0,1]之间生成一个随机数e,如果a(sk,sk *)>e,则接受更新跳变点结果,sk=sk *;反之,拒绝更新跳变点,且保持原有参数不变;
步骤4-1-3,更新阶梯高度h。
从当前m+1个阶梯高度中随机抽取一个阶梯高度hl,k∈{0,1,...,m},该操作只改变该阶梯高度,而保持其它参数不变,候选阶梯高度为hl *,log(hl */hl)服从区间[-1/2,1/2]上的均匀分布;
步骤4-1-4,计算阶梯高度h的接受率其中,在[0,1]之间生成一个随机数e,如果a(hl,hl *)>e,接受更新阶梯高度结果,hl=hl *;反之,拒绝更新阶梯高度,且保持原有参数不变。
步骤4-2,增加或删除阶梯,如图4所示,具体步骤如下:
步骤4-2-1,设置阈值参数,包括:设置跳变点个数最大值Tm,最大迭代次数T;
步骤4-2-2,判断执行增加阶梯操作或执行删除阶梯操作:设选择增加阶梯操作的概率和选择删除阶梯的概率分别为bm和dm,当m=1时,dm=0;当m=Tm时,bm=0;当m≠1且m≠Tm时,bm=dm;在[0,1]之间生成一个随机数u,如果u>0.5,则执行增加阶梯操作,转至步骤4-2-3;否则执行删除阶梯操作,转至步骤4-2-4;
步骤4-2-3,执行增加阶梯操作:
首先选取候选跳变点s*,s*均匀分布在区间[0,L]上,并以概率1位于区间(sk,sk+1)上。如果该候选跳变点被接受,则将其标记为sk+1 *,sk+1,sk+2,...,sm重新标号为sk+2 *,sk+3 *,...,sm+1 *,与之对应阶梯高度的标号做同样的操作。对于子区间(sk,sk+1 *)和(sk+1 *,sk+2 *)生成候选阶梯高度hl *和hl+1 *。
设跳变前后变化的区间及阶梯高度满足加权几何平均关系:并假定hl *和hl+1 *间存在微小的扰动,使得hl+1 */hl *=(1-v)/v,其中,v为均匀分布于[0,1]区间的随机数。生成阶梯的接受概率为:ab=min{1,Rb},其中,
在[0,1]之间生成一个随机数e,如果ab>e,接受增加阶梯操作结果,则s*={s1,...,sk,sk+1 *,sk+2 *,sk+3 *,...,sm+1 *},h*={h0,h1,...,hl *,hl+1 *,...,hm *}。
步骤4-2-4,执行删除阶梯操作:
在m个跳变点中随机抽取一个跳变点sk+1,如果删除该跳变点被接受,顺序改变后跳变点的标号为sk+1 *,sk+2 *,...,sm *,使得(sk *,sk+1 *)=(sk,sk+2),其概率为1/m+1。子区间(sk *,sk+1 *)对应的新高度为hl *,并满足如下关系:
删除阶梯操作是生成阶梯操作的对偶操作,因此其接受概率为:在[0,1]之间生成一个随机数e,如果ad>e,接受删除阶梯操作结果,则s*={s1,...,sk,sk+1 *,...,sm *},h*={h0,h1,...,hl-1,hl *,...,hm *}。
图5为增加阶梯操作的实例,其中s *、hl *和hl+1 *分别为新增加的跳变点以及增加跳变点后新生成的两个阶梯高度。
步骤4-3,当迭代指次数达到最大迭代次数T时,停止迭代,反之,重复步骤4-1,继续迭代。
步骤4-4,构成阶梯的两个跳变点之间的中点为该阶梯对应的地物目标的相对高程,阶梯的高度为高程点对应的有效散射单元个数的均值,有效散射单元个数的均值与高程点个数的乘积为有效散射单元个数,有效散射单元个数与根据遥感影像得到的实际地物目标的面积之比为地物的散射系数,即散射特性。
全波LiDAR波形的最后一个返回波为地面的返回波,设其相对高程为0,得到其它地物目标的相对高程,当前阶梯高度即有效散射单元个数的均值与高程点个数的乘积为有效散射单元个数。根据仿真实验得到的有效散射单元个数与根据遥感影像得到的实际地物目标的面积之比为地物的散射系数。若两个地物目标的散射特性相近则认为两个地物目标相同(材质或类型相同)。
通过实验,进一步说明本发明的有效性:
1.实验环境:在CPU为Core i5-34703.20GHz、内存4GB、Windows 7旗舰版系统上使用MATLAB7.12.0软件编程实现仿真。
2.选择4个实验区域进行波形分解的仿真实验,并将分解结果与原始数据进行叠加,以及将实验得到地物目标的相对高程和散射系数与遥感影像及实地测量得到的对应地物相比较,对实验结果进行评价。
3.如图6所示,图6(a1)至(d1)为4个实验区域((a1)至(d1)分别区域I、区域II、区域III、区域IV)的原始波形数据,图6(a2)至(d2)为4个实验区域((a1)至(d1)分别区域I、区域II、区域III、区域IV)对应的遥感影像,其中不同灰度代表不同目标地物。
4.如图7所示,为每个实验区域由原始数据得到每个采样点对应有效散射单元个数与由本发明得到阶梯函数的叠加图,其中,图7(a)至(d)分别为图6(a1)至(d1)的分解结果。可以看出,由本发明得到的阶梯函数跳变点的位置能够和原始数据重合,阶梯函数高度与数据得到的每个采样点对应的散射单元个数趋势相同;根据图6(a2)至(d2)的遥感影像可以看出本发明的仿真实验得到阶梯函数的每个阶梯对应一个地物目标,但是地物有高低起伏变化时会被分成多个阶梯,如图7(a)和(b)的最后三个阶梯对应的是有坡度的地面。
5.如图8(a)~(d)所示,为每个实验区域由数据计算得到的每个阶梯对应散射单元个数与本发明实验结果得到的每个阶梯对应散射单元个数的叠加图。可以看出,本发明得到的每个阶梯对应散射单元的个数与数据每个回波对应有效散射单元个数基本相同。
6.如图9所示,为每个实验区域10000次迭代过程中后验概率取对数后的变化情况,(a)~(d)分别为不同实验区域后验概率取对数后的变化情况。可以看出,本发明方法在迭代过程中后验概率很快收敛,并趋于稳定。
7.为了精确验证本发明方法有效性,根据实验结果计算每类地物目标对应的相对高程以及散射系数。计算结果如表1所示,
表1每类地物目标对应的相对高程以及散射系数
可以看出,本发明方法得到的相对高程与经过实地勘测得到的实际地物高程基本一致;并且相同材质的目标地物散射系数基本相同,区域I和区域II白色地物的沥青屋顶以及黑色地物石灰地面,区域III浅灰地物和区域IV白色地物的水泥屋顶,如图6所示。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何本技术领域人士在本发明所揭露的技术范围内,可理解到的替换,都应涵盖在本发明的包含范围之内,例如基于有效散射单元的波形分解、基于散射系数的波形分解、基于阶梯函数的波形分解等。因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (3)
1.基于地物散射的全波LiDAR波形分解方法,其特征在于,包括:
步骤1、激光雷达发射信号到地物目标上反馈回来的信号叠加形成全波LiDAR波形,包括分布在有效高程区域内的采样点及各采样点处的散射能量;
步骤2、采用有效散射单元来描述全波LiDAR波形:将全波LiDAR波形转化为各采样点对应的具有相同散射能量的有效散射单元;
步骤3、采用阶梯函数对全波LiDAR波形中的有效散射单元个数建模:用阶梯的高度表示各个高程点对应的有效散射单元个数的均值,阶梯的位置表示地物目标的相对高程;
步骤4、利用RJMCMC算法对全波LiDAR波形进行分解,即求解阶梯函数,得到地物目标相对高程以及散射特性。
2.根据权利要求1所述的基于地物散射的全波LiDAR波形分解方法,其特征在于,步骤2的具体步骤如下:
步骤2-1、设各采样点处的散射能量中的散射能量最小值c= min(y i )为一个有效散射单元;
步骤2-2、计算每个采样点对应的有效散射单元个数z i = y i /c, i =1, 2, …, n;即采样点处的散射能量y i 除以有效散射单元的散射能量c;每个采样点对应的有效散射单元个数最多为z max;
步骤2-3、设有效高程区域内每两个采样点间有z max个高程点,且z i 个有效散射单元随机分布在z max个高程点上,则每个高程点对应的有效散射单元个数w j 为0或1,其中,j= 1, 2,…, t,j 为高程点索引,t为高程点个数,即:t = z max×n。
3.根据权利要求1所述的基于地物散射的全波LiDAR波形分解方法,其特征在于,步骤3的具体步骤如下:
步骤3-1、将全波LiDAR波形分解为m+1个阶梯,包括m个跳变点和m+1个阶梯高度;
步骤3-2、定义全波LiDAR波形或然率服从均值是阶梯函数的泊松分布,建立全波LiDAR波形或然率模型,即参数集F= {s, h, m}条件下高程点对应有效散射单元个数的条件概率,其中,s为跳变点,h为阶梯高度,m为跳变点个数;
步骤3-3、根据贝叶斯定理构建已知高程点对应的有效散射单元个数的条件下参数集F= {s, h, m}的联合后验概率密度函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510564627.4A CN105137411B (zh) | 2015-09-07 | 2015-09-07 | 基于地物散射的全波LiDAR波形分解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510564627.4A CN105137411B (zh) | 2015-09-07 | 2015-09-07 | 基于地物散射的全波LiDAR波形分解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105137411A CN105137411A (zh) | 2015-12-09 |
CN105137411B true CN105137411B (zh) | 2017-05-31 |
Family
ID=54722817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510564627.4A Expired - Fee Related CN105137411B (zh) | 2015-09-07 | 2015-09-07 | 基于地物散射的全波LiDAR波形分解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137411B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599808B (zh) * | 2016-12-01 | 2020-12-15 | 中国科学院光电研究院 | 一种基于全波形激光雷达数据的隐蔽目标提取方法 |
CN110673109B (zh) * | 2019-11-01 | 2021-09-21 | 自然资源部国土卫星遥感应用中心 | 一种星载大光斑激光雷达全波形数据分解方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012172526A1 (en) * | 2011-06-17 | 2012-12-20 | Leddartech Inc. | System and method for traffic side detection and characterization |
CN103926575A (zh) * | 2014-03-05 | 2014-07-16 | 南京理工大学 | 一种激光雷达全波形信号分析方法 |
CN104280725A (zh) * | 2013-09-25 | 2015-01-14 | 中国科学院光电研究院 | 全波形激光雷达数据波形分解方法 |
CN104849722A (zh) * | 2015-05-08 | 2015-08-19 | 南京林业大学 | 一种LiDAR波形综合特征的单木识别方法 |
-
2015
- 2015-09-07 CN CN201510564627.4A patent/CN105137411B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012172526A1 (en) * | 2011-06-17 | 2012-12-20 | Leddartech Inc. | System and method for traffic side detection and characterization |
CN104280725A (zh) * | 2013-09-25 | 2015-01-14 | 中国科学院光电研究院 | 全波形激光雷达数据波形分解方法 |
CN103926575A (zh) * | 2014-03-05 | 2014-07-16 | 南京理工大学 | 一种激光雷达全波形信号分析方法 |
CN104849722A (zh) * | 2015-05-08 | 2015-08-19 | 南京林业大学 | 一种LiDAR波形综合特征的单木识别方法 |
Non-Patent Citations (2)
Title |
---|
利用非参数回归跳变检测模型提取LiDAR建筑边缘;夏春林 等;《测绘通报》;20140831(第8期);17-20 * |
模拟回火马尔可夫链蒙特卡罗全波形分析方法;尹文也 等;《物理学报》;20141231;第63卷(第16期);1-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN105137411A (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101196562B (zh) | 一种基于改进的em算法的激光雷达波形数据分解的方法 | |
CN105354860B (zh) | 基于箱粒子滤波的扩展目标CBMeMBer跟踪方法 | |
Jin et al. | A point-based fully convolutional neural network for airborne lidar ground point filtering in forested environments | |
CN107016464B (zh) | 基于动态贝叶斯网络的威胁估计方法 | |
CN107392875A (zh) | 一种基于k近邻域划分的点云数据去噪方法 | |
CN104021289B (zh) | 一种非高斯非稳态噪声建模方法 | |
CN105574929A (zh) | 一种基于地面LiDAR点云数据的单株植被三维建模方法 | |
CN102608620B (zh) | 一种依据反射强度与地形形态的激光扫描点云植被过滤方法 | |
Lin et al. | Machine learning templates for QCD factorization in the search for physics beyond the standard model | |
CN103971115A (zh) | 一种基于NDVI和PanTex指数的高分辨率遥感影像新增建设用地图斑自动提取方法 | |
CN103759732B (zh) | 一种角度信息辅助的集中式多传感器多假设跟踪方法 | |
CN110514567B (zh) | 基于信息熵的气体源搜索方法 | |
CN105137411B (zh) | 基于地物散射的全波LiDAR波形分解方法 | |
CN112017288A (zh) | 一种基于地形图数据的dem综合与简化方法 | |
CN105389821A (zh) | 一种基于云模型和图割相结合的医学图像分割方法 | |
CN103268632B (zh) | 一种机载激光雷达扫描生成地形信息的方法 | |
CN113610139A (zh) | 一种多视角强化图像聚类方法 | |
CN106021289A (zh) | 基于节点用户的概率矩阵分解模型的构建方法 | |
Wang et al. | Roof pressure prediction in coal mine based on grey neural network | |
CN102521655A (zh) | 基于非支配邻域免疫算法的动态网络社区检测方法 | |
Casleton et al. | A local structure model for network analysis | |
CN106373129A (zh) | 一种基于双隶属度的fcm遥感图像分割方法 | |
CN115173421B (zh) | 一种基于渐进拉丁超立方采样的概率最优潮流计算方法 | |
Duca et al. | An overview of non-Gaussian state-space models for wind speed data | |
CN105096379A (zh) | 一种基于k-邻域的三角网格曲面重建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 123000 Zhonghua Road, Liaoning, Fuxin, China, No. 47 Applicant after: LIAONING TECHNICAL University Address before: 110043 Zhonghua Road, Liaoning, Fuxin, China, No. 47 Applicant before: Liaoning Technical University |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170531 Termination date: 20210907 |
|
CF01 | Termination of patent right due to non-payment of annual fee |