We review three data stacking methods developed for attaining the splitting parameters of normal modes of the Earth's free oscillation, namely, SHS (spherical harmonic stacking), OSE (optimal sequence estimation), and MSE (multistation experiment), that take advantage of the spatial orthogonality of the (scalar and vector) spherical harmonic functions. We further extend these methods to broader applicability as appropriate, in particular SHS extended to matrix‐SHS (MSHS) to better accounting for the nonglobal distribution of recording stations, and OSE and MSE extended to transverse components and arbitrary harmonic degrees. We conduct synthetic numerical experiments and conclude that OSE and MSHS stackings yield superior results. Based on the lessons learned, we apply the methods on real global seismic records after the 2004 Sumatra earthquake and report the complete resolution of all singlets for 13 selected multiplets: 8 are mantle‐sensitive (1S2, 0S4, 1S4, 0S5, 0S6, 1T2, 0T4, and 0T6) and 5 inner‐core sensitive and of anomalous splitting (13S2, 10S2, 2S3, 3S2, and 11S1), among them the singlet resolution of 1T2, 0T4, 10S2, 2S3, 3S2, and 11S1 are reported for the first time.