A method includes obtaining multiple local activation times (LATs) at different respective measurement locations on an anatomical surface of a heart. The method further includes computing respective directions of electrical propagation at one or more sampling locations on the anatomical surface, by, for each sampling location, selecting a respective subset of the measurement locations for the sampling location, constructing a set of vectors, each of at least some of the vectors including, for a different respective measurement location in the subset, three position values derived from respective position coordinates of the measurement location and an LAT value derived from the LAT at the measurement location, and computing the direction of electrical propagation at the sampling location based on a Principal Component Analysis (PCA) of a 4×4 covariance matrix for the set of vectors. The method further includes indicating the directions of electrical propagation on a display.