The order estimation algorithm relies on multidimensional folding combined with ESTimation ERror criterion.
Multidimensional folding and unfolding construct two matrix slices that satisfy the shift invariance equality.
The ESTER criterion estimates number of harmonics based on shift invariance equality.
The algorithm is able to identify significantly more harmonics.
The algorithm is robust against colored noise.