We propose a novel framework to predict pacing sites in the left ventricle (LV) of a heart and its result can be used to assist pacemaker implantation and programming in cardiac resynchronization therapy (CRT), a widely adopted therapy for heart failure patients. In a traditional CRT device deployment, pacing sites are selected without quantitative prediction. That runs the risk of suboptimal benefits. In this work, the spherical harmonic (SPHARM) description is employed to model the ventricular surfaces and a novel SPHARM-based surface correspondence approach is proposed to capture the ventricular wall motion. A hierarchical agglomerative clustering technique is applied to the time series of regional wall thickness to identify candidate pacing sites. Using clinical MRI data in our experiments, we demonstrate that the proposed framework can not only effectively identify suitable pacing sites, but also distinguish patients from normal subjects perfectly to help medical diagnosis and prognosis.