This only works if the polygons reasonably approximate the underlying surface.
A better approximation can be found using a clustering analysis of the normals on the unit sphere.