一种用于圆对称参数分布的双倍分辨率图像重建方法
【技术领域】
[0001] 本发明设及一种用于圆对称参数分布的双倍分辨率图像重建方法，具体设及在参 数满足圆对称分布特性时，能够有效简化系统结构并提高成像分辨率的单一视角投影图像 重建方法。
【背景技术】
[0002] 圆对称参数分布特性广泛存在于工业应用领域。例如，平焰燃烧炉可W产生均匀、 稳定的层流火焰，常用于研究燃烧、传热等过程。一般的实验条件下，由于冷边界层的存在 W及火焰与外界空气的热交换，从火焰中屯、到火焰外部温度逐渐降低，进而导致平焰燃烧 炉火焰区域内温度和燃烧产物气体浓度呈现圆对称分布特性，而非均匀分布。因此，为检 验平焰燃烧炉产生火焰的稳定性和均匀性，需要对平焰燃烧炉的火焰温度分布进行在线监 巧。。又如，航空喷气式发动机尾焰是锥状的高温、高速的燃气射流。研究表明，尾焰温度场 和气体组分浓度也呈现圆对称分布特性。为评估尾焰红外福射性能，优化推进剂配比方案， 需要对发动机尾焰温度分布进行测量。
[0003] 采用多视角投影可W对圆对称参数分布进行图像重建，但该样设计的结构较为复 杂。利用被测区域的圆对称参数分布特性，可W使用同屯、圆对其进行划分并建立方程，利用 Abel反变换求解方程，最终得到圆对称参数分布图像。区别与传统的单视角投影方法，本发 明通过优化探测阵列排布，在投影线数量不变的情况下，实现双倍分辨率重建参数分布图 像。
【发明内容】
[0004] 针对具有圆对称参数分布特性的被测区域，为进一步简化系统结构，本发明提供 一种用于圆对称参数分布的图像重建方法。
[0005] 本方法包括W下主要步骤。
[0006] 步骤一，获取单一视角的投影数据具有圆对称参数分布特性的被测区域中屯、 0为原点，建立一个平面直角坐标系xOy，X轴与y轴正交于0;两组平行于y轴的投影线分 别穿过被测区域，其中一组排列在y轴左侧，另一组排列在y轴右侧，每一组投影线中两条 相邻投影线的间距均为Ar;所述排列在y轴左侧的投影线与y轴的最小距离为Ar/2,所 述排列在y轴右侧的投影线与y轴的最小距离为Ar;所述两组投影线的总条数为N，将所 述两组投影线中的每一条投影线按照与y轴距离的大小，从小到大依次编号为1，2,…，N， 所述排列在y轴左侧的投影线编号为从1开始的连续奇数列，所述排列在y轴右侧的投影 线编号为从2开始的连续偶数列；在所述投影线的一端布置与所述投影线数量相同的接收 器，所述接收器编号与所对应的投影线编号相同；接收器的输出数据即为所对应投影线的 投影数据，所述投影数据的编号与所对应投影线的编号相同；
[0007] 步骤二，依据投影线穿过的区域和所对应的投影数据构建用于求解被测参数的方 程组；将被测区域围绕其中屯、划分为N个同屯、圆环，并假设每个圆环内被测参数取值相同， 若用变量i表示投影线和投影数据的编号，i= 1，2,…，N，则第i个投影数据可表示为被测 参数沿第i条投影线路径上的积分，即
[0008]
巧
[0009] 其中，A(Xi)为被测参数沿第i条投影线路径上的积分，a(rp为第j个圆环内的 被测参数值，f是积分变元，将i分别等于1，2,…，N时获得的方程依次排列，得到W被测参 数为未知数的方程组；
[0010] 步骤S，求解被测参数并重建图像；所述方程组可W通过Abel反变换求解，即
[0011]
巧
[001引其中，5为反变换积分变元，A'（Xj+S)为A(xp在Xj附近的立点Abel反卷积近 似，即
[0013]
， (3)
[0014] 根据公式（2)和公式（3)求得被测参数后，可重建被测参数的分布图像，所得图像 分辨率为Ar/2,该样，由两组间距均为Ar的投影线获得了分辨率为Ar/2的被测参数的 重建图像。
【附图说明】
[0015] 图1是基于可调谐激光吸收光谱技术的测量系统结构原理图。
[0016] 图2是激光器和探测器排列分布示意图。
[0017] 图3是激光穿过待测截面吸收投影示意图。
[0018] 图4是由仿真计算和实测重建得到的温度分布特性曲线。
[0019] 图5是由仿真计算和实测重建得到的水蒸气浓度分布特性曲线。
【具体实施方式】
[0020] 在本实施例中，W平焰燃烧炉产生具有圆对称温度和气体组分浓度分布的火焰参 数为例，利用基于可调谐激光吸收光谱技术的单一视角平行激光束投影对圆对称温度分布 进行图像重建。下面结合附图对本发明作进一步的说明：
[0021] 步骤一，选取谱线中屯、频率分别为7185. 6畑14和7444. 4cnTi的两条谱线进行实验 研究，该是两条水蒸气的特征谱线。选取炉盘直径为60mm的平面燃烧炉为实验载体，燃料 注入当量比为0. 69,此情况下燃烧过程稳定，选择直径为70mm的区域为实验研究对象。
[0022] 步骤二，两个可调谐半导体激光器1和2的中屯、波长分别为7185. 6畑14和 7444. 4畑14,激光控制器LDC-3900通过两个独立的通道对两个激光器进行温度和电流控 审IJ，进而实现对两个激光器输出激光频率的调谐。
[0023] 步骤=，如图1所示，两个激光器输出的两路激光首先通过2X4的光纤禪合-分 路器，其输出的一路激光由准直器准直并经固体标准具后被一个光电探测器接收并转变为 电信号，电信号输入到数据采集器，根据法布里波罗干设原理，该路信号岁激光频率的线性 变化产生一个周期震荡，因此可用于标定激光频率变化；另=路激光分别经过=个1X4的 光纤分路器，由光纤准直器准直，平行穿过被测区域。如图2所示，建立一个原点位于待测 区域中屯、的平面直角坐标系，在y轴左右两侧，各排列6路激光，左右两侧激光束间距均为 Ar= 5.6mm。其中，左侧1号光束y轴距离为Ar/2,右侧2号光束与y轴距离为Ar。在 投影区域布置12个探测器，接收W上12路激光的光强信号，经过转换电路作为最终的投影 值。该里需要强调的是，由于准直器、探测器等体积限制，无法将12路激光束布置在y轴同 侧。
[0024] 步骤四，如图3所示，被测区域是平面燃烧炉火焰的一个横截面，将其围绕其中屯、 划分为半径差等于Ar/2同屯、圆环，区域中屯、为一个半径为Ar/2的圆。该截面与炉盘距 离为40mm，若在原有平面直角坐标系的基础上建立一个空间直角坐标系，那么Z= 40mm。 假设划分后的每个区域内各类分布参数为常值，它们可W表达为W中屯、距r为自变量的函 数。该里，压强记为P(r)(本例平面燃烧炉燃烧过程在开放环境下进行，P与环境气压相等， 为常值latm)，温度记为T(r)，水蒸气浓度记为X(r)。根据吸收光谱相关理论，激光经过待 测区域的吸收密度函数可表示为：
[0025] a(r) =P(r) ?S[T(r) ] ?X(r), (4)
[0026] 其中，S(T)是线强度函数，可化围过光谱数据库查询到。
[0027] 步骤五，根据上面叙述，探测器接受的投影值即一束激光的总吸收强度，因此，第i 个投影值可表示为第i条投影线路径上W吸收密度为权重的积分：
[0028]
巧
[0029] 其中，A(Xi)为投影值，a(rp是第j个圆环上的参数值，r是求解半径。
[0030] 步骤六，通过Abel反变换求解吸收密度，
[0031]
货）
[003引该里，A(Xi)可W在Xi附近通过二次多项式近似，
[0033]
。 （7)
[0034] 步骤走，根据双线波长扫描测温方法中计算温度的方法，第i个圆环内的温度可 由下式计算：
[00对
做
[0036] 而浓度分布则可W在求解温度后由下式计算：
[0037]
(9)
[003引将计算得到的分布参数用于重建，所得温度和水蒸气浓度分布图像分辨率为 Ar/2,该样，由两组间距均为Ar的投影线获得了分辨率为Ar/2的被测参数的重建图像。
[0039] 步骤走，使用与前述相同的几何参数及物理参数进行计算流体力学仿真，可W得 到理想状态下温度和水蒸气浓度的分布数据。
[0040] 通过此方法重建的圆对称温度场分布与仿真得到的温度场分布对比如图4所示， 通过此方法重建的圆对称水蒸气浓度分布与仿真得到的水蒸气浓度分布对比如图5所示。 由图可知，不论温度分布还是水蒸气浓度，重建的分布参数曲线和仿真得到的分布参数曲 线是吻合的，证明该方法是有效的。
[0041] W上对本发明及其实施方式的描述，并不局限于此，附图中所示仅是本发明的实 施方式之一。在不脱离本发明创造宗旨的情况下，不经创造地设计出与该技术方案类似的 结构或实施例，均属本发明保护范围。
【主权项】
1. 一种用于圆对称参数分布的双倍分辨率图像重建方法，包括以下步骤： 步骤一，获取单一视角的投影数据：以具有圆对称参数分布特性的被测区域中心O为 原点，建立一个平面直角坐标系x〇y，X轴与y轴正交于O ;两组平行于y轴的投影线分别穿 过被测区域，其中一组排列在y轴左侧，另一组排列在y轴右侧，每一组投影线中两条相邻 投影线的间距均为ΔΓ;所述排列在y轴左侧的投影线与y轴的最小距离为Λ Γ/2,所述排 列在y轴右侧的投影线与y轴的最小距离为ΔΓ ;所述两组投影线的总条数为Ν，将所述两 组投影线中的每一条投影线按照与y轴距离的大小，从小到大依次编号为1，2,…，Ν，所述 排列在y轴左侧的投影线编号为从1开始的连续奇数列，所述排列在y轴右侧的投影线编 号为从2开始的连续偶数列；在所述投影线的一端布置与所述投影线数量相同的接收器， 所述接收器编号与所对应的投影线编号相同；接收器的输出数据即为所对应投影线的投影 数据，所述投影数据的编号与所对应投影线的编号相同； 步骤二，依据投影线穿过的区域和所对应的投影数据构建用于求解被测参数的方程 组：将被测区域围绕其中心划分为N个同心圆环，并假设每个圆环内被测参数取值相同，若 用变量i表示投影线和投影数据的编号，i = 1，2,…，N，则第i个投影数据可表示为被测参 数沿第i条投影线路径上的积分，即 其中，A(Xi)为被测参数沿第i条投影线路径上的积分，a(rp为第j个圆环内的被测 参数值，F是积分变元，将i分别等于1，2,…，N时获得的方程依次排列，得到以被测参数为 未知数的方程组； 步骤三，求解被测参数并重建图像：所述方程组可以通过Abel反变换求解，即其中，S为反变换积分变元，Α'（Χ」+δ)为A(X j)在\附近的三点Abel反卷积近似，即恨畑U; 氺侍攸测麥数厄，円里姓攸测麥数tfr刀、弗囹诼，所得图像分辨 率为ΛΓ/2,这样，由两组间距均为Ar的投影线获得了分辨率为ΛΓ/2的被测参数的重建 图像。
【专利摘要】本发明提供一种用于圆对称参数分布的双倍分辨率图像重建方法，包括以下步骤：利用覆盖被测参数分布场的单一视角投影，得到对应的投影数据；根据投影关系和投影数据建立用于求解被测参数的方程组；利用Abel反变换求解分被测参数，进而得到圆对称参数双倍分辨率图像。本发明的效果是系统结构简单，可靠性高，成本低，具有广泛的应用前景。
【IPC分类】G06T11/00
【公开号】CN104881887
【申请号】CN201510329261
【发明人】徐立军, 曹章, 曲前伟, 刘畅
【申请人】北京航空航天大学
【公开日】2015年9月2日
【申请日】2015年6月15日