CN111751328B - 一种快速测量高反光空间目标材质的方法 - Google Patents

一种快速测量高反光空间目标材质的方法 Download PDF

Info

Publication number
CN111751328B
CN111751328B CN202010651854.1A CN202010651854A CN111751328B CN 111751328 B CN111751328 B CN 111751328B CN 202010651854 A CN202010651854 A CN 202010651854A CN 111751328 B CN111751328 B CN 111751328B
Authority
CN
China
Prior art keywords
sample
distribution function
light
laser
detector
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
Application number
CN202010651854.1A
Other languages
English (en)
Other versions
CN111751328A (zh
Inventor
代京京
赵思思
王智勇
张景豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Beijing Institute of Space Research Mechanical and Electricity
Original Assignee
Beijing University of Technology
Beijing Institute of Space Research Mechanical and Electricity
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology, Beijing Institute of Space Research Mechanical and Electricity filed Critical Beijing University of Technology
Priority to CN202010651854.1A priority Critical patent/CN111751328B/zh
Publication of CN111751328A publication Critical patent/CN111751328A/zh
Application granted granted Critical
Publication of CN111751328B publication Critical patent/CN111751328B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/47—Scattering, i.e. diffuse reflection
    • G01N2021/4704—Angular selective
    • G01N2021/4711—Multiangle measurement
    • G01N2021/4714—Continuous plural angles

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种快速测量高反光空间目标材质的方法，包括：超连续谱激光光源经过实验系统光路准直后入射到待测样品上；超连续谱激光经待测样品散射后，散射光信号被光电探测器接收；光电探测器输出信号至信号处理系统进行处理，得到待测样品的激光雷达散射截面；基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数；将获得的双向反射分布函数在已有的双向反射分布函数特征库中进行匹配，识别出待测样品的表面材料。本发明可以快速计算太阳光照射下被测材质的双向反射分布函数，基于双向反射分布函数特征库进行匹配分析，可识别出待测样品的表面材料。

Description

一种快速测量高反光空间目标材质的方法
技术领域
本发明涉及光学技术领域，具体涉及一种快速测量高反光空间目标材质的方法。
背景技术
激光雷达散射截面(LRCS)是用来描述目标对照射到它上面的激光光束的散射能力的物理量；可以定义为：在入射光波为均匀平面波时，激光光束光斑完全覆盖目标表面且待测目标表面材料为各项同性材料时，探测器接收到总的激光散射光光功率与入射激光光功率之比为：
Figure BDA0002575273260000011
式中，Ei为目标处入射光波的电场强度振幅，Er为激光雷达雷达接收机接收到的来自目标的散射光波的电场强度振幅，R则为目标和接收机之间的距离。由于探测方向不同时，在该方向上接收到的散射光功率不同，则计算得到的激光雷达散射截面在各个方向上也不同；即在不同的散射方向上，相同目标的激光雷达散射截面一般不相同；对于单站式的激光雷达散射截面称为激光雷达的后向散射截面。
目前常用的BRDF测量装置有两种：一种是采用双圆弧轨道实现半球范围内方位角和入射角度的探测，如专利《一种测量物体表面双向反射分布的装置》，专利号：2009200033317.X，这种测量方法的局限性在于圆弧轨道遮光影响其测量结果以及受圆弧轨道加工精度不高，直接会影响到测量角度的定位精度不高，进而导致整个系统测量可靠性降低。
另一种装置采用光源固定，样品使用四维机械手固定，光电探测器可以安装在垂直面内做±180°旋转的方式完成对入射角
Figure BDA0002575273260000012
和反射角
Figure BDA0002575273260000013
的测量。如专利《一种新型双向反射分布函数测量装置》，专利号：201210075733.2，但是这种方式的局限性在于光电探测器接收面小，要求反射光束扩散角很小，这就意味着只能进行镜面物体表面的散射，如果是粗糙表面则很难实现散射到整个半球空间的各个光信号的采集，导致信号误差。
传统的太阳光照射下样品的双向反射分布函数的实验装置复杂，精密度要求高，通常需要采集上万个实验数据才可以进行模拟计算得到双向反射分布函数。
发明内容
针对现有技术中存在的上述问题，本发明提供一种快速测量高反光空间目标材质的方法。
本发明公开了一种快速测量高反光空间目标材质的方法，包括：
激光光源经过实验系统光路准直后入射到待测样品上；
激光经待测样品散射后，散射光信号被光电探测器接收；
光电探测器输出信号至信号处理系统进行处理，得到待测样品的激光雷达散射截面；
基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数；
将获得的双向反射分布函数在已有的双向反射分布函数特征库中进行匹配，识别出待测样品的表面材料。
作为本发明的进一步改进，激光光源经准直后的光束，采用双光路进行测量。
作为本发明的进一步改进，使用的激光为超连续谱激光，激光波长范围为：400nm-2200nm，该波长的激光可以模拟太阳光和近红外光辐照下的空间背景光。
作为本发明的进一步改进，根据先对超连续谱波长得到的雷达散射截面进行波长信号取均值，再基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数。
作为本发明的进一步改进，待测样品安装在五维位移台上且处于转台中心轴位置，入射激光方向保持不变，探测器所在的转臂保持静止，样品台上的待测样品随着转台转动。
与现有技术相比，本发明的有益效果为：
本发明采用几百个的实验数据就可以分析和模拟太阳光照射下材料的双向放射分布函数；实验装置简单，计算出样品的激光雷达散射截面和双向反射分布函数，基于双向反射分布函数特征库进行匹配分析，可识别出待测样品的表面材料。
附图说明
图1为本发明一种实施例公开的实现快速测量高反光空间目标材质的方法的流程图；
图2为本发明一种实施例公开的实现快速测量高反光空间目标材质的方法的装置结构图；
图3为本发明一种实施例公开的后向散射光线方向示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合附图对本发明做进一步的详细描述：
在实际测量中，入射光波的电场强度振幅不便于直接测量，应该变换成能够直接方便测量的物理量。由于光能流密度正比于光波的电场强度振幅的模的平方，即I∝|E|2，所以激光雷达散射截面可以表示为：
Figure BDA0002575273260000031
式中，R为目标到探测器的距离，Is为散射回波在探测器方向处的光能流密度，Ii为入射到目标表面上的光能流密度。由于公式中各个参量测量简单，因此常用在实验中计算激光雷达散射截面。
在实际应用中，激光雷达散射截面σ还可以定义为：
σ＝ρA0G
式中，ρ为目标表面半球反射率，A0为目标实际的投影面积，G为目标增益。
本发明采用双光路方法，设探测器D1的输出电压为ν1，探测器D2(固定在测量转臂末端)的输出电压为ν2，他们都正比于激光输出功率，而且比例系数均为常数。所以ν2/ν1＝k0是一个与输出功率无关的常量，但v2是与目标的散射截面有关，因为由雷达方程可知在一定条件下，散射截面与接收功率成正比，在此即是与比值成正比。目标散射截面探测定标采用标准体，其中面目标材料定标采用标准白板，三维体目标的散射截面测量定标采用朗伯球。
假设由探测器D2接收到白板散射回来的信号电压为ν1，同时由D1探测器接收到白板的参考电压为ν2，样片的信号电压为ν3，样片的参考电压为ν4，标准白板的电压比值为k0，样片的电压比值为ks，则：
Figure BDA0002575273260000041
设白板单位面积激光雷达散射截面为
Figure BDA0002575273260000042
白板双向反射分布函数BRDF为
Figure BDA0002575273260000043
样片单位面积激光雷达散射截面为
Figure BDA0002575273260000044
样片的双向反射分布函数BRDF为
Figure BDA0002575273260000045
在相同条件下，白板和样片的照射面积相同，所以：
Figure BDA0002575273260000046
双向反射分布函数和单位面积激光散射雷达截面之间的函数关系为：
Figure BDA0002575273260000047
Figure BDA0002575273260000048
联合式(1)、式(2)和式(3)得到：
Figure BDA0002575273260000049
标准白板是朗伯面，所以白板的双向反射分布函数为
Figure BDA00025752732600000410
将其代入得：
Figure BDA00025752732600000411
由式(4)即可通过实验获取白板的单位面积散射截面，将实验结果与理论计算值进行对比，分析测量系统的可靠性以及改进方法等。
在实际测量中，BRDF双向反射分布函数的测定要求照射面元的面积保持一致，因此：
Figure BDA00025752732600000412
其中θ为入射角。
基于上述原理，如图1、2所示，本发明提供一种快速测量高反光空间目标材质的方法，实现该方法的装置包括：测试光源、转台与样品架、信号检测与处理系统；具体包括：激光光源1、斩波器电源2、斩波器3、分束镜A 4、透镜5、分束镜B 6、探测器A 7、探测器B 8、锁相放大器9、微型计算机10、转台11、样品架12和步进电机13；其中，信号检测与处理系统包括锁相放大器和微型计算机；所有光电探测器均与锁相放大器相连，锁相放大器与微型计算机相连。该方法包括：
S1、超连续谱激光光源经过实验系统光路准直后入射到待测样品上；其中，
测试光源使用的激光为超连续谱激光，激光波长范围为：400nm-2200nm，该波长的激光可以模拟太阳光和近红外光辐照下的空间背景光。激光光源经准直后的光束，采用双光路的方法，减小激光强度不稳定引起的误差，提高测量精度；光束入射至转台与样品架的待测样品上。
S2、超连续谱激光经待测样品散射后，散射光信号被光电探测器接收；其中，
待测样品安装在五维位移台上且处于转台中心轴位置，入射激光方向保持不变，探测器所在的转臂保持静止，样品台上的待测样品随着转台转动；样品架为待测样品表面单站后向测量样品架，待测样品处在不同方位时，入射角随之改变，待测样品随着转台的轴线转动；
S3、光电探测器输出信号至信号处理系统进行处理，得到待测样品的激光雷达散射截面；
S4、基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数；具体为：
根据先对超连续谱波长得到的雷达散射截面进行波长信号取均值，再基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数；
S5、将获得的双向反射分布函数在已有的双向反射分布函数特征库中进行匹配，识别出待测样品的表面材料。
具体的：
本发明在空间目标的LRCS测量前，需进行实验系统噪声测量和稳定性测量；测量时，样品平面表面单站后向BRDF及单位面积LRCS测量；首先对平面样品表面进行LRCS测量，测量时该表面应处在转轴上。如图3所示，后向测量时入射激光方向
Figure BDA0002575273260000061
保持不变，探测器所在的转臂保持静止，样品台上的样品随转台B转动，样品处在不同方位时，入射角θi随之改变。对于后向散射，接收方向
Figure BDA0002575273260000062
散射角θs＝θi。
Figure BDA0002575273260000063
为激光入射方向，
Figure BDA0002575273260000064
为接收方向(下同)，样品绕着转台的轴线转动。
空间目标LRCS仿真计算与分析：
由于实际材料表面的复杂性，因此利用实验测量获取典型材料的双向反射分布函数数据，采用经验和半经验的方法建立BRDF统计模型，进而分析和计算目标LRCS是十分必要的。
对于几何简单体可以根据目标表面的BRDF和LRCS之间的关系获得解析形式的σ0。根据BRDF和LRCS的定义式，可得二者关系为：
Figure BDA0002575273260000065
式中：
Figure BDA0002575273260000066
正好是面元材料的BRDF。将上式沿目标表面积分，即得目标的LRCS：
σ＝∫∫Adσ＝∫∫A4πf(θi,θs)cosθicosθsdA (4-2)
Figure BDA0002575273260000067
通常可上述利用公示计算得到在各种入射-散射条件下的LRCS。
换算成BRDF：
Figure BDA0002575273260000068
本发明的优点为：
本发明采用几百个的实验数据就可以分析和模拟太阳光照射下材料的双向放射分布函数；实验装置简单，计算出样品的激光雷达散射截面和双向反射分布函数，基于双向反射分布函数特征库进行匹配分析，可识别出待测样品的表面材料。
以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种快速测量高反光空间目标材质的方法，其特征在于，基于激光雷达散射截面（LRCS）快速测量高反光空间目标材质，包括：
超连续谱激光光源经过实验系统光路准直后入射到待测样品上；其中，激光光源经准直后的光束，采用双光路进行测量；具体包括：沿光线方向依次设置的超连续谱激光光源、斩波器、分束镜A、透镜、分束镜B和样品架，所述斩波器与斩波器电源相连、所述斩波器电源与锁相放大器相连，所述分束镜A的入射光的反射侧设有探测器A、所述探测器A与锁相放大器相连，从分束镜A透射出的光经透镜到达分束镜B，透过分束镜B的光照射至待测样品上；所述分束镜B接收待测样品的散射光的反射侧设有探测器B、所述探测器B与锁相放大器相连，所述样品架安装在转台上，所述转台与步进电机相连，所述锁相放大器和步进电机均连接在微型计算机上；所述双光路为探测器A和探测器B的两个探测光路；其中，探测器A和探测器B为光电探测器；
光电探测器输出信号至信号处理系统进行处理，得到待测样品的激光雷达散射截面；
基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数；其中，根据先对超连续谱激光波长得到的雷达散射截面进行波长信号取均值，再基于激光雷达散射截面和双向反射分布函数之间的函数关系，获得待测样品表面的双向反射分布函数；
将获得的双向反射分布函数在已有的双向反射分布函数特征库中进行匹配，识别出待测样品的表面材料。
2.如权利要求1所述的方法，其特征在于，超连续谱激光的波长范围为：400nm-2200nm，该波长的激光可以模拟太阳光和近红外光辐照下的空间背景光。
3.如权利要求1所述的方法，其特征在于，待测样品安装在五维位移台上且处于转台中心轴位置，入射激光方向保持不变，探测器所在的转臂保持静止，样品台上的待测样品随着转台转动。
CN202010651854.1A 2020-07-08 2020-07-08 一种快速测量高反光空间目标材质的方法 Active CN111751328B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010651854.1A CN111751328B (zh) 2020-07-08 2020-07-08 一种快速测量高反光空间目标材质的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010651854.1A CN111751328B (zh) 2020-07-08 2020-07-08 一种快速测量高反光空间目标材质的方法

Publications (2)

Publication Number Publication Date
CN111751328A CN111751328A (zh) 2020-10-09
CN111751328B true CN111751328B (zh) 2021-04-09

Family

ID=72709955

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010651854.1A Active CN111751328B (zh) 2020-07-08 2020-07-08 一种快速测量高反光空间目标材质的方法

Country Status (1)

Country Link
CN (1) CN111751328B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112835013A (zh) * 2020-12-29 2021-05-25 北京环境特性研究所 一种基于连续激光的目标lrcs模拟测量系统及方法
CN113607694B (zh) * 2021-08-02 2023-09-19 上海无线电设备研究所 一种低散射涂层双向反射分布函数高精度测量装置
CN115201793A (zh) * 2022-06-27 2022-10-18 中国计量科学研究院 定标模组及其使用方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102095688B (zh) * 2011-03-08 2012-07-18 中国人民解放军总装备部工程兵科研一所 一种用于材料激光性能的测量设备
CN103699810A (zh) * 2014-01-14 2014-04-02 西安电子科技大学 一种粗糙面微波段双向反射分布函数的建模方法
CN106770045A (zh) * 2016-11-28 2017-05-31 北京理工大学 一种局部照射下目标激光散射特性的简便计算方法
CN107193897A (zh) * 2017-05-09 2017-09-22 北京理工大学 一种测试及建立航天航空包覆材料brdf数据库的方法和系统
CN109470656A (zh) * 2018-12-26 2019-03-15 北京工业大学 一种新型双向反射分布函数快速测试系统及方法
CN109490253A (zh) * 2018-12-26 2019-03-19 北京工业大学 一种新型模拟自然光双向反射分布函数测试装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102095688B (zh) * 2011-03-08 2012-07-18 中国人民解放军总装备部工程兵科研一所 一种用于材料激光性能的测量设备
CN103699810A (zh) * 2014-01-14 2014-04-02 西安电子科技大学 一种粗糙面微波段双向反射分布函数的建模方法
CN106770045A (zh) * 2016-11-28 2017-05-31 北京理工大学 一种局部照射下目标激光散射特性的简便计算方法
CN107193897A (zh) * 2017-05-09 2017-09-22 北京理工大学 一种测试及建立航天航空包覆材料brdf数据库的方法和系统
CN109470656A (zh) * 2018-12-26 2019-03-15 北京工业大学 一种新型双向反射分布函数快速测试系统及方法
CN109490253A (zh) * 2018-12-26 2019-03-19 北京工业大学 一种新型模拟自然光双向反射分布函数测试装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
复杂目标激光雷达散射截面一体化计算研究;张昊等;《激光与光电子学进展》;20141231;第021404-1至021404-6页 *

Also Published As

Publication number Publication date
CN111751328A (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN111751328B (zh) 一种快速测量高反光空间目标材质的方法
Gayet et al. A new airborne polar Nephelometer for the measurements of optical and microphysical cloud properties. Part I: Theoretical design
US5308981A (en) Method and device for infrared analysis, especially with regard to food
US7813470B2 (en) Three-dimensional contents determination method using transmitted x-ray
US9804087B2 (en) Hemispherical scanning optical scatterometer
WO2019091038A1 (zh) 透反射式集成装置及光谱仪系统
US20030133121A1 (en) Directional reflectometer
CN109490253B (zh) 一种新型模拟自然光双向反射分布函数测试装置
CN109470656A (zh) 一种新型双向反射分布函数快速测试系统及方法
CN109444082B (zh) 漫反射光谱测量装置及测量方法
CN209086170U (zh) 一种高反射镜表面疵病参数表征装置
CN208125613U (zh) 一种反射率测定装置
RU2395073C2 (ru) Способ и оборудование для количественного анализа растворов и дисперсий с помощью ближней инфракрасной спектроскопии
CN112414965A (zh) 太赫兹三维层析成像系统样品位置自动调节装置及方法
CN109297986A (zh) 激光陀螺高反射镜表面疵病参数表征装置和检测方法
CN110196229B (zh) 光纤探头装置
CN216350331U (zh) 一种用于太赫兹无损检测的可视化垂直反射测量装置
CN207908363U (zh) 一种双向散射分布测量装置
CN110646379A (zh) 一种用于对航空煤油浊度进行连续角度散射信号测量的装置
CN209513614U (zh) 一种新型双向反射分布函数快速测试系统
CN209513615U (zh) 一种新型模拟自然光双向反射分布函数测试装置
EP3830553B1 (en) Diffuse reflectance apparatus
CN209745812U (zh) 漫反射光谱测量装置
CN106706557A (zh) 一种基于球面漫反射的便携式近红外检测装置及方法
CN107607194B (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