CN110827215A - ERT image reconstruction artifact removing method based on fuzzy clustering - Google Patents

ERT image reconstruction artifact removing method based on fuzzy clustering Download PDF

Info

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
Application number
CN201910998118.0A
Other languages
Chinese (zh)
Inventor
李奇
岳士弘
高晓峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201910998118.0A priority Critical patent/CN110827215A/en
Publication of CN110827215A publication Critical patent/CN110827215A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/77Retouching; Inpainting; Scratch removal
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic 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

ERT image reconstruction artifact removing method based on fuzzy clustering
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):
Figure BDA0002240398200000022
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
Figure BDA0002240398200000031
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 say
Figure BDA0002240398200000032
The gray value of the middle pixel point is assigned asThe average gray value of the middle pixel point is recorded as
Figure BDA0002240398200000034
I.e. the final gray vector is:
Figure BDA0002240398200000035
6) using the formula (4)
Figure BDA0002240398200000036
Image inversion is performed as shown in fig. 3. The artifact region in fig. 3 is effectively reduced, closer to the original object field distribution.

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.
CN201910998118.0A 2019-10-21 2019-10-21 ERT image reconstruction artifact removing method based on fuzzy clustering Pending CN110827215A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (13)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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