340 likes | 516 Views
Applying Constraints to the Electrocardiographic Inverse Problem. Rob MacLeod Dana Brooks. Electrocardiography. Electrocardiographic Mapping. Bioelectric Potentials Goals Higher spatial density Imaging modality Measurements Body surface Heart surfaces Heart volume.
E N D
Applying Constraints to the Electrocardiographic Inverse Problem Rob MacLeodDana Brooks
Electrocardiographic Mapping • Bioelectric Potentials • Goals • Higher spatial density • Imaging modality • Measurements • Body surface • Heart surfaces • Heart volume
Body Surface Potential Mapping Taccardi et al, Circ., 1963
Cardiac Mapping • Coverage • Sampling Density • Surface or volume
A Forward B Inverse Inverse Problems in Electrocardiography + - + - + - C
Definition Estimate sources from remote measurements Motivation Noninvasive detection of abnormalities Spatial smoothing and attenuation Forward Inverse Epicardial Inverse Problem
Forward problem Epicardial/Endocardial Activation Time Geometric Model Body Surface Potentials Inverse problem Forward/Inverse Problem Thom Oostendorp, Univ. of Nijmegen
Sample Problem: Activation Times Measured Thom Oostendorp, Univ. of Nijmegen Computed
Elements of the Inverse Problem • Components • Source description • Geometry/conductivity • Forward solution • “Inversion” method (regularization) • Challenges • Inverse is ill-posed • Solution ill-conditioned
Inverse Problem Research • Role of geometry/conductivity • Numerical methods • Improving accuracy to clinical levels • Regularization • A priori constraints versus fidelityto measurements
Regularization • Current questions • Choice of constraints/weights • Effects of errors • Reliability • Contemporary approaches • Multiple Constraints • Time Varying Constraints • Novel constraints (e.g., Spatial Covariance) • Tuned constraints
Problem formulation Constraint Solution Tikhonov Approach
For k constraints with solution Multiple Constraints
For two spatial constraints: Dual Spatial Constraints Note: two regularization factors required
Redefine y, h, A: And write a new minimization equation: Joint Time-Space Constraints
General solution: For a single space and time constraint: Note: two regularization factors and implicit temporal factor Joint Time-Space Constraints
Determining Weights • Based on a posteriori information • Ad hoc schemes • CRESO: composite residual and smooth operator • BNC: bounded norm constraint • AIC: Akaike information criterion • L-curve: residual norm vs. solution seminorm
Rhl Thl Ahl- y L-Surface • Natural extension of single constraint approach • “Knee” point becomes a region
Joint Regularization Results with Fixed Laplacian Parameter RMSError Energy Regularization Parameter with Fixed Energy Parameter RMSError Laplacian Regularization Parameter
Admissible Solution Approach Admissible Solution Region Constraint 3 (differentiable) Constraint 1 (non-differentiable but convex) Constraint 4 (differentiable) Constraint 2 (non-differentiable but convex)
Define f(x) s.t. with the constraint such that that satisfies the convex condition Single Constraint
Define multiple constraints fi(x) so that the set of these represents the intersection of all constraints. When they satisfy the joint condition Then the resulting x is theadmissible solution Multiple Constraints
Examples of Constraints • Residual contsraint • Regularization contstraints • Tikhonov constraints • Spatiotemporal contraints • Weighted constraints • Novel constraints
Subgradients Ci Ci+1 + + Ei Ei+1 Constraint set Normal hyperplane Ellipsoid Algorithm
Admissible Solution Results AdmissibleSolution Original Regularized
New Opportunity • Catheter mapping • provides source information • limited sites
Endocardial Epicardial Venous Catheter Mapping
New Opportunity • Catheter mapping • provides source information • limited sites • Problem • how to include this information in the inverse solution • where to look for “best” information • Solutions? • Admissible solutions, Tikhonov? • Statistical estimation techniques
CVRTI Bruno Taccardi Rich Kuenzler Bob Lux Phil Ershler Yonild Lian CDSP Dana Brooks Ghandi Ahmad Agknowledgements