Angular correlations of identified particles measured in ultrarelativistic proton-proton (pp) and heavy-ion collisions exhibit a number of features which depend on the collision system and particle type under consideration. Those features are produced by mechanisms, such as (mini)jets, elliptic flow, resonance decays, and conservation laws. In addition, of particular importance are those related to the quantum statistics (QS) and final-state interactions (FSIs). In this paper we show how to unfold the QS and FSI contributions in angular correlation functions by employing a Monte Carlo approach and using momentum correlations (femtoscopy), focusing on pp reactions. We validate the proposed method with PYTHIA 8 Monte Carlo simulations of pp collisions at $\sqrt{s}=13$ TeV coupled to calculations of QS and FSI effects with the Lednick\'y and Lyuboshitz formalism and provide predictions for the unfolded effects. In particular, we show how those effects modify the shape of the angular correlation function with emphasis on pions and protons. Most importantly, specific structures observed in the near-side region for both baryon-baryon and baryon-antibaryon pairs, originating from the strong interaction, are unveiled with the proposed method.