CN107091858A - A kind of method mapped by Current Voltage map construction voltage x current - Google Patents
A kind of method mapped by Current Voltage map construction voltage x current Download PDFInfo
- Publication number
- CN107091858A CN107091858A CN201710187323.XA CN201710187323A CN107091858A CN 107091858 A CN107091858 A CN 107091858A CN 201710187323 A CN201710187323 A CN 201710187323A CN 107091858 A CN107091858 A CN 107091858A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- msup
- msubsup
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
Abstract
The present invention relates to a kind of method mapped by Current Voltage map construction voltage x current, for having the sensor of N number of electrode in electricity tomography, the parameters such as characteristic value and characteristic vector by calculating current voltage matrix, by the formula derived, voltage x current matrix uniquely is calculated, voltage x current mapping is constructed.The present invention gives the direct building method of voltage x current mapping, it can be applied to electricity chromatography imaging field and directly reconstruct in algorithm, its explicit physical meaning, and it is simple and easy to apply.
Description
Technical field
Imaging field is chromatographed the present invention relates to electricity, more particularly to one kind is reflected by current-voltage map construction voltage-to-current
The method penetrated.
Background technology
Electricity tomography (Electrical Tomography, abbreviation ET) technology is one kind of chromatography imaging technique.Its
Encouraged by applying to testee, and detect the change of its boundary value, rebuild using specific algorithm for reconstructing in measurand
The distribution of portion's electrical characteristic parameter, so as to obtain the distribution situation of interior of articles.Compared with other chromatography imaging techniques, electricity chromatography
Imaging is with the advantage such as radiationless, Noninvasive, portability, fast response time, cheap.
Electricity tomographic reconstruction algorithm generally can be divided into two classes:Algorithm for reconstructing based on sensitivity matrix and directly reconstruct
Algorithm.Using the former it is generally necessary to solve ill linear equation, this means that must be while rebuild all pixels of measured zone
Value.The latter realizes that the gray value of each pixel can be straight by calculating current-voltage mapping or voltage-to-current mapping
Connect, independent calculating.
The building method of voltage-to-current mapping is the important component that electricity tomography directly reconstructs algorithm.For N
The sensor of individual electrode, the data that voltage-to-current mapping is generally obtained by adjacent actuators measurement pattern are rebuild, because this
Measurement pattern is easier the realization on hardware, and the voltage-to-current map construction method based on adjacent actuators measurement pattern can give
Go out the direct physical significance of voltage-to-current mapping.But there is presently no construct voltage-to-current according to current-voltage mapping
The method of mapping.
The content of the invention
It is an object of the invention to propose a kind of direct structure that voltage-to-current mapping is constructed according to current-voltage mapping
Make method.
The technical scheme is that:
Step 1: the voltage on the electrode of the sensor with N number of electrode has general expression V=RAj, write as matrix
Form has:
V=ARJ (1)
Wherein J is current density vector, and A is the internal surface area of each electrode, and R is resistor matrix, and V is voltage vector, table
It is up to formula:
jsIt is the current density on s (1≤s≤N) individual electrode inner surface, VtIt is the magnitude of voltage on t-th of electrode,Table
Show that unitary current flows into tested field domain from i-th electrode, during j-th of electrode outflow, i-th and j electrode and s and t electrode
Between mutual impedance.
Then have:
As number of poles is N, and dielectric constant is distributed as current-voltage mapping matrix during ε (z).
It can be seen from Kirchhoff's law, loop current sum is zero, i.e., current density vector J order is N-1, therefore
Order:
Then:
WhereinAs number of poles is N, and dielectric constant is distributed as voltage-to-current mapping matrix during ε (z).
Definition:
So have:
BTSubscript T represent transposition to matrix, B can be provedTJ preceding N-1 characteristic value is -1, and last
Characteristic value is 0.By calculating BTJ characteristic value and characteristic vector, can be obtained:
BTJ=P { diag ([- 1-1 ...-1 0])N×N}PT (8)
Wherein diag ()N×NRepresent the diagonal matrix of N ranks, PTP=I, P=[p1 p2 … pN-1 pN].Defined feature to
Measure piIt is P i-th (1≤i≤N-1) row, i.e. BTJ·pi=(- 1) pi, pNIt is the corresponding characteristic vector of characteristic value 0, i.e. BTJ·pN
=(0) pN.Because average value often capable in equation (7) is all 0, thereforeSo BTJ can be with
It is written as:
Wherein IN×NIt is N rank unit matrixs.
Similar, by calculating characteristic value and characteristic vector,Order be N-1, can be written as following form:
∑R=diag ([λ1 λ2 … λN-1 0])N×NBe byN number of eigenvalue cluster into diagonal matrix, Q is by phase
The matrix of corresponding characteristic vector composition.QTQ=I, Q=[q1 q2 … qN-1 qN], characteristic vector qiRepresent Q i-th (1≤i≤
N-1) arrange, i.e.,qNIt is the corresponding characteristic vector of characteristic value 0, i.e.,
Therefore have:
Step 2: BTOrder be N, formula (8) both members simultaneously premultiplication (BT)-1:
J=(BT)-1P{diag([-1 -1 … -1 0])N×N}PT (12)
Formula (12) is substituted into formula (8) to obtain:
Both members multiply P, premultiplication B in the right side simultaneouslyT, then premultiplication PT:
If being apparent from matrix that a Matrix Multiplication obtains with diagonal matrix still for diagonal matrix, the matrix is also to angular moment
Battle array.Then understand:
Both members are while premultiplication P, the right side multiplies PT:
Formula (8) and formula (10) are substituted into formula (16) to obtain:
Both members are while premultiplication (BT)-1, the right side multiplies BT:
And BTJ=JBT, then:
Both members are while premultiplication QT, the right side multiplies Q and obtained:
Formula (12) is substituted into formula (20) to obtain:
Then have:
Both members are while premultiplication Q, the right side multiplies QT:
Thus formula is understood, matrixIt can be now uniquely determined.
Further, with reference to this method, any orthogonal set of current excitation pattern is shared in calculating voltage-to-current mapping.
Brief description of the drawings
Fig. 1 is implementing procedure figure.
Fig. 2 is embodiment isoboles.
Embodiment
Participate in Fig. 1, a kind of building method algorithm block diagram mapped according to current-voltage map construction voltage-to-current.To scheme
Exemplified by 16 end to end ring resistance networks shown in 2, illustrate the embodiment of this method.
It the described method comprises the following steps:
Step 1: the voltage on the electrode of the sensor with N=16 electrode has general expression V=RAj, write as
Matrix form has:
V=ARJ (24)
Wherein J is current density vector, and A is the internal surface area of each electrode, and R is resistor matrix, and V is voltage vector, table
It is up to formula:
jsIt is the current density on the individual electrode inner surfaces of s (1≤s≤16),Represent unitary current from i-th of electrode
Tested field domain is flowed into, during j-th of electrode outflow, the mutual impedance between i-th and j electrode and s and t electrode.
Then have:
As number of poles is 16, and dielectric constant is distributed as current-voltage mapping matrix during ε (z).
It can be seen from Kirchhoff's law, loop current sum is zero, i.e., current density vector J order is 15, therefore order:
Then:
WhereinAs number of poles is 16, and dielectric constant is distributed as voltage-to-current mapping matrix during ε (z).
Order:
So:
B can be provedTJ preceding 15 characteristic values are -1, and last characteristic value is 0.By calculating BTJ feature
Value and characteristic vector, can be obtained:
BTJ=P { diag [- 1-1 ...-1 0]16×16}PT (31)
Wherein diag ()16×16Represent the diagonal matrix of 16 ranks, PTP=I, P=[p1 p2 … p15 p16].Defined feature
Vectorial piIt is P i-th (1≤i≤15) row, i.e. BTJ·pi=(- 1) pi, p16It is the corresponding characteristic vector of characteristic value 0, i.e. BTJ·
p16=(0) p16.Because average value often capable in equation (30) is all 0, thereforeSo BTJ can
To be written as:
Wherein I16×16It is 16 rank unit matrixs.
Similar, by calculating characteristic value and characteristic vector,Order be 15, can be written as following form:
∑R=diag ([λ1 λ2 … λ15 0])16×16Be by16 eigenvalue clusters into diagonal matrix, Q be by
The matrix of corresponding characteristic vector composition.QTQ=I, Q=[q1 q2 … q15 q16], characteristic vector qiRepresent Q i-th (1≤
I≤15) arrange, i.e.,q16It is the corresponding characteristic vector of characteristic value 0, i.e.,
In ring resistance network, voltage matrix can be calculated:
It follows that current-voltage mapping matrix is:
Calculated according to formula (33):
Obtain:
Thus formula is understood, matrixIt can be now uniquely determined.
To prove this conclusion, following checking is done:Under neighboring modes, current density matrix is calculated:
I with setting is consistent, thus can verify that this method.
A kind of described method that voltage-to-current mapping is constructed by current-voltage, gives the one of voltage-to-current mapping
Computational methods are planted, this method explicit physical meaning is simple and easy to apply.With reference to this method, any orthogonal set of measurement pattern is encouraged
It may be used to calculate voltage-to-current mapping.
Above to the description of the present invention and embodiments thereof, it is not limited to which this, is only the reality of the present invention shown in accompanying drawing
Apply one of mode.Without departing from the spirit of the invention, it is similar with the technical scheme without designing with creating
Structure or embodiment, belong to the scope of the present invention.
Claims (1)
1. a kind of direct building method this method that voltage-to-current mapping is constructed according to current-voltage mapping includes following steps
Suddenly:
Step 1: the voltage on the electrode of the sensor with N number of electrode has general expression V=RAj, dielectric constant distribution
Current-voltage mapping matrix during for ε (z):
<mrow>
<msubsup>
<mi>R</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mi>A</mi>
<mi>R</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
It can be seen from Kirchhoff's law, loop current sum is zero, i.e., current density vector J order is N-1, therefore is made:
Then:
<mrow>
<mi>J</mi>
<mo>=</mo>
<msubsup>
<mi>&Lambda;</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<mi>V</mi>
<mo>=</mo>
<msubsup>
<mi>&Lambda;</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<mi>J</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
Definition:
So have:
Calculate BTJ characteristic value and characteristic vector, can be obtained:
BTJ=P { diag ([- 1-1 ...-1 0])N×N}PT (6)
Wherein diag ()N×NRepresent the diagonal matrix of N ranks, PTP=I, P=[p1 p2 … pN-1 pN], defined feature vector pi
It is P i-th (1≤i≤N-1) row, i.e. BTJ·pi=(- 1) pi, pNIt is the corresponding characteristic vector of characteristic value 0, i.e. BTJ·pN=
(0)pN,BTJ can be written as:
<mrow>
<msup>
<mi>B</mi>
<mi>T</mi>
</msup>
<mi>J</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>N</mi>
</msub>
<msubsup>
<mi>p</mi>
<mi>N</mi>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
Similar, by calculating characteristic value and characteristic vector,Order be N-1, can be written as following form:
<mrow>
<msubsup>
<mi>R</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>Q&Sigma;</mi>
<mi>R</mi>
</msub>
<msup>
<mi>Q</mi>
<mi>T</mi>
</msup>
<mo>=</mo>
<mi>Q</mi>
<mo>{</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msub>
<mrow>
<mo>(</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>)</mo>
</mrow>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>}</mo>
<msup>
<mi>Q</mi>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
∑R=diag ([λ1 λ2 … λN-1 0])N×NBe byN number of eigenvalue cluster into diagonal matrix, Q is by corresponding
The matrix of characteristic vector composition, QTQ=I, Q=[q1 q2 … qN-1 qN], characteristic vector qiQ i-th (1≤i≤N-1) row are represented,
I.e.qNIt is the corresponding characteristic vector of characteristic value 0, i.e.,
Therefore have:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<mi>Q</mi>
<mi>T</mi>
</msup>
<msup>
<mi>B</mi>
<mi>T</mi>
</msup>
<mi>J</mi>
<mi>Q</mi>
<mo>=</mo>
<msup>
<mi>Q</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>p</mi>
<mi>N</mi>
</msub>
<msubsup>
<mi>p</mi>
<mi>N</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>Q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>Q</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>N</mi>
</msub>
<msubsup>
<mi>q</mi>
<mi>N</mi>
<mi>T</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mi>Q</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msub>
<mrow>
<mo>(</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>)</mo>
</mrow>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
Step 2: BTOrder be N, formula (6) both members simultaneously premultiplication (BT)-1:
J=(BT)-1P{diag([-1 -1 … -1 0])N×N}PT (10)
Formula (10) is substituted into formula (3) to obtain:
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>P</mi>
<mo>{</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msub>
<mrow>
<mo>(</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>)</mo>
</mrow>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>}</mo>
<msup>
<mi>P</mi>
<mi>T</mi>
</msup>
<mo>=</mo>
<msubsup>
<mi>&Lambda;</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<msubsup>
<mi>R</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>B</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>P</mi>
<mo>{</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msub>
<mrow>
<mo>(</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>)</mo>
</mrow>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>}</mo>
<msup>
<mi>P</mi>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
So, voltage-to-current maps corresponding matrixIt can be written as:
<mrow>
<msubsup>
<mi>&Lambda;</mi>
<mi>&epsiv;</mi>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msubsup>
<mo>=</mo>
<mi>Q</mi>
<mo>{</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<msub>
<mrow>
<mo>(</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
</mfrac>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>)</mo>
</mrow>
<mrow>
<mi>N</mi>
<mo>&times;</mo>
<mi>N</mi>
</mrow>
</msub>
<mo>}</mo>
<msup>
<mi>Q</mi>
<mi>T</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
Thus formula is understood, matrixIt can be now uniquely determined.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710187323.XA CN107091858B (en) | 2017-03-27 | 2017-03-27 | Method for constructing voltage-current mapping from current-voltage mapping |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710187323.XA CN107091858B (en) | 2017-03-27 | 2017-03-27 | Method for constructing voltage-current mapping from current-voltage mapping |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107091858A true CN107091858A (en) | 2017-08-25 |
CN107091858B CN107091858B (en) | 2020-02-14 |
Family
ID=59646351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710187323.XA Active CN107091858B (en) | 2017-03-27 | 2017-03-27 | Method for constructing voltage-current mapping from current-voltage mapping |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107091858B (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005022138A1 (en) * | 2003-08-28 | 2005-03-10 | University Of Leeds | Eit data processing system and method |
CN101520478A (en) * | 2009-03-13 | 2009-09-02 | 北京航空航天大学 | Direct image reconstruction method based on capacitance tomography of round sensor |
CN103116101A (en) * | 2013-02-01 | 2013-05-22 | 北京航空航天大学 | Adjacent excitation measurement mode based voltage-current mapping construction method |
CN103149472A (en) * | 2013-02-01 | 2013-06-12 | 北京航空航天大学 | Voltage-current mapping constructing method based on double-terminal excitation measuring model |
CN103163404A (en) * | 2013-02-01 | 2013-06-19 | 北京航空航天大学 | Current-voltage mapping construction method based on adjacent stimulus measurement model |
CN103439375A (en) * | 2013-08-23 | 2013-12-11 | 华北电力大学 | Integrated type capacitance-ultrasound tomography sensor |
CN103776873A (en) * | 2014-02-18 | 2014-05-07 | 北京航空航天大学 | Method for constructing current-voltage mapping by virtue of voltage-current mapping |
CN105717145A (en) * | 2016-02-03 | 2016-06-29 | 北京航空航天大学 | Multi-mounted three-dimensional cone beam computer tomography method and device |
-
2017
- 2017-03-27 CN CN201710187323.XA patent/CN107091858B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005022138A1 (en) * | 2003-08-28 | 2005-03-10 | University Of Leeds | Eit data processing system and method |
CN101520478A (en) * | 2009-03-13 | 2009-09-02 | 北京航空航天大学 | Direct image reconstruction method based on capacitance tomography of round sensor |
CN103116101A (en) * | 2013-02-01 | 2013-05-22 | 北京航空航天大学 | Adjacent excitation measurement mode based voltage-current mapping construction method |
CN103149472A (en) * | 2013-02-01 | 2013-06-12 | 北京航空航天大学 | Voltage-current mapping constructing method based on double-terminal excitation measuring model |
CN103163404A (en) * | 2013-02-01 | 2013-06-19 | 北京航空航天大学 | Current-voltage mapping construction method based on adjacent stimulus measurement model |
CN103439375A (en) * | 2013-08-23 | 2013-12-11 | 华北电力大学 | Integrated type capacitance-ultrasound tomography sensor |
CN103439375B (en) * | 2013-08-23 | 2015-10-21 | 华北电力大学 | A kind of integrated capacitive-ultrasound tomography sensor |
CN103776873A (en) * | 2014-02-18 | 2014-05-07 | 北京航空航天大学 | Method for constructing current-voltage mapping by virtue of voltage-current mapping |
CN103776873B (en) * | 2014-02-18 | 2016-04-06 | 北京航空航天大学 | A kind of method mapped by voltage-to-current map construction current-voltage |
CN105717145A (en) * | 2016-02-03 | 2016-06-29 | 北京航空航天大学 | Multi-mounted three-dimensional cone beam computer tomography method and device |
Non-Patent Citations (7)
Title |
---|
FANG-MING YU 等: "A rotative electrical impedance tomography reconstruction system", 《JOURNAL OF PHYSICS:CONFERENCE SERIES》 * |
ZHOU HAILI 等: "Image reconstruction for invasive ERT in vertical oil well logging", 《CHINESE JOURNAL OF CHEMICAL ENGINEERING》 * |
张立峰 等: "一种修正的电阻层析成像Landweber迭代算法", 《计量学报》 * |
曹章: "基于EIT系统的模型、算法研究及低通滤波器优化设计", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
曹章: "电学层析成像系统模型、算法研究", 《中国博士学位论文全文数据库》 * |
王超 等: "基于L1-L1范数的电学层析成像静态成像算法", 《中国医疗设备》 * |
赵佳: "电阻层析成像稀疏重建方法研究", 《中国博士论文全文数据库》 * |
Also Published As
Publication number | Publication date |
---|---|
CN107091858B (en) | 2020-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Graham et al. | Objective selection of hyperparameter for EIT | |
Graham et al. | Electrode placement configurations for 3D EIT | |
CN109598768A (en) | Electricity tomographic image method for reconstructing based on convolutional neural networks | |
CN111281385B (en) | Electrical impedance imaging method based on tissue space distribution characteristics and impedance variation characteristics along with frequency | |
CN109919844A (en) | A kind of high-resolution electricity tomography distribution of conductivity method for reconstructing | |
CN107369187A (en) | The electricity tomography regularization reconstruction method for the sum that is deteriorated based on adjoint point | |
CN109934885A (en) | A kind of electrical resistance tomography image rebuilding method that sharpened edge is kept | |
Bera et al. | A MATLAB-based boundary data simulator for studying the resistivity reconstruction using neighbouring current pattern | |
CN109035352A (en) | L1-L2 spatially adaptive electricity tomography regularization reconstruction method | |
CN109646000A (en) | A kind of anti-imaging method of node-resistance based on the thin subdivision in part | |
CN109758149A (en) | A kind of electrical impedance imaging method of the imaging of the combination time difference and frequency difference imaging | |
CN105701824B (en) | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method | |
CN103065286A (en) | Image reconstruction method in quasi-static electrical impedance imaging | |
CN107091858A (en) | A kind of method mapped by Current Voltage map construction voltage x current | |
CN103149472B (en) | Voltage-current mapping constructing method based on double-terminal excitation measuring model | |
CN103163404B (en) | Current-voltage mapping construction method based on adjacent stimulus measurement model | |
Wtorek et al. | Impedance mammograph 3D phantom studies | |
CN107845119A (en) | A kind of electricity tomography mixed method | |
CN103116101B (en) | Adjacent excitation measurement mode based voltage-current mapping construction method | |
CN103776873B (en) | A kind of method mapped by voltage-to-current map construction current-voltage | |
Chen et al. | New post-processing method for interpretation of through casing resistivity (TCR) measurements | |
Xu et al. | A fast iterative updated thresholding algorithm with sparsity constrains for electrical resistance tomography | |
CN108646307B (en) | A kind of four-dimensional resistivity inversion methods based on dynamic adjusting data weighted value | |
Nejatali et al. | An iterative algorithm for electrical impedance imaging using neural networks | |
CN108888268B (en) | CT constrained electrical impedance medical imaging method |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |