POD based higher order spectral estimation