A direct substitution, equation error technique for solving the thermographic tomography problem

A new technique for solving the combined state and parameter estimation problem in thermographic tomography is presented. The technique involves the direct substitution of known skin temperatures into the finite difference form of the bio-heat transfer equation as formulated for solving an initial value problem with a convection boundary condition at the skin surface. These equations are then used to solve the inverse bio-heat transfer problem for the unknown subcutaneous tissue temperatures and physiological parameters. For a small number of nodal points, closed form algebraic solutions are obtained. For larger sets of equations, a hybrid technique is used in which the problem is initially posed as an unconstrained optimization problem in which the model equation error is minimized using the conjugate gradient descent technique to get close to a solution. Then a generalized Newton-Raphson technique was used to solve the equations. A numerical simulation of a one-dimensional problem is investigated both with and without noise superimposed on the input (transient) skin temperature data. The results show that the technique gives very accurate results if the skin temperature data contains little noise. It is also shown that if the physical properties of the tissue and the metabolism are known, that a given set of proper transient skin temperature inputs yields a unique solution for the unknown internal temperatures and blood perfusion rates. However, the similar problem with known blood perfusion rates and unknown metabolisms does not yield a unique solution for the internal temperatures and metabolisms.