3D shape classification plays a fundamental role in geometric big data analysis. This work proposes a novel method for shape classification based on optimal mass transportation theory. The Riemann surfaces are mapped onto canonical domains conformally based on uniformization theorem. The conformal factor function is treated as probability distribution on the canonical domain. For each pair of probability distributions, the optimal mass transportation map is computed by solving Monge-Amper´e equation. The transportation cost is theWasserstein distance between two distributions. By using this distance, geometric classification based on clustering can be performed. The method is applied to 3D human facial expression recognition, which demonstrates the efficiency and efficacy of the method.