一种三维电容层析成像的图像重建方法技术

技术编号:12976814 阅读:99 留言:0更新日期:2016-03-04 00:16
本发明专利技术公开了一种三维电容层析成像的图像重建方法,该方法是采用三维电容层析成像传感器获取的测量电容值作为图像重建的投影数据。在迭代重建过程中采用阈值滤波方式对所获得的重建图像进行伪迹抑制,滤波阈值定期自适应调整,滤波阈值的选取是使当前图像模糊度量最小。该方法有效地减少了重建图像伪迹,提高了复杂分布的三维成像效果,并实现了真正的三维电容层析成像图像重建。

【技术实现步骤摘要】

本专利技术涉及一种图像重建新方法,特别是涉及一种三维电容层析成像的图像重建新方法,属于电容层析成像

技术介绍
电容层析成像(electricalcapacitancetomography,简称ECT)技术是一种基于电容敏感原理的过程层析成像技术。它通过安置在管道或封闭容器周围一组特殊设计的电容传感器阵列,获取非导电被测物场在不同观测角度下的投影数据(电容值),采用适当的重建算法来反演出被测物场内介质(介电常数)分布,并以图像的形式给出介质分布结果。具有非侵入、响应快、低成本、安全无辐射等优点。在两相流和多相流检测中应用最为广泛,除此之外还被应用在头骨模型温度分布成像、仓储粮食水分监测、冻土层测量、滑动轴承润滑油膜测量等多个方面,涉及石油、化工、电力、冶金、建材、医学等多个国民经济和工业领域。因此该技术具有广阔的应用前景和发展潜力。目前针对ECT系统的测量和重建主要集中在二维重建上,即以图像的形式给出某一截面上的介电常数分布情况。但多数过程都发生在三维空间,且介质分布的几何形状、相对位置、体积等三维信息无法直接从二维图像中获取。因此直接的三维ECT成像已成为国内外学者的研究热点。但是与二维重建相比,三维重建图像像素数目是投影数据的几十倍、几百倍甚至上千倍,且投影数据信噪比更低,敏感场受被测介质分布的影响及软场误差更为严重。导致重建图像伪迹(即与真实介质分布无关的影像)严重,重建图像质量差。在二维重建过程中通常采用0-1滤波伪迹抑制方法(即将像素单元灰度值小于0的置0,灰度值大于1的置1),对重建图像中部分伪迹起到抑制作用,提高了重建图像质量。但是与二维电容层析成像相比,三维重建图像伪迹更为严重,在图像重建过程中对伪迹抑制能力的要求也更高。而0-1滤波伪迹抑制方法对三维ECT图像重建的伪迹抑制能力较弱,三维重建图像效果不佳。因此,具有较强伪迹抑制能力图像重建方法的研究具有重要的理论研究和实际应用价值。
技术实现思路
专利技术目的本专利技术提供一种三维电容层析成像的图像重建方法,克服现有方法在图像伪迹抑制上的不足,采用自适应的阈值滤波方法进行伪迹抑制。其目的是为了减少重建图像伪迹,提高对于复杂分布的重建能力。技术方案:本专利技术是通过以下技术方案来实现的:一种三维电容层析成像的图像重建方法,其特征在于:该方法步骤如下:步骤一:预先将J个极板分多层布置在被测区域周围,形成J极板电容层析成像传感器;将成像区域剖分为M个像素单元;测量各极板对间的电容值,获得N个电容值,N=J(J-1)/2,形成N维的测量电容向量C,N<<M;预先计算出图像重建所需的N×M维的归一化灵敏度矩阵S和M×N维的归一化灵敏度矩阵的转置阵ST,以及迭代步长α,α=2/λmax,λmax为STS的最大特征值,当α<1时,取α=1.5;定义滤波算子P(x),如式(1)所示,令η=0,由式(2)计算出图像初始值G0(G0=[G0(1),G0(2),…,G0(M)]);并令k=1;P(x)=0x≤ηxη<x≤11x>1---(1);]]>G0=P(STC)(2);步骤二:计算:Gk=Gk-1+αST(C-SGk-1)(3);步骤三:如果((k))q≠1,((k))q表示k对q取余数,q为预先设定好的一个正整数,则转到步骤四;否则按以下方法求η;将0~1的灰度值均匀地划分为0~L-1即L个灰度级,L为预先设定好的一个正整数,且为2的幂指数;令gmax、gmin分别代表图像的最大、最小灰度级,h(g)代表灰度级g在图像中出现的次数;对于给定阈值t,用式(4)~(8)定义阈值t所对应的图像模糊度量E(t)为:E(t)=1Mln2[Σg=gmintF(|g-g‾l|)h(g)+Σg=t+1gmaxF(|g-g‾h|)h(g)]---(4);]]>F(g)=-μ(g)ln[μ(g)]-[1-μ(g)]ln[1-μ(g)](5);μ(g)=11+g/(gmax-gmin),g=0,1,...,gmax-1---(6);]]>g‾l=int[Σg=0tgh(g)Σg=0th(g)]---(7);]]>g‾h=int[Σg=t+1L-1gh(g)Σg=t+1L-1h(g)]---(8);]]>其中,int[x]表示对实数x四舍五入取整数值;令E=1、t=gmin、topt=gmin;i)用式(4)~(8)计算阈值t所对应的图像模糊度量E(t),如E(t)<E,则令E=E(t),topt=t;ii)如果t=gmax-1则转到(iii);如果t<gmax-1,则令t=t+1转到(i);iii)令η=toptL-1;]]>步骤四:用式(9)对图像Gk进行强化伪迹抑制;Gk=P(Gk)(9);步骤五:判断是否满足迭代终止条件;如果满足终止条件则结束迭代,否则令k=k+1,转到步骤二。对每次迭代后的结果都进行一次基于阈值滤波的强化伪迹抑制,但每迭代q次才更新一次阈值,阈值的选取是使当前图像模糊度量最小。迭代终止条件是迭代次数达到指定次数。优点及效果:电容层析成像技术在两相流或多相流的参数检测中应用最多,同时在诸如头骨模型温度分布成像、仓储粮食水分监测、冻土层测量、滑动轴承润滑油膜测量等多个方面也有应用,涉及多个国民经济和工业领域。相对于传统的二维成像技术,三维电容层析成像技术难度更大。而目前电容层析成像技术所使用的重建方法中普遍采用0-1滤波伪迹抑制方法,该伪迹抑制方法对于三维电容层析成像重建来说,重建图像伪迹多,对复杂的介质分布的成像效果差。本专利技术提出一种三维电容层析成像的图像重建方法,该方法使用自适应阈值滤波方式进行伪迹抑制,并定期调整滤波阈值,可有效减少重建图像伪迹,具有方法简单、易于编程实现等优点。三维电容层析成像技术能够给出封闭容器内介质分布的三维图像,可获得介质分布的几何形状、相对位置、体积等信息,这是传统二维电容层析成像技术所获得的二维图像所不能比拟的。因此本专利技术提出的重建方法可以有效地抑制三维重建图像伪迹,提高重建图像质量,对于三维ECT技术在石油、化工、电力、冶金、建材等工业中的应用具有重要的实用价值,蕴涵良好的经济效益和应用前景。附图说明图1为12个极板的三维传感器结构示意图,本文档来自技高网
...
一种三维电容层析成像的图像重建方法

【技术保护点】
一种三维电容层析成像的图像重建方法,其特征在于:该方法步骤如下:步骤一:预先将J个极板分多层布置在被测区域周围,形成J极板电容层析成像传感器;将成像区域剖分为M个像素单元;测量各极板对间的电容值,获得N个电容值,N=J(J‑1)/2,形成N维的测量电容向量C,N<<M;预先计算出图像重建所需的N×M维的归一化灵敏度矩阵S和M×N维的归一化灵敏度矩阵的转置阵ST,以及迭代步长α,α=2/λmax,λmax为STS的最大特征值,当α<1时,取α=1.5;定义滤波算子P(x),如式(1)所示,令η=0,由式(2)计算出图像初始值G0(G0=[G0(1),G0(2),…,G0(M)]);并令k=1;P(x)=0x≤ηxη<x≤11x>1---(1);]]>G0=P(STC)          (2);步骤二:计算:Gk=Gk‑1+αST(C‑SGk‑1)           (3);步骤三:如果((k))q≠1,((k))q表示k对q取余数,q为预先设定好的一个正整数,则转到步骤四;否则按以下方法求η;将0~1的灰度值均匀地划分为0~L‑1即L个灰度级,L为预先设定好的一个正整数,且为2的幂指数;令gmax、gmin分别代表图像的最大、最小灰度级,h(g)代表灰度级g在图像中出现的次数;对于给定阈值t,用式(4)~(8)定义阈值t所对应的图像模糊度量E(t)为:E(t)=1M ln2[Σg=gmintF(|g-g‾l|)h(g)+Σg=t+1gmaxF(|g-g‾h|)h(g)]---(4);]]>F(g)=‑μ(g)ln[μ(g)]‑[1‑μ(g)]ln[1‑μ(g)]          (5);μ(g)=11+g/(gmax-gmin),g=0,1,...,gmax-1---(6);]]>g‾l=int[Σg=0tgh(g)Σg=0th(g)]---(7);]]>g‾h=int[Σg=t+1L-1gh(g)Σg=t+1L-1h(g)]---(8);]]>其中,int[x]表示对实数x四舍五入取整数值;令E=1、t=gmin、topt=gmin;i)用式(4)~(8)计算阈值t所对应的图像模糊度量E(t),如E(t)<E,则令E=E(t),topt=t;ii)如果t=gmax‑1则转到(iii);如果t<gmax‑1,则令t=t+1转到(i);iii)令η=toptL-1;]]>步骤四:用式(9)对图像Gk进行强化伪迹抑制;Gk=P(Gk)              (9);步骤五:判断是否满足迭代终止条件;如果满足终止条件则结束迭代,否则令k=k+1,转到步骤二。...

【技术特征摘要】
1.一种三维电容层析成像的图像重建方法,其特征在于:该方法步骤如下:
步骤一:预先将J个极板分多层布置在被测区域周围,形成J极板电容层析成像传感器;
将成像区域剖分为M个像素单元;测量各极板对间的电容值,获得N个电容值,N=J(J-1)/2,
形成N维的测量电容向量C,N<<M;预先计算出图像重建所需的N×M维的归一化灵敏度矩
阵S和M×N维的归一化灵敏度矩阵的转置阵ST,以及迭代步长α,α=2/λmax,λmax为STS
的最大特征值,当α<1时,取α=1.5;定义滤波算子P(x),如式(1)所示,令η=0,由式(2)
计算出图像初始值G0(G0=[G0(1),G0(2),…,G0(M)]);并令k=1;
P(x)=0x≤ηxη<x≤11x>1---(1);]]>G0=P(STC)(2);
步骤二:计算:
Gk=Gk-1+αST(C-SGk-1)(3);
步骤三:如果((k))q≠1,((k))q表示k对q取余数,q为预先设定好的一个正整数,则转
到步骤四;否则按以下方法求η;
将0~1的灰度值均匀地划分为0~L-1即L个灰度级,L为预先设定好的一个正整数,
且为2的幂指数;令gmax、gmin分别代表图像的最大、最小灰度级,h(g)代表灰度级g在图像
中出现的次数;对于给定阈值t,用式(4)~(8)定义阈值t所对应的图像模糊度量E(t)为:
E(t)=1Mln2[Σg=gmintF(|g-g‾l|)h(g)+Σg=t+1gmaxF(|...

【专利技术属性】
技术研发人员:王伊凡颜华周英钢孙延辉
申请(专利权)人:沈阳工业大学
类型:发明
国别省市:辽宁;21

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1