基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法

文档序号:10651871阅读:364来源:国知局
基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法
【专利摘要】基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置及重建方法,涉及弥散介质光学参数分布的重建探测装置和光学参数分布的重建方法。为了解决传统接触式测量的光学系数重建过程中存在的装置复杂和无法区分各个方向的辐射强度信号的问题。弥散介质光学参数分布的重建探测装置包括激光控制器、激光头、至少一个光场相机和数据采集处理系统;本发明利用光场相机获取调频激光作用下弥散介质边界各个方向上的辐射强度信息,通过模拟弥散介质内的红外辐射传输过程,并结合最优化方法,重建得到介质内部的吸收、散射系数分布图像,从而探测得到弥散介质的内部结构。本发明适用于弥散介质光学参数分布的重建领域。
【专利说明】
基于光场相机与调频激光的弥散介质光学参数分布的重建探 测装置及重建方法
技术领域
[0001] 本发明涉及弥散介质光学参数分布的重建探测装置和光学参数分布的重建方法。
【背景技术】
[0002] 弥散介质内部光学参数分布的重建是通过分析介质边界的测量信号来反演内部 光学参数场,由于介质内部的吸收、散射系数分布与介质内部的结构相关,所以弥散介质内 部吸收、散射系数的重建技术有助于探测介质内部几何结构。作为一种有效的探测技术,弥 散介质光学参数场重建广泛应用于无损探测、光学成像、红外遥感、信息处理和故障诊断等 领域。
[0003] 近红外激光作用于弥散介质,会得到与介质内部吸收、散射系数分布相关的辐射 信号,所以重建介质内部吸收、散射系数首先需要建立模拟激光在介质内的传输模型,以计 算介质边界的出射辐射信号。根据采用激光光源的不同,辐射传输模型可以分为稳态(连续 激光)、时域(脉冲激光)和频域(调频激光)模型,在这三种模型中,采用调频激光的频域辐 射传输模型可以避免时域模型的技术限制,同时能够提供比稳态模型更多的测量信息,成 为在辐射反问题中最具发展前景的计算模型。
[0004] 传统弥散介质光学参数分布的重建过程中,辐射信号的测量主要以接触式的光纤 测量为主,种接触式测量方法需要多光纤通道和光开关来实现多点信号测量,装置复杂;接 触式的光纤测量过程中,进行每次测量时都需要设置多光纤通道,导致该方法的测量效率 较低,并且利用该方法进行光学系数重建时,在测量次数较少的情况下,存在病态特性和吸 收、散射系数之间的串扰问题,所以需要进行多次测量克服这种问题,使得利用光纤进行弥 散介质光学参数分布的重建方法的测量效率较低;同时,利用光纤进行弥散介质光学参数 分布的重建方法所得的信号为辐射热流密度,无法区分各个方向的辐射强度信号。

【发明内容】

[0005] 本发明为了解决传统接触式测量的光学系数重建过程中存在的装置复杂和无法 区分各个方向的辐射强度信号的问题。进而提出一种基于光场相机与调频激光的弥散介质 光学参数分布的重建探测装置和基于光场相机与调频激光的弥散介质光学参数分布的重 建方法。
[0006] 1、基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置,包括:激 光控制器、激光头、至少一个光场相机和数据采集处理系统;
[0007]激光控制器的一端(激光控制信号输出端)连接激光头的激光控制信号输出端,激 光控制器的另一端连接数据采集处理系统;激光头正对弥散介质;光场相机分布在弥散介 质的周围,光场相机的信号输出端连接数据采集处理系统的输入端;数据采集处理系统对 获得的辐射场信号进行处理,并计算得到得到弥散介质的光学参数场的分布。
[0008] 2、利用所述的基于光场相机与调频激光的弥散介质光学参数分布的重建探测装 置进行弥散介质光学参数分布的重建方法,所述方法包括以下步骤:
[0009] 步骤一:开启激光控制器,使激光头发射出的调频激光入射到弥散介质上,然后将 弥散介质旋转η次,并利用激光照射弥散介质,n>0;
[0010] 激光头每发射一次激光则利用光场相机采集一次弥散介质透射或反射出的辐射 场信号,然后将获得的所有辐射场信号发送至数据采集处理系统中,数据采集处理系统分 别对其获得的辐射场信号进行处理,获得弥散介质边界的出射光谱辐射强度值/=作为测 量信号,S代表光源照射序号,d代表探测点位置序号;
[0011]步骤二:假设弥散介质的光学参数场为μ13,将μ13带入频域辐射传输方程,计算得到 介质边界的透反射辐射强度信号与步骤一中的测量信号<1/构成目标函数F(J);
[0012] 步骤三:根据共辄梯度法更新弥散介质光学参数场的分布值:yk=yk4+Ay,
[0013] k = l,2,···;
[0014] 步骤四:根据第k步迭代得到的光学参数分布μ15,通过频域辐射传输方程计算介质 边界的辐射强度信号计算目标函数F(y k),如果目标函数值小于极少值ε,执行步骤六; 否则,执行步骤五;
[0015] 步骤五:如果迭代次数k达到最大迭代次数Κ,使当前的光学参数场分布作为初值, 重新开始迭代,即令Jiyk,即将值赋予μ* 3,执行步骤二;否则,执行步骤三;
[0016] 步骤六:将当前迭代得到的光学参数场作为重建结果,结束反演过程。
[0017]本发明具有以下效果:
[0018] 在目前的弥散介质吸收、散射系数分布同时重建研究中,辐射信号的测量主要以 接触式的光纤测量为主,通过光场相机测量光场信号的技术尚未得到应用。本发明利用光 场成像技术进行弥散介质光学参数分布重建,与传统成像技术相比,获取的信息多出了两 个自由度,因而在信息重建中能获得更加丰富的信息,而且具有系统简单、硬件设备成本 低、结构紧凑,成像效果好,数据处理方便等优点。
[0019] 本发明所述的基于光场相机与调频激光的弥散介质光学参数分布的重建方法能 够分辨出不同方向的辐射强度。
[0020] 本发明的方法无需每次测量时都设置过多的光纤,且测量过程简单;与传统接触 式测量相比,本发明的测量效率能够成倍提升。而且本发明在测量较少的情况下也能够有 效的克服光学系数重建中的病态特性和吸收、散射系数之间的串扰问题。
【附图说明】
[0021 ]图1本发明重建探测装置的结构示意图;
[0022] 图2为基于光场相机与调频激光的弥散介质光学参数分布的重建方法的流程图。
【具体实施方式】
【具体实施方式】 [0023] 一:结合图1说明本实施方式,
[0024] 基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置,包括:激光 控制器1、激光头2、至少一个光场相机4和数据采集处理系统5;
[0025] 激光控制器1的一端(激光控制信号输出端)连接激光头2的激光控制信号输出端, 激光控制器1的另一端连接数据采集处理系统5;激光头2正对弥散介质3;光场相机4分布在 弥散介质3的周围,光场相机4的信号输出端连接数据采集处理系统5的输入端;数据采集处 理系统5对获得的辐射场信号进行处理,并计算得到得到弥散介质的光学参数场的分布。
【具体实施方式】 [0026] 二:结合图1和图2说明本实施方式,
[0027] 利用基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置进行弥 散介质光学参数分布的重建方法,包括以下步骤:
[0028]步骤一:开启激光控制器1,使激光头2发射出的调频激光入射到弥散介质3上,然 后将弥散介质3旋转η次,并利用激光照射弥散介质3,n>0;
[0029] 激光头2每发射一次激光则利用光场相机4采集一次弥散介质透射或反射出的辐 射场信号,然后将获得的所有辐射场信号发送至数据采集处理系统5中,数据采集处理系统 5分别对其获得的辐射场信号进行处理,获得弥散介质边界的出射光谱辐射强度值作为 测量信号,s代表光源照射序号,d代表探测点位置序号;
[0030] 步骤二:假设弥散介质的光学参数场为J,将J带入频域辐射传输方程,计算得到 介质边界的透反射辐射强度信号/&,与步骤一中的测量信号构成目标函数Fb,;
[0031] 步骤三:根据共辄梯度法更新弥散介质光学参数场的分布值: 2,…;
[0032] 步骤四:根据第k步迭代得到的光学参数分布μ15,通过频域辐射传输方程计算介质 边界的辐射强度信号/&,计算目标函数F(y k),如果目标函数值小于极少值ε,执行步骤六; 否则,执行步骤五;
[0033]步骤五:如果迭代次数k达到最大迭代次数Κ,使当前的光学参数场分布作为初值, 重新开始迭代,即令'ιμ15,即将值赋予μ*3,执行步骤二;否则,执行步骤三;
[0034]步骤六:将当前迭代得到的光学参数场作为重建结果,结束反演过程。
[0035]【具体实施方式】三:
[0036] 本实施方式所述的光学参数场μ包括弥散介质的吸收系数ya和散射系数ys,弥散介 质的吸收系数和散射系数y s的分布图像是同时重建的。
[0037] 其它步骤及参数与【具体实施方式】二相同。
[0038]【具体实施方式】四:
[0039] 本实施方式所述步骤二和步骤四中所述的频域辐射传输方程通过离散坐标法求 解,
[0040] 频域辐射传输方程的表达式如下:
[0042]其中,i表示虚数单位,ω为调制频率,c为弥散介质中的光速,Ω为辐射传输方向, ▽表示梯度;ya、ys分别为ya、ys中的元素;r为空间位置,I(r, Ω,ω )为在t时亥I」、位置为r、调 频为ω的辐射强度,Ω 7为辐射入射方向,Ω 7表示立体角;Φ ( Ω、Ω )为弥散介质3的散射 相函数;(1Ω'表示微分。
[0043]其它步骤及参数与【具体实施方式】二或三相同。
【具体实施方式】 [0044] 五:
[0045]本实施方式所述步骤四中目标函数方程的表达式如下:
[0047] 其中,!Κμ)为正则化项,由广义马克尔夫随机场模型构建得到,用以克服反问题的 病态特性。
[0048] 其它步骤及参数与【具体实施方式】二至四之一相同。
【具体实施方式】 [0049] 六:
[0050]本实施方式所述步骤三采用的共辄梯度法具体描述如下:
[0051]
[0052] 其中,ak为第k次迭代的步长,可由一维搜索得到;dk为第k次迭代的下降方向,由当 前目标函数梯度和上一次迭代的下降方向决定:
[0053] dk =-VF + ^ld!iA
[0054] 其中,VP为目标函数关于待重建参数的梯度,妒为第k次迭代的共辄系数。
[0055] 其它步骤及参数与【具体实施方式】二至五之一相同。
【主权项】
1. 基于光场相机与调频激光的弥散介质光学参数分布的重建探测装置,其特征在于所 述装置包括:激光控制器(1)、激光头(2)、至少一个光场相机(4)和数据采集处理系统(5); 激光控制器(1)的一端连接激光头(2)的激光控制信号输出端,激光控制器(1)的另一 端连接数据采集处理系统(5);激光头(2)正对弥散介质(3);光场相机(4)分布在弥散介质 (3)的周围,光场相机(4)的信号输出端连接数据采集处理系统(5)的输入端;数据采集处理 系统(5)对获得的辐射场信号进行处理,并计算得到得到弥散介质的光学参数场的分布。2. 利用权利要求1所述的装置进行弥散介质光学参数分布的重建方法,其特征在于所 述方法包括以下步骤: 步骤一:开启激光控制器(1),使激光头(2)发射出的调频激光入射到弥散介质(3)上, 然后将弥散介质(3)旋转η次,并利用激光照射弥散介质(3),n>0; 激光头(2)每发射一次激光则利用光场相机(4)采集一次弥散介质透射或反射出的辐 射场信号,然后将获得的所有辐射场信号发送至数据采集处理系统(5)中,数据采集处理系 统(5)分别对其获得的辐射场信号进行处理,获得弥散介质边界的出射光谱辐射强度值/? 作为测量信号,s代表光源照射序号,d代表探测点位置序号; 步骤二:假设弥散介质的光学参数场为J,将J带入频域辐射传输方程,计算得到介质 边界的透反射辐射强度信号/&,与步骤一中的测量信号构成目标函数F(J); 步骤三:根据共辄梯度法更新弥散介质光学参数场的分布值:Pk = Pk+1+ ΑμΛ=1,2,…; 步骤四:根据第k步迭代得到的光学参数分布yk,通过频域辐射传输方程计算介质边界 的辐射强度信号 ,计算目标函数F(yk),如果目标函数值小于极少值ε,执行步骤六;否则, 执行步骤五; 步骤五:如果迭代次数k达到最大迭代次数K,使当前的光学参数场分布作为初值,重新 开始迭代,即令J=Pk,即将以喻值赋予Λ执行步骤二;否则,执行步骤三; 步骤六:将当前迭代得到的光学参数场作为重建结果,结束反演过程。3. 根据权利要求2所述的弥散介质光学参数分布的重建方法,其特征在于所述的光学 参数场μ包括弥散介质的吸收系数和散射系数y s,弥散介质的吸收系数和散射系数以5的 分布图像是同时重建的。4. 根据权利要求2或3所述的弥散介质光学参数分布的重建方法,其特征在于步骤二和 步骤四中所述的频域辐射传输方程通过离散坐标法求解, 频域辐射传输方程的表达式如下:其中,i表示虚数单位,ω为调制频率,C为弥散介质中的光速,Ω为辐射传输方向,▽表 示梯度;ya、ys分别为ya、ys中的元素;r为空间位置,I(r, Ω,ω )为在t时亥I」、位置为r、调频为 ω的福射强度,Ω '为福射入射方向,Ω '表不立体角;φ ( Ω ',Ω )为弥散介质(3)的散射相 函数;(1Ω '表不微分。5. 根据权利要求4所述的弥散介质光学参数分布的重建方法,其特征在于步骤四中目 标函数方程的表达式如下:其中,Φ(μ)为正则化项,由广义马克尔夫随机场模型构建得到。6. 根据权利要求5所述的弥散介质光学参数分布的重建方法,其特征在于步骤三采用 的共轭梯度法具体描沭如下:其中,ak为第k次迭代的步长;dk为第k次迭代的下降方向,由当前目标函数梯度和上一 次迭代的下降方向决定:其中,VF为目标函数关于待重建参数的梯度,栌为第k次迭代的共辄系数。7. 根据权利要求2或3所述的弥散介质光学参数分布的重建方法,其特征在于步骤四中 目标函数方程的表达式如下:其中,s和d分别代表光源和探测点的数量,Φ(μ)为正则化项,由广义马克尔夫随机场模 型构建得到。8. 根据权利要求2或3所述的弥散介质光学参数分布的重建方法,其特征在于步骤三采 用的共辄梯度法具体描述如下:其中,ak为第k次迭代的步长;dk为第k次迭代的下降方向,由当前目标函数梯度和上一 次佚代的下隆方向决宙:其中,VF为目标函数关于待重建参数的梯度,栌为第k次迭代的共辄系数。9. 根据权利要求4所述的弥散介质光学参数分布的重建方法,其特征在于步骤三采用 的共辄梯度法具体描述如下:其中,ak为第k次迭代的步长;dk为第k次迭代的下降方向,由当前目标函数梯度和上一 次迭代的下降方向决定:其中,VF为目标函数关于待重建参数的梯度,栌为第k次迭代的共辄系数。
【文档编号】G01N21/17GK106018286SQ201610330924
【公开日】2016年10月12日
【申请日】2016年5月18日
【发明人】齐宏, 乔要宾, 阮世庭, 阮立明, 谈和平, 周鹜, 蔡小舒
【申请人】哈尔滨工业大学
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1