On the basis of a surface shape appearing at a predetermined point moment by moment during rotation of a golf ball having numerous dimples on its surface, a data constellation regarding a parameter dependent on a surface shape of the golf ball is calculated. A preferable parameter is a distance between an axis of the rotation and the surface of the golf ball. Another preferable parameter is a volume of space between a surface of a phantom sphere and the surface of the golf ball. Fourier transformation is performed on the data constellation to obtain a transformed data constellation. On the basis of a peak value and an order of a maximum peak of the transformed data constellation, an aerodynamic characteristic of the golf ball is determined. The peak value and the order of the maximum peak are calculated for each of PH rotation and POP rotation.