The data analysis was performed using the same methodology reported in the previous study . Briefly, the spectra were first processed using asymmetric least square (ALS) baseline correction. Then, baselined spectra were vector normalized and smoothed using Savitzky-Golay filtering (4 th order polynomial, with a frame size of 37). Finally, multivariate analysis techniques and classification algorithms were applied in the spectral range 400-1650 cm -1 . Data processing was carried out using Matlab 2021b.