CN113740873A - 一种基于高斯卷积的海洋激光雷达快速仿真方法 - Google Patents
一种基于高斯卷积的海洋激光雷达快速仿真方法 Download PDFInfo
- Publication number
- CN113740873A CN113740873A CN202111009955.XA CN202111009955A CN113740873A CN 113740873 A CN113740873 A CN 113740873A CN 202111009955 A CN202111009955 A CN 202111009955A CN 113740873 A CN113740873 A CN 113740873A
- Authority
- CN
- China
- Prior art keywords
- scattering
- phase function
- laser radar
- calculating
- water body
- 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
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000004088 simulation Methods 0.000 title claims abstract description 15
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 44
- 238000001514 detection method Methods 0.000 claims abstract description 5
- 238000002834 transmittance Methods 0.000 claims description 13
- 239000013535 sea water Substances 0.000 claims description 11
- 230000003287 optical effect Effects 0.000 claims description 7
- 230000035699 permeability Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
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
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
-
- 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/497—Means for monitoring or calibrating
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
Description
技术领域
本发明属于海洋激光雷达遥感探测技术领域,尤其涉及一种基于高斯卷积的海洋激光雷达快速仿真方法。
背景技术
海洋激光雷达已广泛应用于海洋探测,但目前对水体参数的反演还多使用简化的单次散射激光雷达方程,多忽略了激光在水体传播过程中的多次散射。而因为多次散射的影响,激光雷达回波信号往往要比单次散射的信号强,会对信号的反演带来误差。海洋激光雷达***的设计和回波信号的反演需要多次散射的准确仿真作为依据。
目前对激光雷达多次散射的仿真多采用蒙特卡洛模型,但该模型计算量大、耗时长。本方法根据高斯激光束与水体高斯散射相函数的卷积,提出了一种解析模型,能够快速地计算激光在水体中的多次散射信号。
发明内容
本发明的目的在于获取激光在水体中的多次散射信号,本发明提供一种基于高斯卷积的海洋激光雷达快速仿真方法,该方法通过高斯激光束与水体高斯散射相函数的卷积,获得海洋激光雷达的多次散射回波信号。
本发明的目的通过如下的技术方案来实现:
一种基于高斯卷积的海洋激光雷达快速仿真方法,该方法包括如下步骤:
步骤3:计算单次散射信号强度P1(z);
步骤5:根据多次散射与单次散射比rn计算多次散射信号强度Pn(z);
步骤6:根据多次散射信号强度Pn(z)计算激光回波信号总强度Pt(z)。
作为优选,步骤3所述的单次散射信号强度为P1(z):
其中,η为接收器探测效率;P0为激光器能量;A为探测器接收面积;o为几何重叠因子,TO为接收器光学透过率;Ta为大气透过率;Ts为海表透过率;v为光速;H为激光雷达所在高度;Δt为激光脉宽;n为海水折射率;z为海水深度;pπ(z)/4π为水体180°散射相函数;b(z)为水体散射系数;c(z′)为水体衰减系数。
作为优选,步骤4所述的多次散射与单次散射比rn为:
其中,ρt为激光雷达接收视场角;θ0为激光束发散角。
作为优选,步骤5所述的多次散射信号强度Pn(z)为:Pn(z)=P1(z)rn。
作为优选,步骤6所述的激光回波信号总强度Pt(z)为:Pt(z)=∑n=1Pn(z)。
作为优选,所述的散射相函数为任意水体散射相函数。
本发明的有益效果是:本发明基于高斯卷积,可以快速地计算激光在水体中的多次散射信号,能够适用于任意散射相函数和非均匀水体。
附图说明
图1是本方法的流程图;
图2是实施例1的多次散射信号强度;
图3是实施例2的多次散射信号强度;
图4是实施例3的多次散射信号强度;
图5是实施例4的多次散射信号强度。
具体实施方式
下面根据附图和优选实施例详细描述本发明,本发明的目的和效果将变得更加明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明所述实施例的实施方式相同,区别仅在于采用的激光雷达***参数不同。
本方法采用典型的激光雷达***参数作为实施例1,激光能量为P0=10mJ,接收面积为A=1.76m2,重叠因子O=1,激光雷达高度为H=300m,接收望远镜视场角FOV=10mrad,响应度为η=0.18,接收器光学透射率为To=0.9,典型环境参数海水折射率n=1.33,海表透射率为Ts=0.95,b=0.037m-1,c=0.151m-1,水体相函数采用Henyey-Greenstein相函数ptrue(0)=339.4。
本发明的具体实施方式为:
步骤2所述的散射相函数前向峰γ为:
计算得到的γ=0.5。
步骤3:计算单次散射信号强度P1(z);
步骤3所述的单次散射信号强度为P1(z):
其中,η为接收器探测效率;P0为激光器能量;A为探测器接收面积;O为几何重叠因子,TO为接收器光学透过率;Ta为大气透过率;Ts为海表透过率;v为光速;H为激光雷达所在高度;Δt为激光脉宽;n为海水折射率;z为海水深度;pπ(z)/4π为水体180°散射相函数;b(z)为水体散射系数;c(z′)为水体衰减系数。
步骤4所述的多次散射与单次散射比rn为:
其中,ρt为激光雷达接收视场角;θ0为激光束发散角。
步骤5:根据多次散射与单次散射比rn计算多次散射信号强度Pn(z);
步骤5所述的多次散射信号强度为:
Pn(z)=P1(z)rn。
步骤6:根据多次散射信号强度Pn(z)计算激光回波信号总强度Pt(z);
步骤6所述的激光回波信号总强度为:
Pt(z)=∑n=1Pn(z)。
图2显示实施例1的多次散射和总信号强度。从图2可以看出,激光刚进入水体时,以回波信号强度以单次散射为主。但随着深度的增加,多次散射信号逐渐增强。
实施例2采用的激光雷达***参数为:激光能量为P0=10mJ,接收面积为A=1.76m2,重叠因子O=1,激光雷达高度为H=300m,接收望远镜视场角FOV=10mrad,响应度为η=0.18,接收器光学透射率为To=0.9,海水折射率n=1.33,海表透射率为Ts=0.95,采用二类水体参数b=0.219m-1,c=0.398m-1,水体相函数采用Henyey-Greenstein相函数ptrue(0)=339.4。
实施例2的多次散射信号强度以及总信号强度如图3所示。从图3中可以看出,随着水体浑浊度的增加,海水中的多次散射信号会大幅增加,回波信号主要以多次散射为主。但水体浑浊度的增加会导致水体总的衰减变大,从而导致总回波信号随深度衰减的速度也变得更快,激光在水体中的穿透深度也会减小。
实施例3采用的激光雷达***参数为:激光能量为P0=10mJ,接收面积为A=1.76m2,重叠因子O=1,激光雷达高度为H=300m,接收望远镜视场角FOV=0.1mrad,响应度为η=0.18,接收器光学透射率为To=0.9,海水折射率n=1.33,海表透射率为Ts=0.95,采用二类水体参数b=0.219m-1,c=0.398m-1,水体相函数采用Henyey-Greenstein相函数ptrue(0)=339.4。
实施例3的多次散射信号强度以及总信号强度如图4所示。从图4中可以看出,总信号几乎和单次散射信号重合,说明在视场角非常小的情况下,激光被散射到接收器的视场角以外,从而导致接收到的回波信号以单次散射为主。
实施例4采用的激光雷达***参数为:激光能量为P0=10mJ,接收面积为A=1.76m2,重叠因子O=1,激光雷达高度为H=300m,接收望远镜视场角FOV=10mrad,响应度为η=0.18,接收器光学透射率为To=0.9,海水折射率n=1.33,海表透射率为Ts=0.95,采用二类水体参数b=0.219m-1,c=0.398m-1,水体相函数采用Fournier-Forand相函数ptrue(0)=165.7。
实施例4的多次散射信号强度以及总信号强度如图5所示。从图5中可以看出,在其他参数一致的情况下,不同的水体相函数会导致不同的结果。Fournier-Forand相函数下的水体多次散射有所减少,从而导致总回波信号的减小。在激光雷达仿真过程中需要根据具体的水体环境选择合适的散射相函数。
本领域普通技术人员可以理解,以上所述仅为发明的优选实施例而已,并不用于限制发明,尽管参照前述实施例对发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内,所做的修改、等同替换等均应包含在发明的保护范围之内。
Claims (8)
6.根据权利要求1所述的基于高斯卷积的海洋激光雷达快速仿真方法,其特征在于,步骤5所述的多次散射信号强度Pn(z)为:
Pn(z)=P1(z)rn。
7.根据权利要求1所述的基于高斯卷积的海洋激光雷达快速仿真方法,其特征在于,步骤6所述的激光回波信号总强度Pt(z)为:
Pt(z)=∑n=1Pn(z)。
8.根据权利要求1所述的基于高斯卷积的海洋激光雷达快速仿真方法,其特征在于,所述的散射相函数为任意水体散射相函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111009955.XA CN113740873B (zh) | 2021-08-31 | 2021-08-31 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111009955.XA CN113740873B (zh) | 2021-08-31 | 2021-08-31 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113740873A true CN113740873A (zh) | 2021-12-03 |
CN113740873B CN113740873B (zh) | 2024-05-28 |
Family
ID=78734081
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111009955.XA Active CN113740873B (zh) | 2021-08-31 | 2021-08-31 | 一种基于高斯卷积的海洋激光雷达快速仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113740873B (zh) |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6100976A (en) * | 1998-09-21 | 2000-08-08 | The Board Of Regents For Oklahoma State University | Method and apparatus for fiber optic multiple scattering suppression |
EP1085346A1 (de) * | 1999-09-19 | 2001-03-21 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Messung der Entfernung von Objekten, atmosphärischen Partikeln und dergleichen mittels Lidar- oder Laserradar-Signalen |
US20040130702A1 (en) * | 2001-02-09 | 2004-07-08 | Jupp David L | Lidar system and method |
US20070242269A1 (en) * | 2004-03-06 | 2007-10-18 | Michael Trainer | Methods and apparatus for determining characteristics of particles |
CN105372641A (zh) * | 2015-12-23 | 2016-03-02 | 国家***第一海洋研究所 | 激光雷达水体回波信号仿真方法 |
CN105842698A (zh) * | 2016-05-18 | 2016-08-10 | 西安电子科技大学 | 一种高分辨sar图像的快速仿真方法 |
CN106526611A (zh) * | 2016-10-31 | 2017-03-22 | 鲁东大学 | 一种利用激光雷达***测量海洋固有光学参数的方法 |
CN107831485A (zh) * | 2017-10-19 | 2018-03-23 | 中国科学院海洋研究所 | 船载多视场激光雷达探测多个水体光学特征参数的方法 |
CN107976686A (zh) * | 2017-10-26 | 2018-05-01 | 浙江大学 | 一种多视场角海洋激光雷达及其视场角择优方法 |
US20180188358A1 (en) * | 2017-01-05 | 2018-07-05 | Innovusion Ireland Limited | METHOD AND SYSTEM FOR ENCODING AND DECODING LiDAR |
CN109995427A (zh) * | 2019-03-25 | 2019-07-09 | 西安电子科技大学 | 一种水下上行激光通信的蒙特卡洛仿真方法 |
CN110568449A (zh) * | 2019-10-14 | 2019-12-13 | 自然资源部第二海洋研究所 | 一种风生粗糙海面激光反射、透射矩阵计算方法 |
CN111965608A (zh) * | 2020-07-16 | 2020-11-20 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN111971606A (zh) * | 2018-01-25 | 2020-11-20 | 意大利学院科技基金会 | 具有高空间分辨率的时间分辨成像方法 |
CN112731341A (zh) * | 2020-11-20 | 2021-04-30 | 自然资源部第二海洋研究所 | 基于双级虚像相位阵列的双波长高光谱分辨率激光雷达*** |
CN113011074A (zh) * | 2020-12-30 | 2021-06-22 | 南京先进激光技术研究院 | 水下无线光通信***信道特性建模方法、电子设备及存储介质 |
CN113029342A (zh) * | 2021-04-02 | 2021-06-25 | 西北工业大学 | 基于双向反射理论的水下目标反射偏振光的仿真方法 |
CN113219496A (zh) * | 2021-03-18 | 2021-08-06 | 杭州电子科技大学 | 一种星载海洋大气参数激光雷达探测*** |
-
2021
- 2021-08-31 CN CN202111009955.XA patent/CN113740873B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6100976A (en) * | 1998-09-21 | 2000-08-08 | The Board Of Regents For Oklahoma State University | Method and apparatus for fiber optic multiple scattering suppression |
EP1085346A1 (de) * | 1999-09-19 | 2001-03-21 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Verfahren zur Messung der Entfernung von Objekten, atmosphärischen Partikeln und dergleichen mittels Lidar- oder Laserradar-Signalen |
US20040130702A1 (en) * | 2001-02-09 | 2004-07-08 | Jupp David L | Lidar system and method |
US20070242269A1 (en) * | 2004-03-06 | 2007-10-18 | Michael Trainer | Methods and apparatus for determining characteristics of particles |
CN105372641A (zh) * | 2015-12-23 | 2016-03-02 | 国家***第一海洋研究所 | 激光雷达水体回波信号仿真方法 |
CN105842698A (zh) * | 2016-05-18 | 2016-08-10 | 西安电子科技大学 | 一种高分辨sar图像的快速仿真方法 |
CN106526611A (zh) * | 2016-10-31 | 2017-03-22 | 鲁东大学 | 一种利用激光雷达***测量海洋固有光学参数的方法 |
US20180188358A1 (en) * | 2017-01-05 | 2018-07-05 | Innovusion Ireland Limited | METHOD AND SYSTEM FOR ENCODING AND DECODING LiDAR |
CN107831485A (zh) * | 2017-10-19 | 2018-03-23 | 中国科学院海洋研究所 | 船载多视场激光雷达探测多个水体光学特征参数的方法 |
CN107976686A (zh) * | 2017-10-26 | 2018-05-01 | 浙江大学 | 一种多视场角海洋激光雷达及其视场角择优方法 |
CN111971606A (zh) * | 2018-01-25 | 2020-11-20 | 意大利学院科技基金会 | 具有高空间分辨率的时间分辨成像方法 |
CN109995427A (zh) * | 2019-03-25 | 2019-07-09 | 西安电子科技大学 | 一种水下上行激光通信的蒙特卡洛仿真方法 |
CN110568449A (zh) * | 2019-10-14 | 2019-12-13 | 自然资源部第二海洋研究所 | 一种风生粗糙海面激光反射、透射矩阵计算方法 |
CN111965608A (zh) * | 2020-07-16 | 2020-11-20 | 自然资源部第二海洋研究所 | 一种基于水体叶绿素浓度的星载海洋激光雷达探测能力评估方法 |
CN112731341A (zh) * | 2020-11-20 | 2021-04-30 | 自然资源部第二海洋研究所 | 基于双级虚像相位阵列的双波长高光谱分辨率激光雷达*** |
CN113011074A (zh) * | 2020-12-30 | 2021-06-22 | 南京先进激光技术研究院 | 水下无线光通信***信道特性建模方法、电子设备及存储介质 |
CN113219496A (zh) * | 2021-03-18 | 2021-08-06 | 杭州电子科技大学 | 一种星载海洋大气参数激光雷达探测*** |
CN113029342A (zh) * | 2021-04-02 | 2021-06-25 | 西北工业大学 | 基于双向反射理论的水下目标反射偏振光的仿真方法 |
Non-Patent Citations (7)
Title |
---|
KAIPENG LI: ""A Dual-Wavelength Ocean Lidar for Vertical Profiling of Oceanic Backscatter and Attenuation"", 《REMOTE SENSING》, 31 December 2020 (2020-12-31), pages 1 - 20 * |
刘志鹏: ""海洋激光雷达反演水体光学参数"", 《遥感学报》, vol. 23, no. 5, 31 May 2019 (2019-05-31), pages 944 - 952 * |
崔晓宇: ""采用半解析蒙特卡洛技术模拟星载海洋激光雷达回波信号的软件"", 《红外与激光工程》, vol. 49, no. 2, 29 February 2020 (2020-02-29), pages 1 - 7 * |
徐沛拓: ""海洋激光雷达多次散射回波信号建模与分析"", 《遥感学报》, vol. 24, no. 2, 29 February 2020 (2020-02-29), pages 142 - 149 * |
李晓龙: ""海洋激光雷达水体光学特性参数反演模型对比及分析"", 《光学学报》, vol. 37, no. 10, 31 October 2017 (2017-10-31), pages 1 - 7 * |
王磊: ""蓝绿激光水下成像***的探测灵敏度分析"", 《红外与激光工程》, vol. 41, no. 1, 31 January 2012 (2012-01-31), pages 79 - 85 * |
罗双: ""辐射传输计算中散射相函数截断近似方法对比分析"", 《光学学报》, vol. 37, no. 2, 28 February 2017 (2017-02-28), pages 1 - 7 * |
Also Published As
Publication number | Publication date |
---|---|
CN113740873B (zh) | 2024-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110095784B (zh) | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 | |
CN106154276B (zh) | 基于海底混响和传播损失的深海海底参数反演方法 | |
CN110207695B (zh) | 一种适用于深海auv的无速度辅助单信标定位方法 | |
CN106500671B (zh) | 一种基于lm算法分解激光雷达波形确定海水深度的方法 | |
CN110568449B (zh) | 一种风生粗糙海面激光反射、透射矩阵计算方法 | |
CN107831485A (zh) | 船载多视场激光雷达探测多个水体光学特征参数的方法 | |
CN102012529A (zh) | 基于激光脉冲后向散射的水中目标尾迹探测***及方法 | |
US11237258B1 (en) | Method for optimization of a parametric array shallow profile difference frequency conversion performance | |
CN107976686A (zh) | 一种多视场角海洋激光雷达及其视场角择优方法 | |
CN113671443A (zh) | 基于掠射角声线修正的水声传感器网络深海目标定位方法 | |
CN110135299B (zh) | 一种用于浅水测深的单波段蓝绿激光波形分析方法及*** | |
CN103969694A (zh) | 基于激光后向散射多普勒频移的水中目标尾迹探测方法 | |
CN110376572B (zh) | 一种机载激光雷达水下探测深度模拟装置 | |
CN101464522B (zh) | 基于激光波前变化的水中动态目标尾迹探测***及方法 | |
CN116609758B (zh) | 一种机载激光测深波形旅行时提取方法 | |
Li et al. | Long baseline acoustic localization based on track-before-detect in complex underwater environments | |
CN101533105B (zh) | 基于激光光斑漂移的水中动态目标尾迹探测***及方法 | |
CN113740873A (zh) | 一种基于高斯卷积的海洋激光雷达快速仿真方法 | |
CN110046374A (zh) | 一种基于高斯型声束的声场计算方法 | |
RU2313802C1 (ru) | Способ измерения расстояния до контролируемого объекта | |
CN112434467B (zh) | 一种水下旋转体目标低目标强度外形设计方法 | |
Zong et al. | Research on Laser Dual-Mode Fusion Detection Method of Ship Wake Bubbles | |
CN118094964B (zh) | 基于湍流相位屏水下偏振光传输的矢量蒙特卡洛模拟方法 | |
Jiang et al. | The research on simulation of seabed reverberation model based on multi-beam sounding system | |
CN114791992B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |