CN103868958A - Electrical impedance tomography imaging plant root system architecture in-situ observation method - Google Patents
Electrical impedance tomography imaging plant root system architecture in-situ observation method Download PDFInfo
- Publication number
- CN103868958A CN103868958A CN201410117509.4A CN201410117509A CN103868958A CN 103868958 A CN103868958 A CN 103868958A CN 201410117509 A CN201410117509 A CN 201410117509A CN 103868958 A CN103868958 A CN 103868958A
- Authority
- CN
- China
- Prior art keywords
- resistivity
- unit
- node
- matrix
- finite element
- 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
Images
Abstract
The invention relates to an electrical impedance tomography imaging plant root system architecture in-situ observation method which comprises the steps of 1, evenly distributing N electrodes in the soil containing a plant root system in a circular shape, sequentially applying sinusoidal alternating current (AC) on two of the electrodes, and measuring the voltage between every two adjacent electrodes in the rest electrodes to obtain voltage data; 2, carrying out finite element mesh partition on the cross section of circumference provided with the electrodes in the soil; 3, carrying out system matrix calculation on the voltage data obtained in the step 1 to obtain the resistivities of all the grids obtained in the step 2, wherein each resistivity represents the resistivity of the corresponding position in the solving domain; imaging by utilizing the resistivity data, setting different colors for units with different resistivities according to the resistivity values of the units by a red, green and blue (RGB) color mode, and judging the different resistivities of all the units according to the colors of the image; according to the difference between the root system and the surrounding soil resistivity, determining the shape and the position of the root system to realize the plant root system architecture in-situ observation.
Description
Technical field
The present invention relates to Plant Root Architecture identification field, particularly a kind of electrical impedance tomography Plant Root Architecture in-situ observation method.
Background technology
Existingly mainly contain for realizing Plant Root Architecture in-situ observation identification institute employing technology: mini-rhizotron method, ground penetrating radar (GPR) technology, ray computer tomography (X-CT) technology, Magnetic resonance imaging (MRI) technology.
But common being limited as of these several observation technologies: complete one-shot measurement consuming time long and cannot realize root system of plant Real-Time Monitoring, can only can not distinguish for assessment of root biomass the health status of root system.
They also have separately specifically limitation simultaneously: can not carry out omnibearing observation survey to whole root system configuration as mini-rhizotron method, only can realize the two dimension to root system certain depth, local observation; The radar signal parameter of ground penetrating radar itself is subject to the impact of surrounding environment factor, thereby causes its measuring accuracy poor; Ray computer tomography (imaging with measurand density about) and Magnetic resonance imaging (water cut of imaging and measurand about) thus in have the density of surrounding environment or the water cut intensity profile of image and the overlapped post-processed difficulty that causes of the intensity profile of its surrounding medium when close when root system itself, and their use nucleic or high-intensity magnetic field to have certain impact, need that special messenger operates, detection system is huge and expensive unsuitable outdoor use plant growth.
Summary of the invention
The problem existing for solving above-mentioned prior art, the present invention proposes a kind of electrical impedance tomography Plant Root Architecture in-situ observation method, have equipment simple portable easy to operate, cost is low, the outer interference performance of anti-field domain is strong, image taking speed fast, can realize Real-Time Monitoring, plant health is not formed harm, do not damaging or moving under vegeto-animal prerequisite, function of dominant imaging can provide the feature of the biochemical information that reflects health of root situation.
For achieving the above object, technical scheme of the present invention is:
A kind of electrical impedance tomography Plant Root Architecture in-situ observation method, the method comprises the steps:
Step 1, on the soil that contains root system of plant, be circumference type and be evenly arranged N electrode, on two electrodes, apply sinusoidal ac therein successively, measure the voltage between other adjacent electrodes, obtain voltage data;
Step 2, soil is arranged to the circumferential section of electrode carries out finite element grid division;
Step 3, step 1 gained voltage data is obtained to the resistivity of each grid in step 2 gained grid by system matrix computation, each resistivity reflect the resistivity that solves territory correspondence position, utilize this resistivity data to carry out imaging, then adopt rgb color pattern, according to unit the height of corresponding resistivity, different resistivity unit is arranged to different colours, judge the difference of each cell resistance rate by the color of image; Owing to there are differences between root system and surrounding soil resistivity, and then can determine shape and the position of root system, realize Plant Root Architecture in-situ observation.
Further, in above-mentioned steps one, it is 16 according to the number of electrodes of circumference.
Further, in above-mentioned steps two, grid is divided into, and uses node connection method, is close outside interior dredging by two-dimentional field domain subdivision, and inner sparse 7 layers, the set of 1040 little triangular elements of outside dense 3 layers, the concrete steps that finite element grid is divided are as follows:
(1), determine the each parameter of subdivision field domain: sparse number of turns N1=7, dense number of turns N2=3, section radius R=experimental tank radius, unit number e=1040, nodes en=585, and node, unit topology information defined matrix;
(2), determine that the topology information of each unit, the topology information of unit comprise each summit, this unit numbering, this unit unit information on limit altogether;
(3), determine the coordinate of each node;
(4), determine the topology information of each node, node topology information comprises the unit information of the connected node of this node, public this node;
(5), determine between each node and whether exist directly and be connected and build matrix of coefficients;
(6), realize the drafting of grid;
(7), realize the demonstration of node and element number.
Further, in above-mentioned steps three, the concrete steps of system matrix computing method are:
(1) establishment of finite element unary system matrix number
Appoint and get a triangular element e who solves in the finite element grid of territory, its coordinate is
1(x
i, y
i)
2(x
j, y
j)
3(x
k, y
k), its area is S
e, its initial electrical resistivity is ε
e; The matrix of coefficients of this unit is:
Wherein
(2) establishment of system matrix
System matrix [K] is by the matrix of coefficients [K] of the finite element unit of field domain
eform, system matrix [K] is e
n× e
nrank matrix, wherein, e
nfor the contained nodes in finite element solving territory, the contained each element representation of system matrix is:
Acting as of system matrix: in conjunction with boundary condition, i.e. injecting voltage, utilizes Laplce's finite element equation, [K] [Φ]=0, wherein [Φ] is made up of the each finite element node potential of field domain, obtains the current potential of each finite element node in field domain.
Further, in above-mentioned steps three, the algorithm of imaging comprises the steps: first according to providing arbitrarily the resistivity distribution ρ that solves territory
k, in conjunction with selected excitation and system matrix [K] and then calculate the multi-electrode experiment container boundary voltage V that contains root system of plant
j, utilize and calculate the gained boundary voltage boundary voltage V measured with passing through measuring system
cbuild reconstructed object function, to original resistivity distribution ρ
krevise and obtain new resistivity distribution ρ
k+1; Repeatedly according to new resistivity distribution ρ
k+1solve new boundary voltage V
jby its with record boundary voltage V
crelatively, until the boundary voltage V of gained
jwith record boundary voltage V
cbetween error while being less than fiducial error ε till, last according to new resistivity distribution ρ
k+1rebuild drafting the demonstration of image through row.
Further, in above-mentioned steps three, resistivity cell colors method to set up is: choose the unit that resistivity is the highest and be set to peony: RGB=(255 0 0), choose the unit that resistivity is minimum and be set to mazarine: RGB=(0 0 255), according to unit the height of corresponding resistivity, its color to mazarine gradual change, can by the color of image be judged the difference of each cell resistance rate with this by peony.
With respect to prior art, beneficial effect of the present invention is: the present invention is according to the electrical characteristic parameter difference between root system of plant and surrounding soil, by it is applied to safe exciting current or voltage around, measure the voltage of body surface or current signal simultaneously and know the distribution of interior of articles electrical characteristic parameter, and then reconstruct reflection internal structure of body image, present device is simple portable easy to operate, cost is low, the outer interference performance of anti-field domain is strong, image taking speed is fast, can realize Real-Time Monitoring, plant health is not formed to harm, do not damaging or moving under vegeto-animal prerequisite, functional imaging can provide the biochemical information of root system, measure and realize plant roots in situ mensuration by curtage, there is significant progress, significant in actual production and detection.
Brief description of the drawings
Fig. 1 is imaging algorithm block diagram in the present invention.
Fig. 2 is that finite element grid of the present invention is divided design sketch.
Fig. 3 is root system of plant fault imaging design sketch of the present invention.
Fig. 4 is hardware configuration schematic diagram of the present invention.
Embodiment
Below in conjunction with the drawings and the specific embodiments, the present invention program is described in further details:
A kind of electrical impedance tomography Plant Root Architecture in-situ observation method, adopt hardware as shown in Figure 4, its Computer provides Matlab experiment porch, LCR impedance instrument to have 16 electrodes to contact with experimental subjects for experimental system provides on the excitation selection that also measuring resistance is anti-, multipath high-speed switch control excitation applies and voltage (electric current) potential electrode is right, experiment container.Between computing machine, LCR impedance instrument, multipath high-speed switch three, be connected and communicate by USB/GPIB interface; And input with output line and be connected by signal between LCR impedance instrument and multipath high-speed switch; Multipath high-speed switch is connected with multi-electrode measuring vessel by wire.The method comprises the steps:
Step 1, on the soil that contains root system of plant, be circumference type and be evenly arranged N electrode, on two electrodes, apply sinusoidal ac therein successively, measure the voltage between other adjacent electrodes, obtain voltage data;
Step 2, soil is arranged to the circumferential section of electrode carries out finite element grid division;
Step 3, step 1 gained voltage data is obtained to the resistivity of each grid in step 2 gained grid by system matrix computation, each resistivity reflect the resistivity that solves territory correspondence position, utilize this resistivity data to carry out imaging, then adopt rgb color pattern, according to unit the height of corresponding resistivity, different resistivity unit is arranged to different colours, judge the difference of each cell resistance rate by the color of image; Owing to there are differences between root system and surrounding soil resistivity, and then can determine shape and the position of root system, realize Plant Root Architecture in-situ observation.
Further, in step 1, it is 16 according to the number of electrodes of circumference.
As shown in Figure 2, further, in above-mentioned steps two, grid is divided into, and uses node connection method, is close outside interior dredging by two-dimentional field domain subdivision, inner sparse 7 layers, the set of 1040 little triangular elements of outside dense 3 layers, the concrete steps that finite element grid is divided are as follows:
(1), determine the each parameter of subdivision field domain: sparse number of turns N1=7, dense number of turns N2=3, section radius R=experimental tank radius, unit number e=1040, nodes en=585, and node, unit topology information defined matrix;
(2), determine that the topology information of each unit, the topology information of unit comprise each summit, this unit numbering, this unit unit information on limit altogether;
(3), determine the coordinate of each node;
(4), determine the topology information of each node, node topology information comprises the unit information of the connected node of this node, public this node;
(5), determine between each node and whether exist directly and be connected and build matrix of coefficients;
(6), realize the drafting of grid;
(7), realize the demonstration of node and element number.
Further, in above-mentioned steps three, the concrete steps of system matrix computing method are:
(1) establishment of finite element unary system matrix number
Appoint and get a triangular element e who solves in the finite element grid of territory, its coordinate is
1(x
i, y
i)
2(x
j, y
j)
3(x
k, y
k), its area is S
e, its initial electrical resistivity is ε
e; The matrix of coefficients of this unit is:
Wherein
(2) establishment of system matrix
System matrix [K] is by the matrix of coefficients [K] of the finite element unit of field domain
eform, system matrix [K] is e
n× e
nrank matrix, wherein, e
nfor the contained nodes in finite element solving territory, the contained each element representation of system matrix is:
Acting as of system matrix: in conjunction with boundary condition, i.e. injecting voltage, utilizes Laplce's finite element equation, [K] [Φ]=0, wherein [Φ] is made up of the each finite element node potential of field domain, obtains the current potential of each finite element node in field domain.
Further, as shown in Figure 1, in above-mentioned steps three, the algorithm of imaging comprises the steps: first according to providing arbitrarily the resistivity distribution ρ that solves territory
k, in conjunction with selected excitation and system matrix [K] and then calculate the multi-electrode experiment container boundary voltage V that contains root system of plant
j, utilize and calculate the gained boundary voltage boundary voltage V measured with passing through measuring system
cbuild reconstructed object function, to original resistivity distribution ρ
krevise and obtain new resistivity distribution ρ
k+1; Repeatedly according to new resistivity distribution ρ
k+1solve new boundary voltage V
jby its with record boundary voltage V
crelatively, until the boundary voltage V of gained
jwith record boundary voltage V
cbetween error while being less than fiducial error ε till, last according to new resistivity distribution ρ
k+1rebuild drafting the demonstration of image through row.
Further, in above-mentioned steps three, resistivity cell colors method to set up is: choose the unit that resistivity is the highest and be set to peony: RGB=(255 0 0), choose the unit that resistivity is minimum and be set to mazarine: RGB=(0 0 255), according to unit the height of corresponding resistivity, its color to mazarine gradual change, can by the color of image be judged the difference of each cell resistance rate with this by peony.Income effect as shown in Figure 3.
The above, be only the specific embodiment of the present invention, but protection scope of the present invention is not limited to this, and any variation of expecting without creative work or replacement, within all should being encompassed in protection scope of the present invention.Therefore, protection scope of the present invention should be as the criterion with the protection domain that claims were limited.
Claims (6)
1. an electrical impedance tomography Plant Root Architecture in-situ observation method, is characterized in that, the method comprises the steps:
Step 1, on the soil that contains root system of plant, be circumference type and be evenly arranged N electrode, on two electrodes, apply sinusoidal ac therein successively, measure the voltage between other adjacent electrodes, obtain voltage data;
Step 2, soil is arranged to the circumferential section of electrode carries out finite element grid division;
Step 3, step 1 gained voltage data is obtained to the resistivity of each grid in step 2 gained grid by system matrix computation, each resistivity reflect the resistivity that solves territory correspondence position, utilize this resistivity data to carry out imaging, then adopt rgb color pattern, according to unit the height of corresponding resistivity, different resistivity unit is arranged to different colours, judge the difference of each cell resistance rate by the color of image; Owing to there are differences between root system and surrounding soil resistivity, and then can determine shape and the position of root system, realize Plant Root Architecture in-situ observation.
2. method of testing as claimed in claim 1, is characterized in that, in step 1, is 16 according to the number of electrodes of circumference.
3. method of testing as claimed in claim 2, is characterized in that, in step 2, grid is divided into, and uses node connection method, is close outside interior dredging by two-dimentional field domain subdivision, inner sparse 7 layers, the set of 1040 little triangular elements of outside dense 3 layers, the concrete steps that finite element grid is divided are as follows:
(1), determine the each parameter of subdivision field domain: sparse number of turns N1=7, dense number of turns N2=3, section radius R=experimental tank radius, unit number e=1040, nodes en=585, and node, unit topology information defined matrix;
(2), determine that the topology information of each unit, the topology information of unit comprise each summit, this unit numbering, this unit unit information on limit altogether;
(3), determine the coordinate of each node;
(4), determine the topology information of each node, node topology information comprises the unit information of the connected node of this node, public this node;
(5), determine between each node and whether exist directly and be connected and build matrix of coefficients;
(6), realize the drafting of grid;
(7), realize the demonstration of node and element number.
4. method of testing as claimed in claim 3, is characterized in that, in step 3, the concrete steps of system matrix computing method are:
(1) establishment of finite element unary system matrix number
Appoint and get a triangular element e who solves in the finite element grid of territory, its coordinate is
1(x
i, y
i)
2(x
j, y
j)
3(x
k, y
k), its area is S
e, its initial electrical resistivity is ε
e; The matrix of coefficients of this unit is:
Wherein
(2) establishment of system matrix
System matrix [K] is by the matrix of coefficients [K] of the finite element unit of field domain
eform, system matrix [K] is e
n× e
nrank matrix, wherein, e
nfor the contained nodes in finite element solving territory, the contained each element representation of system matrix is:
all unit, limit e matrix of coefficients and;
Acting as of system matrix: in conjunction with boundary condition, i.e. injecting voltage, utilizes Laplce's finite element equation, [K] [Φ]=0, wherein [Φ] is made up of the each finite element node potential of field domain, obtains the current potential of each finite element node in field domain.
5. method of testing as claimed in claim 4, is characterized in that, in step 3, the algorithm of imaging comprises the steps: first according to providing arbitrarily the resistivity distribution ρ that solves territory
k, in conjunction with selected excitation and system matrix [K] and then calculate the multi-electrode experiment container boundary voltage V that contains root system of plant
j, utilize and calculate the gained boundary voltage boundary voltage V measured with passing through measuring system
cbuild reconstructed object function, to original resistivity distribution ρ
krevise and obtain new resistivity distribution ρ
k+1; Repeatedly according to new resistivity distribution ρ
k+1solve new boundary voltage V
jby its with record boundary voltage V
crelatively, until the boundary voltage V of gained
jwith record boundary voltage V
cbetween error while being less than fiducial error ε till, last according to new resistivity distribution ρ
k+1rebuild drafting the demonstration of image through row.
6. method of testing as claimed in claim 5, it is characterized in that, in step 3, resistivity cell colors method to set up is: choose the unit that resistivity is the highest and be set to peony: RGB=(255 0 0), choose the unit that resistivity is minimum and be set to mazarine: RGB=(0 0255), according to unit the height of corresponding resistivity, its color to mazarine gradual change, can by the color of image be judged the difference of each cell resistance rate with this by peony.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410117509.4A CN103868958A (en) | 2014-03-27 | 2014-03-27 | Electrical impedance tomography imaging plant root system architecture in-situ observation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410117509.4A CN103868958A (en) | 2014-03-27 | 2014-03-27 | Electrical impedance tomography imaging plant root system architecture in-situ observation method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103868958A true CN103868958A (en) | 2014-06-18 |
Family
ID=50907722
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410117509.4A Pending CN103868958A (en) | 2014-03-27 | 2014-03-27 | Electrical impedance tomography imaging plant root system architecture in-situ observation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103868958A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105326502A (en) * | 2014-07-07 | 2016-02-17 | 北京大学深圳研究生院 | Detection device and method for human body surface electrical characteristic distribution |
CN105701825A (en) * | 2016-01-14 | 2016-06-22 | 天津大学 | Circular-arc grid subdivision method used for electrical tomography sensitivity coefficient solving |
CN105701824A (en) * | 2016-01-14 | 2016-06-22 | 天津大学 | Electrical tomographic image reconstruction method based on circular-arc grid subdivision |
CN107455159A (en) * | 2016-06-06 | 2017-12-12 | 松下知识产权经营株式会社 | Plant growth promotes device and method |
CN107981861A (en) * | 2017-12-11 | 2018-05-04 | 中国人民解放军第四军医大学 | A kind of electrical impedance differentiates the imaging method of multiple target |
WO2019114017A1 (en) * | 2017-12-13 | 2019-06-20 | 广州德科投资咨询有限公司 | Greening monitoring method applied to smart city, and intelligent monitoring robot |
CN110186999A (en) * | 2019-05-10 | 2019-08-30 | 天津大学 | Based on LCR table/impedance analyzer capacitor/electromagnetic double-mode state imaging measurement system |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3994916A (en) * | 1975-08-04 | 1976-11-30 | Imc Chemical Group, Inc. | Vinyl oxazolidines |
CN101194834A (en) * | 2006-12-05 | 2008-06-11 | 重庆博恩富克医疗设备有限公司 | Bio-electrical impedance measuring method and apparatus |
CN103065286A (en) * | 2013-01-08 | 2013-04-24 | 中国人民解放军第四军医大学 | Image reconstruction method in quasi-static electrical impedance imaging |
CN103235190A (en) * | 2013-04-19 | 2013-08-07 | 何为 | Impedance test method |
CN103340625A (en) * | 2013-06-18 | 2013-10-09 | 中国人民解放军第四军医大学 | Regularization method of fast optimization in electrical impedance tomography |
-
2014
- 2014-03-27 CN CN201410117509.4A patent/CN103868958A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3994916A (en) * | 1975-08-04 | 1976-11-30 | Imc Chemical Group, Inc. | Vinyl oxazolidines |
CN101194834A (en) * | 2006-12-05 | 2008-06-11 | 重庆博恩富克医疗设备有限公司 | Bio-electrical impedance measuring method and apparatus |
CN103065286A (en) * | 2013-01-08 | 2013-04-24 | 中国人民解放军第四军医大学 | Image reconstruction method in quasi-static electrical impedance imaging |
CN103235190A (en) * | 2013-04-19 | 2013-08-07 | 何为 | Impedance test method |
CN103340625A (en) * | 2013-06-18 | 2013-10-09 | 中国人民解放军第四军医大学 | Regularization method of fast optimization in electrical impedance tomography |
Non-Patent Citations (2)
Title |
---|
TUSHAR KANTI BERA ET AL.: "Studying the resistivity imaging of chicken tissue phantoms with different current patterns in Electrical Impedance Tomography (EIT)", 《MEASUREMENT》 * |
陈姝君: "电阻抗断层成像技术的研究图像重建算法及实现", 《万方数据》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105326502A (en) * | 2014-07-07 | 2016-02-17 | 北京大学深圳研究生院 | Detection device and method for human body surface electrical characteristic distribution |
CN105701825A (en) * | 2016-01-14 | 2016-06-22 | 天津大学 | Circular-arc grid subdivision method used for electrical tomography sensitivity coefficient solving |
CN105701824A (en) * | 2016-01-14 | 2016-06-22 | 天津大学 | Electrical tomographic image reconstruction method based on circular-arc grid subdivision |
CN105701824B (en) * | 2016-01-14 | 2019-02-19 | 天津大学 | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method |
CN105701825B (en) * | 2016-01-14 | 2019-02-19 | 天津大学 | The arc-shaped grid cutting algorithm solved for electricity chromatography image sensitivity coefficient |
CN107455159A (en) * | 2016-06-06 | 2017-12-12 | 松下知识产权经营株式会社 | Plant growth promotes device and method |
CN107455159B (en) * | 2016-06-06 | 2020-04-14 | 松下知识产权经营株式会社 | Plant growth promoting device and method |
CN107981861A (en) * | 2017-12-11 | 2018-05-04 | 中国人民解放军第四军医大学 | A kind of electrical impedance differentiates the imaging method of multiple target |
CN107981861B (en) * | 2017-12-11 | 2020-08-18 | 中国人民解放军第四军医大学 | Electrical impedance resolution multi-target imaging method |
WO2019114017A1 (en) * | 2017-12-13 | 2019-06-20 | 广州德科投资咨询有限公司 | Greening monitoring method applied to smart city, and intelligent monitoring robot |
CN110186999A (en) * | 2019-05-10 | 2019-08-30 | 天津大学 | Based on LCR table/impedance analyzer capacitor/electromagnetic double-mode state imaging measurement system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103868958A (en) | Electrical impedance tomography imaging plant root system architecture in-situ observation method | |
JP2012042463A5 (en) | ||
CN106706715A (en) | Polluted soil detection method based on three-dimensional high-density electrical resistivity method | |
BRPI0915113A2 (en) | method for electromagnetic subsurface inspection using two or more electromagnetic sources activated simultaneously | |
CN206378448U (en) | Polluted Soil detection arrangement of measuring-line structure based on comprehensive geophysical prospecting methods | |
CN101794453B (en) | Reconstruction method of node mapping image based on regression analysis | |
JP2015127701A (en) | Electric survey method | |
Cubbage et al. | A modified Wenner array for efficient use of eight-channel resistivity meters | |
CN108388707A (en) | D.C. magnetic biasing computational methods based on field circuit method under a kind of three-dimensional asymmetric structure soil model | |
CN105136866A (en) | Contaminated soil detection method based on high-density resistivity method | |
Calixto et al. | 3-D soil stratification methodology for geoelectrical prospection | |
CN103544339A (en) | Modeling calculating method of distribution capacitance between transmission lines of underwater XCTD section plotter | |
CN106125143A (en) | DC fields fixed power supply exploitation method | |
CN110118896A (en) | A kind of method and system measuring stratified soil resistivity and dielectric constant frequency dependent characteristic | |
CN105701824B (en) | Electricity based on arc-shaped mesh generation chromatographs image rebuilding method | |
CN104007476B (en) | Borehole-surface electromagnetic prospecting device | |
CN109758149A (en) | A kind of electrical impedance imaging method of the imaging of the combination time difference and frequency difference imaging | |
KR101999107B1 (en) | Modified electrical resistivity survey device and its interpretation method for water leakage detection of water facilities | |
CN103163404B (en) | Current-voltage mapping construction method based on adjacent stimulus measurement model | |
CN103630750A (en) | Electrical impedance tomography-based gel conductivity measurement method | |
CN205139068U (en) | Pollute soil and detect survey line arrangement structure based on high density resistivity method | |
HUANG et al. | A study on FEM modeling of anomalies of 3‐D high‐density E‐SCAN resistivity survey | |
CN206348292U (en) | Polluted Soil detection arrangement of measuring-line structure based on three-dimensional high-density resistivity method | |
CN105572753A (en) | Contour line mapping and high-resistance region identifying method of transient electromagnetic detection data | |
CN113406707A (en) | Magnetotelluric multi-scale and multi-time-period detection method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140618 |