Simultaneous Imaging of Bio- and Non-Conductive Targets by Combining Frequency and Time Difference Imaging Methods in Electrical Impedance Tomography

As a promising medical imaging modality, electrical impedance tomography (EIT) can image the electrical properties within a region of interest using electrical measurements applied at electrodes on the region boundary. This paper proposes to combine frequency and time difference imaging methods in EIT to simultaneously image bio- and non-conductive targets, where the image fusion is accomplished by applying a wavelet-based technique. To enable image fusion, both time and frequency difference imaging methods are investigated regarding the reconstruction of bio- or non-conductive inclusions in the target region at varied excitation frequencies, indicating that none of those two methods can tackle with the scenarios where both bio- and non-conductive inclusions exist. This dilemma can be resolved by fusing the time difference (td) and appropriate frequency difference (fd) EIT images since they are complementary to each other. Through simulation and in vitro experiment, it is demonstrated that the proposed fusion method can reasonably reconstruct both the bio- and non-conductive inclusions within the lung models established to simulate the ventilation process, which is expected to be beneficial for the diagnosis of lung-tissue related diseases by EIT.


Introduction
Electrical impedance tomography (EIT) is a non-invasive imaging technique in which the absolute value or the change of conductivity distributions are reconstructed [1][2][3][4][5][6]. The underlying rotational is based on generating a current density distribution inside the medium and measuring the resultant electromagnetic fields by means of sensors positioned circumferentially on the boundary [2]. Current excitation and voltage acquisition are repeated continuously and rapidly until all the independent electrode combinations are deployed. Attributed to its notable merits of high temporal resolution, non-radiation, low cost, and non-intrusive measurement, EIT has also gained a great deal of interest in monitoring rapid changes in multiphase flow and medical imaging [2][3][4][5][6].
Ever since medical interest in EIT increased in 1980s, EIT has been developed to be a promising medical imaging technique to attain useful diagnostic information where there are obvious conductive changes in the biological tissues and organs with changes in physiological conditions, such as breathing, blood flow, digestive, canceration, and nervous activity [7,8]. Up to now, the medical applications of EIT research have mainly focused on thoracic imaging. Among the organs contained in the thorax, large volumes of air are Amoved into and out of the lungs, leading to obvious conductivity changes during ventilation. Therefore, it is agreed that monitoring of breathing by imaging the lung is the conducted to validate the proposed method by adopting an in vitro lung model made with a white gourd. This paper is structured as follows: In Section 2, the image reconstruction in EIT is briefly reviewed, and the wavelet-based fusion method for td-and fd-EIT images is demonstrated. In Sections 3 and 4, both simulation and experimental results are presented to validate the proposed method with simulated or in vitro lung models. Finally, conclusions are drawn in Section 5.

Principles of EIT
In this investigation the boundary was ∂Ω. Generally, for computational needs, it is assumed that Ω represents the body under the finite element model (FEM), Ω is discretized into small pixels or voxels and assigned a piecewise constant conductivity distribution, i.e., σ = [σ 1 , . . . , σ N ] T ∈ R N is assumed within Ω. The electric potential inside the sensing area induced by the injected AC current can be written as [5,6]: ∇·(σ(x, y)∇u(x, y)) = 0, (x, y) ∈ Ω (1) where u is the electric potential inside Ω. EIT is a nonlinear problem as the excitation current is a function of unknown conductivity distribution [21,22]. Therefore, based on the nonlinear model, the relation between the conductivity inside Ω and the boundary voltages at the boundary, ∂Ω, can be expressed as: where V is the measured voltage, U is the forward operator, and e is the additive noise and measurement error. Regarding the difference EIT imaging, Equation (2) can be linearized given a small conductivity perturbation ∆σ, which can be expressed as (3): where J is the Jacobian matrix, and calculated by solving the EIT forward problem, i.e., Equation (1). Difference EIT calculates a vector of conductivity changes between the measured conductivity distribution and the reference conductivity distribution. One commonly employed method is the single frequency td-EIT imaging. More specifically, the AC current with selected frequency is injected into the sensing area. Hereafter, a measurement at time t 0 is acquired as the reference, and then the voltage change at time t 1 against the reference measurement is employed for image reconstruction, which is expressed as: In most cases, a set of reference measurements is acquired when the EIT sensor is filled with the background medium. To improve its precision, the measured data is typically averaged over many data frames.
Recently, incorporating the advances that have been made on the hardware system, a more attractive method in biomedical applications is the frequency-difference (fd-) EIT imaging. For the fd-EIT imaging, the current source is required to output a pair of complementary currents with two different frequency components, i.e., f 1 and f 2 , which are injected to the sensor simultaneously or separately. At each frequency component, the boundary voltages are measured. The difference between the two sets of measured voltages is adopted for image reconstruction, which is calculated by After attaining the voltage difference, to solve the inverse or image reconstruction problem in EIT, a variety of algorithms have been proposed in the literature [23]. In general, Biosensors 2021, 11, 176 4 of 16 using the regularization technique is the predominant way. The conductivity distribution can be obtained by solving the optimization problem: arg min 1 2 J∆σ − ∆V 2 2 +λR(∆σ) (6) where R is the regularization function incorporating a priori knowledge and λ is the regularization parameter. In this case, for lung EIT, traditional and proprietary algorithms have been an obstacle to interpretation of EIT images as the reconstructed images cannot reveal the character of the lung [2]. GREIT (Graz consensus reconstruction algorithm for EIT) is a type of regularized imaging approach for 2D linear EIT reconstruction of the lung, which is based on the experience with a large number of linear regularized EIT reconstructions. One of the advantages of GREIT is that it conforms to the performance requirements for lung application, while other reconstruction algorithms do not. One can refer to Adler [2] for the detailed description of GREIT.

Wavelet-Based Image Fusion
Image fusion has three levels: pixel level, feature level, and decision level [20]. Pixellevel fusion transcends other image fusion levels as it retains more detailed image information and can achieve higher accuracy. Among the pixel-level fusion methods, the wavelet-based method is recommended because of its non-redundancy and directionality [24][25][26][27]. Figure 1 shows the image fusion process based on the wavelet transform. After attaining the voltage difference, to solve the inverse or image reconstruction problem in EIT, a variety of algorithms have been proposed in the literature [23]. In general, using the regularization technique is the predominant way. The conductivity distribution can be obtained by solving the optimization problem: where is the regularization function incorporating a priori knowledge and is the regularization parameter. In this case, for lung EIT, traditional and proprietary algorithms have been an obstacle to interpretation of EIT images as the reconstructed images cannot reveal the character of the lung [2]. GREIT (Graz consensus reconstruction algorithm for EIT) is a type of regularized imaging approach for 2D linear EIT reconstruction of the lung, which is based on the experience with a large number of linear regularized EIT reconstructions. One of the advantages of GREIT is that it conforms to the performance requirements for lung application, while other reconstruction algorithms do not. One can refer to Adler [2] for the detailed description of GREIT.

Wavelet-Based Image Fusion
Image fusion has three levels: pixel level, feature level, and decision level [20]. Pixellevel fusion transcends other image fusion levels as it retains more detailed image information and can achieve higher accuracy. Among the pixel-level fusion methods, the wavelet-based method is recommended because of its non-redundancy and directionality [24][25][26][27]. Figure 1 shows the image fusion process based on the wavelet transform. The fusion process is performed as follows: Firstly, the source images were decomposed by applying the wavelet transform with an appropriate wavelet basis function. In this procedure, reconstructed images by fd-and td-EIT can be assumed as a continuous change of two-dimensional signal ( , ) ∈ ( ), and are its abscissa and ordinate, respectively. The wavelet transform of ( , ) is expressed as: , , ( , ) = 1 ( − , − ) where ( , ) is a two-dimensional wavelet basis function, is the scaling factor, and is the displacement factor. By applying the basis function, the low-frequency and highfrequency components of the reconstructed image can be obtained, respectively. The lowfrequency components of the image give approximate features, while the high-frequency components provide detailed features [27]. Secondly, appropriate fusion rules are employed to fuse the high frequency and low frequency components. The choices of fusion rules are critical in this procedure. There are many fusion rules, such as larger absolute The fusion process is performed as follows: Firstly, the source images were decomposed by applying the wavelet transform with an appropriate wavelet basis function. In this procedure, reconstructed images by fd-and td-EIT can be assumed as a continuous change of two-dimensional signal f (x 1 , x 2 ) ∈ L 2 R 2 , x 1 and x 2 are its abscissa and ordinate, respectively. The wavelet transform of f (x 1 , x 2 ) is expressed as: where ψ(x 1 , x 2 ) is a two-dimensional wavelet basis function, a is the scaling factor, and b is the displacement factor. By applying the basis function, the low-frequency and highfrequency components of the reconstructed image can be obtained, respectively. The lowfrequency components of the image give approximate features, while the high-frequency components provide detailed features [27]. Secondly, appropriate fusion rules are employed to fuse the high frequency and low frequency components. The choices of fusion rules are critical in this procedure. There are many fusion rules, such as larger absolute value of coefficients, weighted average method, and local variance criterion [23]. According to previous studies, it is suggested that the fusion rules should adapt to the features of the source images. In this study, the fused image needs to maintain the position, size, edge, and shape of the target inclusions in the reconstructed images by fd-and td-EIT. In view of these requirements, the high frequency components were fused by taking the maximum absolute value in the source images. By combining the pixels correlation and regional variance in the image area, a hybrid fusion method was adopted for the fusion of low frequency components, as introduced by Wang et al. [27].
In the hybrid fusion rule, assuming that L(X) represents the low frequency coefficient matrix of the image X , let G(X, p) denote the regional variance of the region Q centred at p(m, n) in the low-frequency coefficient matrix of image X, which can be expressed as [27]: where L(X, p) represents the value of the element at (m, n) in the low-frequency component coefficient matrix, and u(X, p) represents the average value of the elements inside Q. w(q) is the weight coefficient, which would increase if it approached p(m, n). M 2 (p) is the matching degree of the regional variance regarding the low-frequency coefficient matrix of the source images A and B at p(m, n), which is defined as: It is noted that the value of M 2 (p) varies from 0 to 1. The smaller the M 2 (p), the lower the correlation between the low frequency coefficient matrices. Let t 2 (0.5 < t 2 < 1) be the threshold of the matching degree. If M 2 (p) < t 2 , the option fusion strategy is adopted: Otherwise, the average fusion strategy is used:

Image Quality Evaluation
Quantitative measures, e.g., relative coverage ratio (RCR), relative image error (RIE), and structural similarity (SSIM) index are defined to evaluate the quality of reconstructed images, which are the differences and correlations between the true image and the reconstructed image [28]: where CR denotes the coverage ratio defined as the ratio of the size of the inclusions to the total size of the sensing area. Correspondingly, CR True is CR of the true target distribution. g andĝ are the true and reconstructed images, respectively, µ x and µ y are the mean grey values in the true and reconstructed images, respectively, with σ x and σ y being the corresponding standard deviations, respectively. Quantitative evaluation of the reconstructed images of biological objects under different frequencies is challenging because the conductivity changes with frequencies and the ground truth is unknown, especially for the fd-EIT imaging where different images can be obtained with different frequency intervals. Therefore, the popular criteria, such as image errors and correlation coefficients, cannot be employed in this case.

Characterization of td-and fd-EIT Imaging by Experiment
A series of static experiments were conducted to characterize the td-and fd-EIT imaging and validate the feasibility of the proposed method.

Experimental Setup
An experimental system with a cylindrical tank sensor was established to characterize the fd-and td-EIT imaging, as shown in Figure 2. The diameter of the tank was 38 cm with four layers of EIT electrodes mounted axially. In each layer, 32 metallic electrodes (stainless steel screws, radius = 5 mm) were uniformly mounted on the inner surface of the tank. Sixteen of them (choosing one in the other two electrodes) in the second layer were employed for EIT imaging. A National Instruments (NI) PXI-based EIT system developed by Beihang University was employed in the experiment. The system mainly comprises three modules: the AC current source, the data collection system, and the PC with system control functions for image reconstruction. In the experiment, the cylindrical tank was filled with tap water as the background medium with a conductivity of 0.034 S/m. The pairwise injected current had a peak-to-peak magnitude of around 1 mA, while the frequency was set to 1 KHz, 5 KHz, 10 kHz, 20 kHz, and 50 KHz. With the adjacent excitation strategy, the excitation current was injected into a pair of adjacent electrodes. By multiplexing, the potential differences between each possible pair of adjacent electrodes were measured by the data collection system. The measurements were then sent to the PC for image reconstruction via a USB interface. Finally, an image of the target inclusions was reconstructed in MATLAB using the data received. For a single measurement, 20 frames of data were averaged to mitigate the measurement noise. mean grey values in the true and reconstructed images, respectively, with and being the corresponding standard deviations, respectively. Quantitative evaluation of the reconstructed images of biological objects under different frequencies is challenging because the conductivity changes with frequencies and the ground truth is unknown, especially for the fd-EIT imaging where different images can be obtained with different frequency intervals. Therefore, the popular criteria, such as image errors and correlation coefficients, cannot be employed in this case.

Characterization of td-and fd-EIT Imaging by Experiment
A series of static experiments were conducted to characterize the td-and fd-EIT imaging and validate the feasibility of the proposed method.

Experimental Setup
An experimental system with a cylindrical tank sensor was established to characterize the fd-and td-EIT imaging, as shown in Figure 2. The diameter of the tank was 38 cm with four layers of EIT electrodes mounted axially. In each layer, 32 metallic electrodes (stainless steel screws, radius = 5 mm) were uniformly mounted on the inner surface of the tank. Sixteen of them (choosing one in the other two electrodes) in the second layer were employed for EIT imaging. A National Instruments (NI) PXI-based EIT system developed by Beihang University was employed in the experiment. The system mainly comprises three modules: the AC current source, the data collection system, and the PC with system control functions for image reconstruction. In the experiment, the cylindrical tank was filled with tap water as the background medium with a conductivity of 0.034 S/m. The pairwise injected current had a peak-to-peak magnitude of around 1 mA, while the frequency was set to 1 KHz, 5 KHz, 10 kHz, 20 kHz, and 50 KHz. With the adjacent excitation strategy, the excitation current was injected into a pair of adjacent electrodes. By multiplexing, the potential differences between each possible pair of adjacent electrodes were measured by the data collection system. The measurements were then sent to the PC for image reconstruction via a USB interface. Finally, an image of the target inclusions was reconstructed in MATLAB using the data received. For a single measurement, 20 frames of data were averaged to mitigate the measurement noise. Td-and fd-EIT were used to image three different phantoms, which were set up by inserting a conductive inclusion (a peeled radish in phantom 1), a non-conductive inclusion (a nylon rod in phantom 2), and both of them (in phantom 3) at a certain location in the tank. The diameters of the peeled radish and nylon rod were 60 mm and 55 mm, respectively. Their length was sufficiently long to extend beyond the water surface so that all tests were conducted with target inclusions homogeneous in the vertical direction.  Td-and fd-EIT were used to image three different phantoms, which were set up by inserting a conductive inclusion (a peeled radish in phantom 1), a non-conductive inclusion (a nylon rod in phantom 2), and both of them (in phantom 3) at a certain location in the tank. The diameters of the peeled radish and nylon rod were 60 mm and 55 mm, respectively. Their length was sufficiently long to extend beyond the water surface so that all tests were conducted with target inclusions homogeneous in the vertical direction.  According to previous work [29,30], electrical impedance spectroscopy (EIS) can reveal the frequency response of electrical impedance of bio-tissues, which can be employed to demonstrate the frequency-dependencies of the conductivities of the peeled radish and nylon rod. Figure 4 shows the EIS measurements of the peeled radish and nylon rod by an impedance analyzer (MFLI Lock-in Amplifier-500 kHz/5 MHz, Zurich Instruments, Zurich, Switzerland) over the selected frequency range, i.e., 0-100 kHz. The results of EIS indicated that the conductivity of the radish increased linearly with the increase in frequency, and the linear fitting formula was: Its RMSE (root mean squared error) was 0.001834, and the value of R-square (coefficient of determination) was 0.9877.
In particular, it is noted that the conductivity of the radish at 10 kHz was about 0.03 S/m, which is nearly the same as the conductivity of the background medium. With respect to the nylon rod, the EIS data illustrate that its conductivity stayed almost unchanged with frequency. This is consistent with the previous analysis and verifies the feasibility of the use of the radish and nylon rod for the dedicated study on td-and fd-EIT imaging.

Experimental Results and Analysis
By applying the GREIT algorithms for image reconstruction, the td-and fd-EIT imaging results of the three phantoms were attained, and the size of the images w 64 64 pixels. In fd-EIT, 1 kHz was selected as the reference frequency. The voltage differences between the measured voltages at all others frequencies and at the reference frequency were procured for image reconstruction. In td-EIT, the measured voltage when the EIT sensor was filled with the background medium was taken as the reference, i.e., the meas- According to previous work [29,30], electrical impedance spectroscopy (EIS) can reveal the frequency response of electrical impedance of bio-tissues, which can be employed to demonstrate the frequency-dependencies of the conductivities of the peeled radish and nylon rod. Figure 4 shows the EIS measurements of the peeled radish and nylon rod by an impedance analyzer (MFLI Lock-in Amplifier-500 kHz/5 MHz, Zurich Instruments, Zurich, Switzerland) over the selected frequency range, i.e., 0-100 kHz. The results of EIS indicated that the conductivity of the radish increased linearly with the increase in frequency, and the linear fitting formula was: According to previous work [29,30], electrical impedance spectroscopy (EIS) can reveal the frequency response of electrical impedance of bio-tissues, which can be employed to demonstrate the frequency-dependencies of the conductivities of the peeled radish and nylon rod. Figure 4 shows the EIS measurements of the peeled radish and nylon rod by an impedance analyzer (MFLI Lock-in Amplifier-500 kHz/5 MHz, Zurich Instruments, Zurich, Switzerland) over the selected frequency range, i.e., 0-100 kHz. The results of EIS indicated that the conductivity of the radish increased linearly with the increase in frequency, and the linear fitting formula was: Its RMSE (root mean squared error) was 0.001834, and the value of R-square (coefficient of determination) was 0.9877.
In particular, it is noted that the conductivity of the radish at 10 kHz was about 0.03 S/m, which is nearly the same as the conductivity of the background medium. With respect to the nylon rod, the EIS data illustrate that its conductivity stayed almost unchanged with frequency. This is consistent with the previous analysis and verifies the feasibility of the use of the radish and nylon rod for the dedicated study on td-and fd-EIT imaging.

Experimental Results and Analysis
By applying the GREIT algorithms for image reconstruction, the td-and fd-EIT imaging results of the three phantoms were attained, and the size of the images w 64 64 pixels. In fd-EIT, 1 kHz was selected as the reference frequency. The voltage differences between the measured voltages at all others frequencies and at the reference frequency were procured for image reconstruction. In td-EIT, the measured voltage when the EIT sensor was filled with the background medium was taken as the reference, i.e., the meas- Its RMSE (root mean squared error) was 0.001834, and the value of R-square (coefficient of determination) was 0.9877.
In particular, it is noted that the conductivity of the radish at 10 kHz was about 0.03 S/m, which is nearly the same as the conductivity of the background medium. With respect to the nylon rod, the EIS data illustrate that its conductivity stayed almost unchanged with frequency. This is consistent with the previous analysis and verifies the feasibility of the use of the radish and nylon rod for the dedicated study on td-and fd-EIT imaging.

Experimental Results and Analysis
By applying the GREIT algorithms for image reconstruction, the td-and fd-EIT imaging results of the three phantoms were attained, and the size of the images w 64 × 64 pixels. In fd-EIT, 1 kHz was selected as the reference frequency. The voltage differences between the measured voltages at all others frequencies and at the reference frequency were procured for image reconstruction. In td-EIT, the measured voltage when the EIT sensor was filled with the background medium was taken as the reference, i.e., the measurement at t 1 , and the voltage differences were obtained when perturbations were present in the background medium under the tested frequencies. The td-and fd-EIT imaging results are shown in  conductivity of peeled radish was lower than the tap water when the frequency was lower than 10 kHz, and larger than the tap water at frequencies >10 kHz. On the other hand, since the conductivity of the radish was related to frequency, the image quality was sensitive to the frequency intervals. Besides, size overestimation and image distortion are more significant with increase in the frequency intervals, which is attributed to the decrease of SNR with the increase in frequency. The image distortion was the most serious at 10 kHz due to the conductivity of peeled radish being very close to the tap water and the measured voltages had a much smaller signal-to-noise ratio (SNR), which would lead to wrong imaging result since the EIT reconstruction was ill-conditioned.
fd-EIT    For further demonstration, imaging results of phantom 3 by td-and fd-EIT are shown Figure 7. In Figure 7, the position and shape of the nylon rod can be well reconstructed by td-EIT, while due to the existence of the non-conductive nylon rod, the radish can hardly be seen as it causes much smaller responses at the sensor boundary. However, regarding fd-EIT imaging, only the peeled radish is shown in the reconstructed images. The reason is that the conductivity of the radish varies with excitation frequency while that of the nylon rod does not. These findings are verified by the evaluation metrics in Table 3.    For further demonstration, imaging results of phantom 3 by td-and fd-EIT are shown Figure 7. In Figure 7, the position and shape of the nylon rod can be well reconstructed by td-EIT, while due to the existence of the non-conductive nylon rod, the radish can hardly be seen as it causes much smaller responses at the sensor boundary. However, regarding fd-EIT imaging, only the peeled radish is shown in the reconstructed images. The reason is that the conductivity of the radish varies with excitation frequency while that of the nylon rod does not. These findings are verified by the evaluation metrics in Table 3.    6 show the td-and fd-imaging results of phantoms 1 and 2 at 1 kHz, 5 kHz, 10 kHz, 20 kHz, and 50 kHz. Figures 5-7 have the same color bars, and the gray value range was 0-255, corresponding to the conductivity value of 0-0.08 S/m, respectively. The same applied to the relationship between the gray value and the conductivity in Figure 8, despite of the use of different colors for presentation. As indicated, td-EIT can reconstruct the position and shape of the single radish or nylon rob clearly (except for at 10 kHz where the conductivity of radish was close to the background medium as indicated by td-EIT imaging results). Fd-EIT could image the radish in phantom 1, but it failed to indicate the accurate position of the nylon rod in the reconstructed image since the conductivity of non-conductive nylon rod was not frequency-dependent. Regarding the phantoms 1 and 2, td-EIT provided better imaging results of the inclusions, as indicated in Tables 1 and 2 for quantitative evaluation, where td-EIT imaging results have lower RIE, as well as RCR and SSIM closer to 1. In Figure 5, td-EIT imaging results indicate that the color scale of the reconstructed inclusion changed with the increase in excitation frequency, which was due to the frequency-dependency of the conductivity of peeled radish. More specifically, the conductivity of peeled radish was lower than the tap water when the frequency was lower than 10 kHz, and larger than the tap water at frequencies >10 kHz. On the other hand, since the conductivity of the radish was related to frequency, the image quality was sensitive to the frequency intervals. Besides, size overestimation and image distortion are more significant with increase in the frequency intervals, which is attributed to the decrease of SNR with the increase in frequency. The image distortion was the most serious at 10 kHz due to the conductivity of peeled radish being very close to the tap water and the measured voltages had a much smaller signal-to-noise ratio (SNR), which would lead to wrong imaging result since the EIT reconstruction was ill-conditioned.    For further demonstration, imaging results of phantom 3 by td-and fd-EIT are shown Figure 7. In Figure 7, the position and shape of the nylon rod can be well reconstructed by td-EIT, while due to the existence of the non-conductive nylon rod, the radish can hardly be seen as it causes much smaller responses at the sensor boundary. However, regarding fd-EIT imaging, only the peeled radish is shown in the reconstructed images. The reason is that the conductivity of the radish varies with excitation frequency while that of the nylon rod does not. These findings are verified by the evaluation metrics in Table 3. Table 3. RCR, RIE, and SSIM regarding td-and fd-EIT imaging results of phantom 3.

Image Fusion with Respect to td-and fd-EIT Imaging
To evaluate the performance of the above-mentioned fusion strategy, the reconstructed images of phantom 3 by td-and fd-EIT at 5 kHz were taken as the source images. Meanwhile, to investigate the impact of wavelet basis function, various basis functions available in MATLAB were employed in the wavelet transform, including "coif3", "fk14", "bior3.7", "dmey", "sym6", and "db8".
As shown in Figure 8, the source images by td-and fd-EIT cannot either demonstrate the conductive inclusion (i.e., the radish) or the non-conductive inclusion (i.e., the nylon rod). By applying the wavelet-based image fusion strategy, the fused image can successfully recover both the conductive and non-conductive inclusions, which resolves the issues of tdand fd-EIT imaging. It is noted that different basis functions have significant impact on the quality of the fused images, especially for the radish. By employing the "coif3" or "dmey" basis function in the wavelet transform, the shapes of the nylon rob in the fused images are obviously distorted besides the shapes of the radish. Overall, the "fk14", "bior3.7", "sym6", and "db8" basis functions can provide fused images of reasonable quality, which are much better than those reconstructed by either td-EIT or fd-EIT.

Investigation for Improving Lung Ventilation Imaging with Thorax Model
To investigate the performance of the proposed method in medical EIT imaging of lung ventilation, simulations and in vitro experiments were conducted. A real-shaped thorax model, as shown in Figure 9, was built for FEM simulation. The thorax domain and the sub-domains of the lungs were created according to the build-in function of the thorax model in EIDORS [31]. It is well known that the whole respiratory cycle contains two stages: inspiration and expiration. In the simulation, the lung was simulated as varying states by setting different air contents in the given 2-D lung model. The model is assumed ideal as it does not consider the modelling errors regarding the non-homogenous background, the lung shape movement and the change of the chest due to the breathing in the real clinical situations. The following two test cases were designed: Case 1: Lung imaging in the inspiration (Figure 9b) stage, during which the breathing would increase the air content, and air would almost fully fill the lung in the end of inspiration. For this reason, it is suggested to set up a sufficiently large amount of air content inside the lung model to simulate the inspiration.
Case 2: Lung imaging in the expiration (Figure 9c) stage, during which the air content decreases with the breathing and becomes negligible in the end of expiration. In a similar fashion, a small amount of air content was set up inside the lung model to simulate the expiration.
In the simulation, lung conductivity was considered as a function of the excitation frequency, according to literature [32,33]. The maximum and minimum conductivity values were set to 0.079 S/m and 0.11 S/m, respectively. The conductivity of background tissues was set to 0.2 S/m, as suggested by Liu et al. [6]. Sixteen electrodes were attached on the boundary of the thorax domain. The same current injection and measurement strategy employed in the previous experimental studies was also adopted in the simulation. The simulation was conducted with 2-D FEM approximation in COMSOL. To produce a reference for fd-and td-EIT imaging, the reference frequency was chosen to be 1 kHz, and the reference distribution for td-EIT imaging was assumed to be the thorax model only containing the background medium and lung tissues, as shown in Figure 9a. Simulation results are shown in Figures 10 and 11. In the reconstructed images, whose size is 84 128 pixels, the outermost dark solid line and the innermost solid dark line denote the true lung shape and air region, respectively. In the simulated cases 1 and 2, the fd-and td-EIT imaging can successfully reconstruct the lung-shaped inclusion or the air inclusion, respectively, but not both of them. In addition, the air inclusion is overestimated by td-EIT, while the lung-shaped inclusion is underestimated by fd-EIT. In fact, the end of expiration would reduce the content of air in the lung, and hence result in a higher RCR (the air domain is significantly overestimated) in the case 2, as shown in Tables 4 and 5. Due to the limitation of the reconstruction algorithm, the sharp edge and accurate shape of the inclusions were not well preserved in the reconstructed images. Case 2: Lung imaging in the expiration (Figure 9c) stage, during which the air content decreases with the breathing and becomes negligible in the end of expiration. In a similar fashion, a small amount of air content was set up inside the lung model to simulate the expiration.
In the simulation, lung conductivity was considered as a function of the excitation frequency, according to literature [32,33]. The maximum and minimum conductivity values were set to 0.079 S/m and 0.11 S/m, respectively. The conductivity of background tissues was set to 0.2 S/m, as suggested by Liu et al. [6]. Sixteen electrodes were attached on the boundary of the thorax domain. The same current injection and measurement strategy employed in the previous experimental studies was also adopted in the simulation. The simulation was conducted with 2-D FEM approximation in COMSOL. To produce a reference for fd-and td-EIT imaging, the reference frequency was chosen to be 1 kHz, and the reference distribution for td-EIT imaging was assumed to be the thorax model only containing the background medium and lung tissues, as shown in Figure 9a.
Simulation results are shown in Figures 10 and 11. In the reconstructed images, whose size is 84 × 128 pixels, the outermost dark solid line and the innermost solid dark line denote the true lung shape and air region, respectively. In the simulated cases 1 and 2, the fd-and td-EIT imaging can successfully reconstruct the lung-shaped inclusion or the air inclusion, respectively, but not both of them. In addition, the air inclusion is overestimated by td-EIT, while the lung-shaped inclusion is underestimated by fd-EIT. In fact, the end of expiration would reduce the content of air in the lung, and hence result in a higher RCR (the air domain is significantly overestimated) in the case 2, as shown in Tables 4 and 5. Due to the limitation of the reconstruction algorithm, the sharp edge and accurate shape of the inclusions were not well preserved in the reconstructed images. Note that the fd-EIT images at 1 kHz cannot be obtained since it is the chosen reference frequency. Note that the fd-EIT images at 1 kHz cannot be obtained since it is the chosen reference frequency.
Frequency 1 kHz 5 kHz 10 kHz 20 kHz 50 kHz td-EIT fd-EIT fd-EIT  The reconstructed images by td-and fd-EIT at 5 kHz and 50 kHz were used to demonstrate the improvement of the lung imaging by the proposed fusion method. As shown in Figures 12 and 13, the image fusion for lung imaging was fairly successful. Compared with the imaging results by either td-EIT or fd-EIT, the fused images cannot only provide the rational changes of air content, but also the position and shape of the lungshaped inclusion. It is not trivial in practice to obtain the change of air content and lung shape information simultaneously since the latter can serve as a reference and there may, in addition to the air current, be tissue changes related to certain diseases such as pulmonary edema. Meanwhile, the size and shape of the air inclusion or lung-shaped inclusion in the fused image agree well with the actual ones. It is noted that the fused images at 50 fd-EIT  The reconstructed images by td-and fd-EIT at 5 kHz and 50 kHz were used to demonstrate the improvement of the lung imaging by the proposed fusion method. As shown in Figures 12 and 13, the image fusion for lung imaging was fairly successful. Compared with the imaging results by either td-EIT or fd-EIT, the fused images cannot only provide the rational changes of air content, but also the position and shape of the lungshaped inclusion. It is not trivial in practice to obtain the change of air content and lung shape information simultaneously since the latter can serve as a reference and there may, in addition to the air current, be tissue changes related to certain diseases such as pulmonary edema. Meanwhile, the size and shape of the air inclusion or lung-shaped inclusion in the fused image agree well with the actual ones. It is noted that the fused images at 50  The reconstructed images by td-and fd-EIT at 5 kHz and 50 kHz were used to demonstrate the improvement of the lung imaging by the proposed fusion method. As shown in Figures 12 and 13, the image fusion for lung imaging was fairly successful. Compared with the imaging results by either td-EIT or fd-EIT, the fused images cannot only provide the rational changes of air content, but also the position and shape of the lung-shaped inclusion. It is not trivial in practice to obtain the change of air content and lung shape information simultaneously since the latter can serve as a reference and there may, in addition to the air current, be tissue changes related to certain diseases such as pulmonary edema. Meanwhile, the size and shape of the air inclusion or lung-shaped inclusion in the fused image agree well with the actual ones. It is noted that the fused images at 50 kHz indicate a much larger image contrast in terms of the lung shape. This is in accordance with the finding that the conductivity variation of the lung tissues between the measurement frequency (5 kHz and 50 kHz) and the reference frequency (1 kHz) increases with frequency. kHz indicate a much larger image contrast in terms of the lung shape. This is in accordance with the finding that the conductivity variation of the lung tissues between the measurement frequency (5 kHz and 50 kHz) and the reference frequency (1 kHz) increases with frequency. Frequency To further validate the feasibility and effectiveness of the proposed fusion method, static in vitro experiments were conducted. These experiment was conducted in the same tank as in the previous experiments. The tank was also filled with tap water as the background medium, and two lung-shaped inclusions were made of white gourd. Two pairs of air-containing inclusions were inserted into the two lung-shaped inclusions to simulate the lung-filling conditions in the expiration and inspiration stages, respectively, as in the two simulation cases. Figures 14 and 15(a1,b1) show the phantom pictures. To quantify the air content inside the lung-shaped inclusion, the volume ratio was measured. The volume ratio was 3/20 and 6/20 for the expiration and inspiration setup, respectively.  kHz indicate a much larger image contrast in terms of the lung shape. This is in accordance with the finding that the conductivity variation of the lung tissues between the measurement frequency (5 kHz and 50 kHz) and the reference frequency (1 kHz) increases with frequency. Frequency To further validate the feasibility and effectiveness of the proposed fusion method, static in vitro experiments were conducted. These experiment was conducted in the same tank as in the previous experiments. The tank was also filled with tap water as the background medium, and two lung-shaped inclusions were made of white gourd. Two pairs of air-containing inclusions were inserted into the two lung-shaped inclusions to simulate the lung-filling conditions in the expiration and inspiration stages, respectively, as in the two simulation cases. Figures 14 and 15(a1,b1) show the phantom pictures. To quantify the air content inside the lung-shaped inclusion, the volume ratio was measured. The volume ratio was 3/20 and 6/20 for the expiration and inspiration setup, respectively. To further validate the feasibility and effectiveness of the proposed fusion method, static in vitro experiments were conducted. These experiment was conducted in the same tank as in the previous experiments. The tank was also filled with tap water as the background medium, and two lung-shaped inclusions were made of white gourd. Two pairs of air-containing inclusions were inserted into the two lung-shaped inclusions to simulate the lung-filling conditions in the expiration and inspiration stages, respectively, as in the two simulation cases. Figures 14 and 15(a1,b1) show the phantom pictures. To quantify the air content inside the lung-shaped inclusion, the volume ratio was measured. The volume ratio was 3/20 and 6/20 for the expiration and inspiration setup, respectively.
Based on the visual inspection of Figures 14 and 15, the fused images provide a reasonable estimation of the shape and location of the lung-shaped inclusion as well as the air inclusion inside it at different frequencies. The td-EIT performs well in reconstructing the air inclusion, while the fd-EIT is better in imaging the lung-shaped inclusion. It is noted that the contrast variation caused by the change of air content is clearly shown. Therefore, it can be concluded that the lung imaging based on the proposed fusion method is feasible and effective.    Based on the visual inspection of Figures 14 and 15, the fused images provide a reasonable estimation of the shape and location of the lung-shaped inclusion as well as the air inclusion inside it at different frequencies. The td-EIT performs well in reconstructing the air inclusion, while the fd-EIT is better in imaging the lung-shaped inclusion. It is Based on the visual inspection of Figures 14 and 15, the fused images provide a reasonable estimation of the shape and location of the lung-shaped inclusion as well as the air inclusion inside it at different frequencies. The td-EIT performs well in reconstructing the air inclusion, while the fd-EIT is better in imaging the lung-shaped inclusion. It is (b2,b3) are fd-EIT and td-EIT images, respectively; (b4-b9) are fused images using different wavelet bases, respectively.

Conclusions
This paper characterizes and compares fd-and td-EIT for the imaging of target distributions composed of bio-material inclusions and non-conductive inclusions, and the lung ventilation imaging is a potentially applicable scenario. Through an experimental study with a peeled radish or a nylon rod as targets, it was found that td-EIT reconstructions can reconstruct the bio-material inclusion (the radish) and the non-conductive inclusion (the nylon rod) if either the former or the latter exists in the sensing domain. Fd-EIT could only image the bio-materials inclusion (the radish) as the conductivity of non-conductive materials is not frequency-dependent, but has significant distortions with an appropriate frequency interval. However, both the td-EIT and fd-EIT failed to simultaneously reconstruct the radish and nylon rod if both of them exist in the sensing domain. This is because the response to the non-conductive rod (i.e., a much larger conductivity contrast to the background medium) would overwhelm that to the radish for td-EIT imaging, which would occur when imaging certain lung diseases such as pulmonary edema. To tackle with this dilemma, it is proposed to fuse the td-and fd-EIT imaging results using a dedicated wavelet-based fusion strategy, which can integrate the merits of td-and fd-EIT. Based on the simulation and experiment with lung-shaped inclusions or in vitro lung models, it was found that the proposed method has much better reconstruction accuracy than either td-or fd-EIT, which can simultaneously reconstruct both the lung-shaped inclusions and air inclusions in the fused image. Meanwhile, the proposed method can reduce the image distortions by td-or fd-EIT, and provide the relative shapes, sizes, and positions of those inclusions. In summary, the proposed method can improve the performance of EIT in lung ventilation imaging, so that the diagnosis functions of EIT can be expanded to deal with more lung-tissue-related diseases. Although such a fd-and td-EIT system certainly requires the hardware to accommodate to both the EIT modes, which can be readily achieved in practical applications, this paper aims to determine how much improvement can be expected by applying such a fusion method and understanding and identifying the compromises involved in fd-and td-EIT imaging. Future work will focus on optimizing the image reconstruction and fusion methods to further enhance the image quality. On other hand, the proposed method should be tested in more realistic cases of lung imaging.