A new method for incorporating weighted temporal and spatial smoothing in the inverse problem of electrocardiography

In this paper, we present a method for incorporating temporal smoothing (TS) into the estimate of epicardial potentials from body surface potential data. Our algorithm employs a different spatial smoothing parameter, chosen by the composite residual error and smoothing operator criteria, at each time step in the sequence. The total spatial smoothing term is then simply partitioned between temporal and spatial smoothing. The algorithm appears to be quite robust with regard to this partitioning. The new method was evaluated in the setting of additive Gaussian noise, but otherwise realistic conditions of body geometry and reference epicardial potentials. In examining the match between estimated and measured electrograms, or the match between estimated isopotential maps and measured isopotential maps, the estimates constructed using the new TS algorithm produced consistently smaller relative errors than those constructed using a quasi-static (QS) algorithm or those constructed by postprocessing the QS estimate with a moving average filter.