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
- gray
- image reconstruction
- cluster clusters
- fuzzy clustering
- method based
- 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 15
- 239000013598 vector Substances 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims abstract description 13
- 230000005284 excitation Effects 0.000 description 6
- 238000003325 tomography Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 238000007621 cluster analysis Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 230000005855 radiation 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
The invention relates to an ERT image reconstruction artifact removing method based on fuzzy clustering, which comprises the following steps: 1) respectively calculating the gray level vector of the measured field by utilizing at least two image reconstruction algorithms; 2) combining the gray vectors in the step 1) to form a gray matrix as a characteristic vector in a clustering algorithm; 3) applying a fuzzy clustering algorithm to the gray matrix, and performing clustering analysis on the gray matrix; 4) finding out cluster clusters respectively representing a target, a background and an artifact according to the statistical characteristics expressed by the cluster clusters; 5) cluster clusters representing artifacts are classified into cluster clusters representing background.
Description
Technical Field
The invention belongs to the field of electrical tomography image reconstruction algorithms, and particularly relates to an ERT image reconstruction artifact removing method based on fuzzy clustering.
Background
An Electrical Tomography (ET) technique is a technique for reconstructing the internal Electrical property distribution of a measured object field, and the technique adopts simple hardware equipment structure and low cost, and the detection process has the characteristics of no radiation and non-invasiveness. The technology has wide application prospect in the fields of medical clinical monitoring, industrial measurement and the like.
In Electrical Resistance Tomography (ERT), current excitation is applied to a measured object field by using array electrodes, boundary voltage data are measured, and conductivity distribution in a measured field domain is reconstructed. Taking a 16-electrode electrical resistance tomography sensor as an example, a method of adjacent current excitation and adjacent voltage measurement is adopted. At each excitation, the excitation electrode was removed for a total of 13 measurements. Such as: the voltage signals of the electrode pairs 3, 4, 5, 6, …, 15 and 16 are measured by taking the No. 1 and No. 2 electrodes as excitation electrodes. There are 16 excitations in total, i.e. 13 × 16 — 208 measurements. In order to achieve both temporal resolution and spatial resolution, the sensitive field to be measured is generally divided into 812 pixels. The problem of reconstructing the electrical resistance tomography image is to solve the gray value of the 812 pixel points.
The spatial resolution of the reconstructed image is poor due to the "soft-field" effect of ERT and the inherent ill-qualification of the inverse problem. The traditional linear back projection algorithm (LBP) is to perform back projection of voltage data according to the relationship between the boundary voltage value and the projection area of the measured field, but due to the limitation of hardware equipment, the back projection data is far less than the gray value to be solved, so that the reconstructed image has larger artifacts, the boundary of the target and the background is fuzzy, and the spatial resolution is poor.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides an ERT image reconstruction artifact removing method based on fuzzy clustering, which can effectively reduce the influence of artifacts on image resolution. All pixel points of the reconstructed image can be classified into 3 types: objects, background, and artifacts. The artifacts influence the imaging resolution to a great extent, and the quality of the reconstructed image can be improved to a certain extent by effectively processing the artifacts. Considering that the result obtained by the above algorithm is a gray value, i.e. numerical data, the gray value can be divided by using a clustering algorithm, and then the gray value is further processed according to the division result, so as to obtain a more ideal imaging effect.
The technical scheme of the invention is as follows:
an ERT image reconstruction artifact removing method based on fuzzy clustering comprises the following steps:
1) respectively calculating the gray level vector of the measured field by utilizing at least two image reconstruction algorithms;
2) combining the gray vectors in the step 1) to form a gray matrix as a characteristic vector in a clustering algorithm;
3) applying a fuzzy clustering algorithm to the gray matrix, and performing clustering analysis on the gray matrix;
4) finding out cluster clusters respectively representing a target, a background and an artifact according to the statistical characteristics expressed by the cluster clusters;
5) cluster clusters representing artifacts are classified into cluster clusters representing background.
Due to the adoption of the technical scheme, the invention has the following advantages:
1. compared with the traditional electrical tomography algorithm, the algorithm has the advantages that a plurality of classical imaging algorithms are fused, the advantages of the classical imaging algorithms are complemented, the gray value is further subjected to cluster analysis on the basis, the influence of artifacts is reduced, the boundary between the target and the background is clearer, and the spatial resolution of imaging is improved.
2. The clustering algorithm adopts a fuzzy clustering algorithm instead of a traditional hard division algorithm, so that the target edge is smoother and closer to the actual edge.
Drawings
FIG. 1 is a diagram of an original object field model in an embodiment of the invention;
FIG. 2 is a result of image reconstruction using the LBP algorithm for the simulation model shown in FIG. 1;
FIG. 3 is a result of image reconstruction using the algorithm of the present invention for the simulation model shown in FIG. 1.
Detailed Description
The invention aims to overcome the defect of low spatial resolution of the conventional image reconstruction algorithm, and provides an ERT image reconstruction artifact removing method based on fuzzy clustering. The fuzzy clustering object is a set of data, and the data are divided by extracting the characteristic vector of the data, so that the intra-class difference is small, and the inter-class difference is large. By carrying out cluster analysis on the gray matrix, the influence of the artifact on the image resolution is reduced.
The method for removing the ERT image reconstruction artifact based on fuzzy clustering according to the present invention is described in detail below with reference to the accompanying drawings and embodiments. The method comprises the following specific steps:
1) fig. 1 shows a flow pattern containing a discrete target medium, and for the measured object field shown in fig. 1, if the regularization (TK) algorithm and LBP algorithm are used to calculate the gray vector, the gray value of the reconstructed image can be shown in formula (1).
GT={gt1,gt2,...,gt812},GL={gl1,gl2,...,gl812} (1)
Wherein G isTRepresenting the gray vector, G, obtained by the TK algorithmLRepresenting the gray vector resulting from the LBP algorithm.
An image reconstructed by the LBP algorithm is shown in fig. 2, an area with a small upper circular-like area represents a target, an area with a largest lower area is a background area, and a gray area between the upper circular-like area and the lower circular-like area is an artifact area.
2) Combining the gray level vectors in 1) to form a gray level matrix as the feature vector in the cluster, i.e. combining the gray level vectors GT,GLForming a gray matrix and calculating a gray mean vectorAs shown in equation (2):
3) and (3) carrying out clustering analysis on the gray matrix in the step 2) by using a fuzzy clustering algorithm. The fuzzy clustering algorithm needs to give a clustering number, and as can be seen from the meaning of the inverse problem, the clustering number can be set to 3, that is, data is divided into three types: target GoBackground GbAnd an artifact GtI.e. by
4) And calculating the statistical characteristics of each cluster. The 3 clusters in step 3) have different statistical characteristics: the class with the largest mean is the target, in contrast to the class with the smallest mean being the background, the class with the largest variance represents the artifact. According to the characteristics, the meaning of each cluster in the clustering result can be clarified.
5) Assign clusters representing artifacts to clusters representing background, i.e. to sayThe gray value of the middle pixel point is assigned asThe average gray value of the middle pixel point is recorded asI.e. the final gray vector is:
Claims (1)
1. An ERT image reconstruction artifact removing method based on fuzzy clustering comprises the following steps:
1) and respectively calculating the gray level vector of the measured field by utilizing at least two image reconstruction algorithms.
2) Combining the gray vectors in the step 1) to form a gray matrix as a characteristic vector in a clustering algorithm;
3) applying a fuzzy clustering algorithm to the gray matrix, and performing clustering analysis on the gray matrix;
4) finding out cluster clusters respectively representing a target, a background and an artifact according to the statistical characteristics expressed by the cluster clusters;
5) cluster clusters representing artifacts are classified into cluster clusters representing background.
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 | 西安电子科技大学 | Multi-agent genetic clustering algorithm-based image segmentation method |
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 kind of stable imaging method of coplanar array capacitor |
-
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 | 西安电子科技大学 | Multi-agent genetic clustering algorithm-based image segmentation method |
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 kind of stable imaging method of coplanar array capacitor |
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 |
---|---|---|
Collewet et al. | Influence of MRI acquisition protocols and image intensity normalization methods on texture classification | |
Zhu et al. | Multicontext fuzzy clustering for separation of brain tissues in magnetic resonance images | |
CN112200306B (en) | Electrical impedance imaging method based on deep learning | |
CN111047663B (en) | Method for reconstructing electrical tomography artifact inhibition image | |
Yang et al. | Image reconstruction for electrical impedance tomography using enhanced adaptive group sparsity with total variation | |
CN111616708A (en) | Image reconstruction method for accurately identifying cerebral apoplexy intracranial lesion area | |
CN108717694B (en) | Electrical impedance tomography image quality evaluation method based on fuzzy C-means clustering | |
CN111047662A (en) | Self-adaptive non-convex mixing total variation regularization industrial resistance tomography method | |
CN108830875A (en) | One kind being based on the smallest electrical impedance tomography image partition method of residual error | |
CN110910466B (en) | Novel multi-frequency differential electrical impedance tomography reconstruction algorithm | |
CN112798654A (en) | Rapid gradient method and self-adaptive Jacobian matrix reconstruction method for electrical impedance tomography | |
Pan et al. | An optimization method of planar array capacitance imaging | |
CN110208605A (en) | A method of inhibit the electrical resistance tomography distribution of conductivity of alias to rebuild | |
Liu et al. | Noise suppression in brain magnetic resonance imaging based on non-local means filter and fuzzy cluster | |
CN114052701A (en) | Capacitive coupling resistance tomography image reconstruction method | |
Feng et al. | Multimodal medical image fusion based on visual saliency map and multichannel dynamic threshold neural P systems in sub-window variance filter domain | |
Godtliebsen et al. | A visual display device for significant features in complicated signals | |
CN110827215A (en) | ERT image reconstruction artifact removing method based on fuzzy clustering | |
Viani et al. | Where Bayes tweaks Gauss: Conditionally Gaussian priors for stable multi-dipole estimation | |
CN108122229B (en) | Full-reference electrical impedance imaging image quality evaluation method | |
CN113034635B (en) | Image reconstruction method for inhibiting industrial imaging step artifacts | |
Aishwarya et al. | An image fusion framework using morphology and sparse representation | |
Shams et al. | Automated iterative reclustering framework for determining hierarchical functional networks in resting state f MRI | |
CN108175407B (en) | Method for selecting local optimal regularization parameters of craniocerebral EIT | |
CN113034632A (en) | Image reconstruction method for detecting industrial two-phase flow |
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 |
Application publication date: 20200221 |
|
RJ01 | Rejection of invention patent application after publication |