CN110934586B - 一种灰度值矩阵快速分解重建的正则化方法 - Google Patents
一种灰度值矩阵快速分解重建的正则化方法 Download PDFInfo
- Publication number
- CN110934586B CN110934586B CN201911302022.2A CN201911302022A CN110934586B CN 110934586 B CN110934586 B CN 110934586B CN 201911302022 A CN201911302022 A CN 201911302022A CN 110934586 B CN110934586 B CN 110934586B
- Authority
- CN
- China
- Prior art keywords
- gray value
- updating
- reconstruction
- steps
- variable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 239000011159 matrix material Substances 0.000 title claims abstract description 32
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 18
- 238000005259 measurement Methods 0.000 claims abstract description 28
- 238000003384 imaging method Methods 0.000 claims abstract description 11
- 230000035945 sensitivity Effects 0.000 claims abstract description 11
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 230000009191 jumping Effects 0.000 claims abstract description 3
- 238000003325 tomography Methods 0.000 claims description 19
- 238000009826 distribution Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- 230000005284 excitation Effects 0.000 claims description 8
- 230000009977 dual effect Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000002593 electrical impedance tomography Methods 0.000 claims description 5
- 125000004122 cyclic group Chemical group 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 229910052802 copper Inorganic materials 0.000 claims description 2
- 239000010949 copper Substances 0.000 claims description 2
- 230000036039 immunity Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012804 iterative process Methods 0.000 description 2
- 238000004321 preservation Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/053—Measuring electrical impedance or conductance of a portion of the body
- A61B5/0536—Impedance imaging, e.g. by tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种灰度值矩阵快速分解重建的正则化方法,包含如下步骤:(1)根据被测场域获取重建所需的边界测量电压b和灵敏度矩阵A;(2)设置初始化参数;(3)更新辅助变量pk+1;(4)更新辅助变量qk+1;(5)更新变量wk+1;(6)更新变量vk+1;(7)更新灰度值gk+1;(8)设置优化参数;(9)更新灰度值uk+1;(10)更新变量κk+1;(11)更新变量mk+1;(12)判断迭代是否符合迭代终止条件uk+1是否收敛,若是则迭代终止,进行下一步操作;若否,设置k=k+1并跳回第(9)步,继续迭代求解;(13)根据最终求解所得灰度值进行成像。本发明能够有效地降低阶梯伪影、提升背景清晰及提高抗噪性。
Description
技术领域
本发明属于电阻抗层析成像图像重建方法技术领域,具体涉及一种灰度值矩阵快速分解重建的正则化方法。
背景技术
电学层析成像(Electrical Tomography,ET)是层析成像技术的一种,通过对被测物体施加电流激励获得测量电压数据值,利用特定数学手段逆推被测物体内部的电特性参数分布(电导率/介电系数/复导纳/磁导率),从而实现场域内部可视化的过程。与磁共振(MRI),CT等成像技术相比,电学层析成像具有非侵入、无辐射、实时性等优势。例如电阻抗层析成像(EIT)技术,目前已广泛应用于工业过程监测、地球物理地下探测和生物组织诊断等领域。值得注意的是,在医学成像领域上,EIT已被用于乳腺肿瘤的检查、肺通气和脑部成像。
从数学上看,电学层析成像逆问题(即图像重建过程)求解过程是一个高度的非线性过程。一般是通过线性化处理将非线性问题转化为线性逆问题进行求解。由于获取的场域边界电压数量远小于求解场域的像素值,会造成求解逆问题时的不适定性,所以一般选取正则化方法处理,例如Tikhonov正则化方法。这种方法能够很好地恢复边缘光滑的物体,而当恢复边缘锐利的物体时则会出现过度的平滑。为了减弱Tikhonov正则化方法带来的过度光滑性,研究了全变分(Total variation regularization method,TV)正则化方法重建尖锐不连续物体的方法。然而,TV正则化方法在重建图像过程中会产生分段不变的图像,称为“阶梯伪影”。为了抑制“阶梯伪影”,对TV正则化方法的改进进行了大量的工作,提出了总广义变分(Total generalized variation,TGV),在一定程度上缓解了阶梯伪影。
在现有的研究中,这些正则化算法都局限在对噪声敏感的L2范数保真度项和在一定程度上产生阶梯伪影的惩罚项。这使得成像结果容易受到噪声的影响,并且无法很好的去除阶梯伪影,导致图像分辨率低。因此,上述几种变分方法无法取得很好的成像效果,使得很难在实际运用中推广使用。
发明内容
本发明的目的在于提出一种灰度值矩阵快速分解重建的正则化方法,该方法可以有效地降低阶梯伪影、提升背景清晰及提高抗噪性。
本发明为实现上述目的采用如下技术方案,一种灰度值矩阵快速分解重建的正则化方法,其特征在于具体过程为:将电学层析成像看作一个线性不适定问题Ag=b,其中,A为灵敏度矩阵,b为边界测量电压差值,g为灰度值;设计最小化的目标函数为
为了解决上述式子的L1最小化问题,在原始对偶算法的基础上,在最小化目标函数中引入了辅助变量I:p和辅助变量II:q,将最小化目标函数的对偶形式根据Legendre-Fenchel对偶形式导出
由此可推导出
式中σ,τ代表着参数;
为了进一步提高图像空间分辨率,抑制阶梯伪影,针对步骤五迭代的灰度值,本发明提出了一种全新的惩罚项用于图像重建过程。设计了灰度值优化迭代过程:对于优化后灰度值u来说,因为YTY=I,所以可以得到u=YTYu,Y表示灰度值快速分解系数矩阵,YT是Y的逆变换,YT表示灰度值快速重建系数矩阵。可以利用惩罚项具有沉余性和多分辨率结构的特点,结合L0范数,对灰度值进行快速分解重建,这样可以稀疏的逼近图像的分段光滑函数,降低阶梯伪影,提高图像空间分辨率。其目标函数定义为:式中,/>是保真项,C是可以自适应选择图像不同区域的微分算子,g是灰度值,||λ·κ||0是惩罚项,λ是正则化参数,κ是优化系数。
为了解决上述L0范数最小化问题,我们采用拉格朗日算法求解最小化问题。上述式子可转化为
式中m是辅助变量,γ,μ是参数。
通过迭代可得到
由于方程的每一个式子都有一个封闭的形式解,所以L0极小化问题可以改写为
式中,算子H是一个硬阈值函数。
最终,根据求解优化后灰度值u完成成像。
本发明与现有技术相比具有以下有益效果:本发明的灰度值矩阵快速分解重建的正则化方法在传统TGV正则化的基础上,通过再次迭代灰度值来提高背景清晰度和抑制产生的伪影。另外,为了解决所提出的L0极小化问题,采用拉格朗日算法进行求解。相比于TV正则化与TGV正则化,本发明提出的灰度值矩阵快速分解重建的正则化方法在降低阶梯伪影和提升图像清晰度等方面有着很好的效果,同时自适应参数的选择增强了参数选择的客观性和简易性。
附图说明
图1为本发明的灰度值矩阵快速分解重建的正则化方法的流程框图。
图2为本发明的电阻层析成像系统环形单截面被测场域,激励电流和测量电压的模式以及电极分布。
图3为本发明的实施例在选取三种模型的真实分布时,TV正则化算法,TGV正则化算法和伪影抑制图像重建方法的图像重建结果示意图。
图4为三种方法在不同噪声水平下对模型2的重建图像。
图5为三种方法在不同噪声水平下模型2的图像相对误差和相关系数。
图中:1-被测场域,2-电极,3-激励电流,4-测量电压。
具体实施方式
结合附图和实施例对本发明的一种灰度值矩阵快速分解重建的正则化方法加以说明。
本发明提出的一种灰度值矩阵快速分解重建的正则化方法,针对传统正则化算法产生的阶梯伪影问题和背景不清晰问题,以TGV正则化为基础,结合提出的灰度值优化方法,采用L2范数作为保真度项,L0范数作为惩罚项,通过自适应选择正则化参数和权重因子来选取最优值,以一种高效、快速的拉格朗日算法来求解本发明提出的目标函数,完成最终的逆问题求解。
如图1所示,为本发明的一种灰度值矩阵快速分解重建的正则化方法流程图。
如图2所示,为电学层析成像之一的生物电阻抗层析成像系统环形单截面被测场域1,激励电流3和测量电压4的模式以及电极2分布,采用16电极2均匀分布在被测场域1外壁。
选取三种典型的介质模型为实施例,场域内物体真实分布如图3左侧一竖列所示,其他三列从左到右分别表示为TV正则化算法,TGV正则化算法和本发明提出的正则化方法。为了较好地体现本发明中的方法与其它两种方法的不同,分别给出了三种典型模型在这三种正则化算法下的求解结果。
实施例包括如下具体步骤:一种灰度值矩阵快速分解重建的正则化方法,该方法将电学层析成像看作一个线性不适定问题Ag=b,其中,A为灵敏度矩阵,b为边界测量电压差值,g代表着灰度值,待求解的逆问题可以用最小二乘优化形式的目标函数表示:其中,f(g)为目标函数。正则化方法为提高解的稳定性提供了一种新的选择,其形式可以描述为:/>其中,λ是正则化参数。
TV正则化算法在保持锐利边缘方面表现出良好的效果。其一般形式可以表示为:其中,/>表示TV正则化算法中的惩罚项。在边缘保持方面,TV正则化方法有着明显的优势。然而TV正则化分段逼近,在图像重建中可以观察到明显的阶梯伪影。
为了克服TV正则化在图像重建过程中产生的阶梯伪影问题,提高重建图像的质量,提出了TGV正则化,其一般形式可以表示为:其中/>表示TGV正则化算法中的惩罚项。在一定程度上抑制了阶梯伪影,但是在边缘保持方面表现不好。
本发明针对传统正则化算法产生的阶梯伪影和边缘过度光滑问题,以TGV正则化为基础,结合本发明提出的灰度值优化方法,通过自适应选择正则化参数和权重因子来选取最优值,以一种高效、快速的拉格朗日算法求解本发明提出的目标函数。
步骤一:针对三种典型模型,分别获取各自重建所需的灵敏度矩阵A和测量电压b:本发明是将被测对象置于电学层析成像测量系统中,采用绝对成像的方式,被测场域外边界处均匀分布16个铜电极,采用经典的电流激励电压测量方式,持续采集循环激励循环测量下的边界电压,共获得208个测量值,104个独立测量值;b为不含夹杂物的空场边界电压bmeas1和含有夹杂物的有物场的边界测量电压bmeas2之差(即相对边界测量值b=bmeas2-bmeas1)。
灵敏度矩阵是根据不含内含物的空场的边界测量电压,结合灵敏度理论,计算灵敏度矩阵,计算公式为:其中,Aij是第j个电极对对第i个电极对的灵敏度系数,/>分别为第i个电极对及第j个电极对在激励电流为Ii,Ij时场域电势分布。
步骤二:将非线性电阻层析成像的逆问题转化为线性问题;边界测量电压与电导率分布的关系是非线性的f(σ)=b,式中σ是电导率。对于变化不大的电导分布,f(σ)=b可将边界测量电压的变化简化为线性形式式中Δσ是电导率分布的扰动,Δb是电导率的变化引起的边界测量电压的变化。最后导出电阻层析成像逆问题的线性化形式Ag=b,式中g为所求成像的灰度值。
步骤三:设计目标函数为:式中,/>是保真项,λTGV是正则化参数,A代表着雅可比矩阵,g代表着灰度值,/>是惩罚项,α1,α0是自适应权重因子,/>为梯度算子,/>表示对称梯度算子,v是对称场量。在电阻层析成像图像重建中,通过将目标函数极小化来求解最优灰度值,其模型式表示为:/>
步骤五:根据步骤四推导的对偶形式的算法如下
(5)更新变量vk+1:vk+1=vk+τ(pk+1+divεqk+1)。
步骤六:为了进一步提高图像空间分辨率,抑制阶梯伪影,针对步骤五迭代的灰度值,本发明提出了一种全新的惩罚项用于图像重建过程。设计了灰度值优化迭代过程:对于优化后灰度值u来说,因为YTY=I,所以可以得到u=YTYu,Y表示灰度值快速分解系数矩阵,YT是Y的逆变换,YT表示灰度值快速重建系数矩阵。可以利用惩罚项具有沉余性和多分辨率结构的特点,结合L0范数,对灰度值进行快速分解重建,这样可以稀疏的逼近图像的分段光滑函数,降低阶梯伪影,提高图像空间分辨率。其目标函数定义为:式中,/>是保真项,C是可以自适应选择图像不同区域的微分算子,g是灰度值,||λ·κ||0是惩罚项,λ是正则化参数,κ是优化系数。
步骤七:利用拉格朗日算法推导出目标函数的模型式
式中,m是辅助变量,μ,γ是设置的参数。
步骤八:根据步骤七推导的模型式的算法如下
(2)更新优化后灰度值uk+1:uk+1=(CTC+(μ+γ)I)-1(CTg+γuk+μYT(κk-mk))。
(3)更新变量κk+1:κk+1=Hλ,μ,γ(Yuk+1+mk,κk)。
(4)更新变量mk+1:mk+1=mk+Yuk+1-κk+1。
(5)判断迭代是否符合迭代终止条件uk+1是否收敛,若是则迭代终止,进行下一步操作;若否,设置k=k+1并跳回步骤八第(2)步,继续迭代求解。
步骤九:根据最终求解优化后灰度值u,进行成像。
图3为本发明的实施例选取三种典型的介质模型,TV算法,TGV算法和本发明的方法的图像重建结果示意图。可以看出,三种典型模型中,采用TV正则化方法时,背景中存在明显的阶梯伪影,严重影响了图像恢复的质量。与TV相比,TGV有效地降低了阶梯伪影。然而在边缘保持方面表现不好。相比于前面两种正则化算法,本发明提供的灰度值矩阵快速分解重建的正则化方法很好的抑制了阶梯,并提升了图像背景的清晰程度,同时有着更高的图像空间分辨率。在实际应用中,电学层析成像系统容易受到噪声的影响。为了验证本发明提出的一种灰度值矩阵快速分解重建的正则化方法的抗噪声性能,图4比较了三种方法在不同噪声水平下对模型B的重建图像。可以看出,当噪声水平逐渐增大时,本发明提出的方法在抑制阶梯伪影优于其它两种方法。同时,在电学层析成像中,通常采用图像相对误差(Relative Error,RE)和相关系数(Correlation Coefficient,CC)评价算法来定量图像重建质量,表达式如①、②所示,图像相对误差越小,相关系数越大,表明图像重建质量越好。
图5给出了这三种方法不同噪声水平下模型B的图像相对误差和相关系数,可以看出,本发明提出的一种灰度值矩阵快速分解重建的正则化方法与TV算法,TGV算法相比,具有最低的相对误差和最高的相关系数,能够准确的重建出被测场域内分布,明显地提高了电学层析成像逆问题求解精度和抗噪声能力。
以上所述仅为本发明的较佳实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种灰度值矩阵快速分解重建的正则化方法,其特征在于具体过程为:
步骤一:针对三种典型模型,分别获取各自重建所需的灵敏度矩阵A和测量电压b:将被测对象置于电学层析成像测量系统中,采用绝对成像的方式,被测场域外边界处均匀分布16个铜电极,采用经典的电流激励电压测量方式,持续采集循环激励循环测量下的边界电压,共获得208个测量值,104个独立测量值;b为不含夹杂物的空场边界电压bmeas1和含有夹杂物的有物场的边界测量电压bmeas2之差即相对边界测量值b=bmeas2-bmeas1;
灵敏度矩阵是根据空场状态下的边界测量电压数据得到的,结合灵敏度理论通过公式计算得到灵敏度矩阵,计算公式为:其中,Aij是第j个电极对对第i个电极对的灵敏度系数,/>分别为第i个电极对及第j个电极对在激励电流为Ii,Ij时场域电势分布;
步骤二:电阻抗层析成像的逆问题是一个高度的非线性问题,需要转化为线性问题进行求解:边界测量电压与电导率分布的关系是非线性的f(σ)=b,式中σ是电导率,对于变化不大的电导分布,f(σ)=b可将边界测量电压的变化简化为线性形式式中Δσ是电导率分布的扰动,Δb是电导率的变化引起的边界测量电压的变化,最后导出电阻层析成像逆问题的线性化形式Ag=b,式中g为所求成像的灰度值;
步骤三:设计目标函数为:式中,/>是保真项,λTGV是正则化参数,A代表着雅可比矩阵,g代表着灰度值,b是测量电压,是惩罚项,α1,α0是自适应权重因子,▽为梯度算子,ε(v)=(▽v+▽vT)/2表示对称梯度算子,v是对称场量,在电阻抗层析成像图像重建中,通过将目标函数极小化来求解最优灰度值,其模型式表示为:/>
步骤五:根据步骤四推导的对偶形式的迭代算法如下:
(5)更新变量vk+1:vk+1=vk+τ(pk+1+divεqk+1);
(6)更新灰度值gk+1:gk+1=gk+τ(div▽pk+1-ATwk+1);
步骤六:为了进一步提高图像空间分辨率,抑制阶梯伪影,针对步骤五迭代的灰度值,设计了灰度值优化迭代过程:对于优化后灰度值u来说,因为YTY=I,所以得到u=YTYu,Y表示灰度值快速分解系数矩阵,YT是Y的逆变换,YT表示灰度值快速重建系数矩阵;利用惩罚项具有沉余性和多分辨率结构的特点,结合L0范数,对灰度值进行快速分解重建,这样能够稀疏的逼近图像的分段光滑函数,降低阶梯伪影,提高图像空间分辨率,其目标函数定义为:式中,/>是保真项,C是可以自适应选择图像不同区域的微分算子,g是灰度值,||λ·κ||0是惩罚项,λ是正则化参数,κ是优化系数;
步骤七:推导出目标函数的拉格朗日形式:
式中,m是辅助变量,μ,γ是设置的参数;
步骤八:根据步骤七拉格朗日形式的迭代算法如下:
(2)更新优化后灰度值uk+1:uk+1=(CTC+(μ+γ)I)-1(CTg+γuk+μT(κYk-mk));
(3)更新变量κk+11:κk+1=Hλ,μ,γ(Yuk+1+mk,κk);
(4)更新变量mk+1:mk+1=mk+Yuk+1-κk+1;
(5)判断迭代是否符合迭代终止条件uk+1是否收敛,若是则迭代终止,进行下一步操作;若否,设置k=k+1并跳回步骤八的第(2)步,继续迭代求解;
步骤九:根据最终求解优化后灰度值u进行成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911302022.2A CN110934586B (zh) | 2019-12-17 | 2019-12-17 | 一种灰度值矩阵快速分解重建的正则化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911302022.2A CN110934586B (zh) | 2019-12-17 | 2019-12-17 | 一种灰度值矩阵快速分解重建的正则化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110934586A CN110934586A (zh) | 2020-03-31 |
CN110934586B true CN110934586B (zh) | 2023-06-20 |
Family
ID=69911574
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911302022.2A Active CN110934586B (zh) | 2019-12-17 | 2019-12-17 | 一种灰度值矩阵快速分解重建的正则化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110934586B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113034633B (zh) * | 2021-02-18 | 2022-11-01 | 河南师范大学 | 一种目标函数准确限定支配的图像重建方法 |
CN113012250B (zh) * | 2021-03-04 | 2022-08-26 | 贵州润源医信智能有限公司 | 一种提高肺部成像质量的图像重建方法 |
CN113034635B (zh) * | 2021-03-04 | 2022-08-23 | 重庆不贰科技(集团)有限公司 | 一种抑制工业成像阶梯伪影的图像重建方法 |
CN113052927B (zh) * | 2021-03-04 | 2024-02-06 | 河南师范大学 | 一种提高空间分辨率的成像检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018099321A1 (zh) * | 2016-11-30 | 2018-06-07 | 华南理工大学 | 一种基于广义树稀疏的权重核范数磁共振成像重建方法 |
CN109934885A (zh) * | 2019-02-28 | 2019-06-25 | 河南师范大学 | 一种锐利边缘保持的电阻层析成像图像重建方法 |
-
2019
- 2019-12-17 CN CN201911302022.2A patent/CN110934586B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018099321A1 (zh) * | 2016-11-30 | 2018-06-07 | 华南理工大学 | 一种基于广义树稀疏的权重核范数磁共振成像重建方法 |
CN109934885A (zh) * | 2019-02-28 | 2019-06-25 | 河南师范大学 | 一种锐利边缘保持的电阻层析成像图像重建方法 |
Non-Patent Citations (1)
Title |
---|
优化Tikhonov迭代法在电容层析成像中的应用;高鹤明;许传龙;王式民;;东南大学学报(自然科学版)(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110934586A (zh) | 2020-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110934586B (zh) | 一种灰度值矩阵快速分解重建的正则化方法 | |
CN108961237B (zh) | 一种基于卷积神经网络的低剂量ct图像分解方法 | |
CN111047663B (zh) | 电学层析成像伪影抑制图像重建方法 | |
CN109934885B (zh) | 一种锐利边缘保持的电阻层析成像图像重建方法 | |
CN111047662B (zh) | 自适应非凸混合全变分正则化工业电阻层析成像方法 | |
CN111616708A (zh) | 一种精准识别脑卒中颅内病变区域的图像重建方法 | |
CN108830875B (zh) | 一种基于残差最小的电阻抗层析成像图像分割方法 | |
CN110720915A (zh) | 一种基于gan的脑部电阻抗断层成像方法 | |
US20210027461A1 (en) | Systems and methods for determining a fluid and tissue volume estimations using electrical property tomography | |
Shi et al. | Reduction of staircase effect with total generalized variation regularization for electrical impedance tomography | |
CN112798654B (zh) | 用于电阻抗层析成像的快速梯度法和自适应雅可比矩阵重构方法 | |
CN110208605B (zh) | 一种抑制阶梯效应的电阻层析成像电导率分布重建的方法 | |
CN111062999B (zh) | 有效保留锐利边缘的生物医学电阻抗层析成像方法 | |
Shi et al. | A non-convex L₁-norm penalty-based total generalized variation model for reconstruction of conductivity distribution | |
Li et al. | SAR-CGAN: Improved generative adversarial network for EIT reconstruction of lung diseases | |
CN112914543A (zh) | 一种用于检测人体肺部肿瘤的电阻抗层析成像方法 | |
CN114052701A (zh) | 一种电容耦合电阻层析成像图像重建方法 | |
CN111161182A (zh) | Mr结构信息约束的非局部均值引导的pet图像部分容积校正方法 | |
CN110992385B (zh) | 抑制伪影与保护边缘的颅内图像重建方法 | |
Liu et al. | Cascade resunet with noise power spectrum loss for low dose ct imaging | |
CN110223360B (zh) | 一种应用于减弱图像重建中阶梯效应的混合阶正则化方法 | |
CN113034635B (zh) | 一种抑制工业成像阶梯伪影的图像重建方法 | |
Tian et al. | A wavelet frame constrained total generalized variation model for imaging conductivity distribution. | |
CN113012250B (zh) | 一种提高肺部成像质量的图像重建方法 | |
CN113052927B (zh) | 一种提高空间分辨率的成像检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |