CN111616708A - 一种精准识别脑卒中颅内病变区域的图像重建方法 - Google Patents
一种精准识别脑卒中颅内病变区域的图像重建方法 Download PDFInfo
- Publication number
- CN111616708A CN111616708A CN202010449581.2A CN202010449581A CN111616708A CN 111616708 A CN111616708 A CN 111616708A CN 202010449581 A CN202010449581 A CN 202010449581A CN 111616708 A CN111616708 A CN 111616708A
- Authority
- CN
- China
- Prior art keywords
- iteration
- solution
- objective function
- expressed
- matrix
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 208000006011 Stroke Diseases 0.000 title claims abstract description 22
- 238000007917 intracranial administration Methods 0.000 title claims abstract description 20
- 230000003902 lesion Effects 0.000 title claims abstract description 20
- 206010008190 Cerebrovascular accident Diseases 0.000 title description 10
- 230000002490 cerebral effect Effects 0.000 title description 10
- 238000002593 electrical impedance tomography Methods 0.000 claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims abstract description 19
- 230000000694 effects Effects 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 230000035945 sensitivity Effects 0.000 claims abstract description 13
- 238000003384 imaging method Methods 0.000 claims abstract description 9
- 230000005284 excitation Effects 0.000 claims description 9
- 238000003325 tomography Methods 0.000 claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 5
- 230000003190 augmentative effect Effects 0.000 claims description 4
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000009191 jumping Effects 0.000 claims description 2
- 230000004069 differentiation Effects 0.000 claims 2
- 239000006227 byproduct Substances 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 34
- 208000016988 Hemorrhagic Stroke Diseases 0.000 abstract description 7
- 208000020658 intracerebral hemorrhage Diseases 0.000 abstract description 7
- 238000009499 grossing Methods 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 19
- 238000002591 computed tomography Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 230000000740 bleeding effect Effects 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 210000005013 brain tissue Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002008 hemorrhagic effect Effects 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000003925 brain function Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000000622 irritating effect Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 230000004199 lung function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
- A61B5/4064—Evaluating the brain
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Neurology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Psychology (AREA)
- Neurosurgery (AREA)
- Physiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种精准识别脑卒中颅内病变区域的图像重建方法,用于电阻抗层析成像重建出血性脑卒中的电导率分布。该方法将电阻抗层析成像问题看作求解一个线性不适定问题,确定目标函数;重建方法的流程为:获取物场信息,求解相对边界测量值向量和灵敏度矩阵;设置初始化参数;使用加权矩阵转换数据保真项;利用迭代方式求解目标函数;判断迭代是否结束;根据求得的灰度值,进行成像。本发明有效的克服了边缘过光滑效应和抑制了阶梯伪影,提高了算法的稳定性、适用性和重建图像的空间分辨率,扩展了总变差正则化算法在医学检测领域的应用范围。
Description
技术领域
本发明属于电阻抗层析成像技术领域,具体涉及一种精准识别脑卒中颅内病变区域的图像重建方法。
背景技术
目前,出血性脑卒中被认为是一种严重的脑疾病,它会使脑组织暴露于血液的刺激作用下,导致患者的残疾或死亡。在出血性脑卒中的诊断中,计算机断层扫描(ComputedTomography,CT)和磁共振成像(Magnetic Resonance Imaging,MRI)是两种常用的影像学技术,具有较高的空间分辨率。然而,这两种技术都无法检测出血性中风的早期过程,阻碍了及时治疗。另外,CT和MRI的时间分辨率很低,成像时间较长,设备体积庞大,价格相对昂贵。
电阻抗层析成像(Electrical Impedance Tomography,EIT)技术是一种可以重建检测区域中的电导率分布技术,目前引起了人们的广泛关注。生物组织的生理或病理变化会引起电导率的变化,这可以通过用EIT获得的重建图像来反映。同时,与CT和MRI相比,EIT技术便携、安全、经济且无创。此外,还可以获得较高的时间分辨率。EIT技术由于其优越性,已被用于肺功能、乳腺癌、肝脂肪和脑功能的医学成像。特别是当出血性脑卒中发生时,由于血液的电导率高于正常脑组织的电导率,电导率的分布会发生变化,使得EIT甚至在出血性脑卒中的早期监测成为可能。然而,从EIT的测量数据恢复电导率分布通常是一个不适定反问题,不足以保证解的存在性、唯一性和稳定性,从而导致了重建图像的空间分辨率相对较低。
对于不适定问题,为了得到有效稳定的逆问题的解,通常选用正则化方法来逼近真实解。目前,典型的正则化方法有基于L2范数的Tikhonov正则化(TikhonovRegularization,TR)方法和基于L1范数的总变差(Total Variation,TV)正则化方法。Tikhonov正则化算法由于其简单性得到了广泛的应用,但图像重建质量依赖于正则化参数的选取,如M.Rezghi等人2009年发表于《计算和应用数学杂志》(Journal ofComputational and Applied Mathematics)第231卷,第914-924页,题为《基于Tikhonov正则化的L-曲线法的新的变形方法》(A new variant of L-curve for Tikhonovregularization)。Tikhonov正则化算法对于被测介质连续分布具有良好的性能。然而,当被测介质不连续分布时,在边界上施加了过度的光滑性,从而降低了重建图像的分辨率。为了保持锐利的边缘,总变差正则化算法提供了一种解决方案。该方法保留了边界的不连续性,并允许重建锐利的边缘,以产生更清晰的图像,在EIT领域受到了广泛的关注,如G González等人2017年发表于《计算机和数学及其应用》(Computers&Mathematics withApplications)第74卷,第564-576页,题为《电阻抗层析成像中的各向同性和各向异性总变差正则化》(Isotropic and anisotropic total variation regularization inelectrical impedance tomography)。然而,由于总变差方法往往会产生分段常数解,使得该方法在重建锐利边缘的同时也会在光滑区域产生“阶梯效应”,降低了重建图像的整体分辨率,严重限制了总变差方法在EIT中的应用。
为了减轻总变差正则化方法的阶梯效应,提高其重建图像的空间分辨率,本发明提出一种精准识别脑卒中颅内病变区域的图像重建方法来求解电阻抗层析成像的逆问题。通过比较该方法在被测区域内绝对重建和相对重建的图像效果,验证本发明所提算法的稳定性和适用性。
发明内容
本发明解决的技术问题是提供了一种具有较好稳定性和适用性的精准识别脑卒中颅内病变区域的图像重建方法,用于EIT重建出血性脑卒中的电导率分布。为了降低了求解目标函数的难度,提出一种加权矩阵对数据保真项进行转换;为了进一步去除伪影,提高鲁棒性,引入了软阈值算子;为了求解最小化的目标函数,采用了一种有效的迭代算法;为了验证所提方法重建图像的分辨率,使用数值模拟研究了绝对重建和相对重建的图像效果。与Tikhonov正则化方法和总变差正则化方法相比,本发明提出的精准识别脑卒中颅内病变区域的图像重建方法在电阻抗层析成像重建图像方面显示出了其有效抑制阶梯伪影和精准识别物场介质的优势。
本发明为解决上述技术问题采用如下技术方案,一种精准识别脑卒中颅内病变区域的图像重建方法,该方法将电阻抗层析成像作为一个线性不适定问题Ag=bmeas,式中,A为灵敏度矩阵,bmeas为相对边界测量值向量,g为所求成像灰度值,所建立的最小化目标函数为:
其中
式中,S(g)是目标函数;λ是正则化参数,为一个正标量,用来保持保真度项和混合惩罚项之间的平衡;ω∈[0,1]是控制一阶TV和高阶TV惩罚项之间权衡的加权因子;||·||p和||·||1分别是向量Lp范数和L1范数;表示梯度;lx和ly分别表示在x方向和y方向上的有限元剖分;和分别是一阶和二阶离散微分。
重建算法的流程如下:
本发明的有益效果是:本发明所述的精准识别脑卒中颅内病变区域的图像重建方法在计算目标函数时,为了保证算法的稀疏性和降低求解目标函数的难度,引入一种加权矩阵对目标函数进行转换;为了进一步去除伪影,提高鲁棒性,引入了软阈值算子;最后,为了求解最小化的目标函数,采用了一种有效地迭代算法。本发明提高了算法的稳定性和适用性,扩展了总变差正则化算法的应用范围,有效的克服了传统Tikhonov正则化算法的边缘过光滑效应和总变差正则化算法图像阶梯伪影严重的特点,在提高电阻抗层析成像技术重建图像的空间分辨率方面有很好的效果。
附图说明
图1为本发明的一种精准识别脑卒中颅内病变区域的图像重建方法的流程框图;
图2为本发明的电阻层析成像系统圆形单截面被测场域,激励电流和测量电压的模式以及电极分布;
图3为本发明的有限元剖分网格(a)正问题有限元剖分(b)逆问题有限元剖分;
图4为本发明在选取六种模型(模型1-6)的真实分布时,典型的TR算法,TV算法和一种精准识别脑卒中颅内病变区域的图像重建方法,绝对重建的图像结果示意图;
图5为本发明在选取两组模型(模型7-12)的真实分布时,典型的TR算法,TV算法和一种精准识别脑卒中颅内病变区域的图像重建方法,相对重建的图像结果示意图;
图6为相对重建模型空间分辨率的定量评价。
图中:1-被测场域,2-电极,3-激励电流,4-测量电压。
具体实施方式
结合附图和实施例对本发明的一种精准识别脑卒中颅内病变区域的图像重建方法加以说明。
本发明的一种精准识别脑卒中颅内病变区域的图像重建方法,用于EIT重建出血性脑卒中的电导率分布。为了降低了求解目标函数的难度,提出一种加权矩阵对数据保真项进行转换;同时,为了进一步去除伪影,提高鲁棒性,引入了软阈值算子;最后,利用一种有效的迭代算法完成最终的逆问题求解。
如图1所示,为本发明的一种精准识别脑卒中颅内病变区域的图像重建方法流程图。
如图2所示为电学层析成像之一的电阻层析成像系统圆形单截面被测场域,激励电流3和测量电压4的模式以及电极2分布,采用16电极2均匀分布在被测场域1外壁。
采取有限元剖分的方式来对EIT进行求解,如图3所示。其中,对于EIT正问题采用三角形网格剖分,如图3(a)所示;EIT逆问题采用正方形网格剖分,如图3(b)所示;为了避免图像重建过程中逆向违规对成像质量的影响,正问题的剖分网格数要大于逆问题的剖分网格数。
为了较好地体现本发明一种精准识别脑卒中颅内病变区域的图像重建方法与典型TR算法和TV正则化算法的不同,分别给出了六种绝对重建模型和两组相对重建模型在这三种正则化算法下的求解结果。
实施例包括如下具体步骤:一种精准识别脑卒中颅内病变区域的图像重建方法,该方法将电学层析成像作为一个线性不适定问题Ag=bmeas。式中,A为灵敏度矩阵,bmeas为相对边界测量值向量,g为成像灰度值。
逆问题可以用最小二乘优化形式的目标函数表示:式中,f(g)为目标函数。为了解决逆问题的不定性问题,正则化方法为提高解的稳定性提供了一种新的方法。在最小二乘的目标函数中增加正则化项,也称为惩罚项。正则化方法的一般形式可以描述为:其中,λ是一个正则化因子,它控制最小二乘项(通常称为数据保真度项)和正则化项R(g)之间的权衡。
TR算法可将正则化项R(g)替换为正则化项可以描述为:然而,TR算法在被测介质不连续分布时,在边界上施加了过度的光滑性,从而降低了重建图像的分辨率。为了保持锐利的边缘,TV正则化算法提供了一种解决方案。TV正则化算法的形式可以描述为:虽然重建图像保持了边缘上不连续性,但会产生分段常数解,使得该方法在重建锐利边缘的同时也会在光滑区域产生“阶梯效应”,从而降低了重建图像的整体分辨率。
为了抑制TV正则化方法的阶梯效应和提高重建图像的分辨率,本发明针对数据保真项和正则化项,提出了一种精准识别脑卒中颅内病变区域的图像重建方法来求解电阻抗层析成像的逆问题。其最小化的目标函数为:
其中
式中,S(g)是目标函数;λ是正则化参数,为一个正标量,用来保持保真度项和混合惩罚项之间的平衡;ω∈[0,1]是控制一阶TV和高阶TV惩罚项之间权衡的加权因子;||·||p和||·||1分别是向量Lp范数和L1范数;表示梯度;lx和ly分别表示在x方向和y方向上的有限元剖分;和分别是一阶和二阶离散微分。重建步骤如下:
(1)获取物场信息:
针对重建模型,分别获取各自重建所需的边界测量值和灵敏度矩阵。边界测量值是将被测对象置于电学层析成像测量系统中,被测场域外均匀分布16个电极(如图2所示),采用电流激励电压测量且激励电极不测量的模式,采集循环激励循环测量下的边界电压,共获得208个测量值;灵敏度矩阵是根据不含内含物的空场的边界测量电压,结合灵敏度理论,计算灵敏度矩阵,计算公式为:
式中Aij是第j个电极对对第i个电极对的灵敏度系数,φi,φj分别为第i个电极对及第j个电极对在激励电流为Ii,Ij时场域电势分布。
(2)设置初始化参数:
正则化参数λ=75,加权因子ω=0.3,范数p=0.9,正常数ε=10-6,迭代终止θmin=10-4,最大迭代次数kmax=50,惩罚参数β1=100,β2=6,β3=6,软阈值系数κ=0.08,拉格朗日乘子和初始解g0=0n×1。
(3)使用加权矩阵转换数据保真项:
为了降低了求解目标函数的难度,使用一种加权矩阵对数据保真项进行转换,则最小化目标函数可以转变为:
其中
式中,Wk是在上一次迭代中计算并在每次迭代中更新的加权矩阵,gk+1和gk是第(k+1)次解的列向量和第k次解的列向量,ε是用于防止分母为零的正小常数,选择为10-6。
(4)对目标函数进行迭代求解:
通过引入多个辅助变量,将复杂约束问题分解为一系列简单约束子问题。同时,每个子问题都有闭式解,以便用现有的数值方法求解。为此,引入了三个辅助变量f、q和r,可表示:
其增广拉格朗日函数可以表示为:
式中,γ1,γ2和γ3表示增广拉格朗日乘子,β1,β2,β3>0是惩罚参数。因此,辅助变量和灰度值的解可以描述为:
对于f子问题可表示为:
对于q子问题可表示为:
对于r子问题可表示为:
对于灰度值g的子问题可以表示为:
转换可得:
为了进一步减少阶梯效应,引入了一个感兴趣的区域,定义:
式中,R指感兴趣的区域,κ是区分背景和介质的阈值系数。因此,灰度值的解g可转化为:
gk+1=R.*gk+1
式中,.表示的是元素与元素之间的对应相乘;
最后,拉格朗日乘子的计算和更新方式可以表示为:
(5)迭代终止准则:
(6)输出结果:
根据最终求解所得灰度值g,进行成像。
图4为本发明选取六种典型的介质模型,传统的TR算法,TV正则化算法和本方法,绝对重建图像的结果示意图。从重建图像中可以看出,用TR算法重建图像时,在重建图像的边界上观察到过大的平滑度,难以准确识别图像的大小和形状。与TR算法相比,TV方法在重建过程中避免了TR算法出现过多的平滑度,具有更好的重建效果。然而,在重建图像的背景中却发现了严重的阶梯效应。与TR和TV方法相比,本方法在重建过程中很好地消除了过度平滑和阶梯效应。重建图像的质量最高,背景清晰,边界清晰,能准确地检测出颅脑中出血组织的大小和位置。
在临床出血性脑卒中研究中,随着出血过程的发展,重建电导率分布的相对变化具有重要意义。相对重建的目的是获得两个时间间隔之间电导率分布情况的变化(模型9和模型12),还描绘了模型7、模型8、模型10和模型11的绝对重建,以便清楚地了解算法的性能,如图5所示。很明显,TR方法在绝对重建中很难区分介质的边界,而在相对重建中则在重建图像的边界上观察到模糊现象。对于TV方法,在绝对重建和相对重建中,可以很好地识别边界,同时观察到严重的阶梯效应。另外,对于TR和TV两种方法,重建的介质都大于实际的变化。在绝对和相对重建图像中,采用本发明所提方法可以有效地抑制模糊,消除阶梯效应。此外,重构后的介质尺寸与实际尺寸几乎相同。同时,采用空间分辨率(SpatialResolution,SRE)来定量评价图像重建质量,其表达式可以表示为:式中,其中,其中Pd和P0分别表示介质和整个检测区域的面积,Pi是第i个像素的面积,[gd]i=1(如果gi大于g的最大振幅的四分之一),否则选择[gd]i=0。
从图6可以看出,一种精准识别脑卒中颅内病变区域的图像重建方法与典型TR算法和TV正则化算法相比,具有最小的空间分辨率,能够有效的抑制边缘模糊,消除阶梯效应,并准确的描绘出被测场域内部介质分布,明显地提高了电阻抗层析成像逆问题求解精度和图像重建质量。
以上所述仅为本发明的较佳实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种精准识别脑卒中颅内病变区域的图像重建方法,其特征在于:将电阻抗层析成像作为一个线性不适定逆问题Ag=bmeas,式中,A为灵敏度矩阵,bmeas为相对边界测量值向量,g为成像灰度值,建立最小化目标函数为:
其中
式中,S(g)是目标函数;λ是正则化参数,为一个正标量,用来保持保真度项和混合惩罚项之间的平衡;ω∈[0,1]是控制一阶TV和高阶TV惩罚项之间权衡的加权因子;||·||p和||·||1分别是向量Lp范数和L1范数;表示梯度;lx和ly分别表示在x方向和y方向上的有限元剖分;和分别是一阶和二阶离散微分,具体重建步骤如下:
(1)获取物场信息:
针对重建模型分别获取重建所需的边界测量值和灵敏度矩阵,边界测量值是将被测对象置于电学层析成像测量系统中,被测场域外均匀分布16个电极,采用电流激励电压测量且激励电极不测量的模式,采集循环激励循环测量下的边界电压,共获得208个测量值;灵敏度矩阵是根据不含介质的空场的边界测量电压,结合灵敏度理论,计算灵敏度矩阵,计算公式为:式中,Aij是第j个电极对对第i个电极对的灵敏度系数,φi,φj分别为第i个电极对及第j个电极对在激励电流为Ii,Ij时场域电势分布;
(2)设置初始化参数:
正则化参数λ=75,加权因子ω=0.3,范数p=0.9,正常数ε=10-6,迭代终止θmin=10-4,最大迭代次数kmax=50,惩罚参数β1=100,β2=6,β3=6,软阈值系数κ=0.08,拉格朗日乘子和初始解g0=0n×1;
(3)使用加权矩阵转换数据保真项:
为了降低求解目标函数的难度,使用一种加权矩阵对数据保真项进行转换,则最小化目标函数转变为:其中,式中,Wk是在上一次迭代中计算并在每次迭代中更新的加权矩阵,gk+1和gk是第(k+1)次解的列向量和第k次解的列向量,ε是用于防止分母为零的正小常数,选择为10-6;
(4)对目标函数进行迭代求解:
通过引入多个辅助变量,将复杂约束问题分解为一系列简单约束子问题,同时每个子问题都有闭式解,以便用现有的数值方法求解,为此,引入了三个辅助变量f、q和r,具体表示为:
其增广拉格朗日函数表示为:
式中,γ1,γ2和γ3是增广拉格朗日乘子,β1,β2,β3>0是惩罚参数;
辅助变量和灰度值的解描述为:
对于f子问题表示为:式中“°”是逐点积,“sgn”表示signum的函数;对于q子问题表示为:对于r子问题表示为:对于灰度值g的子问题表示为:转换得到:其中,上标T表示转置运算符,和ATA为循环矩阵;
为了进一步减少阶梯效应,引入一个感兴趣的区域,定义:式中,R指感兴趣的区域,κ是区分背景和介质的阈值系数,灰度值的解g转化为:gk+1=R.*gk+1,其中,.表示的是元素与元素之间的对应相乘;拉格朗日乘子的计算和更新方式表示为:
(5)迭代终止准则:
判断迭代是否符合迭代终止条件θ=||gk+1-gk||2/||gk||2≤θmin或者是否达到最大迭代次数k≤kmax,若是则迭代终止,进行下一步操作;若否,设置k=k+1并跳回第(3)步继续迭代求解;
(6)输出结果:
根据最终求解所得灰度值g进行成像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010449581.2A CN111616708A (zh) | 2020-05-25 | 2020-05-25 | 一种精准识别脑卒中颅内病变区域的图像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010449581.2A CN111616708A (zh) | 2020-05-25 | 2020-05-25 | 一种精准识别脑卒中颅内病变区域的图像重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111616708A true CN111616708A (zh) | 2020-09-04 |
Family
ID=72255138
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010449581.2A Pending CN111616708A (zh) | 2020-05-25 | 2020-05-25 | 一种精准识别脑卒中颅内病变区域的图像重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111616708A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112798654A (zh) * | 2021-01-26 | 2021-05-14 | 河南师范大学 | 用于电阻抗层析成像的快速梯度法和自适应雅可比矩阵重构方法 |
CN113012250A (zh) * | 2021-03-04 | 2021-06-22 | 施成成 | 一种提高肺部成像质量的图像重建方法 |
CN113034635A (zh) * | 2021-03-04 | 2021-06-25 | 施成成 | 一种抑制工业成像阶梯伪影的图像重建方法 |
CN113066146A (zh) * | 2021-04-07 | 2021-07-02 | 太原理工大学 | 一种基于Tikhonov技术的改进ECT图像重建方法 |
CN113384257A (zh) * | 2021-06-15 | 2021-09-14 | 河南师范大学 | 针对脑出血患者病变部位长期高精度eit检测的电极排布方法 |
CN114373025A (zh) * | 2021-12-06 | 2022-04-19 | 山东师范大学 | Cbct图像重建方法、系统、存储介质及设备 |
CN115444392A (zh) * | 2022-08-31 | 2022-12-09 | 河南师范大学 | 一种基于电阻抗层析成像的非线性脑卒中分析方法 |
CN115830156A (zh) * | 2022-05-10 | 2023-03-21 | 深圳市元甪科技有限公司 | 精准电阻抗层析成像方法、装置、系统、介质及设备 |
CN117649503A (zh) * | 2024-01-29 | 2024-03-05 | 杭州永川科技有限公司 | 图像重建方法、装置、计算机设备、存储介质和程序产品 |
-
2020
- 2020-05-25 CN CN202010449581.2A patent/CN111616708A/zh active Pending
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112798654B (zh) * | 2021-01-26 | 2023-06-20 | 河南师范大学 | 用于电阻抗层析成像的快速梯度法和自适应雅可比矩阵重构方法 |
CN112798654A (zh) * | 2021-01-26 | 2021-05-14 | 河南师范大学 | 用于电阻抗层析成像的快速梯度法和自适应雅可比矩阵重构方法 |
CN113034635B (zh) * | 2021-03-04 | 2022-08-23 | 重庆不贰科技(集团)有限公司 | 一种抑制工业成像阶梯伪影的图像重建方法 |
CN113012250A (zh) * | 2021-03-04 | 2021-06-22 | 施成成 | 一种提高肺部成像质量的图像重建方法 |
CN113034635A (zh) * | 2021-03-04 | 2021-06-25 | 施成成 | 一种抑制工业成像阶梯伪影的图像重建方法 |
CN113012250B (zh) * | 2021-03-04 | 2022-08-26 | 贵州润源医信智能有限公司 | 一种提高肺部成像质量的图像重建方法 |
CN113066146A (zh) * | 2021-04-07 | 2021-07-02 | 太原理工大学 | 一种基于Tikhonov技术的改进ECT图像重建方法 |
CN113066146B (zh) * | 2021-04-07 | 2023-06-27 | 太原理工大学 | 一种基于Tikhonov技术的改进ECT图像重建方法 |
CN113384257A (zh) * | 2021-06-15 | 2021-09-14 | 河南师范大学 | 针对脑出血患者病变部位长期高精度eit检测的电极排布方法 |
CN113384257B (zh) * | 2021-06-15 | 2024-05-14 | 河南师范大学 | 针对脑出血患者病变部位长期高精度eit检测的电极排布方法 |
CN114373025A (zh) * | 2021-12-06 | 2022-04-19 | 山东师范大学 | Cbct图像重建方法、系统、存储介质及设备 |
CN115830156A (zh) * | 2022-05-10 | 2023-03-21 | 深圳市元甪科技有限公司 | 精准电阻抗层析成像方法、装置、系统、介质及设备 |
CN115830156B (zh) * | 2022-05-10 | 2024-06-07 | 深圳市元甪科技有限公司 | 精准电阻抗层析成像方法、装置、系统、介质及设备 |
CN115444392A (zh) * | 2022-08-31 | 2022-12-09 | 河南师范大学 | 一种基于电阻抗层析成像的非线性脑卒中分析方法 |
CN115444392B (zh) * | 2022-08-31 | 2024-05-14 | 河南师范大学 | 一种基于电阻抗层析成像的非线性脑卒中分析方法 |
CN117649503B (zh) * | 2024-01-29 | 2024-05-28 | 杭州永川科技有限公司 | 图像重建方法、装置、计算机设备、存储介质和程序产品 |
CN117649503A (zh) * | 2024-01-29 | 2024-03-05 | 杭州永川科技有限公司 | 图像重建方法、装置、计算机设备、存储介质和程序产品 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111616708A (zh) | 一种精准识别脑卒中颅内病变区域的图像重建方法 | |
CN101564294B (zh) | 一种结构信息融合的电阻抗断层成像方法 | |
CN112200306B (zh) | 一种基于深度学习的电阻抗成像方法 | |
Miller et al. | Diffeomorphometry and geodesic positioning systems for human anatomy | |
Liu et al. | Deep iterative reconstruction estimation (DIRE): approximate iterative reconstruction estimation for low dose CT imaging | |
WO2010062220A1 (ru) | Способ неинвазивного электро физиологического исследования сердца | |
US11605162B2 (en) | Systems and methods for determining a fluid and tissue volume estimations using electrical property tomography | |
CN110934586B (zh) | 一种灰度值矩阵快速分解重建的正则化方法 | |
CN111047663A (zh) | 电学层析成像伪影抑制图像重建方法 | |
CN110720915A (zh) | 一种基于gan的脑部电阻抗断层成像方法 | |
CN113902825A (zh) | 基于VDD-Net的肺部电阻抗成像方法 | |
CN114708350A (zh) | 一种用于脑部电阻抗层析成像的电导率可视化方法 | |
CN110811596B (zh) | 基于低秩与稀疏约束和非局部全变分的无创心脏电位重建方法 | |
Shi et al. | Sparse image reconstruction of intracerebral hemorrhage with electrical impedance tomography | |
CN116869504A (zh) | 一种用于脑缺血电导率分布重建的数据补偿方法 | |
Zhang et al. | DREAM-Net: Deep residual error iterative minimization network for sparse-view CT reconstruction | |
CN102599907B (zh) | 基于网格位移模型的电阻抗断层成像方法 | |
CN112914543A (zh) | 一种用于检测人体肺部肿瘤的电阻抗层析成像方法 | |
Chen et al. | Temporal feature prior-aided separated reconstruction method for low-dose dynamic myocardial perfusion computed tomography | |
Huang et al. | Finite element modeling of human thorax based on MRI images for EIT image reconstruction | |
CN111192209B (zh) | 一种基于深度学习的磁感应断层成像优化方法 | |
Yu et al. | Deep learning single view computed tomography guided by FBP algorithm | |
CN113052927B (zh) | 一种提高空间分辨率的成像检测方法 | |
CN113012250B (zh) | 一种提高肺部成像质量的图像重建方法 | |
Chen et al. | Global inhomogeneity index evaluation of a DCT-based EIT lung imaging |
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 |