In this paper, a novel methodology is presented for the construction of global effective rotation-vibration spectroscopic models from potential energy and transition moment ab initio surfaces. Non-empirical effective Hamiltonians are obtained via the block-diagonalization of selected variationally-computed eigenvector matrices, with the possibility of matching observation without any fit. For the first time, the derivation of effective transition moment operators is also carried out in a systematic way. This general approach can be implemented quite easily in most of the variational computer codes and turns out to be a clear alternative to the rather involved Van Vleck perturbation method. Symmetry is exploited at all stages to translate first-principles calculations into a set of spectroscopic parameters to be further refined on experiment. Using the strategy proposed in this work, we demonstrate on H2CO, PH3, CH4, C2H4 and SF6 how the analysis of infrared spectra can be drastically simplified compared to a more traditional empirical approach and why completeness is no longer an issue.