Macro velocity model plays a key role in seismic imaging and inversion. The performance of traditional velocity analysis methods is degraded by multiples and amplitude-versus-offset (AVO) anomalies. Local event slopes, containing the subsurface velocity information, have been widely used to accomplish common time-domain seismic processing, imaging and velocity estimation. In this paper, we propose a method for velocity analysis with probability density function (PDF) related to local event slopes. We first estimate local event slopes with phase information in the Fourier domain. An adaptive filter is applied to improve the performance of slopes estimator in the low signal-to-noise ratio (SNR) situation. Second, the PDF is approximated with the histogram function, which is related to attributes derived from local event slopes. As a graphical representation of the data distribution, the histogram function can be computed efficiently. By locating the ray path of the first arrival on the semblance image with straight-ray segments assumption, automatic velocity picking is carried out to establish velocity model. Unlike local event slopes based velocity estimation strategies such as averaging filters and image warping, the proposed method does not make the assumption that the errors of mapped velocity values are symmetrically distributed or that the variation of amplitude along the offset is slight. Extension of the method to prestack time-domain migration velocity estimation is also given. With synthetic and field examples, we demonstrate that our method can achieve high resolution, even in the presence of multiples, strong amplitude variations and polarity reversals.