We present a novel procedure for analyzing the lattice-QCD spectrum via the finite-volume formalism to obtain constraints on multi-hadron scattering amplitudes at both real and complex energies. This approach combines a Bayesian reconstruction of the scattering amplitude on the real axis with Nevanlinna interpolation for analytic continuation to complex-valued energies. The method is non-parametric, inherently accounting for parametrization dependence within the uncertainty. We demonstrate the applicability of this approach using real lattice-QCD data in resonant systems from the HadSpec and BaSc collaborations.