一种激光成像光谱雷达遥感探测方法及设备

文档序号:26588759发布日期:2021-09-10 20:11阅读:81来源:国知局
一种激光成像光谱雷达遥感探测方法及设备

1.本技术涉及激光技术应用领域,尤其涉及一种激光成像光谱雷达遥感探测方法及设备。


背景技术:

2.大气状态对于飞行器大气数据系统、有毒气体泄漏遥感、气体湍流与流场特性测量与诊断等领域至关重要,为测得准确的大气数据,fabry

perot干涉仪、michelson干涉仪、doppler差分干涉仪等相继被用于大气参数的探测中,而激光外差光谱技术因其分辨率高、灵敏度高、体积小等优势也受到了广泛的关注。但是现有的探测遥感方法受探测范围影响,无法进行高精度探测,且测量周期长,数据处理效率低。


技术实现要素:

3.本技术实施例提供一种激光成像光谱雷达遥感探测方法及设备,用以解决现有的遥感探测方法探测精度低,且测量周期长、数据处理效率低的技术问题。
4.一方面,本技术实施例提供了一种激光成像光谱雷达遥感探测方法,包括:基于图像传感器,按帧获取不同波长的调谐激光后向散射光对应的图像数据;针对所述图像数据中的各像元位置,根据该像元在各帧图像数据中的时间序列描述相应激光传输路径的散射光强变化,进行基线重构,得到对应的吸收谱线;对所述吸收谱线进行线形函数拟合,得到对应的吸收谱线参数,以计算所述后向散射光的散射区域内的气体浓度、大气温度、气体压强和气流速度;根据所述图像数据中各像元位置对应的像元数据序列,计算相应空间位置处的空气密度。
5.另一方面,本技术实施例还提供了一种激光成像光谱雷达遥感探测设备,设备包括:处理器;及存储器,其上存储有可执行代码,当可执行代码被执行时,使得处理器执行如上述的一种激光成像光谱雷达遥感探测方法。
6.本技术实施例提供的一种激光成像光谱雷达遥感探测方法及设备,至少具备以下有益效果:
7.将成像技术与光谱技术结合在一起,在探测物体空间特征的同时对每个空间像元吸收光谱进行探测,以较高帧率获取目标的光路图像。由于图像传感器高光谱分辨率的巨大优势,在获取空间散射信号的同时获取气体分子吸收光谱信息,可通过图像传感器获取的信息反演空气的密度、氧气浓度、空气温度、压力和速度等多个参量,实现飞行器所需要的大气数据的多种参数的同步传感和测量,提高了高光谱数据处理的自动化和智能化水平,在航空、航天器中有广泛的应用。
附图说明
8.此处所说明的附图用来提供对本技术的进一步理解,构成本技术的一部分,本技术的示意性实施例及其说明用于解释本技术,并不构成对本技术的不当限定。在附图中:
9.图1为本技术实施例提供的一种激光成像光谱雷达遥感探测方法流程图;
10.图2为本技术实施例提供的一种激光成像光谱雷达遥感探测设备结构示意图。
具体实施方式
11.为使本技术的目的、技术方案和优点更加清楚,下面将结合本技术具体实施例及相应的附图对本技术技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本技术一部分实施例,而不是全部的实施例。基于本技术中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本技术保护的范围。
12.图1为本技术实施例提供的一种激光成像光谱雷达遥感探测方法流程图,具体包括以下步骤:
13.s101、基于图像传感器,按帧获取不同波长的调谐激光后向散射光对应的图像数据。
14.不同波长的调谐激光对空气进行照射,空气中的气溶胶和气体分子对调谐激光发生散射和吸收作用,可得到对应的后向散射光。图像传感器通过光电转换可将采集到的后向散射光转换为对应的图像数据,并逐帧获取图像数据。某一帧图像数据所记录的是某时刻调谐激光在不同空间位置和角度的散射光,该时刻对应了某个特定激光波长,各帧图像数据对应的激光波长不同。其中,图像传感器可以是增强电荷耦合器件(intensified charge coupled device,iccd)相机、ccd相机、cmos相机等任意能够覆盖所使用波长的设备。
15.在一个实施例中,氧气、二氧化碳、甲烷等大气中的主要气体成分,在海拔100km高空以内,其体积比浓度变化不大,且具有特征吸收作用,适合用于进行大气数据的测量。因此,本技术可设定用于进行测量的气体分子,并确定与其对应的激光信号的中心波长。比如,采用760nm的激光信号,或者其他与氧气分子能够发生作用的谱段,使氧气分子发生特征吸收。
16.在一个实施例中,在发射激光时,根据预设周期循环发射不同波长的激光形成周期性的发射。激光照射到运动物体上产生反射或散射,当激光发射到空气中后,气溶胶和气体分子会对调谐激光进行散射或吸收,通过滤光片对经过气体分子散射及吸收后得到的后向散射光中的杂散光进行滤除,滤除杂散光之后的后向散射光进入图像传感器后会被转换为相应的图像数据。通过优化的调谐范围和调谐周期设计、充分利用图像传感器的数据采集帧率,可避免数据采样率对光谱分辨率的限制。
17.s102、针对所述图像数据中的各像元位置,根据该像元位置在各帧图像数据中的时间序列描述相应激光传输路径的散射光强变化,进行基线重构,得到对应的吸收谱线。
18.图像数据中某像元位置的时间序列,代表了激光沿某一路径散射光强的变化。对所有帧的同一像元位置的时间序列进行分析,可得到同一空间位置、角度和不同波长的散射光强变化序列,其中携带了散射光所历经光程的气体分子吸收,包含了光谱信息,对其进行基线拟合,可得到对应的吸收光谱线形。
19.具体地,根据像元中携带的吸收信息,采用描述多普勒展宽线形的高斯函数和描述碰撞展宽线形的洛伦兹函数之间的卷积,得到采用福伊特函数描述的气体吸收谱线展宽的线形。
20.分子吸收具有选择性,只有当入射光波数(或波长)与气体分子的某个跃迁共振时才发生吸收。预设气体分子为氧气分子时,则吸收作用取决于氧气分子的特征吸收光谱。气体分子的吸收遵循朗伯比尔定律,气体吸收光谱取决于分子结构、对应跃迁能级、所处环境条件,受所处环境条件影响而呈现不同的线形。在低压环境或高温环境中,气体吸收谱线增宽的一个主要因素是多普勒展宽,该展宽机制源于气体分子的热运动。分子热运动规律遵循麦克斯韦

玻尔兹曼统计分布,温度越高分子运动越剧烈,多普勒展宽效应增强,可通过高斯函数描述。气体分子间碰撞会引起分子吸收谱线增宽,可通过洛伦兹函数描述。在很多情况下,会同时存在多普勒展宽和碰撞展宽,气体谱线展宽由这二者共同决定。因此,可通过高斯函数和洛伦兹函数的卷积,得到福伊特函数voigt描述的气体谱线展宽的线形。具体如下:
[0021][0022]
其中,φ
d
表示通过高斯函数描述的谱线展宽,φ
l
表示通过洛伦兹函数描述的碰撞展宽,u表示高斯线形的波数,v表示福伊特线形的波数,v0表示高斯线形的中心波数,a
v
和w为无量纲数值,分别表示为:
[0023][0024]
s103、对所述吸收谱线进行线形函数拟合,得到对应的吸收谱线参数,以计算所述后向散射光的散射区域内的气体浓度、大气温度、气体压强和气流速度。
[0025]
通过ai

lm拟合算法,对吸收谱线线形进行voigt线形函数拟合,根据当前迭代过程的迭代拟合结果和对应的拟合参数,预测下一迭代过程中初值赋值的对应参量,通过反馈迭代,确定标准迭代初值,作为最优迭代初值,实现了多参量非线性函数的有效拟合,以得到对应的吸收谱线参数。
[0026]
该方法不会受到信号信噪比、干扰等影响,有效解决了拟合误差大、拟合失败的问题。voigt线形拟合依据的洛伦兹、高斯两种线形的卷积,在一定程度上能够有效去除一些干扰和误差,提高测量结果准确性,所以测量标准差相对较小。并且,voigt线性拟合结果更加稳定,可用于后续的温度等大气数据的计算。
[0027]
在一个实施例中,根据吸收谱线计算后向散射光的散射区域内的气体浓度、大气温度、气体压强和气流速度,具体通过以下方式:
[0028]
首先,通过以下公式,计算得到积分吸光度:
[0029][0030]
其中,a(v)表示气体分子的积分吸光度,i
t
(v)表示出射光强度,i0(v)表示入射光强度,α(v)表示气体分子吸收系数,p表示气体系统总压强,c表示测量路径上气体平均浓度,s(t)表示气体分子吸收线线强,l表示测量路径的光程长度,φ(v)为吸收线线形函数,该函数为归一化函数,其函数积分为1。
[0031]
在计算气体浓度时,根据积分吸光度、气体压强、测量光程长度、吸收线线强,可计
算得到气体浓度。具体通过以下公式计算:其中参数的含义参照公式(4)。
[0032]
在计算大气温度时,由于热力学平衡态下,分子能级的粒子数布局满足玻尔兹曼分布,吸收光谱的线强与对应能级跃迁的粒子数目以及跃迁概率有关,就特定吸收谱线而言,其光谱线强的大小只与温度有关。因此利用同一气体分子的两个吸收谱线,就可以在一定温度范围内反映所处环境的温度。由于吸收光谱的积分吸光度与压力、吸收谱线的线强、物质的摩尔浓度和光程长度密切相关,两条吸收谱线的积分吸光度在同一环境中测得,并且认为物质的摩尔浓度,压力,温度和光程长度相同,则两个吸收峰的积分吸光度之比可以简化为线强之比,即:
[0033][0034][0035]
其中,a1、a2分别是两个吸收谱线的积分吸光度,e
″1、e
″2分别为两个谱线的低跃迁能级的能量,s1(t0)、s2(t0)分别是两个谱线在参考温度t0时的谱线线强,k
b
为玻尔兹曼常数,c为真空中光速。在实践中,参考温度t0时的谱线线强s(t0)可以通过光谱数据库(hitran或hitemp)查询得到,或者由实验测量获得。
[0036]
在计算大气压力时,根据测量得到的气体的吸光度和谱线的洛伦兹展宽,可以计算气体压强,具体通过以下公式实现:
[0037][0038]
δv
c
=p∑
j
x
j
γ
j
(t)=px
o2
γ
o2
+px
air
γ
air
ꢀꢀꢀ
公式(9)
[0039]
其中,δv
c
表示洛伦兹线宽,a表示气体的积分吸光度。
[0040]
在计算速度时,由于空气的气体流动,散射光的中心频率和吸收光谱的中心频率都发生多普勒频移,气体分子运动速度v与多普勒频移的关系为:其中,λ为激光中心波长,δf
d
为多普勒频移,θ是激光束与粒子运动方向的夹角。
[0041]
在空气中发生的散射和分子吸收都有多普勒效应,可以分别用于空气速度测量,通过吸收光谱计算气流速度时无须鉴频器。并且通过同时测速、并以二者的测速结果做比对,可相互印证,获得更为精准的测量结果。二者同时应用时,还可以用于对复杂流场的湍流结构进行测量和诊断。
[0042]
s104、根据所述图像数据中各像元位置对应的像元数据序列,计算相应空间位置处的空气密度。
[0043]
同一帧图像数据中不同的像元位置分别对应了不同角度和空间位置的空间散射
位置。通过对各像元位置在各帧图像数据中的像元数据序列进行分析,可计算对应空间位置的散射(或散射)强度,进而反演出相应的米散射光强和瑞利散射光强,求得大气中气溶胶和大气分子的密度,计算出该帧图像数据所在空间位置处的空气密度。
[0044]
在一个实施例中,调谐激光在空气中经过散射后形成了相应的图像数据,分别针对各像元计算相应空间散射位置的散射光强,以测量该帧图像数据所在空间位置的空气密度。
[0045]
具体地,确定图像数据中各像元对应的空间散射位置后,根据各像元携带的散射信息和吸收信息,计算其对应的散射光强。
[0046]
第一,后向散射光中包括气体分子的瑞利散射和气溶胶的米散射,根据各像元携带的散射信息,以及气体分子的直径和调谐激光的波长之间的关系,分别采用米散射模型或瑞利散射模型,对像元散射能量进行描述。调谐激光与空气粒子之间的散射作用取决于激光波长与粒子直径的大小关系,若粒子直径与波长相近或大于波长,通过米散射模型对像元的散射能量进行描述,若粒子直径远小于波长,通过瑞利散射模型对像元的散射能量进行描述。
[0047]
第二,根据像元中携带的吸收信息,采用描述多普勒展宽线形的高斯函数和描述碰撞展宽线形的洛伦兹函数之间的卷积,得到采用福伊特函数描述的气体吸收谱线展宽的线形。具体过程在s102中已详细说明,对此不再进行赘述。
[0048]
第三,根据像元散射能量的描述和气体谱线展宽的线形的描述,计算各像元对应的散射光强。
[0049]
具体地,调谐激光在空气中同时产生米散射、瑞利散射和分子吸收时,通过以下公式(4)可计算出像元对应的散射光强:
[0050][0051]
其中,i0表示入射光强度,v表示入射激光的波数,λ表示波长,r为空间散射位置到观察点的距离,n为折射率,w表示质量浓度,ρ表示粒子密度,d表示粒子平均粒径,i1(θ)和i2(θ)分别为垂直和平行于散射平面的偏振散射强度函数,是与粒子粒径d、折射率n、波长λ有关的函数,α(v)表示气体分子吸收系数,吸收系数表示单位浓度以及单位光程长度下气体的吸光度,代表了气体分子对不同波数光的吸收能力,c表示测量路径上的气体平均浓度,l表示测量路径的光程长度。公式(4)中包含了散射和吸收的多重过程,信号极为复杂,但是每种信号的特征并不相同。米散射与瑞利散射具有不同的角度分布特征,且具有不明显的波长特征,即表现为宽带特征,而吸收表现为随波长的窄带特征。通过对接收到的后向散射信号进行空间角度、波长特征的测量和解析,获得不同信号的分离,进而可计算信号中携带的微粒密度、分子密度、压力、温度和速度等数据。
[0052]
在一个实施例中,在计算得出各像元对应的散射光强后,可根据不同像元位置对应的亮度数据序列对空气中的气体分子密度和气溶胶密度进行估计,从而计算出该帧图像数据相应空间位置处的空气密度。
[0053]
具体地,根据图像数据中各像元的灰度值,分别计算对应的散射光能量。本技术实施例采用iccd作为图像传感器。散射光能量是通过iccd进行拍照,采集对应位置的像素灰度值。通过以下公式将测得的灰度值转换为散射光能量:
[0054]
n
photon
=s[cnts]/sens[cnts/ph.el]/qe
ꢀꢀꢀ
公式(11)
[0055]
其中,s[cnts]表示测得的像素灰度值;sens[cnts/ph.el]表示灵敏度,即每个光子代表的灰度值大小,与增益数值有关,增益值为100%时,sens[cnts/ph.el]为620cnts/ph.el;qe为量子效率,通过qe[%]=pr[ma/w]
×
124/λ(nm)计算得到,在760nm波长下,qe=3%。在测量过程中,在iccd镜头前放置滤光片,透光率为τ,最终iccd探测到的散射光通量其中,p为激光功率,η为转换效率,α为修正因子,取1。
[0056]
进一步地,探测获得的散射光中包括了气体分子的瑞利散射和气溶胶的米散射,由于激光光路较长,一般在十几米的范围内,该高度下的大气中包含气体分子和较多的气溶胶粒子,因此在描述散射光能量时,需要将米散射和瑞利散射分开,分别计算出米散射光强、瑞利散射光强与散射光能量的占比。
[0057]
首先,针对不同的像元,分别计算出相应的米散射光强和瑞利散射光强。空气中的米散射来源于空气中的气溶胶,米散射强度可以由麦克斯韦方程组推导出来,给出均质球形粒子在电磁场中对平面波散射的精确解。波长为λ、光强度为i0的单色光,照射到折射率为n的粒子上发生散射,p设为观察点。通过公式(12),计算米散射的光强:
[0058][0059]
其中,r为空间散射位置到观察点的距离,n
mie
为散射区域内的气溶胶粒子数;r和z轴组成的平面为散射面,θ为散射角。i1(θ)和i2(θ)分别为垂直和平行于散射平面的偏振散射强度函数,是与微粒粒径d、折射率m、波长λ有关的函数。随着粒子直径增大,米散射效率因子随波长基本不变,具有几乎相同强度的散射,米散射的强度跟波长无关。
[0060]
激光在空气中的瑞利散射,主要来自于空气中的气体分子和极细微气溶胶,散射光强正比于粒子密度。分子热运动会引起多普勒频移以及谱线加宽,因此利用瑞利散射可以测量空气的密度、温度和速度等重要参数。瑞利散射是弹性散射,不改变散射光的波长。通过公式(13),计算瑞利散射的光强:
[0061][0062]
其中,i0是入射光强度;为散射距离与y轴之间的夹角;为散射距离与z轴之间的夹角;d为散射分子半径;n
r
为粒子的折射率;r为散射距离;n
ray
为散射区域内的气体分子数。水平方向的散射光强呈对称性变化,垂直方向的散射光强分布与散射的角度无关。
[0063]
其次,计算米散射光强、瑞利散射光强与散射光能量的占比。将公式(10)表示为以下形式:
[0064]
i(v)=(m+r)i0(v)
ꢀꢀꢀ
公式(14)
[0065]
其中,m表示米散射在整个散射能量中所占比重,r表示瑞利散射在整个散射能量中所占比重。通过对不同散射角度以及散射距离下,m和r的计算,可以得到在760nm激光波长下,米散射和瑞利散射的比例。由此得到大气密度。
[0066]
[0067][0068]
进一步地,激光经过散射后会形成相应的散射区,需要求得各像元对应的散射区体积。对于激光,其光束宽度很窄,因此散射区可近似看作一个圆台型,散射区体积为:
[0069][0070]
其中,h1=r1+r2φ2/2,d1=h1φ1/2,分别为大圆锥的高度和底面半径;h2=r1‑
r2φ2/2,d2=h2φ1/2,分别为小圆锥的高度和底面半径。
[0071]
更进一步地,由公式(15)和公式(16)可知,在已知散射距离r和散射角度θ下,米散射的贡献m可以视为n
mie
的线性函数,其余参数可以计算得到。同理,r可以表示为n
ray
的线性函数。通过二元线性回归法对n'
mie
和n'
ray
进行估计:
[0072]
i
n
=(μ0+n'
mie
·
v
n
·
x
1n
+n'
ray
·
v
n
·
x
2n
)
·
i
0n
ꢀꢀꢀꢀꢀ
公式(18)
[0073]
其中,n
mie
和n
ray
表示散射区内的气溶胶数和气体分子数,与散射体积大小有关,n'
mie
和n'
ray
表示每平方米中所含气溶胶粒子数和分子数,μ0为回归系数,v
n
为每个光斑散射区域大小,单位为平方米,区域大小,单位为平方米,
[0074]
在一定高度下,单位体积内的气溶胶数和气体分子数保持不变,因此,可采用最小二乘法对n'
mie
和n'
mie
进行估计,并对回归方差进行分析,从而确定出最佳回归结果,得到相应的气溶胶密度和气体分子密度。
[0075]
更进一步地,根据气溶胶密度和气体分子密度,及单个气溶胶粒子和的单个气体分子的质量,计算出像元所在图像数据的空间位置处的空气密度。
[0076]
在一个实施例中,通过s101~s104可根据不同波长的调谐激光,计算出不同光路对应的散射区域内的大气数据。针对多条光路测得的多组空气密度、气体浓度、大气温度、气体压强和气流速度,对多组数据均进行平均处理,求多组大气数据的平均值,并将处理结果作为标准大气数据,这样可提高大气数据的准确性,减少误差。
[0077]
在一个实施例中,在三维坐标系的x、y、z轴三个方向上均设置有图像传感器,飞行器在飞行过程中可从三个方向上分别测得相应的空气密度、气体浓度、大气温度、气体压强和气流速度,并作出对应的三维廓线,为飞行器的安全飞行提供数据支持。
[0078]
本技术将成像技术与光谱技术结合在一起,在探测物体空间特征的同时对每个空间像元吸收光谱进行探测,以较高帧率获取目标的光路图像。由于图像传感器高光谱分辨率的巨大优势,在获取空间散射信号的同时获取气体分子吸收光谱信息,可通过图像传感器获取的信息反演空气的密度、氧气浓度、空气温度、压力和速度等多个参量,实现飞行器所需要的大气数据的多种参数的同步传感和测量,提高了高光谱数据处理的自动化和智能化水平,在航空、航天器中有广泛的应用。
[0079]
以上为本技术提出的方法实施例。基于同样的发明构思,本技术实施例还提供了一种激光成像光谱雷达遥感探测设备,其结构如图2所示。
[0080]
图2为本技术实施例提供的一种激光成像光谱雷达遥感探测设备结构示意图。如图2所示,设备包括处理器201、及存储器202,其上存储有可执行代码,当可执行代码被执行
时,使得处理器201执行如上的一种激光成像光谱雷达遥感探测方法。
[0081]
在本技术的一个实施例中,处理器201用于基于图像传感器,按帧获取不同波长的调谐激光后向散射光对应的图像数据;针对所述图像数据中的各像元位置,根据该像元在各帧图像数据中的时间序列描述相应激光传输路径的散射光强变化,进行基线重构,得到对应的吸收谱线;对所述吸收谱线进行线形函数拟合,得到对应的吸收谱线参数,以计算所述后向散射光的散射区域内的气体浓度、大气温度、气体压强和气流速度;根据所述图像数据中各像元位置对应的像元数据序列,计算相应空间位置处的空气密度。
[0082]
本技术中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于设备实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
[0083]
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个
……”
限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0084]
以上所述仅为本技术的实施例而已,并不用于限制本技术。对于本领域技术人员来说,本技术可以有各种更改和变化。凡在本技术的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本技术的权利要求范围之内。
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1