We show by an extensive method of quasi-discrete multiple-scale approximation that nonlinear multi-dimensional lattice waves subjected to intersite and external on-site potentials are found to be governed by (N +1)-dimensional nonlinear Schrödinger (NLS) equation. In particular, the resonant mode interaction of (2+1)-dimensional NLS equation has been identified and the theory allows the inclusion of transverse effect. We apply the exponential function method to the (2+1)dimensional NLS equation and obtain the class of soliton solutions with a purely algebraic computational method. Notably, we discuss in detail the effects of the external on-site potentials on the explicit form of the soliton solution generated recursively. Under the action of the external on-site potentials, the model presents a rich variety of oscillating multidromion patterns propagating in the system.