WO1998023204A1 - Electrical impedance tomography method and electrode arrangement for use therein - Google Patents

Electrical impedance tomography method and electrode arrangement for use therein Download PDF

Info

Publication number
WO1998023204A1
WO1998023204A1 PCT/CA1997/000876 CA9700876W WO9823204A1 WO 1998023204 A1 WO1998023204 A1 WO 1998023204A1 CA 9700876 W CA9700876 W CA 9700876W WO 9823204 A1 WO9823204 A1 WO 9823204A1
Authority
WO
WIPO (PCT)
Prior art keywords
electrodes
electrode
pair
stimulation
conductivity
Prior art date
Application number
PCT/CA1997/000876
Other languages
French (fr)
Inventor
Philip Maurice Church
Philip Michael Wort
Original Assignee
Computing Devices Canada Ltd.
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 Computing Devices Canada Ltd. filed Critical Computing Devices Canada Ltd.
Priority to AU50447/98A priority Critical patent/AU5044798A/en
Priority to US09/308,521 priority patent/US6501984B1/en
Priority to CA002272429A priority patent/CA2272429C/en
Publication of WO1998023204A1 publication Critical patent/WO1998023204A1/en

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0536Impedance imaging, e.g. by tomography

Definitions

  • the invention relates to a method of performing electrical impedance tomography measurements and electrode arrangements for use therein.
  • the invention is applicable generally to the imaging of electrical impedance variations in conductive media, and especially in volume conductors such as parts of the human body, the ground, pipelines, and so on.
  • EIT Electrical impedance tomography
  • APT applied potential tomography
  • the electrical impedance changes are measured by providing an array of electrodes about the volume to be imaged. A stimulus current is supplied to each pair of electrodes in turn and the resulting potential differences recorded between pairs of the remaining electrodes. The process is repeated until all of the independent combinations of stimulation/recording are exhausted. The measurements are used to determine the transfer impedance changes and construct an image of the volume.
  • the current paths in EIT are functions of an unknown electrical conductivity distribution.
  • the quality of the image depends upon the sensitivity of the measuring system, which typically varies with distance from the electrodes.
  • the electrodes are distributed around the body part to be measured, such as the thorax, difficulties may be encountered in resolving impedance changes which occur deep in the body.
  • Most attempts to improve image quality have focused upon the reconstruction algorithms, as disclosed, for example, in U.S. patents Nos. 4,617,939, 5,381,333 and 5,465,730. Limited attention has been devoted to the electrode arrangement. For example, a study of various electrode configurations was described by Booth et al.
  • Hua et al. described experiments with thirty-two compound electrodes encircling a body to be measured. Each compound electrode comprised an inner electrode surrounded by an outer electrode with an annular space between them. In some cases, Hua et al. short-circuited the inner and outer electrodes and used them for both stimulation and recording. In other cases, Hua et al. injected a spatially-sinusoidal current pattern into the ring of 32 outer electrodes and recorded the corresponding generally sinusoidal voltage distribution around the ring of 32 inner electrodes. This approach requires 32 adjustable current generators which must be adjusted individually to give the required pattern as accurately as possible, making it difficult to obtain repeatable measurements accurately. Also, the equipment is complex and costly.
  • Brown et al. describe positioning 16 electrodes around a body, applying a current to a first pair of the electrodes and recording the potential difference between every other pair of the remaining electrodes.
  • Hua et al. seepra
  • Brown et al. avoid using the same electrode for both stimulation and recording.
  • Brown et al. repeat the procedure, applying current to each pair of electrodes in turn. For each pair, Brown et al. measure voltages at all of the remaining electrodes. In practice, such a procedure would not be entirely satisfactory because many of the voltage measurements would be comparable with noise levels.
  • the present invention seeks to mitigate the disadvantages of these known EIT systems and provide an improved electrical impedance tomography procedure with enhanced sensitivity to changes in the electrical conductivity distribution inside the body of interest.
  • a method of determining electrical impedance tomography of a conductive volume comprises the steps of selecting, in turn, a plurality of pairs of locations upon a surface of the conductive volume, supplying a stimulation current to the surface by way of the selected pair of locations, and recording, for each selected pair of locations, a resulting potential difference between at least one pair of the remaining locations, characterized in that, for a particular pair of locations stimulated, the potential difference is recorded for a number of pairs of said remaining locations that is less than the total of said plurality of locations and also comprises that pair of locations which, if stimulated, would produce an electric field with vectors most closely aligned with the corresponding vectors of the electric field produced by stimulation of said particular pair of locations.
  • the above method may be performed using a plurality of electrode means, each disposed upon or adjacent said surface at a respective one of said plurality of locations, each electrode means being used to apply stimulation current to, or record potential at, said respective one of said locations.
  • Each electrode means may comprise two electrodes, that are closely located spatially, one for stimulation and the other for recording. One electrode may be disposed inside the other in a plane which, in use, will be parallel to the said surface. The stimulation current could then be applied to two outer electrodes and the potential difference measured at a said selected pair comprising the corresponding two inner electrodes. Conversely, the stimulation current could be applied to the two inner electrodes and the potential difference measured at the corresponding outer electrodes.
  • each electrode means may comprise a first electrode and a second electrode spaced apart from each other in a direction which, in use, will be substantially normal to the said surface. The stimulation current could be applied by way of two of said second electrodes and the potential difference measured between two of said first electrodes, or vice versa.
  • a method of measuring electrical conductivity distribution within a conductive volume includes the steps of positioning adjacent a surface of the body an array of electrodes, with a first group of the electrodes closer to the surface than a second group of the electrodes, applying a stimulation current to selected pairs of electrodes in turn and measuring, for each pair stimulated, corresponding electrical potential differences produced between pairs of the remaining electrodes, and processing the measured potential differences to determine electrical conductivity variations within the body.
  • an electrode arrangement for an electrical impedance tomography system comprises a plurality of electrodes and support means for supporting the electrodes adjacent a surface of a conductive body the conductivity of which is to be mapped by the system, the arrangement being such that a first group of said plurality of electrodes will be closer to the surface than a second group of the plurality of electrodes.
  • an electrode arrangement for an electrical impedance tomography system comprises a plurality of electrodes in an array mounted in a support medium for supporting the electrodes adjacent a surface of a volume to be imaged, the support medium having anisotropic conductivity, its conductivity in a direction that, in use, is normal to the surface, being significantly greater than its conductivity in a transverse direction.
  • an electrode arrangement for an electrical impedance tomography system comprises an array of electrodes fixed spatially relative to each other in a rigid support medium, and an interface medium for interfacing the rigid support to a surface of a conductive volume to be measured, the interface medium being conductive so as to connect the electrodes electrically to the surface and pliable so as to conform to variations in relief of the surface.
  • An advantage of this electrode arrangement of such further aspect is that the relative position of the electrodes, and the dimensions of the array, are fixed, which facilitates static image reconstruction of the image.
  • an electrode arrangement for use in electrical impedance tomography comprises a plurality of electrode means each formed by a stimulation electrode and a recording electrode closely located spatially.
  • the electrode means comprises an inner electrode surrounded by an outer electrode with an annular channel between them. Electrically insulating sealing means may be provided in the annular channel for contacting said surface, in use, and insulating the inner electrode electrically from the outer electrode.
  • Figure 1 is a diagram of the generic components that are part of a typical EIT apparatus
  • Figure 2 is a diagram of a first embodiment of the invention comprising an annular EIT electrode array
  • Figure 3 is a three-dimensional view of the annular electrode array of Figure 2;
  • Figure 4 is a cross-sectional detail view taken on the line IV-IV of Figure 2;
  • Figure 5 compares the sensitivity values of the annular EIT electrode array of Figure 2 with those of the EIT system of Figure 1, for radial positions;
  • Figure 6 illustrates the sensitivity theorem used in determining electrical impedance tomography (EIT) of a volume conductor
  • Figure 7 is a flowchart representing a typical measurement procedure using the electrode means of Figure 2;
  • Figure 8 is a plan view of an alternative form of electrode means, comprising an inner electrode and an outer electrode (also referred to as a compound electrode;
  • Figure 9 is a cross-sectional view of the electrode means of Figure 8;
  • Figure 10 illustrates a second embodiment of the invention in the form of a planar rectangular electrode array for a ground probing application;
  • Figure 11 illustrates one cell of the electrode array of Figure 10
  • Figure 12 illustrates a third embodiment in the form of a cylindrical electrode array for three-dimensional imaging of a human torso
  • Figures 13 A to 13J illustrate the image reconstruction obtained from a conventional EIT electrode array, such as that shown in Figure 1 , and an electrode array in which the stimulating electrodes and recording electrodes are co-located, such as that shown in Figure 2.
  • a typical known EIT apparatus comprises an electrode array comprising sixteen electrodes E, to E 16 distributed, equally spaced apart, around a conductive volume 12, to which they are attached directly, if necessary with a suitable conductive gel ensuring good electrical contact.
  • the conductive volume 12 may be, for example, a human thorax.
  • the electrodes are connected by a link 14 to electronic conditioning circuitry 16 which will usually include components such as a current source controlled by a D-to-A-converter, an A-to-D converter, differential amplifier, filter, analog multiplexer, clock, and the digital I/O required to control the circuitry from a computer 18 to which the electronic conditioning circuitry is connected.
  • the computer 18 may be a personal computer equipped with a digital signal processor card used for the image reconstruction process and a suitable display 20 for displaying the images.
  • a first pair of the electrodes E, and E 2 are stimulated by injecting a current of, say, 4mA at 50 kHz, into them and the corresponding potential differences measured between each of the other pairs of electrodes, i.e. (E 3 -E 4 ), (E 4 -E 5 ), (E 5 -E 6 ), and so on up to (E 15 -E 16 ).
  • the potential differences between electrodes E 16 and E E, and E ⁇ and Ej and E 3 are not measured because the voltage drop across the contact resistance due to the injected current would affect the measurements.
  • the next pair of electrodes E 2 and E 3 then are stimulated and the potential differences between pairs of the other electrodes measured. The procedure is repeated until all of the pairs of electrodes have been stimulated in turn.
  • the accumulated measurements are processed by the computer 18 and an image of the electrical impedance variation across the conductive volume is created and displayed on display unit 20.
  • Various algorithms have been disclosed for computing these impedance variation, such as the so-called "back-projection" reconstruction algorithm and the regularized inverse of the sensitivity matrix.
  • Embodiments of the present invention favour high signal-to-noise ratio measurements by limiting the measurements, for a given stimulation pair, to that pair or pairs whose "recording" electric field vectors most closely align with the electric field vectors of the stimulated pair of electrodes.
  • the recording electrodes measure potential difference so their "recording" electric field vectors is really that which would be produced if a stimulation current were injected into them.
  • the present invention aims to keep the stimulating and recording electric field vectors as close to each other as possible. (Ideally, they would be superimposed). Various arrangements for doing so will now be described.
  • the electrode array 22 shown in Figures 2 and 3 differs from that shown in Figure 1 in that it comprises two groups of electrodes arranged in an annular configuration, i.e. in two rings, one within the other, and supported by an annular support.
  • the inner group comprises recording electrodes R, to R 16 and the outer group comprises stimulating electrodes S, to S 16 .
  • the annular support comprises annular flanges 24 and 26 (see Figure 4) bonded to opposite sides of an array of radial separators 28, - 28, 6 and a circumferential band 30 is bonded to their respective outer edges, forming a ring of wedge-shaped compartments or cells 32, to 32 16 which are open at their inner ends 34, to 34, 6 and alternate with separators 28, to 28 16 .
  • the annular flanges 24 and 26, separators 28, and outer band 30 are of insulating material.
  • the wedge-shaped compartments 32, to 32, 6 are equally distributed and each subtends a 15 degree angle. They contain a conductive material whose electrical conductivity is comparable to the expected average electrical conductivity of the imaging medium.
  • the compartments 32, to 32, 6 might contain agar jelly formed by mixing agar powder with saline solution to form a conductive gel, giving a conductivity for the wedge-shaped compartment of about 0.002 mho/cm.
  • the separators 28, to 28, 6 are poorly conductive and each subtends an angle of 7.5 degrees so that the annular support 5 is circular.
  • the separators 28, to 28 16 might be of synthetic plastics material, for example polypropylene.
  • a conductive interfacing medium 36 for example sponge material soaked in saline solution, is provided between the conductive volume or body 12' and the innermost surface of the annular support 22, i.e. the innermost surface of the agar jelly
  • the sponge 36 is pliable and conforms to the shape of the conductive volume surface. It connects each of the electrodes R, to R 16 and S, to
  • the conductivity of the interfacing medium 36 is comparable to the conductivity of the imaging medium 12'.
  • the stimulation electrode is mounted upon the inner face of band 30 and connected through the band to terminal 38 4 .
  • the recording electrode R 4 is suspended in the middle of the compartment 32 4 by a plastics ligature 40 4 , extending from the annular flange 26 and a conductive filament 42 4 extending from the other annular flange 24.
  • the filament is suspended in the middle of the compartment 32 4 by a plastics ligature 40 4 , extending from the annular flange 26 and a conductive filament 42 4 extending from the other annular flange 24.
  • terminal 44 4 connects to terminal 44 4 on the outer surface of flange 24.
  • the terminals 38, to 38, 6 and 44, to 44 4 are not shown in Figures 2 and 3 for clarity.
  • the terminals 38 and 44 are connected to the electronic conditioning circuit 16'. For convenience, the individual connections are not shown in Figures 2 and 3.
  • the circuit 16' permits selection of individual pairs of the outer electrodes S, - S, 6 for
  • Each of the electrodes R, to R, 6 and S, to S, 6 comprises a disc mounted with its plane normal to a radius of the annular support, each pair of electrodes in a cell being aligned radially, i.e. on a line extending normal to the surface of the body 12'.
  • Electrode 30 electrode is about 6 mm in diameter.
  • the electrodes R, to R, 6 and S, to S 16 in their respective compartments are spaced apart radially by a distance sufficient to reduce/avoid shielding effects.
  • Other sizes of electrodes for example from 2 mm to 10 mm diameter, are available. If they were used, the spacing would be adjusted appropriately.
  • the annular support can be considered to have anisotropic conductivity, i.e. high conductivity for currents flowing radially from an electrode S towards the volume 12' and low conductivity for cu ⁇ ents flowing circumferentially from one compartment towards its neighbour.
  • Figure 5 illustrates the average sensitivity, in the radial direction, obtained using a finite element model, for an electrical impedance tomography apparatus using the annular electrode a ⁇ ay of Figure 2, and that for an equivalent conventional electrical impedance tomography apparatus using the electrode array of Figure 1.
  • the conductive volume 12' has a radius of 10 cm, with a uniform conductivity.
  • the annulus in Figure 2 extends from 10 cm to 15 cm and has a conductivity 10 times larger than the conductivity of conductive volume 12'.
  • the medium separating the individual compartments, i.e. of separators 28, to 28, 2 is set to a conductivity value of 1 % of the conductivity of conductivity volume 12'.
  • the average sensitivity curve obtained for the conventional electrode a ⁇ ay of Figure 1 uses configurations of stimulator and recording electrode pairs as disclosed by Brown et al. in U.S. patent No. 4,617,939. As mentioned earlier, in this approach, a pair of neighbouring electrodes is used to stimulate while the electrical potential is measured with the remaining pairs of electrodes. The process is repeated until all the independent combinations are exhausted. The electrodes used for stimulation, however, are not used also for measurement.
  • the average sensitivity curve for the electrode configuration shown in Figure 2 is obtained using pairs of electrodes in the outer group S,-S, 6 to apply stimulation cu ⁇ ent and the co ⁇ esponding pairs of electrode in the inner group R,-R, 6 to measure the resulting electrical potential.
  • the resolution of the reconstructed image depends upon the sensitivity of the apparatus to changes of the electric conductivity in the conductive volume 12'.
  • the sensitivities depicted in Figure 5 were obtained using the sensitivity expression described by D.B. Geselowitz in "An Application of Electrocardiographic Lead Theory to Impedance Plethysmography", IEEE Trans. Biomed. Eng., Vol. BME-18, pp. 38-41, Jan. 1971, and by J. Lehr in "A Vector Derivation Useful in Impedance Plethysmographic field calculation", IEEE Trans. Biomed. Eng. , Vol. BME-19, pp. 156-157, Mar. 1972.
  • is the potential distribution when the cu ⁇ ent I ⁇ is applied to the electrode pair (A,B). is the potential distribution when the cu ⁇ ent I ⁇ is applied to the electrode pair (C,D).
  • ⁇ (x,y,z) is the conductivity distribution when is established for I ⁇ ⁇ (x,y,z,) + ⁇ (x,y,z) is the conductivity distribution when is established for I .
  • ⁇ Zj is the difference in transfer impedance for the ⁇ th combination of a pair of stimulating electrodes and a pair of recording electrodes
  • ⁇ j is the unknown conductivity difference for the region element j
  • Sj j is the sensitivity coefficient term for the stimulating/recording electrode pair combination i, at the element j; and m is the number of domains into which the region has been divided.
  • Sensitivity enhancement is achieved by selecting configurations of stimulating and recording electrode pairs that favour the scalar product in the numerator and using an annulus medium with an anisotropic conductivity that limits the cu ⁇ ent intensities at the denominator.
  • the annular electrode a ⁇ ay of Figure 2 allows electrode configurations that can maximize scalar products of stimulating and recording electric field vectors more so than could be obtained with the electrodes of Figure 1. It permits more combinations of electrodes to be made, with a greater number of different geometries.
  • stimulating electrodes need not be diametrically opposite each other. Referring to Figure 2, using pairs of electrodes in the outer group S, to S 16 to stimulate the body and pairs of electrodes in the inner group R, to R, 6 to measure the resulting electrical potentials tends to maximize the scalar product in the numerator of expression (4).
  • the sensitivity of the apparatus to a change in the electrical conductivity at a position (x,y) in the conductive volume is an inverse function of the product of the total electric cu ⁇ ent intensities that would be required to produce the electric field intensity of the stimulating and recording electrode pairs at position (x,y).
  • Disposing the electrodes in sectors whose boundaries are poorly conductive produces an annulus medium with an anisotropic conductivity, vis. conductive in the radial direction and poorly conductive in the direction. This a ⁇ angement forces the cu ⁇ ent to flow in the medium of interest, i.e. the conductive volume, as opposed to the extended medium. It favours a larger ratio of the electric field to cu ⁇ ent intensity, in the medium of interest.
  • Annulus dimension The size of the annulus used for this evaluation has an outer to inner radius ratio of 1.5. Evaluations have shown that as the ratio increases, and the electrodes are set further away from the imaging volume, the sensitivity decreases. Conversely, as the ratio decreases, and the electrodes get closer to the imaging medium and the sensitivity increases. Although this is true for a discrete point electrode model, the optimal practical ratio will be determined by the size of the electrodes and the minimum separation required to avoid mutual shielding.
  • a minimum of 32 electrodes were considered, arranged in two concentric rings. For this evaluation, each wedge-shaped compartment held two electrodes, aligned along the diameter. One electrode was positioned at the edge of the annulus and the other was positioned at the centre of the annulus. Although this was effective for the model, in practice, of course, the other electrode cannot be any closer to the centre than the inner surface of the annulus. It is expected that standard commercial AgCl disk shaped electrodes would be used for this application.
  • an initial conductivity distribution appropriate to the volume would be assumed (step 7.1).
  • a uniform mean conductivity distribution would be assumed. For example, for landmine detection, this might be a typical conductivity distribution of the particular type of soil whereas, for human body imaging, it might be a typical conductivity distribution for that part of the human body.
  • step 7.2 a series of measurements are made, passing a known cu ⁇ ent through each pair of stimulation electrodes S, to S, 6 in turn and recording the potential difference across the "closest" pair of electrodes, i.e. having the recording electric field vectors most closely aligned with those of the stimulating pair of electrodes. For example, when stimulating S, and S 2 , the potential difference between R, and R 2 only would be measured.
  • stimulation-record pairs (stimulate-record pairs). It should be noted that the stimulating electrodes and the recording electrodes are always identically numbered, and hence closely located. Other stimulation/recording electrode configurations which share this feature will have some benefit.
  • An example of such an alternative, possible with only a single set of electrodes as shown in Figure 1 (denoted E,,E 2 ... E, 6 ) is:
  • This pattern has a reduced number of independent combinations, since it starts with a stimulating electrode gap of 2, for example, but will show some of the benefit of the proposed configuration, and is shown as a simple example of an alternative manifestation of the stimulate/record configuration described herein.
  • the region of interest i.e. the region in which cu ⁇ ent might flow
  • the region of interest is the conductive volume 12', the sponge 36 and the annular a ⁇ ay of compartments 32,-32 16 .
  • Each element has dimensions equal to approximately 10 per cent of the spacing between each pair of the recording electrodes R, to R, 6 .
  • the mesh could be triangular or rectangular, as desired, and the region could be two-dimensional or three-dimensional.
  • the sensitivity matrix S y is calculated for the mesh using the electrode geometry, initial conductivity distribution, and applied cu ⁇ ents, according to equation [4], for the region of interest.
  • step 7.5 an estimate is made of the noise present in the potential difference measurements divided by the applied stimulation cu ⁇ ent.
  • This scaled noise value is used to determine a parameter ⁇ which is the ratio of the scaled noise value to the root mean square (RMS) of the initial conductivity distribution.
  • RMS root mean square
  • step 7.6 the pseudo-inverse S ⁇ of the sensitivity matrix S ;j is calculated according to equations (5) and (6).
  • step 7.7 the differential conductivity distribution ⁇ is calculated according to equation (7).
  • Preferred electrode configuration and compound electrodes are Preferred electrode configuration and compound electrodes.
  • Compound electrodes are so named because they comprise two parts, one part being used for stimulation and the other for recording. It has been found that compound electrodes, in conjunction with the novel measurement procedure described above, provide significantly better performance in physically realistic (noisy) situations. Analysis has shown that this improvement in performance is due both to sensitivity amplitude enhancements (due to alignment of the electric field stimulating and recording electric field vectors) and a better representation (more suitable for reconstruction) of the set of measurements possible by the proposed recording configuration.
  • a potential problem associated with compound electrodes is the possibility of a partial electrical short circuit between the stimulating and measuring electrode, due to the close proximity of these electrodes. Simulations have indicated that satisfaction will be obtained with spacing of the component electrodes of a compound electrode such as to render the voltage associated with this extra (surface) cu ⁇ ent path approximately 0.1 % of that associated with the conventional path (through the region of interest). For a human body this requirement translates into a resistance of (for skin surface alone) approximately 5M ⁇ .
  • An example of a suitable compound electrode is shown in Figures 8 and 9.
  • the compound electrode comprises an annular outer electrode 46 su ⁇ ounding a circular inner electrode 48 with an annular slot 50 between them.
  • the electrodes 46 and 48 are bonded to an insulating planar support 52 and their exposed surfaces are coated with layers of conductive gel 54 and 56, respectively.
  • An "O" ring 58 of insulating material is located in the annular slot 50.
  • the diameter of the "O" ring 58 is approximately equal to the depth of the annular slot 50 so that, when the compound electrode is applied to the surface of the volume under investigation, the "O" ring 58 will prevent gel 54/56 from exuding across the annular slot 50 and causing, effectively, a short-circuit between the inner and outer electrodes. In use, a plurality of such compound electrodes would be distributed about the surface of the body being investigated.
  • Cu ⁇ ent could be applied to the respective outer electrodes of a pair of compound electrodes and the potential difference recorded between the inner electrodes of the same pair. Conversely, the inner electrodes could be stimulated and the potential difference recorded between the outer electrodes. It should be appreciated that other configurations of compound electrode could be used. For example, either or both of the electrodes could be rectangular. It is also envisaged that each compound electrode might comprise an a ⁇ ay or grid of stimulation electrodes interspersed with an a ⁇ ay or grid of recording electrodes.
  • FIG. 10 illustrates a two-dimensional a ⁇ ay of electrode cells which is shown as rectangular but could have other shapes.
  • This electrode array is envisaged for ground imaging applications, such as landmine detection.
  • the array comprises juxtaposed individual electrode cell units 60.
  • each cell comprises a pair of electrodes R, S immersed in a conductive material, for example a jelly formed from agar and saline solution, as before.
  • the cell units 60 lie on a bed 62 ( Figure 10) of conductive material that provides the electrical interface with the ground. This bed 62 would have insulating partitions (not shown) to maintain anisotropic conductivity.
  • Figure 12 illustrates another possible embodiment, where a flexible vest 64 comprises similar individual electrode cells 66 forming a cylindrical a ⁇ ay which, as shown, can be wrapped around a human torso and used for three-dimensional imaging.
  • the compound electrodes would still be in a rigid support, if necessary, with an intervening sponge 36, or the like.
  • the proposed reconstruction technique is based on the construction of a probabilistically regulated, pseudo inverse of the sensitivity matrix introduced above and mentioned in step 7.4.
  • the construction of this pseudo inverse is very straightforward, and follows the following steps :- Calculate the singular value decomposition of the sensitivity matrix.
  • ⁇ * diag( ⁇ *) ⁇
  • RMS root mean square
  • n is the number of independent measurements taken.
  • FIGS 13 A to 13J illustrate a comparison between the novel measurement procedure set out in Table 1, using compound electrodes as per Figure 8, and the measurement procedure disclosed by Brown et al. in U.S. Patent No. 4,617,939, but in both cases using the reconstruction technique described above.
  • Figures 13A to 13J depict the reconstruction of a perturbation covering approximately 10% of the diameter of a circular region of interest at three different locations. As shown in Figures 13A, 13D and 13G, with zero noise, both Brown et al. 's procedure and the procedure of the present invention reproduce the perturbation very accurately.
  • Figure 13A shows a maximum a posteriori (MAP) estimate of a perturbation near the circumference of the region.
  • Figure 13B shows the (MAP) estimate of the differential conductivity distribution ⁇ for a circular a ⁇ ay of 16 compound electrodes of the kind shown in Figure 8.
  • Figure 13C shows the co ⁇ esponding MAP estimate for the conventional electrode configuration of Figure 1 using Brown et al. 's measurement technique.
  • Figures 13D to 13F co ⁇ espond to Figures 13 A to 13C for a perturbation which is closer to the middle of the region and Figures 13G to 13J are the co ⁇ esponding reconstructions for a perturbation at the centre of the region.
  • the electrodes could be in a fixed a ⁇ ay and a separate cu ⁇ ent generator provided for each pair stimulated.
  • a single cu ⁇ ent generator could be used to stimulate each pair of electrodes in turn, using a commutation switching circuit, enabling a more expensive cu ⁇ ent generator to be used.
  • a suitable custom integrated circuit preprogrammed with the reconstruction algorithm software, could be combined with the electronic conditioning circuitry in an interface unit which would connect to a conventional personal computer and/or even directly to a display device.
  • the compound electrode a ⁇ angement may be easier to manufacture.
  • the measurement procedure facilitates simpler implementations for a given sensitivity.

Abstract

An electrode arrangement for an electrical impedance tomography system comprises a plurality of electrodes (S1-S16) in an array (22) mounted in a support medium (24, 26, 28, 30, 32) for supporting the electrodes adjacent a surface of a volume (12) the electrical conductivity distribution of which is to be measured, for example the thorax or another body part, a pipeline, the ground, and so on. The electrodes are disposed in two groups, the arrangement being such that, in use, one group (R1-R16) will be closer to the surface than the other group (S1-S16). In operation, each pair of electrodes will be stimulated in turn by applying a known current to them. For each stimulated pair, a resulting potential difference will be recorded at the pair of the remaining electrodes which, if stimulated, would produce an electric field with vectors of the electric field produced by the stimulated pair of electrodes.

Description

ELECTRICAL IMPEDANCE TOMOGRAPHY METHOD AND ELECTRODE ARRANGEMENT FOR USE THEREIN
DESCRIPTION TECHNICAL FIELD:
The invention relates to a method of performing electrical impedance tomography measurements and electrode arrangements for use therein.
The invention is applicable generally to the imaging of electrical impedance variations in conductive media, and especially in volume conductors such as parts of the human body, the ground, pipelines, and so on.
BACKGROUND ART:
Electrical impedance tomography (EIT), also known as "applied potential tomography" (APT), is used to provide images of spatial variations in electrical impedance within a conductive volume conductor. The electrical impedance changes are measured by providing an array of electrodes about the volume to be imaged. A stimulus current is supplied to each pair of electrodes in turn and the resulting potential differences recorded between pairs of the remaining electrodes. The process is repeated until all of the independent combinations of stimulation/recording are exhausted. The measurements are used to determine the transfer impedance changes and construct an image of the volume. In contrast to X-ray computer tomography applications, where the paths of photons through a body are straight lines, the current paths in EIT are functions of an unknown electrical conductivity distribution. This gives rise to a non-linear image reconstruction problem. Nevertheless, algorithms have been devised for converting the series of measurements into an image of the electrical impedance distribution within the body at a sufficiently high rate that changes associated with respiratory and cardiac functions, for example, can be monitored.
The quality of the image depends upon the sensitivity of the measuring system, which typically varies with distance from the electrodes. In clinical applications, for example, where the electrodes are distributed around the body part to be measured, such as the thorax, difficulties may be encountered in resolving impedance changes which occur deep in the body. Most attempts to improve image quality have focused upon the reconstruction algorithms, as disclosed, for example, in U.S. patents Nos. 4,617,939, 5,381,333 and 5,465,730. Limited attention has been devoted to the electrode arrangement. For example, a study of various electrode configurations was described by Booth et al. in "A Comparison of Three Electrode Configurations for Electrical Impedance Tomography", IEE Colloquium on "Innovation in Instrumentation for Electrical Tomography", pp. 11/1, 11/3, 1995. British patent application number 2,257,530 disclosed electrodes in a circular array or "rosette" placed upon a substantially flat part of the body, such as the chest. This was said to be especially desirable when imaging the heart, since impedance changes occurring in the heart during the cardiac cycle will be larger than those occurring in the region of the skin. Little attention has been paid to the way in which the measurements themselves are made. There are two kinds of procedure for taking the measurements. One involves applying stimulation currents simultaneously to a ring of electrodes to generate a current pattern and simultaneously measuring the corresponding voltage distribution around the same or adjacent electrodes. An example of such an approach is described in an article by P. Hua et al. entitled "An Electrical Impedance Tomograph using Compound Electrodes", presented at the IEEE Engineering in Medicine & Biology Society 11th. Annual International Conference, 1989. Hua et al. described experiments with thirty-two compound electrodes encircling a body to be measured. Each compound electrode comprised an inner electrode surrounded by an outer electrode with an annular space between them. In some cases, Hua et al. short-circuited the inner and outer electrodes and used them for both stimulation and recording. In other cases, Hua et al. injected a spatially-sinusoidal current pattern into the ring of 32 outer electrodes and recorded the corresponding generally sinusoidal voltage distribution around the ring of 32 inner electrodes. This approach requires 32 adjustable current generators which must be adjusted individually to give the required pattern as accurately as possible, making it difficult to obtain repeatable measurements accurately. Also, the equipment is complex and costly.
An example of the other kind of procedure is described in US patent number 4,617,939 (Brown et al.) issued October 1986. Brown et al. describe positioning 16 electrodes around a body, applying a current to a first pair of the electrodes and recording the potential difference between every other pair of the remaining electrodes. As also noted by Hua et al. (supra), significant contact impedance between the electrodes and the body surface militates against measurement sensitivity, so Brown et al. avoid using the same electrode for both stimulation and recording. Brown et al. repeat the procedure, applying current to each pair of electrodes in turn. For each pair, Brown et al. measure voltages at all of the remaining electrodes. In practice, such a procedure would not be entirely satisfactory because many of the voltage measurements would be comparable with noise levels.
The present invention seeks to mitigate the disadvantages of these known EIT systems and provide an improved electrical impedance tomography procedure with enhanced sensitivity to changes in the electrical conductivity distribution inside the body of interest.
DISCLOSURE OF INVENTION:
According to one aspect of the invention, a method of determining electrical impedance tomography of a conductive volume comprises the steps of selecting, in turn, a plurality of pairs of locations upon a surface of the conductive volume, supplying a stimulation current to the surface by way of the selected pair of locations, and recording, for each selected pair of locations, a resulting potential difference between at least one pair of the remaining locations, characterized in that, for a particular pair of locations stimulated, the potential difference is recorded for a number of pairs of said remaining locations that is less than the total of said plurality of locations and also comprises that pair of locations which, if stimulated, would produce an electric field with vectors most closely aligned with the corresponding vectors of the electric field produced by stimulation of said particular pair of locations.
The above method may be performed using a plurality of electrode means, each disposed upon or adjacent said surface at a respective one of said plurality of locations, each electrode means being used to apply stimulation current to, or record potential at, said respective one of said locations.
Each electrode means may comprise two electrodes, that are closely located spatially, one for stimulation and the other for recording. One electrode may be disposed inside the other in a plane which, in use, will be parallel to the said surface. The stimulation current could then be applied to two outer electrodes and the potential difference measured at a said selected pair comprising the corresponding two inner electrodes. Conversely, the stimulation current could be applied to the two inner electrodes and the potential difference measured at the corresponding outer electrodes. Alternatively, each electrode means may comprise a first electrode and a second electrode spaced apart from each other in a direction which, in use, will be substantially normal to the said surface. The stimulation current could be applied by way of two of said second electrodes and the potential difference measured between two of said first electrodes, or vice versa.
According to another embodiment of the invention, a method of measuring electrical conductivity distribution within a conductive volume includes the steps of positioning adjacent a surface of the body an array of electrodes, with a first group of the electrodes closer to the surface than a second group of the electrodes, applying a stimulation current to selected pairs of electrodes in turn and measuring, for each pair stimulated, corresponding electrical potential differences produced between pairs of the remaining electrodes, and processing the measured potential differences to determine electrical conductivity variations within the body.
According to another aspect of the present invention, an electrode arrangement for an electrical impedance tomography system comprises a plurality of electrodes and support means for supporting the electrodes adjacent a surface of a conductive body the conductivity of which is to be mapped by the system, the arrangement being such that a first group of said plurality of electrodes will be closer to the surface than a second group of the plurality of electrodes. According to yet another aspect of the invention, an electrode arrangement for an electrical impedance tomography system comprises a plurality of electrodes in an array mounted in a support medium for supporting the electrodes adjacent a surface of a volume to be imaged, the support medium having anisotropic conductivity, its conductivity in a direction that, in use, is normal to the surface, being significantly greater than its conductivity in a transverse direction.
According to a further embodiment of the invention, an electrode arrangement for an electrical impedance tomography system comprises an array of electrodes fixed spatially relative to each other in a rigid support medium, and an interface medium for interfacing the rigid support to a surface of a conductive volume to be measured, the interface medium being conductive so as to connect the electrodes electrically to the surface and pliable so as to conform to variations in relief of the surface.
An advantage of this electrode arrangement of such further aspect, as compared with electrodes which are simply attached directly to the volume, is that the relative position of the electrodes, and the dimensions of the array, are fixed, which facilitates static image reconstruction of the image.
According to yet a further aspect of the invention, an electrode arrangement for use in electrical impedance tomography comprises a plurality of electrode means each formed by a stimulation electrode and a recording electrode closely located spatially. In one embodiment of this aspect, the electrode means comprises an inner electrode surrounded by an outer electrode with an annular channel between them. Electrically insulating sealing means may be provided in the annular channel for contacting said surface, in use, and insulating the inner electrode electrically from the outer electrode. Various objects, features, aspects and advantages of the present invention will become more apparent from the following detailed description of preferred embodiments of the invention, taken in conjunction with the accompanying drawings.
BRIEF DESCRIPTION OF THE DRAWINGS: Figure 1, labelled PRIOR ART, is a diagram of the generic components that are part of a typical EIT apparatus;
Figure 2 is a diagram of a first embodiment of the invention comprising an annular EIT electrode array;
Figure 3 is a three-dimensional view of the annular electrode array of Figure 2; Figure 4 is a cross-sectional detail view taken on the line IV-IV of Figure 2;
Figure 5 compares the sensitivity values of the annular EIT electrode array of Figure 2 with those of the EIT system of Figure 1, for radial positions;
Figure 6 illustrates the sensitivity theorem used in determining electrical impedance tomography (EIT) of a volume conductor; Figure 7 is a flowchart representing a typical measurement procedure using the electrode means of Figure 2;
Figure 8 is a plan view of an alternative form of electrode means, comprising an inner electrode and an outer electrode (also referred to as a compound electrode; Figure 9 is a cross-sectional view of the electrode means of Figure 8; Figure 10 illustrates a second embodiment of the invention in the form of a planar rectangular electrode array for a ground probing application;
Figure 11 illustrates one cell of the electrode array of Figure 10; Figure 12 illustrates a third embodiment in the form of a cylindrical electrode array for three-dimensional imaging of a human torso; and
Figures 13 A to 13J illustrate the image reconstruction obtained from a conventional EIT electrode array, such as that shown in Figure 1 , and an electrode array in which the stimulating electrodes and recording electrodes are co-located, such as that shown in Figure 2.
BEST MODES FOR CARRYING OUT THE INVENTION:
Referring to Figure 1, a typical known EIT apparatus comprises an electrode array comprising sixteen electrodes E, to E16 distributed, equally spaced apart, around a conductive volume 12, to which they are attached directly, if necessary with a suitable conductive gel ensuring good electrical contact. The conductive volume 12 may be, for example, a human thorax. The electrodes are connected by a link 14 to electronic conditioning circuitry 16 which will usually include components such as a current source controlled by a D-to-A-converter, an A-to-D converter, differential amplifier, filter, analog multiplexer, clock, and the digital I/O required to control the circuitry from a computer 18 to which the electronic conditioning circuitry is connected. The computer 18 may be a personal computer equipped with a digital signal processor card used for the image reconstruction process and a suitable display 20 for displaying the images. Using the procedure disclosed in US 4,617,939, a first pair of the electrodes E, and E2 are stimulated by injecting a current of, say, 4mA at 50 kHz, into them and the corresponding potential differences measured between each of the other pairs of electrodes, i.e. (E3-E4), (E4-E5), (E5-E6), and so on up to (E15-E16). The potential differences between electrodes E16 and E E, and E^ and Ej and E3 are not measured because the voltage drop across the contact resistance due to the injected current would affect the measurements.
The next pair of electrodes E2 and E3 then are stimulated and the potential differences between pairs of the other electrodes measured. The procedure is repeated until all of the pairs of electrodes have been stimulated in turn. The accumulated measurements are processed by the computer 18 and an image of the electrical impedance variation across the conductive volume is created and displayed on display unit 20. Various algorithms have been disclosed for computing these impedance variation, such as the so-called "back-projection" reconstruction algorithm and the regularized inverse of the sensitivity matrix.
While this procedure might appear satisfactory theoretically, in practice, noise may obscure many of the readings and reduce the quality of the reconstruction. Embodiments of the present invention favour high signal-to-noise ratio measurements by limiting the measurements, for a given stimulation pair, to that pair or pairs whose "recording" electric field vectors most closely align with the electric field vectors of the stimulated pair of electrodes. It will be appreciated that the recording electrodes measure potential difference so their "recording" electric field vectors is really that which would be produced if a stimulation current were injected into them. Generally, the present invention aims to keep the stimulating and recording electric field vectors as close to each other as possible. (Ideally, they would be superimposed). Various arrangements for doing so will now be described.
A first embodiment of the invention will now be described with reference to Figures 2 and 3 in which the electronic conditioning circuitry 16', computer 18' and display unit 20' are similar to those of Figure 1 and so have the same reference numerals, but with a prime.
The electrode array 22 shown in Figures 2 and 3 differs from that shown in Figure 1 in that it comprises two groups of electrodes arranged in an annular configuration, i.e. in two rings, one within the other, and supported by an annular support. The inner group comprises recording electrodes R, to R16 and the outer group comprises stimulating electrodes S, to S16. The annular support comprises annular flanges 24 and 26 (see Figure 4) bonded to opposite sides of an array of radial separators 28, - 28,6 and a circumferential band 30 is bonded to their respective outer edges, forming a ring of wedge-shaped compartments or cells 32, to 3216 which are open at their inner ends 34, to 34,6 and alternate with separators 28, to 2816. The annular flanges 24 and 26, separators 28, and outer band 30 are of insulating material. Each pair of an inner electrode R and an outer electrode S mounted in a corresponding one of the plurality of wedge-shaped compartments or cells 32, to 32,6. The wedge-shaped compartments 32, to 32,6 are equally distributed and each subtends a 15 degree angle. They contain a conductive material whose electrical conductivity is comparable to the expected average electrical conductivity of the imaging medium. For example, for measuring a human body part with a conductivity of, say, 0.002 mho/cm, the compartments 32, to 32,6 might contain agar jelly formed by mixing agar powder with saline solution to form a conductive gel, giving a conductivity for the wedge-shaped compartment of about 0.002 mho/cm. The separators 28, to 28,6 are poorly conductive and each subtends an angle of 7.5 degrees so that the annular support 5 is circular. The separators 28, to 2816 might be of synthetic plastics material, for example polypropylene.
In use, a conductive interfacing medium 36, for example sponge material soaked in saline solution, is provided between the conductive volume or body 12' and the innermost surface of the annular support 22, i.e. the innermost surface of the agar jelly
10 filling the compartments 32, to 32,6. The sponge 36 is pliable and conforms to the shape of the conductive volume surface. It connects each of the electrodes R, to R16 and S, to
S,6 electrically to the surface of conductive volume 12'. The conductivity of the interfacing medium 36 is comparable to the conductivity of the imaging medium 12'.
The mounting of the electrodes S, to S,6 in compartments 32, to 32,6 is the same
15 in each case, so only the mounting of electrodes S4 and R4 described - as an example. Thus, the stimulation electrode is mounted upon the inner face of band 30 and connected through the band to terminal 384. The recording electrode R4 is suspended in the middle of the compartment 324 by a plastics ligature 404, extending from the annular flange 26 and a conductive filament 424 extending from the other annular flange 24. The filament
20 424 connects to terminal 444 on the outer surface of flange 24. The terminals 38, to 38,6 and 44, to 444 are not shown in Figures 2 and 3 for clarity.
The terminals 38 and 44 are connected to the electronic conditioning circuit 16'. For convenience, the individual connections are not shown in Figures 2 and 3. The circuit 16' permits selection of individual pairs of the outer electrodes S, - S,6 for
25 stimulation by a suitable current of, say, 0.1 - 1mA at 10 - 100kHz, and selection of pairs of the inner electrodes R, - R,6 for measurement of potential differences across them. Each of the electrodes R, to R,6 and S, to S,6 comprises a disc mounted with its plane normal to a radius of the annular support, each pair of electrodes in a cell being aligned radially, i.e. on a line extending normal to the surface of the body 12'. Each
30 electrode is about 6 mm in diameter. The electrodes R, to R,6 and S, to S16 in their respective compartments are spaced apart radially by a distance sufficient to reduce/avoid shielding effects. Other sizes of electrodes, for example from 2 mm to 10 mm diameter, are available. If they were used, the spacing would be adjusted appropriately. In view of the alternating "high" conductivity of compartments 32, to 32,6 and
"low" conductivity of separators 28, to 28,6, the annular support can be considered to have anisotropic conductivity, i.e. high conductivity for currents flowing radially from an electrode S towards the volume 12' and low conductivity for cuπents flowing circumferentially from one compartment towards its neighbour.
The manner in which this arrangement of electrode geometry and media favours the enhancement of the sensitivity of the apparatus will now be described.
Figure 5 illustrates the average sensitivity, in the radial direction, obtained using a finite element model, for an electrical impedance tomography apparatus using the annular electrode aπay of Figure 2, and that for an equivalent conventional electrical impedance tomography apparatus using the electrode array of Figure 1. In both cases the conductive volume 12' has a radius of 10 cm, with a uniform conductivity. For proof of concept, the annulus in Figure 2 extends from 10 cm to 15 cm and has a conductivity 10 times larger than the conductivity of conductive volume 12'. The medium separating the individual compartments, i.e. of separators 28, to 28,2, is set to a conductivity value of 1 % of the conductivity of conductivity volume 12'.
The average sensitivity curve obtained for the conventional electrode aπay of Figure 1 uses configurations of stimulator and recording electrode pairs as disclosed by Brown et al. in U.S. patent No. 4,617,939. As mentioned earlier, in this approach, a pair of neighbouring electrodes is used to stimulate while the electrical potential is measured with the remaining pairs of electrodes. The process is repeated until all the independent combinations are exhausted. The electrodes used for stimulation, however, are not used also for measurement.
The average sensitivity curve for the electrode configuration shown in Figure 2 is obtained using pairs of electrodes in the outer group S,-S,6 to apply stimulation cuπent and the coπesponding pairs of electrode in the inner group R,-R,6 to measure the resulting electrical potential.
The resolution of the reconstructed image depends upon the sensitivity of the apparatus to changes of the electric conductivity in the conductive volume 12'. The sensitivities depicted in Figure 5 were obtained using the sensitivity expression described by D.B. Geselowitz in "An Application of Electrocardiographic Lead Theory to Impedance Plethysmography", IEEE Trans. Biomed. Eng., Vol. BME-18, pp. 38-41, Jan. 1971, and by J. Lehr in "A Vector Derivation Useful in Impedance Plethysmographic field calculation", IEEE Trans. Biomed. Eng. , Vol. BME-19, pp. 156-157, Mar. 1972.
The theoretical basis for the reconstruction algorithm employed by embodiments of the invention will be described with reference to Figure 6. When the conductivity distribution changes from σ(x,y,z) to σ(x,y,z) + Δσ(x,y,z), the change in the transfer impedance ΔZ for the pair of stimulating electrodes (A,B) and recording electrodes
(C,D) can be given as:
Figure imgf000012_0001
where:
Φ is the potential distribution when the cuπent Iφ is applied to the electrode pair (A,B). is the potential distribution when the cuπent IΨ is applied to the electrode pair (C,D). σ(x,y,z) is the conductivity distribution when is established for Iψ σ(x,y,z,) + Δσ(x,y,z) is the conductivity distribution when is established for I .
The term ~ ψ (σ+Δσ) expression (1) can be expanded with respect to
Δσ. Assuming small Δσ's, the higher order terms can be neglected and expression (1) is expressed as follows:
Figure imgf000012_0002
Since the term ΔZ, the difference between the measured transfer impedance and the assumed transfer impedance, is known, the difference in the conductivity distribution Δσ can be obtained. Dividing the region of interest into small domains over which the conductivity is assumed constant enables expression (2) to be discretized and expressed as the following matrix equation:
ΔZ,. = ∑ 5, Δσ. (3)
where:
ΔZj is the difference in transfer impedance for the \th combination of a pair of stimulating electrodes and a pair of recording electrodes;
Δσj is the unknown conductivity difference for the region element j;
Sjj is the sensitivity coefficient term for the stimulating/recording electrode pair combination i, at the element j; and m is the number of domains into which the region has been divided.
The term S^ is given as follows:
Figure imgf000013_0001
An examination of expression (4) indicates that the sensitivity term can be enhanced by an apparatus that could maximize the integrand. Sensitivity enhancement is achieved by selecting configurations of stimulating and recording electrode pairs that favour the scalar product in the numerator and using an annulus medium with an anisotropic conductivity that limits the cuπent intensities at the denominator.
Electrode configurations. The annular electrode aπay of Figure 2 allows electrode configurations that can maximize scalar products of stimulating and recording electric field vectors more so than could be obtained with the electrodes of Figure 1. It permits more combinations of electrodes to be made, with a greater number of different geometries.
It should be appreciated that the stimulating electrodes need not be diametrically opposite each other. Referring to Figure 2, using pairs of electrodes in the outer group S, to S16 to stimulate the body and pairs of electrodes in the inner group R, to R,6 to measure the resulting electrical potentials tends to maximize the scalar product in the numerator of expression (4).
Annulus conductivity. The sensitivity of the apparatus to a change in the electrical conductivity at a position (x,y) in the conductive volume is an inverse function of the product of the total electric cuπent intensities that would be required to produce the electric field intensity of the stimulating and recording electrode pairs at position (x,y). Disposing the electrodes in sectors whose boundaries are poorly conductive, produces an annulus medium with an anisotropic conductivity, vis. conductive in the radial direction and poorly conductive in the direction. This aπangement forces the cuπent to flow in the medium of interest, i.e. the conductive volume, as opposed to the extended medium. It favours a larger ratio of the electric field to cuπent intensity, in the medium of interest. The sensitivity function used to determine the average sensitivities represented in
Figure 5 was not optimized for the various parameters that define the annulus. Key parameters are expected to be the following:
A. Annulus dimension: The size of the annulus used for this evaluation has an outer to inner radius ratio of 1.5. Evaluations have shown that as the ratio increases, and the electrodes are set further away from the imaging volume, the sensitivity decreases. Conversely, as the ratio decreases, and the electrodes get closer to the imaging medium and the sensitivity increases. Although this is true for a discrete point electrode model, the optimal practical ratio will be determined by the size of the electrodes and the minimum separation required to avoid mutual shielding.
B. Annulus electrical conductivity: Enhanced sensitivity is obtained when the annulus is configured with an anisotropic electrical conductivity distribution, that is, a large ratio of the radial to the azimuthal electrical conductivity. No significant improvement was gained by decreasing the azimuthal electrical conductivity below 1 % of the average imaging medium electric conductivity. C. Electrodes: A minimum of 32 electrodes were considered, arranged in two concentric rings. For this evaluation, each wedge-shaped compartment held two electrodes, aligned along the diameter. One electrode was positioned at the edge of the annulus and the other was positioned at the centre of the annulus. Although this was effective for the model, in practice, of course, the other electrode cannot be any closer to the centre than the inner surface of the annulus. It is expected that standard commercial AgCl disk shaped electrodes would be used for this application.
Measurement procedure. As stated earlier, it is desirable to maximize the scalar product of the electric field vectors for the stimulating pair of electrodes and the "recording" electric field vectors for the recording pair of electrodes. Consequently, for a given pair of stimulating electrodes, it is sufficient to measure the potential difference between that pair or pairs of electrodes which are closest spatially to the stimulating electrodes and not to measure the potential differences between other pairs of electrodes. This measurement procedure will now be described with reference to the flowchart of Figure 7. In order to illustrate the stimulation/measurement strategy explicitly, an example will be given using the geometry depicted in Figure 2. The outer electrode set (S,-S,6) are used as stimulating electrodes and the inner electrode set (R,-R,6) as recording electrodes. It should be appreciated that the outer set could be replaced with a set on the inner surface, as per the comments above.
Following positioning of the electrodes S,-S,6 and R,-R16 upon the conductive volume to be imaged, an initial conductivity distribution appropriate to the volume would be assumed (step 7.1). In the absence of other information from, say, previous measurements of the same or similar volumes, a uniform mean conductivity distribution would be assumed. For example, for landmine detection, this might be a typical conductivity distribution of the particular type of soil whereas, for human body imaging, it might be a typical conductivity distribution for that part of the human body.
In step 7.2, a series of measurements are made, passing a known cuπent through each pair of stimulation electrodes S, to S,6 in turn and recording the potential difference across the "closest" pair of electrodes, i.e. having the recording electric field vectors most closely aligned with those of the stimulating pair of electrodes. For example, when stimulating S, and S2, the potential difference between R, and R2 only would be measured.
Figure imgf000016_0001
Table 1 Stimulate/Measure Configuration This pattern yields 120 independent combinations for 16 compound electrodes
(stimulate-record pairs). It should be noted that the stimulating electrodes and the recording electrodes are always identically numbered, and hence closely located. Other stimulation/recording electrode configurations which share this feature will have some benefit. An example of such an alternative, possible with only a single set of electrodes as shown in Figure 1 (denoted E,,E2 ... E,6) is:
Figure imgf000016_0002
Table 2 Alternative Stimulate/Record Configuration
This pattern has a reduced number of independent combinations, since it starts with a stimulating electrode gap of 2, for example, but will show some of the benefit of the proposed configuration, and is shown as a simple example of an alternative manifestation of the stimulate/record configuration described herein.
Referring again to the flowchart of Figure 7, once the measurements have been made, in step 7.3, the region of interest, i.e. the region in which cuπent might flow, is divided into a mesh of elements. In Figure 2, the region of interest is the conductive volume 12', the sponge 36 and the annular aπay of compartments 32,-3216. Each element has dimensions equal to approximately 10 per cent of the spacing between each pair of the recording electrodes R, to R,6. The mesh could be triangular or rectangular, as desired, and the region could be two-dimensional or three-dimensional. In step 7.4, the sensitivity matrix Sy is calculated for the mesh using the electrode geometry, initial conductivity distribution, and applied cuπents, according to equation [4], for the region of interest.
In step 7.5, an estimate is made of the noise present in the potential difference measurements divided by the applied stimulation cuπent. This scaled noise value is used to determine a parameter μ which is the ratio of the scaled noise value to the root mean square (RMS) of the initial conductivity distribution.
In step 7.6, the pseudo-inverse S^ of the sensitivity matrix S;j is calculated according to equations (5) and (6).
In step 7.7, the differential conductivity distribution Δσ is calculated according to equation (7).
Preferred electrode configuration and compound electrodes.
Compound electrodes are so named because they comprise two parts, one part being used for stimulation and the other for recording. It has been found that compound electrodes, in conjunction with the novel measurement procedure described above, provide significantly better performance in physically realistic (noisy) situations. Analysis has shown that this improvement in performance is due both to sensitivity amplitude enhancements (due to alignment of the electric field stimulating and recording electric field vectors) and a better representation (more suitable for reconstruction) of the set of measurements possible by the proposed recording configuration.
A potential problem associated with compound electrodes is the possibility of a partial electrical short circuit between the stimulating and measuring electrode, due to the close proximity of these electrodes. Simulations have indicated that satisfaction will be obtained with spacing of the component electrodes of a compound electrode such as to render the voltage associated with this extra (surface) cuπent path approximately 0.1 % of that associated with the conventional path (through the region of interest). For a human body this requirement translates into a resistance of (for skin surface alone) approximately 5MΩ. An example of a suitable compound electrode is shown in Figures 8 and 9. The compound electrode comprises an annular outer electrode 46 suπounding a circular inner electrode 48 with an annular slot 50 between them. The electrodes 46 and 48 are bonded to an insulating planar support 52 and their exposed surfaces are coated with layers of conductive gel 54 and 56, respectively. An "O" ring 58 of insulating material is located in the annular slot 50. The diameter of the "O" ring 58 is approximately equal to the depth of the annular slot 50 so that, when the compound electrode is applied to the surface of the volume under investigation, the "O" ring 58 will prevent gel 54/56 from exuding across the annular slot 50 and causing, effectively, a short-circuit between the inner and outer electrodes. In use, a plurality of such compound electrodes would be distributed about the surface of the body being investigated. Cuπent could be applied to the respective outer electrodes of a pair of compound electrodes and the potential difference recorded between the inner electrodes of the same pair. Conversely, the inner electrodes could be stimulated and the potential difference recorded between the outer electrodes. It should be appreciated that other configurations of compound electrode could be used. For example, either or both of the electrodes could be rectangular. It is also envisaged that each compound electrode might comprise an aπay or grid of stimulation electrodes interspersed with an aπay or grid of recording electrodes.
The annular electrode aπay of Figures 2 and 3, whether with the individual electrodes or with the compound electrodes, is well suited for medical diagnostics applications applied to a patient's head, limb or torso. Nevertheless, the enhanced sensitivity approach illustrated with the annular geometry can also be applied to other electrode configurations. Thus, Figure 10 illustrates a two-dimensional aπay of electrode cells which is shown as rectangular but could have other shapes. This electrode array is envisaged for ground imaging applications, such as landmine detection. The array comprises juxtaposed individual electrode cell units 60. As shown in Figure 11, each cell comprises a pair of electrodes R, S immersed in a conductive material, for example a jelly formed from agar and saline solution, as before. The cell units 60 lie on a bed 62 (Figure 10) of conductive material that provides the electrical interface with the ground. This bed 62 would have insulating partitions (not shown) to maintain anisotropic conductivity.
Figure 12 illustrates another possible embodiment, where a flexible vest 64 comprises similar individual electrode cells 66 forming a cylindrical aπay which, as shown, can be wrapped around a human torso and used for three-dimensional imaging.
It will be appreciated that other geometries are also feasible, such as a special cap for head imaging applications or a special brassiere for breast imaging applications.
It should be appreciated that, if the radially-spaced electrodes S/R of Figures 2-4 and 10-12 were replaced by compound electrodes, the agar-fiUed compartments or cells would not be required, since the compound electrodes could lie flat upon the surface.
However, it is envisaged that, to facilitate calculations, the compound electrodes would still be in a rigid support, if necessary, with an intervening sponge 36, or the like.
Reconstruction technique
The proposed reconstruction technique is based on the construction of a probabilistically regulated, pseudo inverse of the sensitivity matrix introduced above and mentioned in step 7.4. The construction of this pseudo inverse is very straightforward, and follows the following steps :- Calculate the singular value decomposition of the sensitivity matrix.
S = U A Vτ (5)
Replace each singular value λ, in Λ = diag(λ,) with λ* =λ. / (λ + μ2),
producing Λ* = diag(λ*) ■ The parameter μ is the ratio of the noise in the measured potential differences to the root mean square (RMS) of the conductivity distribution. Thus, in this sense it is the inverse of the SNR. Form the pseudo inverse
S* = V A« UT (6) which can be used to compute conductivity perturbations from potential difference recordings, i.e. reconstruct the EIT image, since Δ σ-∑ S'ΔZ,, (7)
where n is the number of independent measurements taken.
This approach has been found to respect the inherent low-pass filtering of the EIT process, and offers a rapid and robust option for EIT image reconstruction. Figures 13 A to 13J illustrate a comparison between the novel measurement procedure set out in Table 1, using compound electrodes as per Figure 8, and the measurement procedure disclosed by Brown et al. in U.S. Patent No. 4,617,939, but in both cases using the reconstruction technique described above. Figures 13A to 13J depict the reconstruction of a perturbation covering approximately 10% of the diameter of a circular region of interest at three different locations. As shown in Figures 13A, 13D and 13G, with zero noise, both Brown et al. 's procedure and the procedure of the present invention reproduce the perturbation very accurately. Thus, Figure 13A shows a maximum a posteriori (MAP) estimate of a perturbation near the circumference of the region. Figure 13B shows the (MAP) estimate of the differential conductivity distribution Δσ for a circular aπay of 16 compound electrodes of the kind shown in Figure 8. Figure 13C shows the coπesponding MAP estimate for the conventional electrode configuration of Figure 1 using Brown et al. 's measurement technique.
Figures 13D to 13F coπespond to Figures 13 A to 13C for a perturbation which is closer to the middle of the region and Figures 13G to 13J are the coπesponding reconstructions for a perturbation at the centre of the region.
Although embodiments of the invention have been described and illustrated in detail, it is to be clearly understood that the same are by way of illustration and example only and not to be taken by way of the limitation. The spirit and scope of the present invention embrace various modifications and substitutions. For example, the electrodes could be in a fixed aπay and a separate cuπent generator provided for each pair stimulated. Alternatively, a single cuπent generator could be used to stimulate each pair of electrodes in turn, using a commutation switching circuit, enabling a more expensive cuπent generator to be used.
It is also envisaged that only two stimulating electrodes and two recording electrodes could be used. The various locations would be predefined and the electrodes repositioned at selected locations in turn until all locations had been stimulated as required.
It is envisaged that, instead of providing a personal computer with a DSP card for performing the reconstruction, a suitable custom integrated circuit, preprogrammed with the reconstruction algorithm software, could be combined with the electronic conditioning circuitry in an interface unit which would connect to a conventional personal computer and/or even directly to a display device.
INDUSTRIAL APPLICABILITY The compound electrode aπangement may be easier to manufacture. The measurement procedure facilitates simpler implementations for a given sensitivity.

Claims

CLAIMS:
1. A method of determining electrical impedance tomography of a conductive volume (12) comprising the steps of selecting, in turn, a plurality of pairs of locations (S) upon a surface of the volume (12) and supplying a stimulation cuπent to the surface by way of the selected pair of locations, and recording, for each selected pair, a resulting potential difference between at least one pair (R) of the remaining locations, and processing the recorded potential differences to plot electrical conductivity distribution within said volume characterized in that, for a particular pair of locations stimulated, the potential difference is recorded for a number of pairs of said locations that is less than the total of said plurality of locations and also comprises that pair of locations which, if stimulated, would produce an electric field with vectors most closely aligned with the coπesponding vectors of the electric field produced by stimulation of said particular pair of locations.
2. A method of determining electrical impedance tomography according to claim 1 , using a plurality of electrodes (S,-S16, R,-Ri6) disposed upon or adjacent said each surface at a respective one of said plurality of locations, characterized by the steps of;
(i) applying a stimulation cuπent to said volume by way of a first and a second of said plurality of electrodes at coπesponding first and second locations;
(ii) simultaneously recording the resulting potential difference between that pair of the remaining electrodes at said pair of locations which, if stimulated, would produce electric field vectors most closely aligned with the electric field vectors produced by stimulation of said particular pair of locations;
(iii) applying a stimulation cuπent to said body by way of two electrodes other than those stimulated in a previous step; (iv) simultaneously recording the resulting potential difference between that pair of the remaining electrodes at said locations which, if stimulated, would produce an electric field with vectors most closely aligned with the coπesponding vectors of the electric field produced by the stimulation cuπent applied in step (iii); (v) repeating steps (iii) and (iv) until all of said plurality of electrodes have been stimulated and the coπesponding potential differences recorded.
3. A method according to claim 2, characterized in that stimulation cuπent is applied to a said first and second of said plurality of electrodes with no intervening electrodes.
4. A method according to claim 2, characterized in that said stimulation cuπent is applied to a said first and second of said plurality of electrodes that are spaced apart with one or more other electrodes therebetween.
5. A method according to claim 4, characterized in that . each performance of steps (iii) and (iv) applies said stimulation cuπent to a pair of electrodes that are spaced apart with a different number of intervening electrodes therebetween.
6. A method according to claim 2, characterized in that, in said steps (i) and (ii), said stimulation cuπent is applied to a said first and second of said electrodes that have no electrodes therebetween, and said potential difference recorded, and, in said steps (iii) and (iv), a stimulation cuπent is applied to a said two electrodes that have one or more electrodes therebetween, and the resulting potential difference recorded.
7. A method according to claim 2, characterized in that the stimulating cuπents are applied to selected pairs of a first group of the plurality of electrode means and the resulting potential differences are recorded at selected pairs of a second group of the plurality of electrodes.
8. A method according to claim 7, characterized in that the electrodes of one of the groups are arranged in a ring closer to the surface than the electrodes of the other of the groups, each electrode of the first group being substantially aligned with a coπesponding electrode of the second group in a direction extending substantially normal to said surface and said potential difference is recorded between that pair of said electrodes which are aligned with the pair of electrodes to which the stimulation cuπent is applied.
9. A method according to claim 7, each electrode means comprising a stimulation electrode and a recording electrode closely located to each other and both the same distance from said surface, characterized in that the step of applying stimulation cuπent applies the stimulation cuπent to the stimulation electrodes of two electrode means,
5 respectively, and the coπesponding potential difference is measured between the recording electrodes of the same two electrode means.
10. An electrode aπangement for an electrical impedance tomography system comprising a first group of electrodes (RrR,6) and a second group of electrodes (S,-S,6)
10 and support means (24,26,28,30) for supporting the electrodes adjacent a surface of a volume the electrical impedance of which is to be mapped by the system, characterized in that the first group of electrodes is closer than the second group of electrodes to a part of the support means which, in use, will be adjacent said surface.
15 11. An electrode aπangement according to claim 10, characterized in that each electrode of the first group (R,-R,6) is substantially aligned with an adjacent electrode of the second group (S,-S,6) in a direction extending substantially normal, in use, to said surface.
20 12. An electrode aπangement according to claim 10 or 11, characterized in that the support means comprises an annular structure, the first group of electrodes and the second group of electrodes being aπanged in an inner ring and an outer ring, respectively, about a central part of the annular structure.
25 13. An electrode aπangement according to claim 12, characterized in that the annular structure is circular and the spacing between the inner ring of electrodes and the outer ring of electrodes is about one quarter of the diameter of the inner ring.
14. An electrode aπangement according to claim 10, characterized in that the support
30 means comprises a plurality of compartments (32,-32,6), each comprising a conductive medium for contacting the surface of a volume to be imaged, and lower conductivity partitions (28,-28,6) between adjacent compartments, each compartment housing a pair of electrodes comprising one electrode from said first group and one electrode from said second group.
15. An electrode aπangement according to claim 14, characterized in that the 5 electrodes are supported in a planar aπay (60,62), and each electrode of said first group is spaced from an associated electrode of said second group in a direction substantially perpendicular to the plane of the aπay.
16. An electrode aπangement according to any one of claims 10 to 15, characterized 10 in that the support means is rigid and a conductive interface medium (36) is disposed for juxtaposing between an inner surface of said support means and the surface of the volume (12) to be measured, the interface medium comprising a material that is conformable to contour variations of said surface of the volume.
15 17. An electrode arrangement according to claim 16, characterized in that the interface medium material has a conductivity comparable to that of the volume to be measured.
18. An electrode aπangement according to any one of claims 10 to 17, characterized 20 in that the support means has anisotropic conductivity with its conductivity in said direction normal to said surface, in use, being significantly higher than its conductivity in directions transverse thereto.
19. An electrode aπangement according to claim 18, characterized in that the support 25 means comprises a plurality of higher conductivity sectors alternating with, and defined by, a plurality of lower conductivity regions.
20. An electrode aπangement for an electrical impedance tomography system comprising a plurality of electrodes (S,-S16, R,-Rι6) in an aπay mounted in a support
30 (24,26,28,30) for supporting the electrodes adjacent a surface of a volume (12') to be imaged, characterized in that the support medium has anisotropic conductivity, its conductivity in a direction that, in use, is normal to the surface, being significantly greater than its conductivity in a transverse direction.
21. An electrode aπangement according to claim 20, characterized in that the support comprises a plurality of compartments (32,-32,6), each separated from neighbouring compartments by an insulating partition (28,-28,6), each compartment comprising a material having a conductivity comparable to the conductivity of the volume the
5 conductivity distribution of which is to be measured and containing at least one of the electrodes.
22. An electrode aπangement for an electrical impedance tomography system characterized by an aπay of electrodes fixed spatially relative to each other in a rigid
10 support, and an interface medium for interfacing the rigid support to a surface of a volume the conductivity distribution of which is to be measured, the interface medium being conductive so as to connect the electrodes electrically to the surface and pliable so as to conform to variations in relief of the surface.
15 23. An electrical impedance tomograph comprising an electrode aπangement according to any one of claims 10 to 22, further characterized by means (16,18,20) for stimulating, sequentially, pairs of electrodes and, for each pair, measuring coπesponding potential differences between pairs of the remaining electrodes, and processing the measured potentials to map electrical conductivity distribution of a volume.
PCT/CA1997/000876 1996-11-26 1997-11-25 Electrical impedance tomography method and electrode arrangement for use therein WO1998023204A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
AU50447/98A AU5044798A (en) 1996-11-26 1997-11-25 Electrical impedance tomography method and electrode arrangement for use therein
US09/308,521 US6501984B1 (en) 1996-11-26 1997-11-25 Electrical impedance tomography method and electrode arrangement for use therein
CA002272429A CA2272429C (en) 1996-11-26 1997-11-25 Electrical impedance tomography method and electrode arrangement for use therein

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CA2,191,285 1996-11-26
CA002191285A CA2191285A1 (en) 1996-11-26 1996-11-26 Electrode arrangement for electrical impedance tomography system

Publications (1)

Publication Number Publication Date
WO1998023204A1 true WO1998023204A1 (en) 1998-06-04

Family

ID=4159321

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CA1997/000876 WO1998023204A1 (en) 1996-11-26 1997-11-25 Electrical impedance tomography method and electrode arrangement for use therein

Country Status (4)

Country Link
US (1) US6501984B1 (en)
AU (1) AU5044798A (en)
CA (1) CA2191285A1 (en)
WO (1) WO1998023204A1 (en)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1999004689A1 (en) * 1997-07-25 1999-02-04 Werner Alois Kaiser Device for detecting, characterizing by differential diagnosis, and treating tumours
WO2002008794A2 (en) * 2000-07-26 2002-01-31 Wisys Technology Foundation, Inc. Method and apparatus for producing an electrical property image using a charge correlation matrix
EP1505906A1 (en) * 2002-05-20 2005-02-16 Richard J. Davies System for detecting precancerous and cancerous tissue
GB2396426B (en) * 2002-12-21 2005-08-24 Draeger Medical Ag Artificial respiration system
US6940286B2 (en) 2000-12-30 2005-09-06 University Of Leeds Electrical impedance tomography
WO2009018620A1 (en) * 2007-08-09 2009-02-12 Impedimed Limited Impedance measurement process
EP1952764A3 (en) * 2003-11-25 2009-05-06 University-Industry Cooperation Group of Kyung Hee University System and method for visualizing conductivity and current density distribution in object
US7630759B2 (en) 2002-05-20 2009-12-08 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous breast tissue and epithelium
US7853319B2 (en) 2005-04-21 2010-12-14 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue and epithelium
GB2471873A (en) * 2009-07-15 2011-01-19 Wzvi Ltd Electrode array arranged on a triangular grid for electrical impedance imaging
US8103337B2 (en) 2004-11-26 2012-01-24 Impedimed Limited Weighted gradient method and system for diagnosing disease
US8262575B2 (en) 2002-05-20 2012-09-11 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue
US8332025B2 (en) 2005-12-06 2012-12-11 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue and epithelium
RU2470580C1 (en) * 2011-06-08 2012-12-27 Государственное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный электротехнический университет "ЛЭТИ" им. В.И. Ульянова (Ленина)" Method of determining electric resistance of internal tissues of part of biological object body and rheoanalyser
US8738124B2 (en) 2007-12-11 2014-05-27 Epi-Sci, Llc Electrical bioimpedance analysis as a biomarker of breast density and/or breast cancer risk
US8836345B2 (en) 2007-11-05 2014-09-16 Impedimed Limited Impedance determination
US9149235B2 (en) 2004-06-18 2015-10-06 Impedimed Limited Oedema detection
EP3026430A1 (en) * 2014-11-28 2016-06-01 Commissariat à l'Énergie Atomique et aux Énergies Alternatives Method for imaging a medium by electrical measurements with contact impedance correction
US9392947B2 (en) 2008-02-15 2016-07-19 Impedimed Limited Blood flow assessment of venous insufficiency
US9504406B2 (en) 2006-11-30 2016-11-29 Impedimed Limited Measurement apparatus
US9585593B2 (en) 2009-11-18 2017-03-07 Chung Shing Fan Signal distribution for patient-electrode measurements
US9615767B2 (en) 2009-10-26 2017-04-11 Impedimed Limited Fluid level indicator determination
US9615766B2 (en) 2008-11-28 2017-04-11 Impedimed Limited Impedance measurement process
US9724012B2 (en) 2005-10-11 2017-08-08 Impedimed Limited Hydration status monitoring
US10307074B2 (en) 2007-04-20 2019-06-04 Impedimed Limited Monitoring system and probe
US11660013B2 (en) 2005-07-01 2023-05-30 Impedimed Limited Monitoring system
US11737678B2 (en) 2005-07-01 2023-08-29 Impedimed Limited Monitoring system

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7729756B2 (en) * 2000-01-18 2010-06-01 Siemens Aktiengesellschaft Measurement system for examining a section of tissue on a patient and the use of a measurement system of this type
AU2001283849B2 (en) * 2000-06-09 2005-03-17 Timpel S.A. Method and apparatus for displaying information obtained by electrical impedance tomography data
CA2363821A1 (en) * 2000-11-24 2002-05-24 Alvin Wexler High definition electrical impedance tomography methods for the detection and diagnosis of early stages of breast cancer
US6823212B2 (en) * 2001-06-13 2004-11-23 The Procter & Gamble Company Method and apparatus for measuring properties of a target surface
WO2004049936A2 (en) * 2002-11-29 2004-06-17 Z-Tech (Canada) Inc. Improved breast electrode array and method of analysis for detecting and diagnosing diseases
WO2004062464A2 (en) * 2003-01-10 2004-07-29 Kohwang Foundation, Kohwang Board Of Trustee System and method for three-dimensional visualization of conductivity and current density distribution in electrically conducting object
EP1613212B1 (en) * 2003-03-25 2016-08-24 Fresenius Medical Care Holdings, Inc. Device and method for performing electrical impedance tomography
US20040243019A1 (en) * 2003-03-27 2004-12-02 Z-Tech (Canada) Inc. Weighted gradient method and system for diagnosing disease
DE10315863B4 (en) * 2003-04-08 2013-03-14 Dräger Medical GmbH electrode belt
CN100374076C (en) * 2003-09-23 2008-03-12 中国人民解放军第四军医大学 Electric impedance tomographic imaging method and apparatus for imaging monitoring beside sickbed
US20080076998A1 (en) * 2003-12-01 2008-03-27 Z-Tech (Canada) Inc. Breast electrode array and method of analysis for detecting and diagnosing diseases
US6965836B2 (en) * 2004-04-19 2005-11-15 Battelle Energy Alliance, Llc Method and apparatus for two dimensional surface property analysis based on boundary measurement
GB0516158D0 (en) * 2005-08-05 2005-09-14 Univ Montfort An apparatus and method for `non-contact' electrical impedance imaging
WO2007035887A2 (en) * 2005-09-21 2007-03-29 Beth Israel Deaconess Medical Center, Inc. Electrical impedance myography
ATE512693T1 (en) * 2007-02-28 2011-07-15 St Jude Medical HEART STIMULATOR WITH AN IMPEDANCE MEASUREMENT CIRCUIT FOR DETECTING ISCHEMIA
BRPI0704023A8 (en) * 2007-10-03 2015-11-10 Dixtal Biomedica Ind E Comercio Ltda EQUIPMENT FOR ELECTRICAL IMPEDANCE TOMOGRAPHY
US8102182B2 (en) * 2007-10-11 2012-01-24 University Of Utah Research Foundation Systems and methods for measuring the electrical properties of a microparticle
WO2009118701A1 (en) * 2008-03-27 2009-10-01 Koninklijke Philips Electronics N.V. Method and system for measuring an object of interest
US8386010B2 (en) * 2008-10-23 2013-02-26 Covidien Lp Surgical tissue monitoring system
US8914101B2 (en) * 2008-10-30 2014-12-16 Medtronic, Inc. System and method to localize changes in intrathoracic fluid content using measured impedance in an implantable device
US8963562B2 (en) 2011-08-31 2015-02-24 General Electric Company Transducer configurations and methods for transducer positioning in electrical impedance tomography
US9149225B2 (en) 2011-12-14 2015-10-06 Intesection Medical, Inc. Methods for determining the relative spatial change in subsurface resistivities across frequencies in tissue
WO2014201163A1 (en) 2013-06-11 2014-12-18 Somatis Sensor Solutions LLC Systems and methods for sensing objects
US9804112B2 (en) 2013-11-20 2017-10-31 Transtech Systems, Inc. Selective characterization of material under test (MUT) with electromagnetic impedance tomography and spectroscopy
US10357177B2 (en) 2013-12-13 2019-07-23 General Electric Company Systems and methods for electrical impedance imaging
EP3785624A1 (en) * 2014-02-25 2021-03-03 School Juridical Person The Kitasato Institute Image generation apparatus, conductivity acquisition apparatus, image generation method, and program
DE102015006902B3 (en) * 2015-06-04 2016-06-30 Drägerwerk AG & Co. KGaA Device for processing and visualizing data from an electro-impedance tomography device for the detection and visualization of regional delays in ventilation in the lungs
CN105748072B (en) * 2016-02-01 2019-01-18 上海交通大学 A kind of distributed medical electrical impedance imaging measuring system of high-precision multi-frequency and method
EP3638465A4 (en) 2017-06-15 2021-07-07 OnRobot A/S Systems, devices, and methods for sensing locations and forces
CN108709909B (en) * 2018-06-25 2023-06-30 中国地质大学(武汉) Electrode device based on EIT nondestructive test cement-based flat plate
GB2580164A (en) * 2018-12-21 2020-07-15 Imperial College Sci Tech & Medicine A sensor

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2138148A (en) * 1983-04-13 1984-10-17 Denis Nigel Smith Method and apparatus for deriving currents and potentials representative of the impedance of zones of a body
US4617939A (en) 1982-04-30 1986-10-21 The University Of Sheffield Tomography
WO1993000038A1 (en) * 1991-06-27 1993-01-07 British Technology Group Limited Applied potential tomography
US5381333A (en) 1991-07-23 1995-01-10 Rensselaer Polytechnic Institute Current patterns for electrical impedance tomography
US5465730A (en) 1991-07-26 1995-11-14 British Technology Group Ltd. Electrical impedance tomography
WO1997040743A1 (en) * 1996-04-29 1997-11-06 Dmitry Viktorovich Nikolaev Method for local measurement of bioimpedance and device for implementing the same

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9013177D0 (en) 1990-06-13 1990-08-01 Brown Brian H Real-time imaging, etc.
US5284142A (en) 1991-12-16 1994-02-08 Rensselaer Polytechnic Institute Three-dimensional impedance imaging processes
US5351697A (en) 1991-12-16 1994-10-04 Rensseleaer Polytechnic Institute Three-dimensional impedance imaging process
WO1994012947A1 (en) 1992-11-20 1994-06-09 British Technology Group Limited Image reconstruction
AU699170B2 (en) * 1994-03-11 1998-11-26 British Technology Group Limited Electrical impedance tomography
US5919142A (en) * 1995-06-22 1999-07-06 Btg International Limited Electrical impedance tomography method and apparatus

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4617939A (en) 1982-04-30 1986-10-21 The University Of Sheffield Tomography
GB2138148A (en) * 1983-04-13 1984-10-17 Denis Nigel Smith Method and apparatus for deriving currents and potentials representative of the impedance of zones of a body
WO1993000038A1 (en) * 1991-06-27 1993-01-07 British Technology Group Limited Applied potential tomography
GB2257530A (en) 1991-06-27 1993-01-13 British Tech Group Applied potential tomography
US5381333A (en) 1991-07-23 1995-01-10 Rensselaer Polytechnic Institute Current patterns for electrical impedance tomography
US5465730A (en) 1991-07-26 1995-11-14 British Technology Group Ltd. Electrical impedance tomography
WO1997040743A1 (en) * 1996-04-29 1997-11-06 Dmitry Viktorovich Nikolaev Method for local measurement of bioimpedance and device for implementing the same

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BOOTH ET AL.: "A comparison of Three Electrode Configurations for Electrical Impedance Tomography. IEE Colloquium on "Innovation in Instrumentation for Electrical Tomography.", 1 January 1995 (1995-01-01), pages 11/1 - 3
E.J. WOO ET AL.: "Skin impedance measurements using simple and compound electrodes", MEDICAL & BIOLOGICAL ENGINEERING & COMPUTING, vol. 30, no. 1, January 1992 (1992-01-01), STEVENAGE (GB), pages 97 - 102, XP000246244 *
P. HUA ET AL.: "An electrical Impedance Tomograph using Compound Electrodes - presented at the IEEE Engineering in Medicine & Biology Societey 11th Annual Int. Conference 1989.", 1 January 1989 (1989-01-01)
PING HUA ET AL.: "Using Compound Electrodes in Electrical Impedance Tomography", IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, vol. 40, no. 1, - January 1993 (1993-01-01), NEW YORK (US), pages 29 - 34, XP002062243 *

Cited By (42)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6363275B1 (en) 1997-07-25 2002-03-26 Werner Alois Kaiser Device for detecting, for characterizing by differential diagnosis, and treating tumors
WO1999004689A1 (en) * 1997-07-25 1999-02-04 Werner Alois Kaiser Device for detecting, characterizing by differential diagnosis, and treating tumours
WO2002008794A2 (en) * 2000-07-26 2002-01-31 Wisys Technology Foundation, Inc. Method and apparatus for producing an electrical property image using a charge correlation matrix
WO2002008794A3 (en) * 2000-07-26 2002-04-18 Wisys Technology Found Inc Method and apparatus for producing an electrical property image using a charge correlation matrix
US6763263B2 (en) 2000-07-26 2004-07-13 Wisys Technology Foundation Method and apparatus for producing an electrical property image using a charge correlation matrix
US6940286B2 (en) 2000-12-30 2005-09-06 University Of Leeds Electrical impedance tomography
US9037228B2 (en) 2002-05-20 2015-05-19 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue
US8262575B2 (en) 2002-05-20 2012-09-11 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue
US6922586B2 (en) * 2002-05-20 2005-07-26 Richard J. Davies Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue
US9042977B2 (en) 2002-05-20 2015-05-26 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous breast tissue and epithelium
EP1505906A4 (en) * 2002-05-20 2009-04-08 Richard J Davies System for detecting precancerous and cancerous tissue
US8280504B2 (en) 2002-05-20 2012-10-02 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue
US8275453B2 (en) 2002-05-20 2012-09-25 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous breast tissue and epithelium
EP1505906A1 (en) * 2002-05-20 2005-02-16 Richard J. Davies System for detecting precancerous and cancerous tissue
US7630759B2 (en) 2002-05-20 2009-12-08 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous breast tissue and epithelium
GB2396426B (en) * 2002-12-21 2005-08-24 Draeger Medical Ag Artificial respiration system
EP1952764A3 (en) * 2003-11-25 2009-05-06 University-Industry Cooperation Group of Kyung Hee University System and method for visualizing conductivity and current density distribution in object
US9149235B2 (en) 2004-06-18 2015-10-06 Impedimed Limited Oedema detection
US8103337B2 (en) 2004-11-26 2012-01-24 Impedimed Limited Weighted gradient method and system for diagnosing disease
US7853319B2 (en) 2005-04-21 2010-12-14 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue and epithelium
US11737678B2 (en) 2005-07-01 2023-08-29 Impedimed Limited Monitoring system
US11660013B2 (en) 2005-07-01 2023-05-30 Impedimed Limited Monitoring system
US9724012B2 (en) 2005-10-11 2017-08-08 Impedimed Limited Hydration status monitoring
US11612332B2 (en) 2005-10-11 2023-03-28 Impedimed Limited Hydration status monitoring
US8332025B2 (en) 2005-12-06 2012-12-11 Epi-Sci, Llc Method and system for detecting electrophysiological changes in pre-cancerous and cancerous tissue and epithelium
US9504406B2 (en) 2006-11-30 2016-11-29 Impedimed Limited Measurement apparatus
US10307074B2 (en) 2007-04-20 2019-06-04 Impedimed Limited Monitoring system and probe
WO2009018620A1 (en) * 2007-08-09 2009-02-12 Impedimed Limited Impedance measurement process
AU2008286194B2 (en) * 2007-08-09 2014-05-15 Impedimed Limited Impedance measurement process
US8836345B2 (en) 2007-11-05 2014-09-16 Impedimed Limited Impedance determination
US8738124B2 (en) 2007-12-11 2014-05-27 Epi-Sci, Llc Electrical bioimpedance analysis as a biomarker of breast density and/or breast cancer risk
US9392947B2 (en) 2008-02-15 2016-07-19 Impedimed Limited Blood flow assessment of venous insufficiency
US9615766B2 (en) 2008-11-28 2017-04-11 Impedimed Limited Impedance measurement process
GB2471873B (en) * 2009-07-15 2015-09-23 Wzvi Ltd Electrical impedance imaging
US8928332B2 (en) 2009-07-15 2015-01-06 Wzvi Limited Electrical impedance imaging
GB2471873A (en) * 2009-07-15 2011-01-19 Wzvi Ltd Electrode array arranged on a triangular grid for electrical impedance imaging
US9615767B2 (en) 2009-10-26 2017-04-11 Impedimed Limited Fluid level indicator determination
US9585593B2 (en) 2009-11-18 2017-03-07 Chung Shing Fan Signal distribution for patient-electrode measurements
RU2470580C1 (en) * 2011-06-08 2012-12-27 Государственное образовательное учреждение высшего профессионального образования "Санкт-Петербургский государственный электротехнический университет "ЛЭТИ" им. В.И. Ульянова (Ленина)" Method of determining electric resistance of internal tissues of part of biological object body and rheoanalyser
US9933380B2 (en) 2014-11-28 2018-04-03 Commissariat à l'énergie atomique et aux énergies alternatives Method for imaging a medium through electrical measurements with a contact impedance correction
FR3029287A1 (en) * 2014-11-28 2016-06-03 Commissariat Energie Atomique METHOD FOR IMAGING A MEDIUM BY ELECTRICAL MEASUREMENTS WITH CONTACT IMPEDANCE CORRECTION
EP3026430A1 (en) * 2014-11-28 2016-06-01 Commissariat à l'Énergie Atomique et aux Énergies Alternatives Method for imaging a medium by electrical measurements with contact impedance correction

Also Published As

Publication number Publication date
US6501984B1 (en) 2002-12-31
AU5044798A (en) 1998-06-22
CA2191285A1 (en) 1998-05-26

Similar Documents

Publication Publication Date Title
US6501984B1 (en) Electrical impedance tomography method and electrode arrangement for use therein
EP2073702B1 (en) Data collection for electrical impedance tomography
US4617939A (en) Tomography
McEwan et al. A review of errors in multi-frequency EIT instrumentation
Geselowitz On bioelectric potentials in an inhomogeneous volume conductor
US20100148798A1 (en) Apparatus and Method for Electrical Impedance Imaging
Price Electrical impedance computed tomography (ICT): A new CT imaging technique
McFEE et al. Research in electrocardiography and magnetocardiography
JPH08504632A (en) Tomography
WO1992007509A1 (en) Method and apparatus for imaging electrical activity in a biologocal system
Thevenet et al. The finite element method for a realistic head model of electrical brain activities: preliminary results
US20050054911A1 (en) Current density impedance imaging (CDII)
CA1207083A (en) Applied potential tomography
Jossinet et al. Quantitative technique for bio-electrical spectroscopy
Gracia et al. Multilead measurement system for the time-domain analysis of bioimpedance magnitude
Dimas et al. An efficient point-matching method-of-moments for 2D and 3D electrical impedance tomography using radial basis functions
CA2272429C (en) Electrical impedance tomography method and electrode arrangement for use therein
Kotre Subsurface electrical impedance imaging using orthogonal linear electrode arrays
CN209847158U (en) Electrical impedance imaging apparatus
WO2022258963A1 (en) Impedance tomography
US20140095094A1 (en) Imaging Systems and Methods
Xu et al. The acquisition hardware system with direct digital synthesis and filtered back-projection imaging in electrical impedance tomography
Hyttinen et al. Lead field of ECG leads calculated by a computer thorax model-an application of reciprocity
Mason et al. Non-invasive imaging of neural activity with magnetic detection electrical impedance tomography (MDEIT): a modelling study
Rusho et al. Reconstruction algorithm for pigeon hole imaging (PHI)

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AU CA CZ JP NO US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
ENP Entry into the national phase

Ref document number: 2272429

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 09308521

Country of ref document: US

122 Ep: pct application non-entry in european phase