CN105929380A - 一种卫星激光高度计全波形激光雷达数据去噪方法 - Google Patents
一种卫星激光高度计全波形激光雷达数据去噪方法 Download PDFInfo
- Publication number
- CN105929380A CN105929380A CN201610236447.8A CN201610236447A CN105929380A CN 105929380 A CN105929380 A CN 105929380A CN 201610236447 A CN201610236447 A CN 201610236447A CN 105929380 A CN105929380 A CN 105929380A
- Authority
- CN
- China
- Prior art keywords
- data
- laser radar
- wave shape
- full wave
- shape laser
- 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.)
- Granted
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/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/483—Details of pulse systems
- G01S7/486—Receivers
- G01S7/487—Extracting wanted echo signals, e.g. pulse detection
- G01S7/4876—Extracting wanted echo signals, e.g. pulse detection by removing unwanted signals
-
- 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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开一种卫星激光高度计全波形激光雷达数据去噪方法,该去噪方法步骤如下:(1)确定全波形激光雷达数据滤波窗口和多项式次数变化范围;(2)用最小二乘法对原始全波形激光数据进行滑动拟合,并求出卷积系数。用卷积系数对原始全波形激光数据进行卷积计算,完成SG滤波,得到去噪后的全波形激光数据;(3)重复(1)‑(2),遍历所有可能的全波形激光雷达数据滤波窗口和多项式次数组合;(4)最优滤波结果的提取。本发明所提出的方法可有效去除全波形激光雷达数据的噪声,计算复杂度低,可有效保持原始全波形激光波形数据的相对极大值、相对极小值和脉冲宽度。
Description
技术领域
本发明涉及对地观测领域,特别是涉及一种以SG滤波器的对卫星激光高度计全波形激光雷达数据去噪的方法。特别是涉及一种以最小二乘法原理对原始全波形激光雷达数据进行滑动拟合,将与原始全波形激光雷达数据差值平方和最小的去噪结果作为最终的全波形激光雷达数据去噪结果。
背景技术
全波形激光雷达数据是一种以时间为横坐标,以回波强度为纵坐标的波形数据。地物表面形态直接影响全波形激光雷达的回波波形。通过应用波形分解技术对全波形激光雷达数据进行分解,可以得到足印光斑内地物表面与传感器之间的距离,再结合传感器位置就可以得到足印光斑内地物表面的高度分布状况。因此,美国国家航空航天管理局(NASA)与2003年发射率专门用于极地冰盖、云和陆地的高程卫星ICESAT1(Ice,Cloud,and land Elevation Satellite),ICESAT1上所搭载的全波形激光高度计GLAS(Geoscience Laser AltimeterSystem)获取了全球不同地物的全波形激光雷达数据,并用于反演冰盖、云高、森林生物量、湖面高度、城市建筑物密度等,获得了良好的应用效果,极大地促进了星载主动光学遥感乃至整个遥感技术的发展。而且未来几年内美国国家航空航天管理局还技术发射第二颗用于用于极地冰盖、云和陆地的高程卫星ICESAT2。于此同时,各种机载激光高度计发展也方兴未艾。
波形分解是获取传感器至地物表面距离的关键步骤之一,波形分解的精度直接影响所获取的距离精度。波形分解的精度处理与波形分解技术本身的精度有关之外,还会受到输入的全波形激光数据中噪声的影响,输入数据的噪声越大,波形分解结果误差越大。而在全波形激光雷达测量过程中,传感器所引入的噪声又不可避免的混入所输出的原始全波形激光雷达数据中。因此,在波形分解之前,去除原始全波形激光雷达数据中的噪声是十分必要的。目前,用于全波形激光雷达数据去噪的常用方法主要有基于时域的滑动平均(MovingAverage,MA)和基于频域的快速傅里叶变换(Fast Fourier Transform,FFT)。滑动平均容易将波形数据中的尖峰平滑,将脉冲展宽。通过对原始全波形激光雷达数据进行快速傅里叶变换后,在频率域滤波会导致边缘振荡和振铃。另外,快速傅里叶变换不适用于处理信噪比过低的全波形激光雷达数据。
一般来说,由于要穿过整个大气层,星载全波形激光雷达所接受到的地面回波能量大大减小,其信噪比要低于光程较短的机载全波形激光雷达数据。因此,用于星载全波形激光雷达数据的去噪方法应该在去除噪声的同时,要保证原始的波形趋势、脉冲宽度不发生变化。另外,还应该适用于低信噪比数据的处理。GLAS所获取的全波形激光雷达数据用高斯滤波器进行去噪,高斯滤波器也是一种线性平滑滤波器,高斯滤波器可以较好地消除高斯噪声。但是,高斯滤波器常常会使平滑波形边缘、不能很好地去除噪声中的尖峰。SG滤波器(Savitzky-Golay Filter)是一种在时域内基于滑动全波形激光雷达数据滤波窗口多项式拟合的滤波方法。SG滤波器计算量小,计算速度快,可以很好地保持波形的相对极大值、极小值和脉冲宽度等特性,比较适用于星载全波形激光雷数据的去噪。
发明内容
针对现有技术的空白和缺点。本发明所要解决的技术问题是提供一种运行效率高、且能保持波形相对极大值、极小值和脉冲宽度的全波形激光雷达数据噪声去除方法。
为了解决上述技术问题,本发明提供的一种基于SG滤波器的卫星激光高度计全波形激光雷达数据去噪方法,数据处理步骤如下:
(1)确定全波形激光雷达数据滤波窗口和多项式次数变化范围。根据全波形激光雷达数据的采样特点(采样时间间隔、采样点数量)确定全波形激光雷达数据滤波窗口大小(2M+1)及其变化范围(M1≤M≤M2,M∈¢),多项式次数(r)及其变化范围(r1≤r≤r2,r≤2M+1,)。
(2)用最小二乘法对原始全波形激光数据进行滑动拟合,并求出卷积系数。用卷积系数对原始全波形激光雷达数据进行卷积计算,完成SG滤波,得到去噪后的全波形激光雷达数据。
(3)重复步骤(1)-步骤(2),遍历所有可能的全波形激光雷达数据滤波窗口和多项式次数组合;
(4)最优滤波结果的提取,计算每一个全波形激光雷达数据滤波窗口和多项式次数组合得到的去噪后的全波形激光数据与原始数据的差值平方和,差值平方和的计算方法为:
其中,为滤波后数据;z(t)为原始数;t为采样波形数据中的数据点序数,差值平方和最小的去噪结果作为输出数据,从而得到去噪后的全波形激光雷达数据。
本发明所提出的方法可有效去除全波形激光雷达数据的噪声,计算复杂度低,可有效保持原始全波形激光波形数据的相对极大值、相对极小值和脉冲宽度。
附图说明
图1是全波形激光雷达去噪技术流程图。
具体实施方式
以下对本发明的实施例做了进一步详细描述,但本实施例并不限于本发明,凡是采用本发明的相似方法及其相似变化,均应列入本发明的保护范围。
(1)确定全波形激光雷达数据滤波窗口和多项式次数变化范围。根据全波形激光雷达数据的采样特点(采样时间间隔、采样点数量)确定全波形激光雷达数据滤波窗口大小(2M+1)及其变化范围(M1≤M≤M2,M∈¢),多项式次数(r)及其变化范围(r1≤r≤r2,r≤2M+1,)。
(2)用最小二乘法对原始全波形激光数据进行滑动拟合,并求出卷积系数。用卷积系数对原始全波形激光雷达数据进行卷积计算,完成SG滤波,得到去噪后的全波形激光雷达数据。在2M+1全波形激光雷达数据滤波窗口内,可以用多项式来拟合波形数据,多项式表达为
其中,p(n)为多项式;n为2M+1全波形激光雷达数据滤波窗口内波形数据点的位置,-M≤n≤M;k为多项式次数;ak为多项式系数。用最小二乘法求出系数ak,就可以计算拟合后的全波形激光雷达数据滤波窗口内某点的拟合数值,该数值作为去噪后的数值。按照上述方式,遍历全波形激光雷达数据的所有点,就得到了该条波形数据的SG滤波结果。
(3)重复(1)-(3),遍历所有可能的全波形激光雷达数据滤波窗口和多项式次数组合;
(4)最优滤波结果的提取。计算每种全波形激光雷达数据滤波窗口和多项式次数组合得到的去噪后的全波形激光数据与原始数据的差值平方和,最小者作为最终的输出数据,即去噪后的全波形激光数据。计算原始波形数据与去噪后波形数据之间的差值平方和的方法为
其中,为滤波后数据;z(t)为原始数;t为采样波形数据中的数据点序数。差值平方和最小的去噪结果作为输出数据,从而得到去噪后的全波形激光雷达数据。
Claims (1)
1.一种卫星激光高度计全波形激光雷达数据去噪方法,其特征在于包括如下步骤:
(1)确定全波形激光雷达数据滤波窗口和多项式次数变化范围;全波形激光雷达数据滤波窗口大小为2M+1,其变化范围:M1≤M≤M2,多项式次数r及其变化范围为:r1≤r≤r2,r≤2M+1,
(2)用最小二乘法对原始全波形激光雷达数据进行滑动拟合,并求出卷积系数,用卷积系数对原始全波形激光数据进行卷积运算,完成SG滤波,得到去噪后的全波形激光雷达数据;
(3)重复步骤(1)-步骤(2),遍历所有可能的全波形激光雷达数据滤波窗口和多项式次数组合;
(4)最优滤波结果的提取,计算每一个全波形激光雷达数据滤波窗口和多项式次数组合得到的去噪后的全波形激光数据与原始数据的差值平方和,差值平方和的计算方法为:
其中,为滤波后数据;z(t)为原始数;t为采样波形数据中的数据点序数;差值平方和最小的去噪结果作为输出数据,从而得到去噪后的全波形激光雷达数据。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510866362 | 2015-12-01 | ||
CN2015108663623 | 2015-12-01 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105929380A true CN105929380A (zh) | 2016-09-07 |
CN105929380B CN105929380B (zh) | 2018-06-26 |
Family
ID=56838216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610236447.8A Active CN105929380B (zh) | 2015-12-01 | 2016-04-15 | 一种卫星激光高度计全波形激光雷达数据去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105929380B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111580068A (zh) * | 2020-05-28 | 2020-08-25 | 云南电网有限责任公司电力科学研究院 | 一种基于卫星激光雷达技术的遥感数据处理方法 |
CN112233022A (zh) * | 2019-07-15 | 2021-01-15 | 天津大学 | 一种基于梯度域修正的激光雷达引导滤波去噪方法 |
CN112230198A (zh) * | 2019-07-15 | 2021-01-15 | 天津大学 | 一种基于梯度窗宽权重修正的激光雷达回波波形去噪方法 |
CN115393547A (zh) * | 2022-08-15 | 2022-11-25 | 中国地质大学(北京) | 一种月球卫星重力异常数据的全方向滤波方法及系统 |
CN117630883A (zh) * | 2023-11-30 | 2024-03-01 | 无锡中科光电技术有限公司 | 一种全波形可控的大气颗粒物监测激光雷达信号降噪方法、装置及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278961B1 (en) * | 1997-07-02 | 2001-08-21 | Nonlinear Solutions, Inc. | Signal and pattern detection or classification by estimation of continuous dynamical models |
CN102721968A (zh) * | 2012-06-21 | 2012-10-10 | 中国人民解放军电子工程学院 | 星对星sal成像中跟瞄光束偏移修正方法 |
-
2016
- 2016-04-15 CN CN201610236447.8A patent/CN105929380B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6278961B1 (en) * | 1997-07-02 | 2001-08-21 | Nonlinear Solutions, Inc. | Signal and pattern detection or classification by estimation of continuous dynamical models |
CN102721968A (zh) * | 2012-06-21 | 2012-10-10 | 中国人民解放军电子工程学院 | 星对星sal成像中跟瞄光束偏移修正方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112233022A (zh) * | 2019-07-15 | 2021-01-15 | 天津大学 | 一种基于梯度域修正的激光雷达引导滤波去噪方法 |
CN112230198A (zh) * | 2019-07-15 | 2021-01-15 | 天津大学 | 一种基于梯度窗宽权重修正的激光雷达回波波形去噪方法 |
CN112230198B (zh) * | 2019-07-15 | 2022-10-25 | 天津大学 | 一种基于梯度窗宽权重修正的激光雷达回波波形去噪方法 |
CN111580068A (zh) * | 2020-05-28 | 2020-08-25 | 云南电网有限责任公司电力科学研究院 | 一种基于卫星激光雷达技术的遥感数据处理方法 |
CN115393547A (zh) * | 2022-08-15 | 2022-11-25 | 中国地质大学(北京) | 一种月球卫星重力异常数据的全方向滤波方法及系统 |
CN115393547B (zh) * | 2022-08-15 | 2023-03-10 | 中国地质大学(北京) | 一种月球卫星重力异常数据的全方向滤波方法及系统 |
CN117630883A (zh) * | 2023-11-30 | 2024-03-01 | 无锡中科光电技术有限公司 | 一种全波形可控的大气颗粒物监测激光雷达信号降噪方法、装置及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105929380B (zh) | 2018-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105929380A (zh) | 一种卫星激光高度计全波形激光雷达数据去噪方法 | |
CN108845306B (zh) | 基于变分模态分解的激光雷达回波信号去噪方法 | |
CN111583214A (zh) | 基于rbf神经网络的航海雷达图像反演海面风速方法 | |
CN104360332B (zh) | 基于地基合成孔径雷达干涉的大气相位屏提取方法 | |
CN109459745B (zh) | 一种利用辐射噪声估计运动声源速度的方法 | |
KR101914550B1 (ko) | 레이더의 표적 위치 추적 방법 | |
CN102426700B (zh) | 基于局部和全局区域信息的水平集sar图像分割方法 | |
CN102831436B (zh) | 基于参考图互信息波段选择和emd的高光谱图像分类方法 | |
CN110134976B (zh) | 一种机载激光测深信号提取方法及系统 | |
CN104793253A (zh) | 基于数学形态学的航空电磁数据去噪方法 | |
CN104730519B (zh) | 一种采用误差迭代补偿的高精度相位解缠方法 | |
CN107992676A (zh) | 一种运动目标时域散射回波的快速仿真建模方法 | |
CN108957403B (zh) | 一种基于广义互相关的高斯拟合包络时延估计方法及系统 | |
CN104062642A (zh) | 一种对激光雷达波形数据进行高斯回波分解的方法 | |
CN107329137A (zh) | 基于box样条小波框架的探地雷达地下目标定位方法及系统 | |
CN103809180B (zh) | 用于InSAR地形测量的方位向预滤波处理方法 | |
CN103983969B (zh) | 基于二次拟合距离方程的地面加速运动目标成像方法 | |
CN104020456A (zh) | 一种基于多特显点线阵成像雷达系统幅相误差校正方法 | |
KR101426226B1 (ko) | 레이더의 신호 처리 방법 | |
CN108008416B (zh) | 一种估算斜路径对流层延迟的积分方法 | |
CN106597556A (zh) | 一种机场跑道异物检测系统背景对消的方法 | |
CN103377465B (zh) | 基于素描图和核选择的sar图像降斑方法 | |
CN109116359B (zh) | 一种机载雷达低空风切变风场回波风速的估计方法 | |
CN101871780A (zh) | 微惯性器件信号的虚拟野值降噪方法 | |
CN108646248B (zh) | 一种针对低速运动声源的被动声学测速测距方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |