Sensitivity Analysis for Optimizing Electrical Impedance Tomography Protocols

Claire Onsager,Chulin Wang,Charles Costakis,Can Aygen,Lauren Lang,Suzan van der Lee,Matthew A. Grayson

DOI: https://doi.org/10.48550/arXiv.2111.01397

2021-11-18

Abstract:Electrical impedance tomography (EIT) is a noninvasive imaging method whereby electrical measurements on the boundary of a conductive medium (the data) are taken according to a prescribed protocol set and inverted to map the internal conductivity (the model). This paper introduces a sensitivity analysis method and corresponding inversion and protocol optimization that generalizes the criteria for tomographic inversion to minimize the model-space dimensionality and maximize data importance. Sensitivity vectors, defined as rows of the Jacobian matrix in the linearized forward problem, are used to map targeted conductivity features from model-space to data-space, and a volumetric outer-product of these vectors in model-space called the sensitivity parallelotope volume provides a figure-of-merit for data protocol optimization. Orthonormal basis functions that accurately constrain the model-space to features of interest can be defined from a priori information. By increasing the contact number to expand the number of possible measurements Dmax, and by reducing the model-space to a minimal number M0 of basis functions that describe only the features of interest, the M0 << Dmax sensitivity vectors of greatest length and maximal orthogonality that span this model-space can be identified. The reduction in model-space dimensionality accelerates the inversion by several orders of magnitude, and the enhanced sensitivity can tolerate noise levels up to 1,000 times larger than standard protocols.

Medical Physics

What problem does this paper attempt to address?

The problems that this paper attempts to solve are two long - standing challenges in Electrical Impedance Tomography (EIT):

1. **Long computation time**: The computation time of existing algorithms will increase geometrically when the number of model space parameters is increased. In order to achieve a fast EIT refresh rate, resolution is usually sacrificed.
2. **Sensitivity to experimental noise**: Current EIT measurement protocols include small - signal measurements, which are vulnerable to experimental noise, thus reducing the fidelity of the inverse problem and decreasing the resolution.

To solve these problems, the paper introduces a new sensitivity analysis method to optimize the EIT protocol in the following ways:

- **Reducing the dimension of the model space**: By selecting the minimum number of basis functions that only contain the required resolution or features of interest, the dimension of the model space is reduced, thereby accelerating the computation speed.
- **Improving robustness to noise**: By excluding low - sensitivity measurements from the EIT protocol, the tolerance to noise is enhanced, and a noise level 1,000 times higher than that of the standard protocol can be tolerated.

Specifically, the paper proposes the following key steps:

1. **Selection of the model space**:
- Determine the reference model based on prior knowledge to reduce the number of computational steps required to converge to the final solution.
- Identify the minimum number of basis functions that describe the required features to minimize the computational load.

2. **Assignment of contact points**:
- Determine the number of electrical contact points \( C \) required to provide the required resolution or features, such that the number of independent measurements \( D_I \) far exceeds the number of measurements \( D_0 \) in the protocol.
- It is recommended to use a compliance factor \( c\geq3 \) to ensure that the protocol can be optimized using large - signal measurements.

3. **Definition of the measurement space**:
- Define all possible measurements \( D_{\text{max}} \) as the measurement space and introduce a metric method to evaluate the sensitivity of different measurements.
- By calculating the sensitivity vector of each measurement, select the measurement with the largest volume of the sensitivity parallelepiped as the optimized protocol.

4. **Selection of the data space**:
- Select the most sensitive and linearly independent measurements from the measurement space to form the data space to maximize sensitivity and linear independence.

Through these methods, the paper aims to significantly improve the computational efficiency and robustness to noise of EIT, thereby promoting the application of EIT in fields such as medicine, engineering, and geoscience.