In this paper, we consider an inverse rough surface scattering problem in near-field optical imaging. This problem is actually to reconstruct the scattering surface as well as its impedance coefficient from multifrequency near-field data, and can be reduced into an integral scheme by employing an integral representation. We solve this integral scheme by a non-linear integral equation method, and further develop a fast inversion algorithm for reconstructing both the rough surface and the impedance coefficient. Numerical experiments are presented...