CN105701824B - Electricity based on arc-shaped mesh generation chromatographs image rebuilding method - Google Patents
Electricity based on arc-shaped mesh generation chromatographs image rebuilding method Download PDFInfo
- Publication number
- CN105701824B CN105701824B CN201610024856.1A CN201610024856A CN105701824B CN 105701824 B CN105701824 B CN 105701824B CN 201610024856 A CN201610024856 A CN 201610024856A CN 105701824 B CN105701824 B CN 105701824B
- Authority
- CN
- China
- Prior art keywords
- radius
- subdivision
- arc
- sensitivity
- pixel
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- 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
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
The present invention relates to a kind of, and the electricity based on arc-shaped mesh generation chromatographs image rebuilding method, comprising: according to tested field domain, boundary survey value vector needed for obtaining reconstruction image, and determine the radius of round tested field domain；Radius the number M, M of the concentric circles number N and Equal round that determine subdivision circle field domain are the integral multiple of electricity imaging system number of poles；Obtain the specific coordinate of each pixel of tested region；Arc-shaped pixel based on institute's subdivision solves sensitivity coefficient of the one pair of electrodes relative to other electrodes pair according to conventional method, is rotated to obtain remaining sensitivity coefficient by acquired results, therefore can acquire sensitivity coefficient matrix A；According to resulting sensitivity matrix A above and boundary survey value vector b, reverse temperature intensity is carried out；It is imaged according to gained gray value is solved.Present invention reduces the computation complexities of sensitivity matrix, improve the speed of reconstruction image.
Description
Technical field
The invention belongs to field of electric measurement, and in particular to a kind of electricity chromatography image reconstruction based on arc-shaped mesh generation
Method.
Background technique
Electricity tomography (Electrical Tomography) technology is the one kind to grow up in Theory of Electromagnetic Field
Detection of Process Parameters technology can cause different boundary responses based on different electrical parameter distributions, to rebuild tested field domain
Substance distribution.Compared with other chromatography imaging techniques, electricity tomography has fast response time, radiationless and cheap
Etc. advantages, have tempting development prospect in commercial measurement and medical attendance.
The method for solving of electricity chromatography imaging problem generally uses numerical computation method, needs to carry out grid to tested field domain
Subdivision, to realize that the discretization of pixel calculates.General traditional method is the square net using N × N to tested circle
Field domain carries out subdivision, then removes the central point square net overseas in circular field, remaining is exactly that electricity tomography solves
Used in the process of the pixel that arrives, for example, typically carrying out subdivision to tested round field domain with 32 × 32 square net, finally
It obtains 812 central points and falls in pixel within the circle.Such mesh generation replaces circular boundary by straight line, to make weight
Building image boundary becomes rough, and fraction measurement field domain is not completely covered by the pixel of subdivision.Meanwhile such grid cuts open
Point method also lacks the flexibility of rotation, and then is unable to the solution sensitivity matrix of fast and flexible.
Summary of the invention
The purpose of the present invention is to provide a kind of, and the electricity based on arc-shaped mesh generation chromatographs image rebuilding method, the party
Method can make the pixel of subdivision that measured round field domain be completely covered, and the circular boundary for avoiding reconstruction image is rough
Line segment substituted, while the subdivision of this grid have rotation flexibility, thus by the computation complexity of sensitivity matrix
An order of magnitude is reduced, so that image reconstruction times are significantly faster than that traditional mesh generation mode.
A kind of electricity based on arc-shaped mesh generation provided by the present invention chromatographs image rebuilding method, and technical solution is such as
Under:
Step 1: according to tested field domain, boundary survey value vector b needed for obtaining reconstruction image, and determine round quilt
Survey the radius r of field domainIt surveys；
Step 2: determining the concentric circles number N of subdivision circle field domain and the radius number M of Equal round, then the picture of subdivision can be obtained
Prime number are as follows: N × M, wherein M be electricity imaging system number of poles integral multiple and N number of concentric circles radius must satisfy it is following
Formula:
Wherein, r1,r2,r3,...,rnIt is followed successively by the radius of concentric circles from inside to outside, rIt surveysFor the radius for being tested field domain；
Step 3: the specific coordinate of each pixel of tested region can be obtained, it is in i-th of concentric circles, the picture of j-th of equal subregion
The coordinate representation of element are as follows:
Step 4: the arc-shaped pixel based on institute's subdivision solves sensitivity of the one pair of electrodes relative to other electrodes pair
Coefficient is rotated to obtain remaining sensitivity coefficient by acquired results, acquires sensitivity coefficient matrix A；
Step 5: carrying out reverse temperature intensity according to resulting sensitivity matrix A above and boundary survey value vector b；
Step 6: being imaged according to gained gray value is solved.
Compared with the prior art, the advantages of the present invention are as follows: using the radius of N number of concentric circles and M Equal round to tested
Round field domain carries out subdivision, can make the pixel of subdivision that measured round field domain be completely covered, and reconstruction image boundary is presented
Smooth circle, while the subdivision of this grid has the flexibility of rotation, so that the calculating for reducing sensitivity matrix is complicated
Degree, so that image reconstruction times are significantly faster than that traditional mesh generation mode.
Detailed description of the invention
Fig. 1 is the subdivision result figure of the arc-shaped mesh generation that image reconstruction is chromatographed for electricity of the invention；
Fig. 2 is that electricity imaging system circular cross section is tested field domain and distribution of electrodes in the embodiment of the present invention；
Fig. 3 (a) is the object field actual distribution of measurement model in the embodiment of the present invention；
Fig. 3 (b) is the reconstruction image obtained in the embodiment of the present invention with traditional subdivision method；
Fig. 3 (c) is the reconstruction image obtained in the embodiment of the present invention with circular arc grid cutting algorithm.
In figure: 1, being tested field domain 2, electrode
Specific embodiment
The present invention, i.e., a kind of electricity based on arc-shaped mesh generation chromatograph image rebuilding method, including the following steps:
Step 1: according to tested field domain, boundary survey value vector b needed for obtaining reconstruction image, and determine round quilt
Survey the radius r of field domainIt surveys；
Step 2: determining the concentric circles number N of subdivision circle field domain and the radius number M of Equal round, then the picture of subdivision can be obtained
Prime number are as follows: N × M, wherein M be electricity imaging system number of poles integral multiple and N number of concentric circles radius must satisfy it is following
Formula:
Wherein, r1,r2,r3,...,rnIt is followed successively by the radius of concentric circles from inside to outside, rIt surveysFor the radius for being tested field domain；
Step 3: the specific coordinate of each pixel of tested region can be obtained, it is in i-th of concentric circles, the picture of j-th of equal subregion
The coordinate representation of element are as follows:
Meanwhile round equal part radius number M should be the integral multiple of electricity imaging system number of poles, can rotate solution in this way
Sensitivity coefficient matrix, such as (2,3) electrode can be by (1,2) electrodes to phase to the sensitivity coefficient relative to (4,5) electrode pair
The sensitivity coefficient of (3,4) electrode pair is directly rotated and is obtained, therefore need to only calculate one pair of electrodes relative to other electrodes
Pair sensitivity coefficient, thus remaining sensitivity coefficient can rotate to obtain.Then according to obtained sensitivity matrix and side
Boundary's measurement data can carry out image reconstruction to measured field domain.
Step 4: the arc-shaped pixel based on institute's subdivision solves one pair of electrodes relative to other electricity according to conventional method
Sensitivity coefficient extremely pair, is rotated to obtain remaining sensitivity coefficient by acquired results, therefore can acquire sensitivity coefficient matrix A；
Step 5: carrying out reverse temperature intensity according to resulting sensitivity matrix A above and boundary survey value vector b；
Step 6: being imaged according to gained gray value is solved.
The present invention will be further described with reference to the accompanying drawings and examples:
To a kind of electricity based on arc-shaped mesh generation provided by the present invention in Comsol and Matlab software for calculation
It learns tomographic map method for reconstructing and carries out emulation experiment.Experimental method is as follows:
For the electricity imaging system of typical 16 electrode, carry out 10 × 128 arc-shaped mesh generation, i.e., it is same with 10
(radius ratio is heart circle) and 128 Equal rounds radius to tested circle
Field domain carries out subdivision, obtains 1280 pixels, subdivision result is as shown in figure (1).16 electrode sensors are evenly distributed on round quilt
It surveys outside field domain, such as schemes shown in (2), using four terminal cycle motivation measurement methods of current excitation voltage measurement, can get 208
Boundary survey data.For the measurement model of Fig. 3 (a), with conventional mesh subdivision method and arc-shaped mesh generation of the invention
Method carries out inverting reconstruction image respectively, shown in the reconstruction image that traditional subdivision method obtains such as Fig. 3 (b), present invention circular arc
Shown in the reconstruction image that grid cutting algorithm obtains such as Fig. 3 (c).
The reconstruction image side that circular arc grid cutting algorithm provided by the present invention obtains it can be seen from the simulation experiment result
Smooth circle is presented in boundary, and measured round field domain is completely covered in the pixel of subdivision；Meanwhile image reconstruction institute in embodiment
Sensitivity matrix size is 208 × 1280, using circular arc grid cutting algorithm provided by the present invention, it need to only calculate 13 ×
1280 times, remaining can be obtained by rotating 16 expansions, the computation complexity of sensitivity matrix be reduced, so that image reconstruction
Time is significantly faster than that traditional mesh generation mode.
Emulation experiment demonstrates the good result of circular arc grid cutting algorithm provided by the present invention.
Above to the description of the present invention and embodiments thereof, it is not limited to which this, is only reality of the invention shown in attached drawing
Apply one of mode.Without departing from the spirit of the invention, similar with the technical solution without creatively designing
Structure or embodiment, category the scope of the present invention.
Claims (1)
1. a kind of electricity based on arc-shaped mesh generation chromatographs image rebuilding method, including the following steps:
Step 1: according to tested field domain, boundary survey value vector b needed for obtaining reconstruction image, and determine round tested field
The radius r in domainIt surveys；
Step 2: determining the concentric circles number N of subdivision circle field domain and the radius number M of Equal round, the pixel number of subdivision is obtained are as follows: N
× M, wherein M is the integral multiple of electricity imaging system number of poles and the radius of N number of concentric circles must satisfy following formula:
Wherein, r1,r2,r3,...,rnIt is followed successively by the radius of concentric circles from inside to outside；
Step 3: obtaining the specific coordinate of each pixel of tested region, it is in i-th of concentric circles, the seat of the pixel of j-th of equal subregion
Mark indicates are as follows:
Round equal part radius number M is the integral multiple of electricity imaging system number of poles, to rotate solution sensitivity coefficient matrix,
(2,3) electrode to the sensitivity coefficient relative to (4,5) electrode pair by (1,2) electrode to relative to the sensitive of (3,4) electrode pair
Degree coefficient, which directly rotates, to be obtained；
Step 4: the arc-shaped pixel based on institute's subdivision solves sensitivity system of the one pair of electrodes relative to other electrodes pair
Number, is rotated to obtain remaining sensitivity coefficient by acquired results, acquires sensitivity coefficient matrix A；
Step 5: carrying out reverse temperature intensity according to resulting sensitivity matrix A above and boundary survey value vector b；
Step 6: being imaged according to gained gray value is solved.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610024856.1A CN105701824B (en) | 2016-01-14 | 2016-01-14 | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610024856.1A CN105701824B (en) | 2016-01-14 | 2016-01-14 | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105701824A CN105701824A (en) | 2016-06-22 |
| CN105701824B true CN105701824B (en) | 2019-02-19 |
Family
ID=56227510
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610024856.1A Expired - Fee Related CN105701824B (en) | 2016-01-14 | 2016-01-14 | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105701824B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109086810A (en) * | 2018-07-19 | 2018-12-25 | 天津大学 | A kind of electrical impedance tomography optimization algorithm based on fuzzy clustering |
| CN110464329B (en) * | 2019-07-19 | 2020-06-02 | 东北大学 | Method and device for measuring blood flow velocity distribution |
| CN117649503B (en) * | 2024-01-29 | 2024-05-28 | 杭州永川科技有限公司 | Image reconstruction method, apparatus, computer device, storage medium, and program product |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1745713A (en) * | 2005-08-19 | 2006-03-15 | 清华大学 | X-ray CT multi-phase current inspection based on hereditary algorithm |
| CN102540276A (en) * | 2010-12-29 | 2012-07-04 | 通用电气公司 | Soft field tomography system and method |
| CN103868958A (en) * | 2014-03-27 | 2014-06-18 | 李星恕 | Electrical impedance tomography imaging plant root system architecture in-situ observation method |
-
2016
- 2016-01-14 CN CN201610024856.1A patent/CN105701824B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1745713A (en) * | 2005-08-19 | 2006-03-15 | 清华大学 | X-ray CT multi-phase current inspection based on hereditary algorithm |
| CN102540276A (en) * | 2010-12-29 | 2012-07-04 | 通用电气公司 | Soft field tomography system and method |
| CN103868958A (en) * | 2014-03-27 | 2014-06-18 | 李星恕 | Electrical impedance tomography imaging plant root system architecture in-situ observation method |
Non-Patent Citations (2)
| Title |
|---|
| 多相流电容层析成像系统的仿真研究;张立峰 等;《电力科学与工程》;20051231(第1期);5-7 |
| 电容层析成像系统阵列电极的优化设计;王化祥 等;《天津大学学报》;20030531;第36卷(第3期);307-310 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105701824A (en) | 2016-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109919844A (en) | A kind of high-resolution electricity tomography distribution of conductivity method for reconstructing | |
| CN105701824B (en) | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method | |
| CN102688041B (en) | Three-dimensional electrical impedance tomography method based on crisscross-arranged electrodes | |
| CN101794453B (en) | Reconstruction method of node mapping image based on regression analysis | |
| CN106530367A (en) | Electrical tomography sparse reconstruction method based on Firm threshold iteration | |
| Liu et al. | Time-difference imaging of magnetic induction tomography in a three-layer brain physical phantom | |
| Wang et al. | Image reconstruction algorithm for capacitively coupled electrical resistance tomography | |
| Dong et al. | Image reconstruction for electrical capacitance tomography by using soft-thresholding iterative method with adaptive regulation parameter | |
| CN103040466A (en) | Abnormal electrode connection detecting method for impedance detection | |
| CN109646000B (en) | Node electrical impedance imaging method based on local subdivision | |
| CN103868958A (en) | Electrical impedance tomography imaging plant root system architecture in-situ observation method | |
| CN110852025A (en) | Three-dimensional electromagnetic slow diffusion numerical simulation method based on hyperconvergence interpolation approximation | |
| Kang et al. | A sub-domain based regularization method with prior information for human thorax imaging using electrical impedance tomography | |
| GB2486967A (en) | Soft Field Tomography Iteration Method | |
| Tan et al. | New image reconstruction algorithm for capacitively coupled electrical resistance tomography | |
| CN109758149B (en) | Electrical impedance imaging method combining time difference imaging and frequency difference imaging | |
| Ren et al. | Tomographic wire-mesh imaging of water-air flow based on sparse minimization | |
| CN104574462B (en) | A kind of improvement L-curve electricity tomographic reconstruction method based on curvature estimation | |
| CN104535294B (en) | Corrected L-curve electrical tomography reconstruction method based on second-order differential | |
| CN109118553A (en) | Electrical impedance tomography content Boundary Reconstruction method based on geometric constraints | |
| CN105701825B (en) | The arc-shaped grid cutting algorithm solved for electricity chromatography image sensitivity coefficient | |
| CN110501587A (en) | A kind of Radon inverse transformation capacitance chromatography imaging method based on electric force lines distribution | |
| CN103630750A (en) | Electrical impedance tomography-based gel conductivity measurement method | |
| Xu et al. | A new regularization algorithm based on the neighborhood method for electrical impedance tomography | |
| Cao et al. | 2D image reconstruction of a human chest by using Calderon's method and the adjacent current pattern |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190219 Termination date: 20220114 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |