METHOD OF MEASURING BY ELECTRICAL IMPEDANCE TOMOGRAPHY

Abstract

An electrical impedance tomography method for the measurement of a body comprising a cylindrical part containing a fluid, the method comprising arranging a number of electrodes around a periphery of the cylindrical part of the body, simultaneously exciting each of the electrodes, each electrode being excited by a potential of a selected form, measuring the electrical properties of the body using electrodes, and processing the data from measuring step so as to obtain a signed data matrix representative of an image.

Claims

1. An electrical-impedance-tomography measuring method of a body comprising a cylindrical portion containing a fluid, the method comprising: arranging a number n.sub.e of electrodes around a periphery of the cylindrical portion of the body, simultaneously exciting each of the n.sub.e electrodes, each electrode being excited by a potential V.sub.n.sup.exc having a form: $V_n^{\mathrm{exc}}\left(t\right)=A\underset{m=1}{\overset{n_e}{.Math.}}{}_{\cos }\left(2?f_{m}t\right)\left[{?}_{m^{\cos }}^{?}\left(m{?}_n\right)+{?}_{m^{\sin }}^{?}\left(\frac{m{?}_n}{2}\right)\right]$ where A is a signal amplitude, ?.sub.n is an angular position of electrode n, f.sub.m=m*f.sub.0 is an oscillation frequency, f.sub.0 is a fundamental frequency chosen such that f.sub.m is less than a Nyquist frequency of the system for all m, ? is the Kronecker delta, and
={2k+1:k?} is the set of odd integers
={2k:k?*} is the set of non-zero even integers. measuring electrical properties V.sub.n.sup.meas of the body using the electrodes, and processing data generated in the measuring step, by: a) for each electrode E.sub.n, computing data points M.sub.n defined by: $M_n\left(k\right)=\frac{1}{\mathrm{RP}}.Math.\underset{p=0}{\overset{P-1}{.Math.}}{V}_n^{\mathrm{meas}}\left(p\right)e^{\mathrm{ik}{?}_p}.Math.$ where R is a resistance of a resistor used to measure V.sub.n.sup.meas with V.sub.n.sup.meas=R I.sub.n across terminals of the resistor, P is a number of points in a discrete sequence of measurement of the current I.sub.n, p is a discrete time, k is a Fourier coefficient comprised between 1 and (n.sub.e?1) and ?.sub.p=(2?p/P), b) constructing a data matrix D from the data points M.sub.n(k) for all n and for all k, according to the equation: $D=\left(\begin{array}{c}\left\{M_n\left(1\right)\right\}\\ \left\{M_n\left(2\right)\right\}\\ \left\{M_n\left(3\right)\right\}\\ .Math.\\ \left\{M_n\left(n_e-1\right)\right\}\end{array}\right)\left[[.]\right],$ and c) constructing a signed data matrix, the elements of which are defined by the following equation when a phase shift ?.sub.n,l(k) between t an he excitation potential at electrode l and a current measured at electrode n is less than ?/2: ${\tilde{D}}_n^m=\frac{\sin \left({?}_{n,l}\left(k\right)\right)}{.Math.\sin \left({?}_{n,l}\left(k\right)\right).Math.}{D}_n^m\left[[.]\right]$ and the elements of which are defined by the following equation when the phase shift ?.sub.n,l(k) between the excitation potential at electrode l and the current measured at electrode n is greater than or equal to ?/2:
{tilde over (D)}.sub.n.sup.m=?.sub.n.sup.m where ? is a sign matrix defined such that an i-th element of a j-th row of ? is the sign of cosine([2?/([j 1]/2)]*(i?1)/n.sub.e) for odd j and a sign of sine([2?/(j/2)]*(i?1)/n.sub.e) for even j.

2. The method as claimed in claim 1, wherein the set of potentials V.sub.n.sup.exc meeting satisfy the condition: $\underset{n=1}{\overset{n_e}{.Math.}}{V}_n^{\mathrm{exc}}\left(t\right)=0.$

3. The method as claimed in claim 2, wherein an image is produced using a one-step iterative least-squares reconstruction algorithm applied to the signed data matrix.

4. The method as claimed in claim 1, wherein the electrodes are angularly distributed in a regular manner around a periphery of the body.

5. The method as claimed in claim 1, further comprising performing tomographic measurement of a two-phase flow, the body being a pipe of a nuclear installation.

6. A computer program product comprising a non-transitory computer-readable medium storing instructions that are readable by a processor so that, when executed, the instructions cause an acquiring system to be controlled in order to implement the measuring method as claimed in claim 1.

7. A device for implementing the method as claimed in claim 1, the device comprising: an acquiring system comprising at least one programmable logic array, a module for generating analog signals, and a module for measuring analog signals; a computer configured to control the acquiring system; and a plurality of electrodes connected to the acquiring system.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0053] FIG. 1 shows a device for implementing the measuring method according to the invention.

[0054] FIG. 2 illustrates spatial cosine patterns.

[0055] FIG. 3 shows one portion of an electronic circuit allowing an electrode to be excited.

[0056] FIG. 4 illustrates a method for generating excitation signals.

[0057] FIG. 5 illustrates a method for measuring the signals generated by the electrodes.

[0058] FIG. 6 shows, in graph form, excitation signals and certain of their properties.

[0059] FIG. 7 illustrates a sign matrix for a device comprising 16 electrodes.

[0060] FIG. 8 illustrates a sign matrix for a device comprising 32 electrodes.

DETAILED DESCRIPTION

[0061] FIG. 1 illustrates a device 1 for implementing the EIT measuring method according to the invention.

[0062] The method aims to take, by EIT, measurements of a body 6 containing a fluid.

[0063] The device 1 comprises 16 electrodes (referenced 2) arranged non-intrusively on a circular periphery of the body 6. The electrodes 2 are preferably angularly distributed in a regular manner around the periphery of the body 6.

[0064] The electrodes 2 are connected to a printed circuit board 3, itself connected to a data-acquiring system 4. A screen 5 allows the data and images produced from these data to be viewed. The data-acquiring system 4 contains a Linux operating system (HOST), which controls a programmable logic array (FPGA) also contained in the data-acquiring system 4.

[0065] The acquiring system 4 allows analog excitation signals to be generated and the analog measurement signals delivered by the electrodes 2 to be measured.

[0066] The system 4 for example comprises a cRIO-9039 controller from the manufacturer National Instruments, which includes a programmable logic network, NI-9262 modules from the manufacturer National Instruments for generating analog excitation signals and NI-9223 modules from the manufacturer National Instruments for measuring the analog signals delivered by the electrodes 2.

[0067] Implementation of the measuring method according to the invention with a device such as shown in FIG. 1 will now be described.

[0068] In a first step of the measuring method, a number ne of electrodes 2 is placed around a periphery of a cylindrical portion of the body 6.

Exciting the Electrodes

[0069] In a second step of the measuring method, the electrodes are simultaneously excited by a potential having a well-chosen form.

[0070] The electrodes 2 form a set of linearly independent electrodes. They are used to create electrical excitations on the surface of the body, and to measure the electrical properties thereof.

[0071] For a number n.sub.e of electrodes, there are (n.sub.e?1) linearly independent excitation patterns. To describe these linearly independent patterns it is natural to choose Fourier basis functions, according to the equation:

[00006]$\begin{array}{cc}\frac{1}{\sqrt{2?}}{e}^{\mathrm{in}?}& \left[\mathrm{Math}2\right]\end{array}$

[0072] The method according to the invention consists in simultaneously exciting all of the electrodes using excitations of trigonometric form.

[0073] This set of simultaneous excitations is broken down into spatial oscillations in the Fourier basis and into temporal oscillations.

[0074] Various frequencies are employed in order to establish a distinction between the various frequency-multiplexed trigonometric signals.

[0075] For each trigonometric excitation pattern, each electrode E.sub.n is associated with one static voltage V.sub.n.sup.sta.

[0076] The set of n.sub.e static voltages V.sub.n.sup.sta forms a set of sinusoidal and cosinusoidal functions having various spatial frequencies m.

[0077] FIG. 2 shows the spatial cosinusoidal patterns for m varying from 1 to 5. The sinusoidal patterns have not been shown.

[0078] The n.sub.e electrodes are arranged on the body's periphery ??, which has been represented by the dashed lines in FIG. 2. The continuous lines represent the excitation potential applied to the electrodes.

[0079] The static potential V.sub.n,m.sup.sta associated with each electrode E.sub.n is defined by the following equation, for n comprised between 1 and n.sub.e:

[00007]$\begin{array}{cc}{V}_{n,m}^{\mathrm{sta}}=A\left[{?}_m^{?}\cos \left(m{?}_n\right)+{?}_m^{?}\sin \left(\frac{m{?}_n}{2}\right)\right],& \left[\mathrm{Math}3\right]\end{array}$

where:

[Math 3a] [0080] [0081] where m?1, . . . , (n.sub.e?1) represents spatial frequency, [0082] ={2k+1:k?} is the set of odd integers, [0083] ={2k:k?N*} is the set of non-zero even integers, [0084] ? is the Kronecker delta with =1 if m? and =0 otherwise, [0085] ?.sub.n=2?i.sub.e is the angular coordinate of electrode E.sub.n, [0086] A is the amplitude of the applied voltage.

[0087] For a given set of n.sub.e electrodes, all the independent excitation patterns are fully described with (n.sub.e?1) different spatial frequencies.

[0088] For the particular case where ?.sub.n=2?n/n.sub.e, the electrodes are regularly distributed over the periphery of the body.

[0089] Each spatial frequency m is associated with a temporal frequency f.sub.m and is applied simultaneously to each of the electrodes.

[0090] Thus, each simultaneous excitation potential V.sub.n.sup.exc comprises a superposition of (n.sub.e?1) trigonometric functions, each function oscillating at a particular frequency f.sub.m. The excitation signal of an electrode n is V.sub.n.sup.exc, as defined by the following equation:

[00008]$\begin{array}{cc}\begin{array}{c}{V}_n^{\mathrm{exc}}\left(t\right)=\underset{m=1}{\overset{n}{.Math.}}\cos \left(2?f_{m}t\right)V_{n,m}^{\mathrm{sta}}\\ =A\underset{m=1}{\overset{n_e}{.Math.}}\cos \left(2?f_{m}t\right)\left[{?}_m^{?}\cos \left(m{?}_n\right)+{?}_m^{?}\sin \left(\frac{m{?}_n}{2}\right)\right]\end{array}.& \left[\mathrm{Math}4\right]\end{array}$

[0091] The method according to the invention comprises a step consisting in simultaneously applying each potential V.sub.n.sup.exc(t) thus defined to each electrode E.sub.n.

[0092] Care is taken to ensure that the sum of the excitation voltages of the n.sub.e electrodes is zero whatever the time t. This is expressed by the condition:

[00009]$\begin{array}{cc}\underset{n=1}{\overset{n_e}{.Math.}}{V}_n^{\mathrm{exc}}\left(t\right)=0.& \left[\mathrm{Math}5\right]\end{array}$

Measuring the Electrical Properties of the Body

[0093] In a third step of the measuring method, the electrical properties of the body 6 are measured using the electrodes 2.

[0094] A computer controls the programmable logic network of the acquiring system 4 so as to generate 16 excitation signals having the described properties. These 16 digital signals are converted into analog signals by the NI-9262 modules and transmitted by coaxial cables to the electrodes 2.

[0095] The printed circuit board 3 comprises one exciting circuit for each electrode 2, each of these circuits comprising a resistor R as illustrated in FIG. 3. As may be seen in this figure, the potential V.sub.n.sup.exc is applied to one side of the resistor R, the other side being connected to the electrode E.sub.n.

[0096] The Neumann boundary condition at electrode E.sub.n is the current I.sub.n passing through the exciting circuit. This current is obtained by measuring the voltage V.sub.n.sup.meas=R I.sub.n across the terminals of the resistor. As for V.sub.n.sup.exc, this signal is a sum of trigonometric functions.

Processing the Data

[0097] In a fourth step of the measuring method, the data measured in the third step are processed in order to obtain a signed data matrix representative of an image.

Generating Data Points M.SUB.n

[0098] In a first sub-step of the fourth step of the measuring method, data points M.sub.n are computed for each electrode n.

[0099] The Fourier transform of the measured signal V.sub.n.sup.meas is computed from a P-point current measurement sequence I.sub.n(p), where p is discrete time and respects 0?p?P, namely:

[00010]$\begin{array}{cc}{X}_n\left(k\right)=?\left[I_n\left(p\right)\right]=\frac{?}{R}\underset{p=0}{\overset{P-1}{.Math.}}{V}_n^{\mathrm{meas}}\left(p\right)e^{\mathrm{ik}{?}_p}\mathrm{where}{?}_p=\frac{2?p}{P},i=\sqrt{-1}\mathrm{and}\mathrm{the}\mathrm{normalization}\mathrm{factor}?=\frac{1}{P}.& \left[\mathrm{Math}6\right]\end{array}$

[0100] The Fourier transform may be computed at a frequency f.sub.1, which corresponds to the frequency at which P Fourier coefficients are computed.

[0101] The frequencies f.sub.m of the voltages V.sub.n.sup.exc are chosen to be harmonics of f.sub.1. This makes it possible to distinguish between the measured signals. Thus, each coefficient k is associated with one particular frequency f.sub.m.

[0102] The data are thus generated at a frequency f.sub.1 and the resolution in Fourier space is ?f=f.sub.m+1?f.sub.m=f.sub.1. It should be noted that the highest frequency is chosen to be lower than the Nyquist frequency of the system f.sub.Nyq=??p, where ?p is the sampling time.

[0103] The data points M.sub.n(k) are the moduli of each Fourier coefficient k for each electrode n:

[00011]$\begin{array}{cc}{M}_n\left(k\right)=\frac{1}{\mathrm{RP}}.Math.\underset{p=0}{\overset{P-1}{.Math.}}{V}_n^{\mathrm{meas}}\left(p\right)e^{\mathrm{ik}{?}_p}.Math.& \left[\mathrm{Math}7\right]\end{array}$

[0104] Each data point defines the current of a given trigonometric pattern at a given electrode. The set of data points M.sub.n(k) for all n and all k forms the measurement data.

[0105] The excitation frequencies are determined in the following way. The sampling frequency f.sub.DAQ of the data-acquiring system limits the maximum frequency f.sub.m of the voltages V.sub.n.sup.exec, the Nyquist frequency f.sub.Nyq=f.sub.DAQ/2 being an upper limit.

[0106] For a data-acquiring system such that f.sub.DAQ=1 MS/s, the Nyquist frequency is equal to 500 KHz.

[0107] To avoid small residual-voltage errors resulting from energy stored in the electrode-electrolyte contact impedance, continuous signals must be delivered. To generate continuous signals at various frequencies, the signals are chosen to be harmonics of the lowest generated frequency f.sub.1.

[0108] When application of potential to an electrode is stopped, some of the electrical energy is stored for a few tens of microseconds at the interface between the electrode and the medium. This contact-impedance effect leads to an error in the measurement and implies the need to introduce a dead time between two successive excitations, to allow this energy to dissipate. Generation of continuous signals has the advantage of meaning that the applied voltages never stop, and therefore of avoiding errors related to contact impedance and the need to introduce dead times.

[0109] With 16 electrodes, a set of excitation signals at 15 different frequencies is generated. Given the sampling rate of the acquiring system, for example 1 MS/s, the frequencies may be chosen such that f.sub.i=i*f.sub.0 where f.sub.0 is the fundamental frequency and i is comprised between 1 and 15.

[0110] Furthermore, the discrete Fourier transforms may be performed at P=32 points because only positive results are considered. This results in a rate of acquisition of image data of 1*10.sup.6/32=31 250 images per second. This choice implies that the lowest sinusoidal-signal frequency f.sub.1 must be equal to the frequency of computation of the discrete Fourier transforms. The highest frequency f.sub.15=15*f.sub.0 is 468.875 kHz, which is below the Nyquist limit of 500 kHz for the system under consideration.

[0111] The excitation amplitudes are determined in the following way.

[0112] The voltage generation and acquiring modules operate in an interval of ?10 V. Considering the voltages V.sub.n.sup.exc, the amplitude A of the sinusoids must be significantly lower than the sum of the sinusoids generated due to constructive interference. However, the amplitude A of the signal must be as large as possible to minimize signal-to-noise ratio.

[0113] Another limit to take into consideration is the permitted maximum variation between two successively generated potentials. Real-time control in the acquiring system makes it possible to choose a satisfactory value of A=0.15 V, which gives resonant peaks of ?2.25 V. The rapid transition between positive and negative values of the signal prevents the appearance of electrolytic effects. For example, an electrolysis-related effect appears in water when a DC voltage higher than 1.2 V is applied. This effect does not appear at AC voltages higher than 1.2 V if these voltages vary rapidly enough.

Obtaining a Data Matrix D

[0114] In a second sub-step, a data matrix D is obtained from the data points M.sub.n.

[0115] For the 16-electrode system in question, which employs NI-9223 modules, a fixed-point data format of 20 number-allocated bits, including 5 bits for precision digits, may be used. The index n of the electrodes, which is comprised between 1 and 16, and the Fourier coefficient k, which is comprised between 1 and 15, may be described by binary numbers coded on 4 bits.

[0116] For a 32-electrode system, each data-point modulus M.sub.n(k) is coded with the numbers n and k in the following form:

[00012]$\begin{array}{cc}{M}_n\left(k\right)=\underset{k}{\underset{?}{.Math.+,8,0.Math.}}+\underset{n}{\underset{?}{.Math.+,8,0.Math.}}+\underset{M}{\underset{?}{.Math.?,16,11.Math.}},& \left[\mathrm{Math}8\right]\end{array}$

where the fixed-point format <s, b, p> is used with s: signed/unsigned; b: number of bits allocated; p: number of precision bits. The data M.sub.n(k) are encoded in U32 format.

[0117] The n.sub.e data points measured on the ne electrodes for a given Fourier coefficient k, i.e. for a given frequency f.sub.m, give the following data vector:

{M.sub.n(k)}=(M.sub.1(k)M.sub.2(k)M.sub.3(k) . . . M.sub.n.sub.e(k)) [Math 9]

where the M.sub.n(k) are integers encoded in U32 format.

[0118] The n.sub.e?1 vectors may be concatenated into a data matrix D:

[00013]$\begin{array}{cc}D=\left(\begin{array}{c}\left\{M_n\left(1\right)\right\}\\ \left\{M_n\left(2\right)\right\}\\ \left\{M_2\left(3\right)\right\}\\ .Math.\\ \left\{M_n\left(n_e-1\right)\right\}\end{array}\right).& \left[\mathrm{Math}10\right]\end{array}$

[0119] The size of the data is then S=n.sub.e (n.sub.e?1).

[0120] Only the moduli of the Fourier transforms form part of the data.

[0121] For an image, the size of the data is S*32 bits=4 kB. In comparison, the method described in publications [2], [3] and [4] results in 127 kB of data per image, without additional information on the boundary conditions.

[0122] Thus, the method according to the invention allows a higher image acquisition rate and also allows the size of the data to be reduced by a factor n.sub.e/2 with respect to the method described in publications [2], [3] and [4].

Obtaining a Signed Data Matrix

[0123] In a third sub-step, a signed data matrix representative of an image is obtained.

[0124] The Fourier transforms yield a modulus and a phase. The sign of each data point is estimated from the phase.

[0125] Thus, considering an excitation signal applied to each electrode E.sub.n at a frequency f.sub.m having the form:

[00014]$\begin{array}{cc}{V}_{n,m}^{\mathrm{exc}}\left(t\right)=A\cos \left(2?f_{m}t\right)\left[{?}_m^{?}\cos \left(m{?}_n\right)+{?}_m^{?}\sin \left(\frac{m{?}_n}{2}\right)\right],& \left[\mathrm{Math}11\right]\end{array}$

[0126] The phase of the signal is expressed:

[00015]$\begin{array}{cc}{?}_n^V\left(k\right)=\mathrm{arc}\tan \left(\frac{-{.Math.}_{p=0}^{P-1}{V}_{n,m}^{\mathrm{meas}}\left(p\right)\sin \left(k{?}_p\right)}{{.Math.}_{p=0}^{P-1}{V}_{n,m}^{\mathrm{meas}}\left(p\right)\cos \left(k{?}_p\right)}\right).& \left[\mathrm{Math}12\right]\end{array}$

[0127] The phase of the current I.sub.l.sup.meas(t) measured at electrode E.sub.l at the frequency f.sub.m is:

[00016]$\begin{array}{cc}{?}_l^I\left(k\right)=\mathrm{arc}\tan \left(\frac{-{.Math.}_{p=0}^{P-1}{I}_l^{\mathrm{meas}}\left(p\right)\sin \left(k{?}_p\right)}{{.Math.}_{p=0}^{P-1}{I}_l^{\mathrm{meas}}\left(p\right)\cos \left(k{?}_p\right)}\right).& \left[\mathrm{Math}13\right]\end{array}$

[0128] Assuming synchronous sampling of the input of the analog signal AI and the output of the sampled signal AO, the phase shift between the excitation potential and the measured current is:

?.sub.n,l(k)=?.sub.n.sup.I(k)??.sub.l.sup.V(k). [Math 14]

[0129] The phase shift depends on the design of the EIT sensors and on the nature of the flow within the body. If the phase shift is large, wrap-around effects may make reconstruction of the sign of the data impossible. Specifically, phase is symmetric to a transformation of 2?N, where N is an integer. If the phase is greater than 2?, it wraps around on itself.

[0130] The following two cases have been identified:

|?.sub.n,l(k)|<?/2 and |?.sub.n,l(k)|??2. [Math 15]

[0131] In the first case, the sign of the data points is computed using the preceding equation. The element of the signed data matrix corresponding to the k-th Fourier coefficient and to the n-th electrode is then:

[00017]$\begin{array}{cc}{\tilde{D}}_n^m=\frac{\sin \left({?}_{n,l}\left(k\right)\right)}{.Math.\sin \left({?}_{n,l}\left(k\right)\right).Math.}{D}_n^m.& \left[\mathrm{Math}16\right]\end{array}$

[0132] In the second case, the wrap-around effects prevent estimation of the sign of D.sub.n.sup.k. A sign matrix ? is then introduced to allocate an arbitrary sign to the data.

[0133] The sign matrix ? is estimated based on the sign of the excitation signal at t=0 for a given harmonic at a given electrode, by introducing the sign function:

[00018]$\begin{array}{cc}{.Math.}_{n,m}=\{\begin{array}{cc}1& \mathrm{if}{V}_n^{\mathrm{meas}}>0,\\ 0& \mathrm{if}{V}_n^{\mathrm{meas}}=0,\\ -1& \mathrm{if}{V}_n^{\mathrm{meas}}<0\end{array},& \left[\mathrm{Math}17\right]\end{array}$

with V.sub.n.sup.meas such as defined above.

[0134] The sign matrix ? is more particularly defined so that its rows represent the signs of a cosine function in alternation with those of a sine function having an integer number of periods on each row. The first two rows have a single period and the number of periods increases by one for each pair of subsequent rows.

[0135] In other words, the rows of the sign matrix ? are defined such that the i-th element of the j-th row is the sign of cosine([2?/([j+1]/2)]*(i?1)/n.sub.e) for odd j and sine([2?/(j/2)]*(i?1)/n.sub.e) for even j.

[0136] For example, the first element of the first row of the sign matrix is the sign of cosine(0), i.e. +. For a number of electrodes n.sub.e equal to 16, the sign of the fifth element of the first row is the sign of cosine(2?*[4/16]), i.e. 0.

[0137] Thus, the first row of the sign matrix represents the signs of one period of a cosine function, i.e. the signs of cosine(2?*(i?1)/n.sub.e) for the i-th element in the row. The second row represents the signs of one period of a sine function, i.e. the signs of sine(2?*(i?1)/n.sub.e). The third row represents the signs of two periods of a cosine function, i.e. the signs of cosine([2?/2]*(i?1)/n.sub.e). The fourth row represents the signs of two periods of a sine function, i.e. the signs of sine([2?/2]*(i?1)/n.sub.e). The fifth row represents the signs of three periods of a cosine function, i.e. the signs of cosine([2?/3]*(i?1)/n.sub.e), and so on until the (n.sub.e?1)-th row.

[0138] For example, for a 16-electrode device, the sign matrix takes the form shown in FIG. 7. For a device with 32 electrodes, the sign matrix takes the form illustrated in FIG. 8.

[0139] Use of such a sign matrix makes it possible to optimize the processing of the data to form an image. It allows the sign of each data point to be estimated and an image to be reconstructed.

[0140] For a large phase shift, i.e. larger than or equal to ?/2, the arbitrarily signed amplitudes of the n.sub.e?1 excitation patterns of the n.sub.e electrodes, or in other words the elements of the signed data matrix corresponding to the k-th Fourier coefficient and to the n-th electrode, are given by:

{tilde over (D)}.sub.n.sup.m=?.sub.n.sup.mD.sub.n.sup.m. [Math 18]

Algorithmic Implementation

[0141] The HOST portion of the acquiring system 4 continuously sends the frequency and amplitude parameters to the FPGA of the acquiring system 4.

[0142] FIG. 4 illustrates an algorithm for generating the excitation signals of the 16 electrodes 2.

[0143] With the system considered, the FPGA receives data points 16 by 16 in a loop clocked at 1 MS/s to create 16 analog signals.

[0144] In a first step, the system is initialized.

[0145] To start with the FPGA is empty. The HOST loads the FPGA, then the NI-9262 modules are reset.

[0146] In a second step, an interrupt request is sent and received.

[0147] A hardware interrupt is used to notify the HOST when the FPGA is ready to begin data acquisition. The FPGA waits for HOST validation to begin acquisition.

[0148] In a third step, sampling is checked.

[0149] A sample-pulse generation function is called to begin generating data points. The frequency at which the function is called determines the sampling rate used to generate the data points. In parallel, an I/O status write function is called at the same frequency to check the status of each generated sample.

[0150] In a fourth step, digital excitation signal functions are generated.

[0151] The HOST commands the FPGA to begin generating the excitation signal functions.

[0152] In a fifth step, analog excitation signals are generated. The sixteen excitation signals are sent to the electrodes.

[0153] In a sixth step, the HOST confirms the generation of the signals and reports any errors at the HOST or FPGA level.

[0154] Steps 1 and 2 are performed once, when the algorithm is launched. Steps 3 to 6 are iterated for each output point, at the sampling frequency.

[0155] The sampling frequency may be comprised between 10 kHz and 500 MHZ, and preferably between 500 kHz and 50 MHz.

[0156] FIG. 5 shows an algorithm for processing data received from the electrodes 2.

[0157] Measurement of the voltages across the terminals of the resistors R allows the Neumann boundary conditions to use in the implementation of the image reconstruction algorithm to be deduced. The acquisition rate of 1 MS/s on 16 channels corresponds to a data transfer rate of 320 MB/s.

[0158] Using fast Fourier transforms while only considering Fourier coefficients related to a generated signal allows the size of the data to be reduced without impacting their quality. It also acts as an effective band filter. However, the real-time computation of 16 fast Fourier transforms requires a high computing capacity. The FPGA, which permits real-time and parallel transformation of signals into their Fourier components on a plurality of channels, is a tool suitable for carrying out this task.

[0159] In a first step A, the system is initialized.

[0160] The FPGA resets the NI-9223 analog-signal-acquiring modules.

[0161] In a second step B, the memories are configured.

[0162] The HOST configures and initiates direct access to the memory of the FPGA. The FPGA configures and initiates sixteen first-in-first-out queues to ensure communication between the sixteen measurement channels and their fast-Fourier-transform computations.

[0163] In a third step C, an interrupt request is sent and received.

[0164] A hardware interrupt makes it possible to ensure that the queues and direct memory access are ready.

[0165] In a fourth step D, sampling is checked.

[0166] A sample pulse generation function is called to control the sampling frequency and an I/O status read function is called at the same frequency to check the status of each sample and report any errors to the HOST.

[0167] In a fifth step E, the analog measurement signals are acquired.

[0168] The I/O read function is configured to read out a single sample from each channel of each NI-9223 module. This function is called at 1 MHz and regulated by the sample pulse generation function.

[0169] In a sixth step F, the fast Fourier transform of each channel is computed in a 1 MHz loop. The computation time is determined by the number P of points considered in the Fourier transform. Once the P measured data points have been transferred for the Fourier-transform computations, the function returns the P Fourier coefficients one by one in each iteration of the Fourier-transform loop. Subsequently, the amplitudes of the Fourier coefficients are computed at a frequency of 1 MHz.

[0170] In a seventh step G, the data are addressed.

[0171] The amplitude data points in U32 format form, with the corresponding Fourier coefficient (in U16 format) and the corresponding channel (in U16 format), a data element in U64 format. In each iteration of the Fourier transform loop, 16 elements are written for the 16 channels to the direct-access memory, with a view to being transmitted to the HOST.

[0172] In an eighth step H, the data matrix is constructed.

[0173] The HOST waits for the direct-access memory to collect at least 240 elements, representing a complete data image. The address of the n-th electrode and the Fourier coefficient associated with the amplitude are used to form the data matrix D.

[0174] In a ninth step I, the data matrices are stored.

[0175] The data are used to perform real-time image reconstruction or are stored. The image is based on a one-step iterative least-squares reconstruction algorithm, for example the algorithm described in publication [5].

[0176] The real-time image reconstruction may generate about one hundred images per second.

[0177] In a tenth step J, acquisition of the signals is confirmed and an error check is carried out.

[0178] Synchronization of the analog-signal-generating and -measuring modules is checked and any errors are reported.

[0179] Steps A to C are performed once, when the algorithm is launched. Steps D to G are iterated for each output point, at the sampling frequency. Once a complete data matrix containing n.sub.e*(n.sub.e?1) data points has been acquired, the algorithm passes to step H. Steps H to J are iterated at the image-acquisition frequency.

[0180] The sampling frequency may be comprised between 10 KHz and 500 MHZ, and preferably between 500 KHz and 50 MHz.

[0181] FIG. 6 shows, in the graphs of the left-hand side, the excitation signals of six of the sixteen electrodes 2. Each signal is made up of a sum of 15 sinusoidal functions.

[0182] The continuous lines of the graphs on the right-hand side of FIG. 6 represent, in the equivalent Fourier space, the amplitudes of the voltages V.sup.meas measured across the resistors R shown in FIG. 2, with R=200 ?. The dash-dotted lines represent the Fourier transforms of the voltages of the generated signals.

[0183] The invention is not limited to the examples that have just been described: features of the illustrated examples may in particular be combined together within variants that have not been illustrated.

[0184] Further variants and improvements may be envisioned without however departing from the scope of the invention. In particular, the method according to the invention may be implemented by means of acquiring systems other than the one described.

LIST OF CITED REFERENCES

[0185] [1] Teague, G. (2002). Mass flow measurement of multi-phase mixtures by means of tomographic techniques. University of Cape Town, Faculty of Engineering, Department of Electrical Engineering. [0186] [2] Dupr?, A., Mylvaganam, S. (2017). Simultaneous and Continuous Excitation Strategy for High Speed EIT: the ONE-SHOT method. In Proceedings of the 9th World Congress on Industrial Process Tomography, pages 667-674. [0187] [3] Darnajou, M., Dupr?, A., Dang, C., Ricciardi, G., Bourennane, S., Bellis, C. (2019). On the implementation of simultaneous multi-frequency excitations and measurements for electrical impedance tomography. Sensors, 19(17). [0188] [4] Darnajou, M., Dupr?, A., Dang, C., Ricciardi, G., Bourennane, S., Bellis, C., Mylvaganam, S. (2020). High Speed EIT with Multifrequency Excitation using FPGA and Response Analysis using FDM. IEEE Sensors. [0189] [5] https://www.math.colostate.edu/?siamcsu/files/NOSER.pdf