CN110827215A - ERT image reconstruction artifact removing method based on fuzzy clustering - Google Patents
ERT image reconstruction artifact removing method based on fuzzy clustering Download PDFInfo
- Publication number
- CN110827215A CN110827215A CN201910998118.0A CN201910998118A CN110827215A CN 110827215 A CN110827215 A CN 110827215A CN 201910998118 A CN201910998118 A CN 201910998118A CN 110827215 A CN110827215 A CN 110827215A
- Authority
- CN
- China
- Prior art keywords
- grayscale
- image reconstruction
- fuzzy clustering
- artifact
- clusters
- 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 9
- 239000013598 vector Substances 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 238000007621 cluster analysis Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 description 6
- 238000003325 tomography Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于模糊聚类的ERT图像重建伪影去除方法,包括下列步骤:1)利用至少两种图像重建算法分别计算被测场域的灰度向量;2)将1)中灰度向量结合起来形成一个灰度矩阵作为聚类算法中的特征向量;3)将模糊聚类算法作用于灰度矩阵,对灰度矩阵进行聚类分析;4)按照聚类簇表现出来的统计特征找出分别代表“目标”、“背景”和“伪影”的聚类簇;5)将代表伪影的聚类簇划归到代表背景的聚类簇。
The present invention relates to an ERT image reconstruction artifact removal method based on fuzzy clustering, comprising the following steps: 1) using at least two image reconstruction algorithms to respectively calculate the grayscale vector of the measured field; 2) 1) medium grayscale The vectors are combined to form a gray matrix as the feature vector in the clustering algorithm; 3) The fuzzy clustering algorithm is applied to the gray matrix, and the gray matrix is clustered; 4) According to the statistical characteristics of the clusters Find out the clusters representing "target", "background" and "artifact"respectively; 5) classify the cluster representing the artifact into the cluster representing the background.
Description
技术领域technical field
本发明属于电学层析成像图像重建算法领域,具体涉及一种基于模糊聚类的ERT图像重建伪影去除方法。The invention belongs to the field of electrical tomography image reconstruction algorithms, and in particular relates to an ERT image reconstruction artifact removal method based on fuzzy clustering.
背景技术Background technique
电学层析成像技术(Electrical Tomography,简称ET)是一种重建被测物场内部电学特性分布的技术,该技术采用的硬件设备结构较为简单,成本较为低廉,且检测过程具有无辐射、非侵入性的特点。该技术在医学临床监护和工业测量等领域有着广阔的应用前景。Electrical tomography (ET) is a technology for reconstructing the distribution of electrical properties inside the object field under test. The hardware equipment used in this technology is relatively simple in structure and low in cost. sexual characteristics. The technology has broad application prospects in the fields of medical clinical monitoring and industrial measurement.
在电阻层析成像(ERT)中,利用阵列电极对被测物场施加电流激励,测量边界电压数据,进而重建被测场域内的电导率分布。以16电极的电阻层析成像传感器为例,采用相邻电流激励、相邻电压测量的方式。在每次激励下,除去激励电极,共有13个测量值。如:以1、2号电极作为激励电极,测量3、4,4、5,5、6,…,15、16电极对上的电压信号。共有16次激励,即有13×16=208个测量值。为了兼顾时间分辨率及空间分辨率,一般将被测敏感场域划分为812个像素点。电阻层析成像图像重建问题即为求解这812个像素点的灰度值。In electrical resistance tomography (ERT), current excitation is applied to the measured object field with array electrodes, boundary voltage data is measured, and then the conductivity distribution in the measured field is reconstructed. Taking a 16-electrode resistance tomography sensor as an example, the adjacent current excitation and adjacent voltage measurement methods are used. At each excitation, with the excitation electrodes removed, a total of 13 measurements were made. For example, take the No. 1 and No. 2 electrodes as excitation electrodes, and measure the voltage signals on the 3, 4, 4, 5, 5, 6, ..., 15, 16 electrode pairs. There are 16 excitations in total, that is, 13×16=208 measurements. In order to take into account both the temporal resolution and the spatial resolution, the measured sensitive field is generally divided into 812 pixels. The reconstruction problem of resistance tomography image is to solve the gray value of these 812 pixels.
由于ERT的“软场”效应及逆问题本身的不适定性,导致重建图像的空间分辨率较差。传统的线性反投影算法(LBP)即为根据边界电压值和被测场域投影区的关系进行电压数据的反投,但由于硬件设备的限制,导致反投数据远远少于需要求解的灰度值,使得重建图像有较大的伪影,目标及背景的边界模糊,空间分辨率较差。Due to the "soft field" effect of ERT and the ill-posed nature of the inverse problem, the spatial resolution of the reconstructed image is poor. The traditional linear back-projection algorithm (LBP) is to back-project the voltage data according to the relationship between the boundary voltage value and the projection area of the measured field. However, due to the limitation of hardware equipment, the back-projection data is far less than the grayscale that needs to be solved. If the degree value is high, the reconstructed image has large artifacts, the boundary between the target and the background is blurred, and the spatial resolution is poor.
发明内容SUMMARY OF THE INVENTION
本发明的目的在于克服现有技术的上述不足,提出一种基于模糊聚类的ERT图像重建伪影去除方法,可以有效地减小伪影对图像分辨率造成的影响。重建图像的所有像素点可以被分为3类:目标、背景和伪影。其中,伪影在很大程度上影响了成像的分辨率,对伪影的有效处理可以在一定程度上提高重建图像的质量。考虑到上述算法得到的结果是灰度值,即是数值类数据,因此,可以利用聚类算法对灰度值进行划分,再根据划分结果对灰度值进行进一步的处理,从而得到较为理想的成像效果。The purpose of the present invention is to overcome the above-mentioned shortcomings of the prior art, and to propose a method for removing artifacts in ERT image reconstruction based on fuzzy clustering, which can effectively reduce the impact of artifacts on image resolution. All pixels of the reconstructed image can be divided into 3 categories: object, background and artifact. Among them, artifacts greatly affect the resolution of imaging, and effective processing of artifacts can improve the quality of reconstructed images to a certain extent. Considering that the result obtained by the above algorithm is the gray value, that is, the numerical data, therefore, the gray value can be divided by the clustering algorithm, and then the gray value can be further processed according to the division result, so as to obtain the ideal value. Imaging effect.
本发明的技术方案如下:The technical scheme of the present invention is as follows:
一种基于模糊聚类的ERT图像重建伪影去除方法,包括下列步骤:A method for removing artifacts from ERT image reconstruction based on fuzzy clustering, comprising the following steps:
1)利用至少两种图像重建算法分别计算被测场域的灰度向量;1) Utilize at least two kinds of image reconstruction algorithms to calculate the gray level vector of the measured field respectively;
2)将1)中灰度向量结合起来形成一个灰度矩阵作为聚类算法中的特征向量;2) Combine the grayscale vectors in 1) to form a grayscale matrix as the feature vector in the clustering algorithm;
3)将模糊聚类算法作用于灰度矩阵,对灰度矩阵进行聚类分析;3) Apply the fuzzy clustering algorithm to the grayscale matrix, and perform cluster analysis on the grayscale matrix;
4)按照聚类簇表现出来的统计特征找出分别代表“目标”、“背景”和“伪影”的聚类簇;4) Find out the clusters representing "target", "background" and "artifact" respectively according to the statistical characteristics exhibited by the clusters;
5)将代表伪影的聚类簇划归到代表背景的聚类簇。5) The clusters representing artifacts are classified into clusters representing the background.
本发明由于采取以上技术方案,其具有以下优点:The present invention has the following advantages due to taking the above technical solutions:
1、与传统的电学层析成像算法相比,本算法将多个经典的成像算法融合起来,将其优势互补,并在此基础上对灰度值进行进一步的聚类分析,减小了伪影的影响,使目标与背景之间的界限更加清晰,提高了成像的空间分辨率。1. Compared with the traditional electrical tomography algorithm, this algorithm integrates several classical imaging algorithms to complement each other's advantages. On this basis, the gray value is further clustered and analyzed to reduce the false value. The influence of the shadow can make the boundary between the target and the background clearer and improve the spatial resolution of imaging.
2、聚类算法采用模糊聚类算法而非传统的硬划分算法,使得目标边缘更加平滑,更加贴近实际边缘。2. The clustering algorithm adopts the fuzzy clustering algorithm instead of the traditional hard division algorithm, which makes the target edge smoother and closer to the actual edge.
附图说明Description of drawings
图1为本发明的实施例中原始被测物场模型图;1 is a model diagram of an original object field under test in an embodiment of the present invention;
图2为对图1所示仿真模型用LBP算法进行图像重建的结果;Fig. 2 is the result of image reconstruction using LBP algorithm to the simulation model shown in Fig. 1;
图3为对图1所示仿真模型用本发明的算法进行图像重建的结果。FIG. 3 is the result of image reconstruction using the algorithm of the present invention for the simulation model shown in FIG. 1 .
具体实施方式Detailed ways
本发明的目的是克服现有图像重建算法空间分辨率较低的不足,提出了一种基于模糊聚类的ERT图像重建伪影去除方法。模糊聚类的对象是一组数据的集合,通过提取数据的特征向量来对数据进行划分,使类内差别较小,类间差异较大。通过对灰度矩阵进行聚类分析,减小了伪影对图像分辨率带来的影响。The purpose of the present invention is to overcome the problem of low spatial resolution of the existing image reconstruction algorithms, and propose a method for removing artifacts from ERT image reconstruction based on fuzzy clustering. The object of fuzzy clustering is a set of data, and the data is divided by extracting the feature vector of the data, so that the intra-class difference is small and the inter-class difference is large. By clustering the grayscale matrix, the influence of artifacts on image resolution is reduced.
下面结合附图和实施例对本发明的一种基于模糊聚类的ERT图像重建伪影去除方法进行详细的描述。具体步骤如下:A method for removing artifacts from ERT image reconstruction based on fuzzy clustering of the present invention will be described in detail below with reference to the accompanying drawings and embodiments. Specific steps are as follows:
1)图1表示含有一个离散目标介质的流型,对图1所示被测物场,若采用正则化(TK)算法和LBP算法进行灰度向量的计算,则重建图像的灰度值可以用公式(1)所示。1) Figure 1 shows the flow pattern with a discrete target medium. For the measured object field shown in Figure 1, if the regularization (TK) algorithm and the LBP algorithm are used to calculate the grayscale vector, the grayscale value of the reconstructed image can be It is shown in formula (1).
GT={gt1,gt2,...,gt812},GL={gl1,gl2,...,gl812} (1) GT = {g t1 ,g t2 ,...,g t812 }, GL ={g l1 , g l2 ,...,g l812 } (1)
其中,GT表示由TK算法得到的灰度向量,GL表示由LBP算法得到的灰度向量。Among them, GT represents the grayscale vector obtained by the TK algorithm, and GL represents the grayscale vector obtained by the LBP algorithm.
由LBP算法重建得到的图像如图2所示,上方类圆型面积较小的区域代表目标,下方面积最大的区域为背景区域,介于二者之间的灰白色的区域为伪影区域。The image reconstructed by the LBP algorithm is shown in Figure 2. The area with a small circular area at the top represents the target, the area with the largest area at the bottom is the background area, and the gray-white area in between is the artifact area.
2)将1)中灰度向量结合起来形成一个灰度矩阵,作为聚类中的特征向量,即结合灰度向量GT,GL形成一个灰度矩阵,并计算灰度均值向量如公式(2)所示:2) Combine the grayscale vectors in 1) to form a grayscale matrix, which is used as a feature vector in the clustering, that is, combine the grayscale vectors G T and G L to form a gray matrix, and calculate the gray mean vector As shown in formula (2):
3)利用模糊聚类算法对2)中灰度矩阵进行聚类分析。模糊聚类算法需要给定聚类数,由逆问题本身的意义可知,可以将聚类数设置为3,即将数据分成三类:目标Go、背景Gb和伪影Gt,即3) Use fuzzy clustering algorithm to perform cluster analysis on the gray matrix in 2). The fuzzy clustering algorithm requires a given number of clusters. From the meaning of the inverse problem itself, the number of clusters can be set to 3, that is, the data can be divided into three categories: target G o , background G b and artifact G t , namely
4)计算各个聚类簇的统计特征。步骤3)中的3个聚类簇具有不同的统计特征:具有最大均值的一类为目标,与之相反,具有最小均值的一类为背景,具有最大方差的一类代表伪影。按照上述特征可以明确聚类结果中每个聚类簇代表的意义。4) Calculate the statistical characteristics of each cluster. The three clusters in step 3) have different statistical characteristics: the one with the largest mean is the target, on the contrary, the one with the smallest mean is the background, and the one with the largest variance represents the artifact. According to the above characteristics, the meaning of each cluster in the clustering result can be clarified.
5)将代表伪影的聚类划归到代表背景的聚类,即将中像素点的灰度值赋值为中像素点的平均灰度值,记为即最终的灰度向量为:5) Classify the cluster representing the artifact into the cluster representing the background, namely The gray value of the middle pixel is assigned as The average gray value of the middle pixel, denoted as That is, the final grayscale vector is:
6)利用公式(4)中进行图像反演,如图3所示。图3中伪影区域被有效地减小,更加接近原始物场分布。6) Using formula (4) in Perform image inversion, as shown in Figure 3. The artifact area in Fig. 3 is effectively reduced, closer to the original object field distribution.
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910998118.0A CN110827215A (en) | 2019-10-21 | 2019-10-21 | ERT image reconstruction artifact removing method based on fuzzy clustering |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910998118.0A CN110827215A (en) | 2019-10-21 | 2019-10-21 | ERT image reconstruction artifact removing method based on fuzzy clustering |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110827215A true CN110827215A (en) | 2020-02-21 |
Family
ID=69549796
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910998118.0A Pending CN110827215A (en) | 2019-10-21 | 2019-10-21 | ERT image reconstruction artifact removing method based on fuzzy clustering |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110827215A (en) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2695964Y (en) * | 2003-10-21 | 2005-04-27 | 浙江大学 | Oil-gas two-phase flow investigating device based-on capacitive chromatographic imaging system |
CN101980298A (en) * | 2010-11-04 | 2011-02-23 | 西安电子科技大学 | Image Segmentation Method Based on Multi-agent Genetic Clustering Algorithm |
CN102968779A (en) * | 2012-12-03 | 2013-03-13 | 北方民族大学 | Non-subsample Contourlet domain type MRI (Magnetic Resonance Imaging) image enhancing method based on FCM (Fuzzy C-means) cluster |
CN103654776A (en) * | 2013-11-18 | 2014-03-26 | 中国人民解放军第四军医大学 | Electrical impedance tomography method integrating skull specific resistance non-uniform distribution information |
CN104867150A (en) * | 2015-05-22 | 2015-08-26 | 武汉大学 | Wave band correction change detection method of remote sensing image fuzzy clustering and system thereof |
CN105147310A (en) * | 2015-09-30 | 2015-12-16 | 上海联影医疗科技有限公司 | Breast image obtaining method and device and breast photography system |
CN105741279A (en) * | 2016-01-27 | 2016-07-06 | 西安电子科技大学 | Rough set based image segmentation method for quickly inhibiting fuzzy clustering |
CN106959325A (en) * | 2017-04-07 | 2017-07-18 | 燕山大学 | A kind of planar array electrode capacitance sensor, imaging system and imaging method |
CN107369187A (en) * | 2017-06-20 | 2017-11-21 | 天津大学 | The electricity tomography regularization reconstruction method for the sum that is deteriorated based on adjoint point |
CN107977952A (en) * | 2016-10-21 | 2018-05-01 | 冯原 | Medical image cutting method and device |
CN108717694A (en) * | 2018-04-24 | 2018-10-30 | 天津大学 | Electrical impedance tomography image quality evaluating method based on fuzzy C-means clustering |
CN109086810A (en) * | 2018-07-19 | 2018-12-25 | 天津大学 | A kind of electrical impedance tomography optimization algorithm based on fuzzy clustering |
CN109813772A (en) * | 2019-03-18 | 2019-05-28 | 燕山大学 | A Stable Imaging Method for Coplanar Array Capacitors |
-
2019
- 2019-10-21 CN CN201910998118.0A patent/CN110827215A/en active Pending
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2695964Y (en) * | 2003-10-21 | 2005-04-27 | 浙江大学 | Oil-gas two-phase flow investigating device based-on capacitive chromatographic imaging system |
CN101980298A (en) * | 2010-11-04 | 2011-02-23 | 西安电子科技大学 | Image Segmentation Method Based on Multi-agent Genetic Clustering Algorithm |
CN102968779A (en) * | 2012-12-03 | 2013-03-13 | 北方民族大学 | Non-subsample Contourlet domain type MRI (Magnetic Resonance Imaging) image enhancing method based on FCM (Fuzzy C-means) cluster |
CN103654776A (en) * | 2013-11-18 | 2014-03-26 | 中国人民解放军第四军医大学 | Electrical impedance tomography method integrating skull specific resistance non-uniform distribution information |
CN104867150A (en) * | 2015-05-22 | 2015-08-26 | 武汉大学 | Wave band correction change detection method of remote sensing image fuzzy clustering and system thereof |
CN105147310A (en) * | 2015-09-30 | 2015-12-16 | 上海联影医疗科技有限公司 | Breast image obtaining method and device and breast photography system |
CN105741279A (en) * | 2016-01-27 | 2016-07-06 | 西安电子科技大学 | Rough set based image segmentation method for quickly inhibiting fuzzy clustering |
CN107977952A (en) * | 2016-10-21 | 2018-05-01 | 冯原 | Medical image cutting method and device |
CN106959325A (en) * | 2017-04-07 | 2017-07-18 | 燕山大学 | A kind of planar array electrode capacitance sensor, imaging system and imaging method |
CN107369187A (en) * | 2017-06-20 | 2017-11-21 | 天津大学 | The electricity tomography regularization reconstruction method for the sum that is deteriorated based on adjoint point |
CN108717694A (en) * | 2018-04-24 | 2018-10-30 | 天津大学 | Electrical impedance tomography image quality evaluating method based on fuzzy C-means clustering |
CN109086810A (en) * | 2018-07-19 | 2018-12-25 | 天津大学 | A kind of electrical impedance tomography optimization algorithm based on fuzzy clustering |
CN109813772A (en) * | 2019-03-18 | 2019-05-28 | 燕山大学 | A Stable Imaging Method for Coplanar Array Capacitors |
Non-Patent Citations (3)
Title |
---|
QI LI 等: "The optimization of ET by FCM and Its Variants", 《2018 13TH WORLD CONGRESS ON INTELLIGENT CONTROL AND AUTOMATION (WCICA)》 * |
SHIHONGYUE 等: "Fuzzy clustering based ET image fusion", 《INFORMATION FUSION 14(2013)》 * |
ZEYING WANG 等: "An Unsupervised Method for Evaluating Electrical Impedance Tomography Images", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Efficient multitask structure-aware sparse Bayesian learning for frequency-difference electrical impedance tomography | |
CN112200306B (en) | A deep learning-based electrical impedance imaging method | |
Zhu et al. | Multicontext fuzzy clustering for separation of brain tissues in magnetic resonance images | |
CN109559359B (en) | Artifact removal method for reconstructed image from sparse angle data based on deep learning | |
CN109375125B (en) | Compressed sensing magnetic resonance imaging reconstruction method for correcting regularization parameters | |
CN111616708A (en) | An image reconstruction method for accurate identification of intracranial lesions in stroke | |
CN110613483A (en) | Method and system for detecting fetal craniocerebral abnormality based on machine learning | |
CN112465905A (en) | Characteristic brain region positioning method of magnetic resonance imaging data based on deep learning | |
CN110720915A (en) | Brain electrical impedance tomography method based on GAN | |
CN110298364B (en) | A multi-threshold-based feature selection method for functional brain networks | |
CN108960341B (en) | Brain network-oriented structural feature selection method | |
Venu | Performance and evalution of Guassian kernals for FCM algorithm with mean filtering based denoising for MRI segmentation | |
CN110910466B (en) | A Novel Reconstruction Algorithm for Multi-Frequency Differential Electrical Impedance Tomography | |
CN113052800B (en) | Alzheimer disease image analysis method and device | |
CN114052701A (en) | Capacitive coupling resistance tomography image reconstruction method | |
CN118521672A (en) | Electrical impedance image reconstruction method based on self-normalization neural network | |
CN110827215A (en) | ERT image reconstruction artifact removing method based on fuzzy clustering | |
CN115067920B (en) | An electrical impedance tomography method with high resolution capable of enhancing edge features of reconstructed images | |
CN115251889B (en) | A feature description method of dynamic connection network in fMRI images | |
CN109086810A (en) | A kind of electrical impedance tomography optimization algorithm based on fuzzy clustering | |
CN115670421A (en) | A Denoising Autoencoder Based Electrical Impedance Depth Imaging Method | |
CN111192209B (en) | Magnetic induction tomography optimization method based on deep learning | |
CN112150572A (en) | A method and device for suppressing image contact impedance artifacts for dynamic electrical impedance imaging | |
Yang et al. | Learning sparsifying transforms for image reconstruction in electrical impedance tomography | |
CN106354990B (en) | Method for detecting consistency of EEG and fMRI |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200221 |