With growing demands on product performance and growing complexity of engineering structures, efficient tools for analyzing their dynamic behavior are essential. Linear techniques are well developed and often utilized. However, sometimes the errors due to linearization are too large to be acceptable, making it necessary to take nonlinear effects into account. In many practical applications it is common and reasonable to assume that the nonlinearities are highly local and thus only affect a limited set of spatial coordinates. The purpose of this paper is to present an approach to finding the spatial location of nonlinearities from measurement data, as this may not always be known beforehand. This information can be used to separate the underlying linear system from the nonlinear parts and create mathematical models for efficient parameter estimation and simulation. The presented approach builds on the reverse-path methodology and utilizes the coherence functions to determine the location of nonlinear elements. A systematic search with Multiple Input/Single Output models is conducted in order to find the nonlinear functions that best describe the nonlinear restoring forces. The obtained results indicate that the presented approach works well for identifying the location of local nonlinearities in structures. It is verified by simulation data from a cantilever beam model with two local nonlinearities and experimental data from a T-beam experimental set-up with a single local nonlinearity. A possible drawback is that a relatively large amount of data is needed. Advantages of the approach are that it only needs a single excitation point that response data at varying force amplitudes is not needed and that no prior information about the underlying linear system is needed.