The temporal prior in bioelectromagnetic source imaging problems

Abstract

The multiplicity of temporal priors proposed for regularization of the bioelectromagnetic source imaging problems [e.g., the inverse electrocardiogram (ECG) and inverse electroencephalogram (EEG) problems], is discordant with the fact that fundamental statistical principles sharply limit the choice. Thus, our objective is to derive the form of the prior consistent with the general unavailability of temporal constraints. Writing linear formulations of the inverse ECG and inverse EEG problems as H = FG + N (where the ith columns of matrices H, G, and N, are data, signal, and noise vectors at time step i, and F is the transfer matrix), and using the noninformative principle that features of the spatiotemporal prior not supplied a posteriori should be invariant under temporal transformations, we show that the implied spatiotemporal signal autocovariance matrix (of the vector formed by the entries of G) is given in block matrix form [equation in text] where Cg is a matrix of unit trace proportional to the autocovariance matrix of any column of G (representing supplied information regarding the spatial prior), epsilon[.] denotes expectation, superscript ' indicates transpose, [symbol in text] is the Kronecker product, [symbol in text] is Frobenius norm, and the "matrix scalar product" [symbol in text] indicates the inner product of the two vectors formed by the entries of the two adjacent matrices (i.e., A [symbol in text] B [triple bond] trace[A'B]). This result eliminates some uncertainties and ambiguities that have characterized spatiotemporal regularization methods--including eight methods previously introduced in this transactions. Ultimately, the result derives from an implied symmetry principle under which the form of a nontrivial noninformative temporal component of the prior can be identified. Among other things, separability of the spatiotemporal prior in terms of the above Kronecker product can be thought of as the expression of the lack of "entanglement" of the spatial and temporal contributions (a consequence of noninformativity). The approach is generalized to the important cases of non-Gaussian spatial priors, and signal and noise that are not independent (transfer matrix noise). We also demonstrate a means for computational complexity reduction, related to the application of a particular orthogonal transformation, having features dependent on whether or not the transfer matrix represents a surjective mapping.

Performance of Inverse Problem Regularization Methods for Driver Location during Atrial Fibrillation

“…Temporal correlation can also be included in the problem by using the isotropy assumption [13]. Then, the spatio-temporal covariance matrix can be computed as C X = C t ⊗ C x , where C x is the spatial covariance matrix and C t is the temporal covariance matrix.…”

A Kalman filter-based approach to reduce the effects of geometric errors and the measurement noise in the inverse ECG problem

“…Some of these methods include Tikhonov regularization [5], truncated singular value decomposition [40], Bayesian Maximum A Posteriori estimation [43,36], Laplacian weighted minimum norm [12], and genetic algorithms [19]. Recent studies have employed spatio-temporal approaches such as the multiple constraints method [4], state-space models (Kalman filter) [24,6,3,1], and the isotropy assumption [9]. In most of those studies, only the additive measurement noise is considered; the geometric errors in the forward model relating the epicardial potentials to the BSPMs are neglected.…”

“…In [23], the L-curve, CRESO, and zero crossing methods were compared to determine the regularization parameter for the zeroth order Tikhonov regularization in the presence of geometric errors. Greensite considered geometric noise in the transfer matrix in his formulation, but did not provide results for this case [9]. In [37], truncated total least squares algorithm was used to reduce the effects of geometric errors.…”

Imaging the electrical activity of the heart using a Kalman filter based approach: Comparison of results using different STM's

“…Some of these techniques use only spatial information but in order to represent the electrical behavior of the heart properly, it is more appropriate to employ both spatial and temporal constraints. These techniques are called spatio-temporal solutions [5,6,7].…”