一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法技术

一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法，用于对燃烧火焰剖面重建。该方法基于化学自发光层析成像原理（CTC），利用燃烧火焰多角度同步成像，首先得到火焰剖面上的凸轮廓，进一步结合代数迭代算法（ART）重建火焰剖面。该方法有效避免了在成像角度较少的情况下、采用常规CTC方法对火焰剖面重建容易出现的“星形拖影”；特别当剖面上火焰强度较弱时，拖影现象会对严重影响重建效果。基于该方法对火焰剖面的重建在避免拖影的同时可加快重建进程的结束，减少运算量，且无需增加其他硬件成本。

A Combustion Flame CTC Tomography Method Based on Profile Cam Profile Recognition and ART Algorithms

【技术实现步骤摘要】
一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法
本专利技术涉及一种燃烧火焰层析成像方法，尤其涉及一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法。
技术介绍
近年来化学自发光层析成像技术(CTC)逐渐发展起来并应用于燃烧火焰的可视化检测。CTC技术利用燃烧过程中产生的化学发光，基于辐射传递模型和计算机断层扫描技术经过重建得到燃烧火焰三维分布信息，实现对燃烧过程的三维可视化诊断，为燃料燃烧机理研究提供有效手段。基于CTC技术对燃烧火焰进行三维重建（层析成像）由于算法自身原因，火焰剖面重建易出现“星形拖影”，即剖面上焰体外部区域出现几何亮度拖影，与实际情况不符；检测系统成像角度越少、火焰越微弱拖影造成的干扰越严重；而实际系统由于受空间限制无法实现密集角度成像，因此由于拖影干扰对基于CTC检测原理的燃烧过程可视化诊断精度造成了显著影响。
技术实现思路
本专利技术针对CTC化学自发光层析成像技术中存在的“亮度拖影”干扰问题，公开基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法；基于该方法可对燃烧火焰剖面进行精确重建，提高火焰三维可视化诊断水平；具体方法步骤如下：步骤一、燃烧火焰多角度同步成像与预处理搭建多角度同步成像检测系统采集系统，如图1所示，由系统记录同一时刻燃烧火焰在不同成像角度下的数字图像；每个角度成像由一个成像光路实现，如图2所示，不同成像光路主光轴交于同一点O，各成像光路像平面主轴线重合；设火焰数字图像分辨率(像素元个数)为M×N，则M为图像(像元)矩阵行数，N为图像矩阵列数；因此火焰在不同成像角度下的数字图像由M个像元行向量{P1，P2，P3，···，PM}构成。将系统采集得到的各角度火焰数字图像二值化，以区分图像火焰区与背景区；并根据二值化结果，将各角度火焰数字图像原像中对应背景区的像素元赋值为零，以减少杂散光干扰。二值化的火焰图像所包含像元行向量记为；步骤二、火焰剖面凸轮廓识别如图3所示，垂直于火焰成像主轴线取火焰剖面12，由剖面在火焰主轴线上的相对位置，从各角度火焰二值化图像提取行向量11，行向量11在不同角度火焰图像中均指同一行像元向量，由火焰数字图像分辨率为M×N知，所提取各行向量包含N个像元；故所提取向量可表示为：（1≤i≤R,1≤m≤M），其中i代表成像角度编号，R为系统成像角度总数量，m为行序；同样地，从各角度火焰图像（原图非二值化）提取同一行向量Pim（1≤i≤R,1≤m≤M）；为每个成像角度定义旋转矩阵Ti（1≤i≤R）并赋初值为零，矩阵规模为M×N（M行N列）；依次读取向量各元素值（1≤i≤N），若Pim(j)=1，将该成像角度旋转矩阵Ti第j列元素重新赋值为1，其他元素值不变。按照该方法完成对T1~TR的重新赋值，由重新赋值后的旋转矩阵T1~TR相乘得，则T′即为火焰剖面二值化凸轮廓图像。内部元素值为1的连通区域代表火焰剖面焰体“明亮”区域，元素值为0的部分代表剖面背景区域，“明亮”区域即为火焰剖面凸轮廓；步骤三、基于火焰剖面凸轮廓识别的剖面重建火焰剖面重建采用ART算法结合剖面凸轮廓的识别完成；重建过程为循环迭代运算，具体迭代运算流程如下：定义N×N矩阵R并赋初值为零，由初值矩阵R表示火焰剖面重建图像像元矩阵;依据燃烧火焰在不同角度的成像对初值矩阵R进行“反投影”矫正得到一次矫正矩阵R1,并基于步骤二中对火焰剖面凸轮廓的识别对R1进行进一步修正:，至此完成第一次迭代运算，并将作为下次迭代运算的初值矩阵；启动循环迭代运算，由最后一次迭代得到N次修正矩阵，基于矩阵还原得到火焰剖面图像。附图说明图1为根据本
技术实现思路
提供的多角度同步成像检测系统；图2为根据本
技术实现思路
提供的多角度同步成像检测系统中成像光路像平面及主光轴几何关系的说明；图3为根据本
技术实现思路
提供的火焰剖面与火焰图像行向量间的映射关系示意图；图4为根据本专利技术实施实例提供的具有三个成像光路的多角度同步成像检测系统；图5为根据本专利技术实施实例提供的本生灯层流火焰三个成像角度的数字图像；图6为根据本专利技术实施实例提供的本生灯层流火焰三个成像角度的二值化图像；图7为根据本专利技术实施实例提供的，用于火焰剖面凸轮廓识别的旋转矩阵的构造以及所得火焰剖面凸轮廓；图8为根据本专利技术实施实例提供的，中采用基于本专利技术方法，利用三个成像角度对火焰根部剖面的重建，以及采用传统CTC检测方法对火焰根部剖面的重建对比；附图标记：11—火焰图像行向量；12—火焰剖面；61—图像背景区；62—图像火焰区；63—火焰根部初燃区；71—0°旋转矩阵；72—120°旋转矩阵；73—240°旋转矩阵；74—火焰剖面图轮廓；81—采用本专利技术方法重建所得火焰根部剖面；82—采用传统CTC检测方法重建所得火焰根部剖面。具体实施方式现以本生灯燃烧器产生的层流火焰剖面重建为例，结合附图对本专利技术所公开的一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法具体实施过程做进一步说明；图4为本实施实例搭建的具有三个成像光路的多角度同步成像检测系统，分别从0°、120°、240°三个成像角度同步获取燃烧火焰数字图像，如图5所示；系统利用同步盒发出的脉冲信号作为外部信号源触发三个成像光路中的CCD同步工作，触发信号以脉冲上升沿有效；将各角度火焰图像二值化处理，如图6所示，通过二值化可区分各角度火焰图像中的着火区和背景区，黑色部分61为背景区，白色部分62为着火区；选取火焰根部剖面63进行重建。根据
技术实现思路
步骤二，基于火焰在三个成像角度的二值化图像提取向量并构造旋转矩阵，并基于旋转矩阵得到火焰剖面凸轮廓，如图7所示；图中71、72、73分别为基于本实施实例中0°、120°、240°三个成像角度构造的旋转矩阵，74为基于旋转矩阵得到的火焰根部剖面凸轮廓；根据
技术实现思路
步骤三，基于火焰剖面凸轮廓，采用ART算法对火焰根部初燃区剖面进行重建，迭代次数设置为100；图8中81为采用本专利技术方法重建所得火焰根部剖面，重建结果具有环带亮度分布的特点，符合层流火焰初燃区理论结构，即火焰初燃区外部与空气接触部分为燃烧带，内部不与空气直接接触部分为未燃区，整体呈现环状亮度分布特点。尽管火焰根部亮度较微弱，在仅有三个成像角度的前提下采用本专利技术方法仍基本重建出火焰根部剖面的真实样貌；为增加对比，82给出未采用本专利技术方法，而采用传统CTC检测方法对火焰根部剖面的重建结果，可看出采用传统CTC检测方法对火焰根部剖面的重建具有较显著的“拖影”干扰，重建效果较差。本文档来自技高网...

【技术保护点】
1.一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法，用于燃烧火焰剖面图像重建；该方法通过采集燃烧火焰在多个同步成像角度下的数字图像，基于化学自发光层析成像原理（CTC），实现火焰剖面图像重建，具体步骤为：步骤一、燃烧火焰多角度同步成像采集与预处理；利用同步成像检测系统记录燃烧火焰同一时刻在R个成像角度下的数字图像；系统中每个角度的成像由一个成像光路实现，不同成像光路主光轴交于同一点O，且各光路像平面主轴线（Z轴）重合；设火焰数字图像分辨率(像素元个数)为M×N，则M为图像(像元)矩阵行数，N为图像矩阵列数；因此火焰在不同成像角度下的数字图像由M个像元行向量{

【技术特征摘要】
1.一种基于剖面凸轮廓识别及ART算法的燃烧火焰CTC层析成像方法，用于燃烧火焰剖面图像重建；该方法通过采集燃烧火焰在多个同步成像角度下的数字图像，基于化学自发光层析成像原理（CTC），实现火焰剖面图像重建，具体步骤为：步骤一、燃烧火焰多角度同步成像采集与预处理；利用同步成像检测系统记录燃烧火焰同一时刻在R个成像角度下的数字图像；系统中每个角度的成像由一个成像光路实现，不同成像光路主光轴交于同一点O，且各光路像平面主轴线（Z轴）重合；设火焰数字图像分辨率(像素元个数)为M×N，则M为图像(像元)矩阵行数，N为图像矩阵列数；因此火焰在不同成像角度下的数字图像由M个像元行向量{P1，P2，P3，···，PM}构成；将系统采集得到的各角度火焰数字图像二值化，以区分图像火焰区与背景区；并根据二值化结果，将各角度火焰数字图像原像中对应背景区的像素元赋值为零，以减少杂散光干扰；二值化的火焰图像所包含像元行向量记为；步骤二、火焰剖面凸轮廓识别；垂直于火焰成像主轴线（Z轴）取火焰剖面，根据CTC成像原理，由剖面在火焰主轴线上的相对位置，从各角度火焰二值化图像提取同一行向量（1≤i≤R，1≤m≤M），其中i代表成像角度编号，R为系统成像角度总数量，m为行序；由火焰数字图像分辨率为M×N知，所提取各行向量包含N个像元；同样地，从各角度火焰图像（原图非二值化）提取同一行向量Pim（1≤i≤R，1≤m≤M）；为每个成像角度定义旋转矩阵Ti（1≤i≤R）并赋初值为零，矩阵规模为M×N（M行N列）；依次读取向量各元素值(j)（1≤j≤N），若(j)＝1（1≤j≤N），则将旋转矩阵Ti第j列元素重新赋值为1，其他元素值不变；按照该...

【专利技术属性】
技术研发人员：成艳亭，王建军，孙砚飞，
申请(专利权)人：山东理工大学，
类型：发明
国别省市：山东,37