The proposed technique allows one to efficiently and quickly retrieve the dispersion diagram of a metasurface from its numerically simulated scattering parameters without resorting to a time-consuming numerical eigensolver. It may also be used to find the metasurface structure corresponding to desired surface-wave propagation characteristics.