Estimating three-dimensional temperature fields during hyperthermia: studies of the optimal regularization parameter and time sampling period

During hyperthermia therapy it is desirable to know the entire temperature field in the treatment region. However, accurately inferring this field from the limited number of temperature measurements available is very difficult, and thus state and parameter estimation methods have been used to attempt to solve this inherently ill-posed problem. To compensate for this ill-posedness and to improve the accuracy of this method, Tikhonov regularization of order zero has been used to significantly improve the results of the estimation procedure. It is also shown that the accuracies of the temperature estimates depend upon the value of the regularization parameter, which has an optimal value that is dependent on the perfusion pattern and magnitude. In addition, the transient power-off time sampling period (i.e., the length of time over which transient data is collected and used) influences the accuracy of the estimates, and an optimal sampling period is shown to exist. The effects of additive measurement noise are also investigated, as are the effects of the initial guess of the perfusion values, and the effects of both symmetric and asymmetric blood perfusion patterns. Random perfusion patterns with noisy data are the most difficult cases to evaluate. The cases studied are not a comprehensive set, but continue to show the feasibility of using state and parameter estimation methods to reconstruct the entire temperature field.