一种基于卷积受限玻尔兹曼机的二维重建方法技术

本发明专利技术公开了一种基于卷积受限玻尔兹曼机的二维重建方法，包括如下步骤：根据目标探测气体进行谱线遴选以确定探测所用红外激光器的中心波长；对待测区域进行离散化并获取经待测气体吸收后的不同投影角度的红外激光信号；计算得到积分吸光度并作为训练集输入；应用模拟退火算法对浓度场进行重建，重建结果作为训练集输出；进行网络训练后确定CRBM的模型结构；使用训练好的模型实现对二维气体浓度的检测。本发明专利技术相比于原有的气体浓度二维重建方法，在完成对卷积受限玻尔兹曼机训练集的构建后，能够精确而迅速的重建出待测气体的二维浓度场，具有广泛的应用前景。具有广泛的应用前景。具有广泛的应用前景。

【技术实现步骤摘要】
一种基于卷积受限玻尔兹曼机的二维重建方法


[0001]本专利技术属于排放测试领域，具体涉及一种基于卷积受限玻尔兹曼机的二维重建方法。

技术介绍

[0002]可调谐激光吸收光谱技术(TDLAS)是一种成熟的气体测量技术，特别是用于温度和浓度的测量。但其本质上是视线，限制了其在空间不均匀的气体测量中的应用。随着医学扫描断层成像技术的发展，扫描系统和图像重建得到了快速发展。为了测量复杂燃烧流场的二维分布情况，美国华盛顿大学的P.J.Emmerman首次提出将计算机断层成像(Computed Tomography，CT)技术与TDLAS技术相结合，测量了甲烷浓度的二维分布。这项技术经过发展目前已广泛应用到燃烧诊断领域，同时也有许多重建算法被研究人员应用到重建过程中：常用的有以代数迭代重建算法(ART)为基础的一系列改进算法；以Radon变换为基础的滤波反投影算法；Tiknonov正则化算法；模拟退火(SA)算法等。应用这些算法在重建过程中，要想得到高质量的重建图像，需要通过增加光线的投影数目实现，不仅使得系统的复杂性和成本进一步升高，还将导致重建时间的增加，难以实现在线测量。同时对于模拟退火(SA)算法这种全局优化算法，虽然其保持了较高的重建精度，但也因为计算效率较低而无法实现实时测量。
[0003]基于能量模型的受限玻尔兹曼机RBM(restricted boltzmann machine)以其简单的人工神经网络形式和快速的学习算法，已经广泛应用于数据降维、语音识别和图像处理等多个机器学习领域，其中RBM具有两层结构，即隐藏层和可视层，它用隐藏层单元对可视层单元的分布进行建模，层内无连接，层间全连接。Honglak Lee等人提出了一种卷积受限玻尔兹曼机CRBM(Convolutional Restricted Boltzmann Machine)模型。CRBM是一种将CNN模型与RBM模型相结合的具有平移不变性的层次化生成模型，支持自顶向下和自底向上的概率推理，并且采用概率最大池进行降维和正则化操作，其模型的设计考虑了图像的二维结构，并且有效地使用卷积核。CRBM在医学图像处理领域已表现出十分优秀的性能。

技术实现思路

[0004]本专利技术其中一个专利技术目的在于提供一种基于卷积受限玻尔兹曼机的二维重建方法，以解决现有方法重建精度对于投影光线数目的依赖，同时在保持较高重建精度的同时提高计算效率，减少重建时间。
[0005]为实现上述目的，本专利技术采用以下技术方案，一种基于卷积受限玻尔兹曼机的二维重建方法方法，包括以下步骤：
[0006]步骤1，根据目标探测气体进行谱线遴选确定红外激光器中心波长；利用HITRAN数据库，针对待测气体，选取4条符合条件的吸收谱线。其中选取多条谱线是为了便于使用模拟退火算法构建目标优化函数。
[0007]步骤2，获取不同投影角度的红外激光吸收信号并计算得到积分吸光度；将待测区
域离散成M
×
N个测量网络，并由M
×
N条正交布置的红外激光光路均匀穿过待测区域进行测量。由于采用了时分复用技术，正交光路中每束红外激光都包含步骤1所选4种波长红外激光的光谱吸收率信息。计算积分吸光度时采用Voigt线型函数。
[0008]步骤3，应用模拟退火(SA)算法对浓度场进行重建；将测量得到的积分吸光度分布定义为P
m
(L
j
,λ
i
)，根据假设的气体浓度分布计算得到的值定义为P
c
(L
j
,λ
i
)，定义理想化的目标优化函数为式中i,j分别为波长和光束序号。加入平滑度用来改善目标病态问题的求解，定义平滑度定义平滑系数为γ
T
，γ
P
，待优化目标函数可表述为F(T,P)＝DA(TX)+γ
T
R
T
(T)+γ
P
R
P
(P)。完成目标函数构建进而重建出浓度场。
[0009]步骤4，构建卷积受限玻尔兹曼机(CRBM)训练集；通过改变待测气体浓度分布，可得到不同的积分吸光度，重复以上工作，得到若干组数据，根据生成的数据作为训练集训练神经网络。步骤2计算得到的积分吸光度作为训练集的输入X，以模拟退火算法重建结果作为训练集的输出Y，共采用15000组样本进行网络训练，根据网络的输出结果对网络进行参数调整；
[0010]步骤5，确定CRBM模型结构；卷积受限玻尔兹曼机CRBM模型在结构上只有两层，即输入层v和隐层h，隐层可以有多个卷积核。CRBM模型是概率生成模型。通过模型的能量函数得到模型所表达的数据的联合概率分布，能够从模型中生成符合该分布的样本。定义能量函数
[0011]其中是模型参数。
[0012]用定义的符号表示能量函数为：
[0013]得到模型的联合概率分布：
[0014]由联合概率分布可以得到条件概率分布：由联合概率分布可以得到条件概率分布：其中为sigmoid函数。
[0015]进一步的CRBM模型通过训练确定模型的参数θ，以拟合给定的训练样本，即由该参数确定的CRBM模型表达训练样本的概率分布，CRBM的训练方法采用对比散度算法对CRBM进行训练，以训练集中的样本作为模型的初始值。记训练集样本x＝{x0,x1,...,x
m
}，模型初始值记为v0，用训练样本初始化模型可视层单元的状态，即令v0＝x0，执行多步Gibbis采样。首先是上采样，根据公式计算隐层单元的激活概率，确定隐单元
的状态，即由p(h|v)求得h；然后由上采样得到的h，根据公式进行下采样，重构可视层单元的状态，即由p(v|h)求得v1。不断重复上述的上采样、下采样，如此经过k步采样，得到x0对应的采样结果vk。使用梯度下降法更新CRBM模型参数θ，因此基于CD算法的参数θ的梯度更新公式为：k取值为1，即采用CD
‑
1对CRBM进行训练
[0016]步骤6，使用训练好的模型实现对二维气体浓度的检测。优化后的CRBM有两个卷积层，分别具有8、16个卷积核，卷积核的大小为2
×
2，部长均为1；在第一个卷积层后连接一个下采样层，下采样层大小为2
×
2，步长为1；完全连接后输出浓度分布(M
×
N)
×
1的列，对其进行重排维数，可得到M
×
N的浓度矩阵。
[0017]本专利技术所述的基于卷积受限玻尔兹曼机和TDLAS的二维重建方法具有以下有益效果：
[0018]1.本专利技术所使用的CRBM模型对噪声不敏感，在噪声干扰下，依旧能够完整且准确地反演重建出二维浓度场。
[0019]2.本专利技术相比于原有的模拟退火算法，在完成网络训练后，能够迅速、准确的重建出火焰场的形貌结构，无需构建及改善目标优化函数。
[0020]3.本专利技术的训练集由实测数据经模拟退火算法重建获得，相较于本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于卷积受限玻尔兹曼机的二维重建方法，其特征在于，该方法包括以下具体步骤：步骤1：根据目标探测气体进行谱线遴选确定红外激光器中心波长；步骤2：获取不同投影角度的红外激光吸收信号并计算得到积分吸光度；步骤3：应用模拟退火(SA)算法对浓度场进行重建；步骤4：构建卷积受限玻尔兹曼机(CRBM)训练集；步骤5：确定CRBM模型结构；步骤6：使用训练好的模型实现对二维气体浓度的检测。2.根据权利要求1所述的基于卷积受限玻尔兹曼机和TDLAS的二维重建方法，其特征在于，步骤1针对一种待测气体，选取4条符合条件的吸收谱线，对应4个红外激光器，采用时分复用使得每一路激光都含有4种波长的吸收率信息。3.根据权利要求1所述的基于卷积受限玻尔兹曼机和TDLAS的二维重建方法，其特征在于，步骤2首先将待测区域离散成M
×
N个测量网络，并由M
×
N条正交布置的红外激光光路均匀穿过待测区域进行测量，得到每条光路的吸收光谱信号，采用Voigt线型计算积分吸光度4.根据权利要求1所述的基于卷积受限玻尔兹曼机和TDLAS的二维重建方法，其特征在于，步骤3模拟退火算法中测量得到的积分吸光度分布定义为P
m
(L
j
,λ
i
)，根据假设的气体浓度分布计算得到的值定义为P
c
(L
j
,λ
i
)，定义理想化的目标优化函数为式中i,j分别为波长和光束序号；加入平滑度用来改善目标病态问题的求解，定义平滑度为定义平滑系数为γ
T
，γ
P
，待优化目标函数可表述为F(T,P)＝DA(TX)+γ
T
R
T
(T)+γ
P
R
P
(P)。5.根据权利要求1所述的基于卷积受限玻尔兹曼机和TDLAS的二维重建方法，其特征在于，步骤4构建卷积受限玻尔兹曼机(CRBM)训练集时，步骤2计算得到的积分吸光度作为训练集的输入X，以模拟退火算法重建结果作为训练集的输出Y，共采用15000组样本进行网络训练，根据网...

【专利技术属性】
技术研发人员：郭奥，丰雷，常鹏浩，张李明，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：