Motivated by computing p-values for multiple testing of local maxima in signal and change point detections, we study the height distribution of local maxima of smooth isotropic Gaussian random fields parameterized on Euclidean space or spheres. The obtained formulae hold in general in the sense that there are no restrictions on the covariance function of the field except for smoothness and...