<p>We address the inverse problem of simultaneous reconstruction of the parameters and coupling functions of a discrete single-species metapopulation from measured data of network dynamics. The unknown parameters include growth rates, total migration fractions and the entries of a circulating coupling matrix. From a practical point of view, the parameters are estimated by solving a least squares problem, using input data consisting of measurements of sub-population densities. The method of solution is based on a trust region reflective algorithm which allows the user to provide upper and lower limits on the system parameters. In contrast with many existing methods, our reconstruction procedure does not require any previous linearization, simplification or knowledge of the dynamics behaviour of the measured population densities. Another advantage of our approach is that it enables the simultaneous reconstruction by taking into account the entire metapopulation dynamics. Furthermore, the method is robust as it works well using input data with or without additive noise. The effectiveness of the method has been verified in reconstruction problems involving arbitrarily complex dynamics, such as chaotic or periodic, stationary, synchronous or asynchronous. Several numerical results are presented and discussed.</p>