基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置及其重建方法
【专利摘要】基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置及其重建方法，属于红外光学成像领域。目前缺少通过微透镜阵列测量光场信号的技术方法。基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置，数据采集处理系统(7)连接激光控制器(1)和三个微透镜阵列光场相机；激光控制器(1)连接激光头(2)，激光头(2)和三个微透镜阵列光场相机围成的区域中心设置弥散介质(3)。基于微透镜阵列与脉冲激光的弥散介质光学参数场重建方法，包括介质边界出射辐射测量、近红外脉冲激光在弥散介质的传输计算、光学参数场重建环节。本发明实现对弥散介质光学参数场的重建，为近红外光学成像和红外探测技术提供新的技术手段。
【专利说明】
基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置 及其重建方法
技术领域
[0001] 本发明设及一种通过微透镜阵列测量光场信号的装置和技术方法，具体设及基于 微透镜阵列与脉冲激光的弥散介质光学参数场重建装置及其重建方法。
【背景技术】
[0002] 弥散介质内部光学参数场的重建是通过分析介质边界的测量信号来反演内部光 学参数场，由于介质内部的光学参数分布与介质内部的结构相关，所W弥散介质内部光学 参数场的重建技术有助于探测介质内部几何结构。作为一种有效的探测技术，弥散介质光 学参数场重建广泛应用于无损探测、光学成像、红外遥感、信息处理和故障诊断等领域。
[0003] 超短脉冲激光持续时间非常短，其脉冲宽度与激光在介质内传输的时间相当，达 到皮秒和飞秒量级，此时测量信号(如福射信号)随时间的变化与光子的传播速度相比不可 忽略，其福射传输过程中必须考虑瞬态效应。与稳态福射问题相比，瞬态福射具有大量的随 时间变化的信息可W利用，瞬态光学性质及其时间衰减蕴藏了介质内的一些重要信息，尤 其是关于介质内部结构的信息。
[0004] 光场成像技术与传统成像技术相比，获取的信息多出了两个自由度，因而在信息 重建中能获得更加丰富的信息。微透镜阵列是最常用的光场获取方式，具有光学系统简单、 硬件设备成本低、结构紧凑，标定简单，成像效果好，数据处理方便等优点。利用微透镜阵列 获取福射光场信息，提供比传统接触式测量更为丰富的测量信号，能够有效的克服光学参 数场重建中的病态特性和吸收散射系数之间的串扰问题。在目前的弥散介质光学参数场重 建研究中，福射信号的测量主要W接触式的光纤测量为主，缺少通过微透镜阵列测量光场 信号的技术方法。
【发明内容】
[0005] 本发明的目的是为了解决目前缺少通过微透镜阵列测量光场信号的技术方法的 问题，而提出一种基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置，和一种基 于微透镜阵列与脉冲激光的弥散介质光学参数场的重建方法。
[0006] -种基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置，其组成包括： 激光控制器、激光头、弥散介质、第一微透镜阵列光场相机、第二微透镜阵列光场相机、第= 微透镜阵列光场相机和数据采集处理系统，数据采集处理系统连接激光控制器，激光控制 器连接激光头;数据采集处理系还同时与第一微透镜阵列光场相机、第二微透镜阵列光场 相机、第=微透镜阵列光场相机连接；
[0007] 其中，激光头、第一微透镜阵列光场相机、第二微透镜阵列光场相机、第=微透镜 阵列光场相机的中屯、点在同一平面上;激光头和第二微透镜阵列光场相机中屯、点的连线， 与第一微透镜阵列光场相机和第=微透镜阵列光场相机中屯、点的连线垂直相交，两条中屯、 点连线的交点处分别为各自中点，且在两条中屯、点连线的交点处设置弥散介质。
[0008] -种基于微透镜阵列与脉冲激光的弥散介质光学参数场的重建方法，所述方法通 过W下步骤实现，
[0009] 步骤一、开启激光控制器，使激光头发射的脉冲激光射入弥散介质，然后将弥散介 质W第一微透镜阵列光场相机和第立微透镜阵列光场相机中屯、点的连线为轴旋转90°，使 脉冲激光作用于弥散介质的另外一个边界，直至弥散介质的四个边界都受到一次激光头发 射的脉冲激光；
[0010] 激光头每发射一次脉冲激光时都利用第一微透镜阵列光场相机、第二微透镜阵列 光场相机、第=微透镜阵列光场相机进行一次经弥散介质透射或反射出的福射场信号的采 集，再将获得的所有福射场信号发送至数据采集处理系统中；
[0011] 数据采集处理系统分别对获得的福射场信号进行处理，获得弥散介质各边界射出 的光谱福射强度值;其中，S表示光源位置，i表示探测点位置，M表示测量信号；
[0012] 步骤二、设弥散介质的光学参数场为4*\通过时域福射传输方程计算出弥散介质 边界的透射及反射福射强度信号巧,?，将计算得到的弥散介质边界的透射及反射福射强度信 号^,与步骤一采集的福射场信号构成目标函数。(4<^);
[0013] 其中，透射及反射福射强度信号巧,.中上角标P表示估算信号或模拟信号;y表示待 重建的弥散介质的光学参数场；
[0014] 步骤=、利用共辆梯度法进行迭代计算，更新第k次迭代计算得到的弥散介质光学 参数场的分布值：
[0015]
[0016] 式中，k = l,2,…；Ay表示参数的改变量；
[0017] 步骤四、根据步骤=第k次迭代计算得到的弥散介质光学参数场的分布值yk，计算 弥散介质边界的福射强度信号W及目标函数F(yk);
[0018] 若目标函数F(yk)的值小于极值，则执行步骤六；
[0019] 否则，执行步骤五；
[0020] 步骤五、若迭代次数k达到最大迭代次数K，则使当前的弥散介质光学参数场的分 布值作为初值，重新开始迭代，即令4*^ = 4'，返回至步骤二；
[0021] 否则，返回至步骤S;
[0022] 步骤六、将当前迭代计算得到的光学参数场分布值作为重建结果，结束反演过程。
[0023] 本发明的有益效果为：
[0024] 本发明属于一种近红外光学成像技术，具体地说是一种基于微透镜阵列与近红外 脉冲激光提供的弥散介质边界时域测量信号，对介质内部的未知光学参数场进行重建的方 法。本发明通过介质边界出射福射测量、近红外脉冲激光在弥散介质的传输计算、光学参数 场重建环节实现基于微透镜阵列与脉冲激光的弥散介质光学参数场重建方法。利用具有微 透镜阵列的光场相机获取脉冲激光作用下弥散介质边界各个方向上的福射强度信息，通过 模拟弥散介质内的红外福射传输过程，并结合最优化方法，重建得到介质内部的光学参数 场分布图像，从而探测得到弥散介质的内部结构。通过本发明的仿真计算，可实现对弥散介 质光学参数场的重建，为近红外光学成像和红外探测技术提供新的技术手段。
【附图说明】
[0025] 图1为本发明所采用的硬件的结构示意图；
[0026] 图2为基于微透镜阵列与脉冲激光的弥散介质光学参数场重建方法的流程图。
【具体实施方式】
【具体实施方式】 [0027] 一：
[0028] 本实施方式的基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置，结合 图1其组成包括:激光控制器1、激光头2、弥散介质3、第一微透镜阵列光场相机4、第二微透 镜阵列光场相机5、第=微透镜阵列光场相机6和数据采集处理系统7,数据采集处理系统7 连接激光控制器1，激光控制器1连接激光头2;数据采集处理系统7还同时与第一微透镜阵 列光场相机4、第二微透镜阵列光场相机5、第=微透镜阵列光场相机6连接；
[0029] 其中，激光头2、第一微透镜阵列光场相机4、第二微透镜阵列光场相机5、第=微透 镜阵列光场相机6的中屯、点在同一平面上;激光头2和第二微透镜阵列光场相机5中屯、点的 连线，与第一微透镜阵列光场相机4和第立微透镜阵列光场相机6中屯、点的连线垂直相交， 两条中屯、点连线的交点处分别为各自中点，且在两条中屯、点连线的交点处设置弥散介质3。
【具体实施方式】 [0030] 二：
[0031] 与【具体实施方式】一不同的是，本实施方式的基于微透镜阵列与脉冲激光的弥散介 质光学参数场重建装置，所述弥散介质3是指是指对红外光呈吸收、散射特性的介质，所述 弥散介质3包括为生物组织、火焰或者陶瓷中的一种。
【具体实施方式】 [0032]
[0033] 本实施方式的基于微透镜阵列与脉冲激光的弥散介质光学参数场的重建方法，所 述方法通过W下步骤实现，
[0034] 步骤一、开启激光控制器1，使激光头2发射的脉冲激光射入弥散介质3,然后将弥 散介质3W第一微透镜阵列光场相机4和第=微透镜阵列光场相机6中屯、点的连线为轴旋转 90%使脉冲激光作用于弥散介质3的另外一个边界，直至弥散介质3的四个边界都受到一次 激光头2发射的脉冲激光；
[0035] 激光头2每发射一次脉冲激光时都利用第一微透镜阵列光场相机4、第二微透镜阵 列光场相机5、第=微透镜阵列光场相机6进行一次经弥散介质3透射或反射出的福射场信 号的采集，再将获得的所有福射场信号发送至数据采集处理系统7中；
[0036] 数据采集处理系统7分别对获得的福射场信号进行处理，获得弥散介质3各边界射 出的光谱福射强度值4^其中，3表示光源位置，1表示探测点位置，1表示测量信号，是英文 单词Measured的缩写；
[0037] 步骤二、设弥散介质3的光学参数场为tiD，通过时域嬌射传输方程计算出弥散介质 3边界的透射及反射福射强度信号巧,，将计算得到的弥散介质3边界的透射及反射福射强 度信号巧,与步骤一采集的福射场信号构成目标函数。(4<^);
[0038] 其中，透射及反射福射强度信号/P，,中上角标P表示估算信号或模拟信号，是英文 单词Predicted的缩写;y表示待重建的弥散介质的光学参数场；
[0039] 步骤=、利用共辆梯度法进行迭代计算，更新第k次迭代计算得到的弥散介质光学 参数场的分布値：
[0040]
[0041] 式中，k = l,2,…；Ay表示参数的改变量；
[0042] 步骤四、根据步骤=第k次迭代计算得到的弥散介质光学参数场的分布值yk，计算 弥散介质3边界的福射强度信号/]；,. W及目标函数F(yk);
[0043] 若目标函数F(yk)的值小于极值e，则执行步骤六；
[0044] 否则，执行步骤五；
[0045] 步骤五、若迭代次数k达到最大迭代次数K，则使当前的弥散介质光学参数场的分 布值作为初值，重新开始迭代，即令4*^ = 4"，返回至步骤二；
[0046] 否则，返回至步骤S;
[0047] 步骤六、将当前迭代计算得到的光学参数场分布值作为重建结果，结束反演过程。 [004引【具体实施方式】四：
[0049] 与【具体实施方式】=不同的是，本实施方式的基于微透镜阵列与脉冲激光的弥散介 质光学参数场的重建方法，步骤一所述激光头2发射的脉冲激光射入弥散介质3时，激光头2 发射的脉冲激光射入弥散介质3各边界的中屯、点。
【具体实施方式】 [0050] 五：
[0051] 与【具体实施方式】=或四不同的是，本实施方式的基于微透镜阵列与脉冲激光的弥 散介质光学参数场的重建方法，步骤二所述待重建的弥散介质的光学参数场y包括吸收系 数Wa和散射系数iis两部分参数，且两部分参数场同时进行重建。
【具体实施方式】 [0052] 六：
[0053] 与【具体实施方式】五不同的是，本实施方式的基于微透镜阵列与脉冲激光的弥散介 质光学参数场的重建方法，步骤二所述通过时域福射传输方程计算出弥散介质3边界的透 反射福射强度信号/?的过程为，利用离散坐标法对时域福射传输方程进行求解，其中，瞬 态的时域福射传输方程的表达式为：
[0化4]
[0055]式中，t表示时刻，r表示空间位置，Q表示福射传输方向，Q/表示福射入射方向，C 表示弥散介质3中的光速，▽表示梯度；O ( Q /，Q )表示弥散介质3的散射相函数；
[005引 Kr, Q，t)为在t时刻、位置为r、方向为Q的福射强度。
【具体实施方式】 [0化7] 屯：
[0058]与【具体实施方式】=、四或六不同的是，本实施方式的基于微透镜阵列与脉冲激光 的弥散介质光学参数场的重建方法，步骤二所述目标函数。(4<^)的表达式为：
[0化9]
[0060] 式中，S表示光源，d表示探测点，n表示时层的数量，iKii)为正则化项，由广义马克 尔夫随机场模型构建得到，用W克服反问题的病态特性。
【具体实施方式】 [0061] 八：
[0062] 与【具体实施方式】屯不同的是，本实施方式的基于微透镜阵列与脉冲激光的弥散介 质光学参数场的重建方法，步骤=所述利用共辆梯度法进行迭代计算，更新第k次迭代计算 得到的弥散介质光学参数场的分布值的过程为，所述共辆梯度法具体描述为：
[0063] yk 二 yk-i+曰 kdk
[0064] 式中，ak表示第k次迭代的步长，可由一维捜索得到;dk表示第k次迭代的下降方向， 由当前目标函数梯度和上一次迭代的下降方向决定，且dk = -VF+护dW;式中，▽表示▽表 示梯度;e表示共辆系数。
【具体实施方式】 [0065] 九：
[0066] 与【具体实施方式】一、二、四、六或八不同的是，本实施方式的基于微透镜阵列与脉 冲激光的弥散介质光学参数场的重建方法，所述弥散介质3包括为生物组织、火焰或者陶瓷 中的一种。
[0067] 本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域 技术人员当可根据本发明作出各种相应的改变和变形，但运些相应的改变和变形都应属于 本发明所附的权利要求的保护范围。
【主权项】
1. 一种基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置，其组成包括:激 光控制器（1)、激光头(2)、弥散介质（3)、第一微透镜阵列光场相机(4)、第二微透镜阵列光 场相机(5)、第三微透镜阵列光场相机(6)和数据采集处理系统(7)，其特征在于:数据采集 处理系统(7)连接激光控制器(1)，激光控制器(1)连接激光头(2);数据采集处理系统(7)还 同时与第一微透镜阵列光场相机(4)、第二微透镜阵列光场相机(5)、第三微透镜阵列光场 相机(6)连接； 其中，激光头(2)、第一微透镜阵列光场相机(4)、第二微透镜阵列光场相机(5)、第三微 透镜阵列光场相机(6)的中心点在同一平面上;激光头(2)和第二微透镜阵列光场相机(5) 中心点的连线，与第一微透镜阵列光场相机(4)和第三微透镜阵列光场相机(6)中心点的连 线垂直相交，两条中心点连线的交点处分别为各自中点，且在两条中心点连线的交点处设 置弥散介质(3)。2. 根据权利要求1所述基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置， 其特征在于:所述弥散介质(3)是指对红外光呈吸收、散射特性的介质。3. -种利用基于微透镜阵列与脉冲激光的弥散介质光学参数场重建装置的重建方法， 其特征在于:所述方法通过以下步骤实现， 步骤一、开启激光控制器（1)，使激光头（2)发射的脉冲激光射入弥散介质（3)，然后将 弥散介质(3)以第一微透镜阵列光场相机(4)和第三微透镜阵列光场相机(6)中心点的连线 为轴旋转90°，使脉冲激光作用于弥散介质(3)的另外一个边界，直至弥散介质(3)的四个边 界都受到一次激光头(2)发射的脉冲激光； 激光头(2)每发射一次脉冲激光时都利用第一微透镜阵列光场相机(4)、第二微透镜阵 列光场相机(5)、第三微透镜阵列光场相机(6)进行一次经弥散介质（3)透射或反射出的辐 射场信号的采集，再将获得的所有辐射场信号发送至数据采集处理系统(7)中； 数据采集处理系统(7)分别对获得的辐射场信号进行处理，获得弥散介质（3)各边界射 出的光谱辐射强度值其中，s表示光源位置，i表示探测点位置，M表示所测量信号； 步骤二、设弥散介质(3)的光学参数场为μ'通过时域辐射传输方程计算出弥散介质(3) 边界的透射及反射辐射强度信号/&，将计算得到的弥散介质（3)边界的透射及反射辐射强 度信号与步骤一采集的辐射场信号构成目标函数F(J); 其中，透射及反射辐射强度信号中上角标P表示估算信号或模拟信号;μ表示待重建 的弥散介质的光学参数场； 步骤三、利用共辄梯度法进行迭代计算，更新第k次迭代计算得到的弥散介质光学参数 场的分布值： μ1? = μ1?-1+Λμ 式中，k=l，2，···； Δμ表示参数的改变量； 步骤四、根据步骤三第k次迭代计算得到的弥散介质光学参数场的分布值μ15，计算弥散 介质(3)边界的辐射强度信号以及目标函数F(yk); 若目标函数F(yk)的值小于极值，则执行步骤六； 否则，执行步骤五； 步骤五、若迭代次数k达到最大迭代次数K，则使当前的弥散介质光学参数场的分布值 作为初值，重新开始迭代，即令J=Pk，返回至步骤二； 否则，返回至步骤三； 步骤六、将当前迭代计算得到的光学参数场分布值作为重建结果，结束反演过程。4. 根据权利要求3所述基于微透镜阵列与脉冲激光的弥散介质光学参数场的重建方 法，其特征在于:步骤一所述激光头(2)发射的脉冲激光射入弥散介质（3)时，激光头(2)发 射的脉冲激光射入弥散介质(3)各边界的中心点。5. 根据权利要求3或4所述基于微透镜阵列与脉冲激光的弥散介质光学参数场的重建 方法，其特征在于:步骤二所述待重建的弥散介质的光学参数场μ包括吸收系数~和散射系 数y s两部分参数，且两部分参数场同时进行重建。6. 根据权利要求5所述的基于微透镜阵列与脉冲激光的弥散介质光学参数场的重建方 法，其特征在于，步骤二所述通过时域辐射传输方程计算出弥散介质(3)边界的透反射辐射 强度信号/?的过程为，利用离散坐标法对时域辐射传输方程进行求解，其中，瞬态的时域辐 射传输方程的表达式为：式中，t表不时刻，r表不空间位置，Ω表不福射传输方向，Ω '表不福射入射方向，c表不 弥散介质(3)中的光速，▽表示梯度；Φ(Ω'，Ω)表示弥散介质(3)的散射相函数。7. 根据权利要求3、4或6所述的基于微透镜阵列与脉冲激光的弥散介质光学参数场的 重建方法，其特征在于，步骤二所述目标函数F(J)的表达式为：式中，s表示光源，d表示探测点，η表示时层的数量，Φ(μ)为正则化项，由广义马克尔夫 随机场模型构建得到，用以克服反问题的病态特性。8. 根据权利要求7所述的基于微透镜阵列与脉冲激光的弥散介质光学参数场的重建方 法，其特征在于，步骤三所述利用共辄梯度法进行迭代计算，更新第k次迭代计算得到的弥 散介质光学参数场的分布值的过程为，所述共辄梯度法具体描述为： yk=yk^+akdk 式中，ak表示第k次迭代的步长，可由一维搜索得到;dk表示第k次迭代的下降方向，由当 前目标函数梯度和上一次迭代的下降方向决定，且= 式中，▽表示梯度;β表 示共辄系数。
【文档编号】G06F17/50GK106023082SQ201610348998
【公开日】2016年10月12日
【申请日】2016年5月24日
【发明人】乔要宾, 齐宏, 阮世庭, 阮立明, 谈和平, 许传龙, 张彪
【申请人】哈尔滨工业大学