Historically, using electrical impedance of the human body for physiological study or diagnosis flourished after the introduction of the Tetrapolar Impedance Measurement (TPIM) technique by Nyboer [1], eliminating the skin-electrode contact impedance, which is very large compared to that of the bulk at typical measurement frequencies that had posed a challenge before. In this method, since the current drive and potential measuring electrodes are different, the value obtained is strictly not impedance, rather it is a “transfer Impedance” [2]. In a volume conductor, the transfer impedance very much depends on the geometry of the conductor, placement of the electrodes and the conductivity distribution inside. Therefore, different transfer impedance values will be obtained for different electrode placements, having more sensitivity to certain regions of the volume. TPIM has poor geometrical localization in terms of sensitivity and was used mainly for whole body measurements [3, 4].
A distribution image of impedance within a volume conductor named, “Electrical Impedance Tomography (EIT)” became popular after the pioneering works by several groups [5, 6]. The technique uses multiple TPIM measurements using 16, 32 or more electrodes fixed around a volume conductor and was initially proposed for obtaining 2D tomographic images. This is a complex technique, requiring sophisticated hardware and computer algorithms. Since electrical current spreads out in all directions and does not follow a straight path, 2D images in fact contained significant contributions from 3D [7].
Eventually EIT was extended for 3D imaging using multiple level electrode arrays and is still going through an R&D phase [8]. Obviously, the technology, when ready, would be more complex, contributing to a high cost and sophistication, and may find little use in low and medium income countries (LMIC), which is an area of interest of the present authors.
The authors’ extended group at Dhaka University introduced the Focused Impedance Method (FIM) as a bridge between TPIM and EIT [9, 10, 11]. Using 8, 6 or 4 electrodes it can localize a target region [12]. The complexity is much less than that of EIT and is marginally greater than that of TPIM. FIM has shown success in localized lung ventilation and gastric emptying, giving similar results as obtained using 2D EIT [9, 13]. However, most of the applications attempted so far using FIM targeted measurement from one side of the human body in which contribution of the deep organs were limited, only partial contribution were coming from the regions at shallow depths with respect to the electrodes. In order to probe deeper regions of these organs, some new ideas on the placements of electrodes were suggested [14] using visualizations based on the concepts of equipotential surfaces [15] and of sensitivity equations involving current density vectors due to two separate pairs of electrodes in TPIM [16]. Finite element simulation and phantom measurements using one of these suggested methods targeting the lungs produced promising results [17]. The present work was taken up using another of the novel ideas on electrode placement suggested in the above work [14] that has the possibility of targeting the stomach. This idea was originally suggested for an asymmetrically placed deep object as shown in a cross sectional diagram in Figure 1. Here, two pairs of electrodes indicated by I1-I1’ and V1-V1’ form the first TPIM configuration while I2-I2’ and V2-V2’ indicate the second. The algebraic average of the two TPIM values make the FIM. The dashed lines in the left hand figure show the maximum current lines between respective pairs of electrodes, which also includes the potential measuring electrodes, while the dotted lines on the right hand figure show approximate equipotential lines through the potential measuring electrodes. The left hand diagram suggests that the region between the maximum current lines will have high sensitivity to the overall transfer impedances through both the TPIM configurations, being higher near to the bottom of the diagram, since the electrode separation between the respective I and V electrodes are less on this side of the volume conductor. The right hand diagram indicates that the target object, falling within the four equipotential lines, will have an increased contribution to the overall FIM value, measured using this configuration of electrodes. Of course, this configuration is just a qualitative suggestion, simulation using a finite element method (FEM) will give a more optimized configuration.
Based on the above suggestions, this work carries out a detailed FEM simulation using the COMSOL software package [18] on a simplified model of the human trunk and
the stomach and several electrode configurations, choosing the electrode configurations rather intuitively. Through this, the present work attempts to determine the contribution of the modeled stomach in the measurement of the total transfer impedance to find out if the suggested idea provides outcomes as anticipated. Then changing the electrode configurations a little from the starting configuration, this work will try to find the optimum electrode positions to get a maximum contribution. It will also explore if a simple TPIM with a suitable electrode placement may provide contribution figures close to that from FIM. Furthermore, based on these results, it will try to explore if some more novel electrode configurations may be more effective for the stomach.
As mentioned above, the present work intends to get a reasonably high contribution from the stomach within a human trunk. Following the method suggested above, this measurement will be performed on a transverse plane of the trunk at the level of the stomach. Placing electrodes suitably to opposite sides of this cross-sectional plane it is expected that a reasonably high contribution from the stomach will be obtained in the total transfer impedance measurement.
As a first attempt at simulation for the stomach, the present work chose a very simplified model of the human trunk with an embedded simplified stomach model as shown in Figure 2. The trunk was simulated by a vertical rectangular parallelepiped with rounded edges on all four sides. The horizontal dimensions chosen were 30 cm (representing trunk breadth, x-axis) and 25 cm (representing trunk thickness or depth, y-axis) while the vertical dimension was chosen at 40 cm (representing trunk height, z-axis). The radii of curvature for rounding of the rectangle at the corners were made greater on the assumed back of the body compared to that for the front; the values arbitrarily chosen at 8 cm and 5 cm respectively. This was done since the human trunk has similar features. Figure 3 shows a cross section of the trunk with the stomach model embedded.
In order to model the simplified stomach, a CT scan [19] was used as a guide. For this simulation, the stomach was modeled as a truncated circular cylinder with diameter 16 cm and height 8 cm (Figures 2 and 3). The cylinder was truncated vertically on the right side, 4 cm from the edge of the circle, making the lateral dimension of the modeled stomach as 12 cm. The lower edge of the modeled stomach was placed at a height of 15 cm from the bottom of the trunk (Figure 2). It was placed towards the left side of the trunk, its left edge placed at a distance of 3 cm from the left trunk surface and the front edge placed at a distance of 3 cm from the front trunk surface as shown in Figure 3.
For the simulation, the conductivity σ and the relative permittivity ε of the volume conductor were chosen as appropriate for isotonic 0.9% saline [20], while those for the stomach were arbitrarily chosen, of a higher conductivity. Here, these were chosen as that for seawater, (3.5% saline) [21]. This work was just to determine the contribution of the stomach region in the overall measurements, so was not related to real human tissues or real parameters for any food.
Electrodes were initially modeled as circular cylinders with a cross sectional diameter of 2 cm and length 3 cm each, but embedded along the boundary of the volume conductor such that the portion outside the surface of the trunk has a high conductivity, that for a metal, while the portion inside the surface, within the volume conductor, equals that for the volume conductor. Effectively, the electrodes just touched the outer surface of the volume conductor. Since on moving the conductor at different sites on the periphery changed the part of the electrode inside and outside the trunk volume, this feature was selected in COMSOL through appropriate partitioning of the electrode volume.
The placement of the electrodes chosen for the initial simulation are shown in Figure 4. Here,
The above logic is the result of a visualization of the equation for conventional point sensitivity
where
One can visualize that
COMSOL simulation was made to obtain TPIM1 and TPIM2 and averaging them to obtain FIM for the two configurations as mentioned above. In each case the total transfer impedance value ZT as well as the contribution Zb of the background volume outside the stomach region and the contribution Zst of the stomach volume were determined for each of the measurements (TPIM1, TPIM2, FIM) using the following equations:
where ρb is the resistivity of the background of the trunk volume outside the stomach, assumed uniform throughout, ρst is the resistivity of the stomach region, assumed uniform, and
Obviously, the target would be to make Cst as large as possible. The simulations were performed to obtain Cst for the following two situations.
For both, the conductivity and relative permittivity of the part of the electrodes outside the volume conductor was taken as 6×107 S/m and 1, respectively. As mentioned before, the part of the electrode falling inside the volume conductor was made to have the same electrical parameters as for the background volume of the trunk.
A frequency of 200 kHz was assumed for the simulation. The electrical parameters used are summarized in Table 1. The geometrical parameters appropriate for the chosen model are shown in Table 2. The geometries of the electrode placement configurations are given together with the calculated values of Cst in the results section.
Assumed electrical values for simulation. Frequency = 200 kHz. Current = 1 A.
Electrical Properties | Background trunk volume* | Stomach volume** | Electrode (copper) |
---|---|---|---|
Conductivity (σ) [S/m] | 1.2 | 5 | 6 × 107 |
Relative Permittivity (ϵ) | 80 | 80 | 1 |
* Equivalent to that of 0.9% saline at 200C, also for stomach in Simulation A.
** Equivalent to that of 3.5% Saline (sea water) at 200C, for simulation B.
Geometrical parameters appropriate for the chosen model
Volume of whole Trunk (including stomach) | 28430 cm3 |
Volume of Stomach | 1293 cm3 |
Stomach volume with respect to whole volume | 4.54 % |
Horizontal Cross Sectional Area of Trunk | 711.4 cm2 |
Horizontal Cross Sectional Area of Stomach | 161.4cm2 |
Area of Stomach with respect to Area of Trunk | |
(cross-sectional) | 22.69 % |
While performing the analyses, it was found out that the initial trial configuration of TPIM2 has a lower Cst value compared to that for TPIM1. It was suggested that the electrode placement need some further changes. At this point, further guidance was taken from an equipotential concept [15], that region bounded by the equipotentials passing through the potential electrodes have a high contribution to the transfer impedance. Based on this concept the position of the electrodes for both TPIM1 and TPIM2 were examined qualitatively with the help of Figure 5.
For TPIM1, the equipotential lines bound the stomach region well on the top and bottom. On the other hand, the equipotential lines for Trial-TPIM2 include a large amount of background volume on the left side of the stomach.
Therefore, based on this pictorial information, the lower potential electrode V2’ was shifted to V2N’ as shown for TPIM2 so that the stomach is bounded well between the two equipotential lines. Consequently, it was envisaged that TPIM2 will have a higher contribution from the stomach so that the resulting FIM will also have a higher contribution from the stomach. This was actually observed, and so all the analyses shown in this work use TPIM2, not the Trial-TPIM2.
In order to have an insight into the effect of the separation between the adjacent current and voltage electrodes, the electrode separations near to the stomach (I1-V1 and I2-V2) were increased or decreased by about 20% from the intuitively chosen ones and the results were compared to see which separation gives the maximum contribution from the stomach.
The electrode configuration giving the maximum contribution of the stomach was used to obtain color maps of the point sensitivity in the horizontal plane at the vertical midpoint of the stomach, which is also the plane of the centers of the electrodes. These color maps provide further insight into the contribution of the different regions to the total transfer impedance, which is very helpful (not shown).
Small electrodes, as used in the above work have wide diverging current patterns within the volume conductor. Since the stomach has a vertical height, in order to get more contribution from this 3D stomach, long strip electrodes, as suggested by Rabbani [14] was also used for the electrode configuration giving the maximum contribution of the stomach using the small cylindrical electrodes (2 cm diameter). The comparison would give valuable insight. The strip electrodes were rectangular, 2 cm × 2 cm in cross section, and of variable lengths: 5 cm, 6 cm and 8 cm. These were positioned centrally on the stomach in the vertical direction. In this case, the current paths are expected to diverge less within the vertical extension of the stomach region giving higher current densities in the stomach contributing to a higher contribution in the total transfer impedance.
Another electrode configuration was suggested and tested during the present work. In this configuration, the TPIM1 of the previous configuration was used in the horizontal plane while a similar configuration of electrodes, TPIM1-O was placed orthogonal to TPIM1 as shown in Figure 6. Since the
stomach is nearer to the front surface, to show the nearside electrodes, I1-V1 and the orthogonal ones, this diagram reverses the geometrical representation shown in Figure 2 earlier. Here,
All the results and the comparative analyses are presented in the results section below.
The conducted research is not related to either human or animal use.
For the geometry, the origin (0,0,0) has been taken at the lower left corner of the trunk model as shown in Figure 2, which is on the right side of the back of the trunk, in physical perspective. The co-ordinate (
For the long vertical strip electrodes used in one of the studies, the center of the strips were located at z = 19 cm, while the x and y coordinates are as referenced in Figure 7. Separations between pairs of electrodes were calculated as the direct straight path length between the mid-points of these electrodes touching the trunk volume, not based on their geometrical reference points.
The results of the various simulations performed using COMSOL are given below, together with brief observations wherever feasible. The names of the electrodes are also as in Figure 7. The coordinates of the electrodes are also given in the tables, where necessary. The values in scientific notation follows that of EXCEL where 1.0E6 would mean 1×106.
Table 3 shows the simulated values obtained for TPIM1 with changing the separation of electrodes I1-V1, nearest to the stomach on the front side. The separation of the other electrodes I1’-V1’ remain the same at 15.4 cm. This is the result for Simulation-A as mentioned before, with the stomach and the background having the same electrical parameters as that for isotonic (0.9%) saline. The maximum and minimum sensitivity shown in the table refers to values on the electrode plane (x-y plane) only, while the impedance values refer to the whole volume of the trunk or the stomach as appropriate. The percentage contribution of the stomach for the three separations are shown in the table as well as in Figure 8.
Simulation-A for TPIM1 for varying separation of I1-V1 Electrode Co-ordinates for 3 separations: [I1:(16,25,19), V1:(20,25,19)]; [I1:(14,25,19), V1:(22,25,19)]; [I1:(12,25,19), V1:(24,25,19)]. Fixed [I1’=(16,-3,19); V1’=(29,6.5,19)] Direct straight-line separation for I1'-V1' = 15.4 cm
Separation of I1-V1, cm | |||
---|---|---|---|
Simulated parameters | 4 | 8 | 12 |
Maximum Sensitivity* | 7.71E+06 | 1.27E+06 | 8.99E+05 |
Minimum Sensitivity* | -2.22E+06 | -3.66E+05 | -2.36E+05 |
Total Impedance (ohm) | 3.23 | 1.57 | 1.07 |
Stomach Impedance (ohm) | 0.61 | 0.35 | 0.23 |
* values on the electrode plane (x-y plane) only
Table 4 shows the simulated values (Simulation-A) obtained for TPIM2 with changing the separation of electrodes I2-V2, nearest to the stomach on the left side. The separation of the other electrodes I2’-V2’ remain the same at 24.7 cm. The percentage contribution of the stomach for the three separations are shown in the table as well as in Figure 9.
Simulation-A for TPIM2 for varying separation of I2-V2 Electrode Co-ordinates for 3 separations: [I2:( 30,14,19), V2: (30,18,19)]; [I2:( 30,12,19), V2:( 30,20,19)]; [I2:( 30,10,19), V2: ( 30,22,19)]. Fixed [I1’=(13,-3,19), V1’=(2,21,19)] Direct straight-line separation for I2'-V2' = 24.7 cm
I2-V2, cm | |||
---|---|---|---|
Simulated parameters | 4 | 8 | 12 |
Max Sensitivity* | 9.65E+05 | 5.10E+05 | 4.01E+05 |
Min Sensitivity* | -1.35E+07 | -2.99E+06 | -1.14E+06 |
Total Impedance (ohm) | 3.08 | 1.63 | 1.06 |
Stomach Impedance (ohm) | 0.60 | 0.38 | 0.21 |
* values on the electrode plane (x-y plane) only
It is interesting to note that for both TPIM1 and TPIM2, the separations of the I and V electrodes nearest to the stomach giving the maximum Cst values were the ones that were initially chosen based on qualitative visualisations based on point sensitivity considerations and equipotentials. Increase and decrease of the separation led to reduced Cst values. Table 5 gives a comparison of the best Cst values obtained using TPIM1 and TPIM2 and FIM using these two electrode configurations. Figure 10 shows this result through a histogram.
Simulation-A for FIM with best TPIM configurations Direct straight-line separations: I1-V1 = 8 cm, I1'-V1' = 15.4 cm, I2-V2 = 8 cm, I2'-V2' = 24.7 cm
Simulated parameters | TPIM1 | TPIM2 | FIM |
---|---|---|---|
Max Sensitivity* | 1.27E+06 | 5.10E+05 | 6.40E+05 |
Min Sensitivity* | -3.66E+05 | -2.99E+06 | -1.48E+06 |
Total Impedance (ohm) | 1.57 | 1.63 | 1.60 |
Stomach Impedance (ohm) | 0.35 | 0.38 | 0.36 |
*values on the electrode plane (x-y plane) only
Table 6 gives the results of Simulation-B in which the stomach is assumed to have a higher conductivity, equivalent to that of 3.5% saline (sea water) using the same electrode configurations as for Table 5. The rest of the trunk volume is assumed to have the electrical parameters appropriate for isotonic (0.9%) saline. Figure 11 shows this result through a histogram. It is interesting to note that although the Cst values are almost similar for all the three measurements for Simulation-A, TPIM1 attains a higher value than TPIM2 for Simulation B, with FIM being in-between.
Simulation-B for FIM with best TPIM configurations Direct straight-line separations: I1-V1 = 8 cm, I1'-V1' = 15.4 cm, I2-V2 = 8 cm, I2'-V2' = 24.7 cm
Simulated parameters | TPIM1 | TPIM2 | FIM |
---|---|---|---|
Max Sensitivity* | 8.41E+05 | 4.86E+05 | 4.23E+05 |
Min Sensitivity* | -2.41E+05 | -2.09E+06 | -1.05E+06 |
Total Impedance (ohm) | 1.01 | 1.096 | 1.055 |
Stomach Impedance (ohm) | 0.389 | 0.333 | 0.3611 |
*values on the electrode plane (x-y plane) only
Table 7 shows Simulation-A values obtained for all the three configurations for different lengths of rectangular strip electrodes, for the chosen best placements of electrodes and Figure 12 shows the results graphically.
Simulation-A for Strip Electrodes of different lengths Strip electrode cross section: 2 cm × 2 cm Direct straight-line separations: I1-V1 = 8 cm, I1'-V1' = 15.4 cm, I2-V2 = 8 cm, I2'-V2' = 24.7 cm
Length = 5 cm | Length = 6 cm | Length = 8 cm | |||||||
---|---|---|---|---|---|---|---|---|---|
Simulated Parameters | TPIM1 | TPIM2 | FIM | TPIM1 | TPIM2 | FIM | TPIM1 | TPIM2 | FIM |
Max Sensitivity* | 3.90E+04 | 4.47E+04 | 2.43E+04 | 2.64E+04 | 2.37E+04 | 1.67E+04 | 5.42E+04 | 3.05E+04 | 3.27E+04 |
Min Sensitivity* | -6.32E+04 | -7.53E+04 | -3.60E+04 | -5.14E+04 | -6.73E+04 | -3.26E+04 | -5.77E+04 | -1.57E+05 | -7.80E+04 |
Total Impedance, Ω | 1.49 | 1.57 | 1.53 | 1.40 | 1.48 | 1.44 | 1.47 | 1.55 | 1.51 |
Stomach Impedance, Ω | 0.33 | 0.36 | 0.34 | 0.29 | 0.32 | 0.30 | 0.32 | 0.35 | 0.33 |
*values on the x-y plane at the vertical center of the electrode and stomach
Figure 13 shows percentage contribution of stomach in Simulation-A for 2 cm round electrodes and for 5 cm long strip electrodes of 2 cm × 2 cm cross section, for all the three configurations, TPIM1, TPIM2 and FIM.
Table 8 shows the impedance values for the newly proposed orthogonal configuration, targeting a solid angle into the trunk volume around the stomach. All the three give reasonably good contribution for the stomach.
Simulation-A for Orthogonal Electrode I1-V1 = 8 cm, I1'-V1' = 15.4 cm, I2-V2 = 8 cm, I2'-V2' = 15.4 cm
Simulated parameter | TPIM1 | TPIM1-O | FIM |
---|---|---|---|
Total Impedance (ohm) | 1.58 | 1.40 | 1.49 |
Stomach Impedance (ohm) | 0.35 | 0.29 | 0.32 |
For any electrical impedance measurement probing deep inside a body, the popular way would be using EIT. However, the main motivation of the authors’ group is to develop technologies for the low income countries, which should be simple and robust. Therefore, if simple TPIM or FIM can be used to get information on deep organs with reasonable sensitivity, it may find a useful application, not only for the low income countries, but for the high income countries as well.
The success of using a TPIM or FIM depends on judicial placement of electrodes, where understanding of the basic visualization based on point sensitivity equation and equipotential concepts, both are necessary. This has been illustrated through this work.
The present work has tried small round electrodes and long strip electrodes. Apparently they all gave percentage contribution of the stomach of a similar order, with no preference for either. This suggests that using a small round electrode is preferable in view of practical considerations.
Almost all the methods give a percentage contribution of the stomach a little above 20%. It is interesting to note that the percentage of cross sectional area of the modeled stomach to that of the whole trunk is of similar order (22.7%). Therefore, the values obtained are possibly the maximum that one can get through TPIM or FIM.
From the similar results for TPIM and FIM, one might argue that since TPIM is simpler, it should be preferable. However, it should be recognized that TPIM has more contribution from negative sensitivity zones, which cancel to some extent in FIM [11]. The maximum and minimum sensitivities presented in the tables indicate that these are slightly less in FIM in general. However, these values only pertain to the electrode plane. Usually, these values are expected to be less in 3D.
If the changes only occur in the stomach, then either TPIM or FIM may give useful results. However, if other organs like the duodenum and the intestine change, which is very likely, these changes will affect the result. Therefore, FIM may have an advantage here in reducing the effects of nearby organs. Further studies taking such organs in consideration should be taken up which will show which method is best suited for a particular application.
Therefore, through further studies following the approach presented in this paper, one may get useful physiological and diagnostic information in the deep organs like the stomach using simple techniques like TPIM or FIM.