Noninvasive imaging of cardiac electrical function has begun to move towards clinical adoption. Here we consider one common formulation of the problem, in which the goal is to estimate the spatial distribution of electrical activation times during a cardiac cycle. We address the challenge of understanding the robustness and uncertainty of solutions to this formulation. This formulation poses a non-convex, non-linear least squares optimization problem. We show that it can be relaxed to be convex, at the cost of some degree of physiological realism of the solution set, and that this relaxation can be used as a framework to study model inaccuracy and solution uncertainty. We present two examples, one using data from a healthy human subject and the other synthesized with the ECGSIM software package. In the first case, we consider uncertainty in the initial guess and regularization parameter. In the second case, we mimic the presence of an ischemic zone in the heart in a way which violates a model assumption. We show that the convex relaxation allows understanding of spatial distribution of parameter sensitivity in the first case, and identification of model violation in the second.