Korean Journal of Food Preservation
The Korean Society of Food Preservation
Article

# Hyperspectral imaging technigue for monitoring moisture content of blueberry during the drying process

Ji-Young Choi1, Jiyoon Kim1, Jungsoo Kim1, Saeul Jeong1, Minhyun Kim1, Sanghyeok Park1, Kwang-Deog Moon1,*
1School of Food Biotechnology, Kyungpook National University, Daegu 41566, Korea
*Corresponding author. E-mail：kdmoon@knu.ac.kr, Phone：+82-53-950-5773, Fax：+82-53-950-6772

Copyright © The Korean Society of Food Preservation. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: May 19, 2020; Revised: Oct 28, 2020; Accepted: Nov 16, 2020

Published Online: Jul 31, 2021

## Abstract

Changes in the moisture content (MC) of blueberries during drying was monitored by hyperspectral image analysis, and the degree of drying was determined using the partial least squares (PLS) model. Blueberries (n=820) were dried at 35°C for 0 (control), 3, 6, 9 and 12 days. The PLS discriminant analysis prediction accuracy of smoothing the pre-processed data was the highest. Regression coefficients were high at 706, 790, 827, 868, and 894 nm, corresponding to water molecules and carbohydrates (830-840 nm). To develop a prediction model for blueberry MC, 150 hyperspectral images were obtained from 30 samples per group. The MC of each group was also analyzed. The accuracy of the MC prediction model pretreated by the multiplicative scatter correction method was the highest at 0.9302. As indicated by Pearson’s correlation analysis, the blueberry MC showed a high correlation of 0.95 with the total soluble solid contents, brightness, and total flavonoid contents. These results suggest that hyperspectral imaging techniques can be used to predict and monitor various quality characteristics as well as the MC of blueberries during drying.

Keywords: blueberry; moisture content; partial least-squares regression; Pearson’s correlation analysis; prediction

## Introduction

Blueberries (Ericaceae and Chlorella) are mainly distributed in North America but cultivated throughout the world (Westwood, 1993). The fruits are regarded as one of the world's top 10 super foods by the New York Times (Hwang and Ko, 2010), and can be consumed raw or incorporated into various processed foods, such as jams, jellies, muffins and cakes, as functional ingredients to attract consumer interest (Lee and Lee, 2007; Lee et al., 2008). The anthocyanins, flavonoids and lycopene present in blueberries are known to have excellent antioxidant properties, enhance vision and immunity, and aid in the prevention of stroke (Park et al., 2012). As the domestic production and import of blueberries increase, farmers are spurring the development of processed foods to promote consumption.

Blueberries are soft and prone to mechanical damage during harvesting and transportation (Xu et al., 2015). Due to their high water content, blueberries are also susceptible to food degradation factors, such as microbial growth, biochemical changes, and chemical reactions. Drying is widely used to overcome this issue (Park et al., 2014). Drying reduces the moisture content (MC) of food, which aids in preventing spoilage and degradation by microorganisms and enzymes, thereby contributing to improved shelf life and expanding the options of food product development (Shin et al., 2015).

Dried blueberry production technology has been well studied in the context of drying conditions, such as the drying method and drying temperature, because the MC of dried blueberries is an important indicator that determines the texture and quality of the product. As, all dry products must have a uniform MC non-destructive techniques are desired for measuring the quality of dried blueberries.

Hyperspectral imaging (HSI) is a non-destructive method. This technique combines conventional imaging and spectroscopy techniques and can provide both spatial and spectral information from the target, which is extremely useful for evaluating individual food items (Qin et al., 2013). HSI is an invaluable technique for investigating the local distribution of components and for monitoring particular chemicals using near-infrared rays (Huang et al., 2014). To monitor the quality, a large amount of food must be checked in real time. In conventional quality measurement, the refractometer is often used for total soluble solids (TSS) content and the oven for water content; therefore, only a fraction of the target can be monitored, and the analyzed food is damaged and becomes inedible. Other advantages of non-destructive techniques are that there is no need for pre-treatment of samples, and many foods can be analyzed at once within a short time. In addition, a non-destructive technique is less affected by the proficiency of the researcher. This system is more practical for the food industry because it is suitable for conveyor systems (Iqbal et al., 2014). Gowen et al. (2007) demonstrated the usability of HSI in the visible and near-infrared (NIR) wavelength regions (350-1700 nm) for assessing quality parameters of agricultural (corn, wheat, rye, lentils, soya bean, canola) and horticultural (apple, potato, cucumber, citrus, peach, strawberry, and cherry) crops. Recently, hyperspectral image analysis and chemometrics have been used to predict the MC of various foods (Lin et al., 2019; Wei et al., 2019; Xu et al., 2019). Pu and Sun (2015) visualized the distribution of moisture during microwave vacuum drying of mango slices using Vis-NIR HSI. To non-destructively monitor and control the drying process, Nguyen-Do-Trong et al. (2018) evaluated the quality of dried banana slices using VNIR hyperspectral reflectance imaging. The color and MC of vegetable soybeans during drying were also predicted using hyperspectral imaging technology (Huang et al., 2014)

There have been studies on the prediction of internal bruising, maturity, foreign matter, nutraceutical content and firmness of blueberries using NIR spectroscopy (Fan et al., 2017; Leiva-Valenzuela et al., 2014; Sugiyama et al., 2010; Yang et al., 2014). Although there have been many studies using HSI to analyze raw blueberries, there have been no studies on the quality characteristics of dry-processed products.

Therefore, this study aimed to (i) determine whether the degree of blueberry drying could be predicted through HSI and (ii) develop a prediction model for the MC of dried blueberries during drying using HSI, which is a non- destructive and rapid detection tool.

## Materials and methods

Chemicals

Ethanol (EtOH; Duksan Pure Chemical Company, Ansan, South Korea) was used to extract samples for antioxidant activity analysis. Folin-Ciocalteu reagent (Junsei Chemical Company, Tokyo, Japan) and diethylene glycol (Samchun Pure Chemicals Company, Pyeongtaek-si, South Korea) were used to analyze the content of total phenolic compounds (TPC) and total flavonoid content (TFC), using gallic acid (Sigma Chemical, Co., St Louis, MO, USA) and naringin (Sigma Chemical Co.) as references.

Sample preparation

Blueberries were harvested in May and June 2019 in Korea, and they were purchased from local markets in Daegu. In order to reduce the error of the experiment, only those in good condition and weighing 1.9±0.2 g were selected and used. Shin et al. (2015) reported that the lower the drying temperature, the brighter and more purple the fruit, and the less the change in color; therefore, the blueberry in this study were dried at 35°C-40°C using a climatic chamber (JSRH-150CP, JS Research Inc., Gongju-si, Korea). The blueberries were divided into five groups based on drying time (days): 0 (BB-0), 3 (BB-3), 6 (BB-6), 9 (BB-9) and 12 (BB-12). To accurately determine the degree of drying, the weights were selected to be 1.9±0.2 g (BB-0), 1.5±0.2 g (BB-3), 1.1±0.2 g (BB-6), 0.7±0.2 g (BB-9) and less than 0.5 g (BB-12), respectively. Even if the drying time is the same, the MC may not be constant. The groups were classified by weight to deliberately vary the MC by group. Hyperspectral images were acquired for each group. Altogether, 670 blueberries (134 samples per group) were used for the In partial least-square discriminant analysis (PLS-DA) to determine and predict the degree of drying based on the MC. In partial least-squares regression (PLSR) for predicting MC, 150 (30 samples per group) of the blueberries purchased from the local market were used to develop an MC prediction model.

Determination of quality properties

The MC, pH, hardness and strength, total soluble solids content, color, TPC content, TFC and microbiological content of the blueberries were analyzed. The MC (%) of the dried blueberries was calculated using the following formula:

(1)

where Wc and Wd denote the fresh and dry mass weights of the blueberries, respectively. The TSS content and pH were measured using a refractometer (Master-α, Atago Co., Tokyo, Japan) and a pH meter (Orion 3-Star, Thermo Electron Co. Waltham, MA, USA), respectively, after dilution (1:4, v/v) with distilled water and grinding using a homogenizer (AM-9, Nihonseiki Kashima Co., Ltd., Japan).

The hardness and strength of the blueberries were measured using a rheometer (Compac-100II, SunScientific Co., Tokyo, Japan) by compressing 20% of the sample using a 3- mm-diameter probe.

The color values of blueberries were measured using a colorimeter (CR-400, Minolta Co., Osaka, Japan) calibrated against a standard white plate (L*=97.79, a*=−0.38, b*= 2.05). The L*a*b* color values refer to brightness, redness, and yellowness, respectively.

To analyze the TPC and TFC of blueberries, a sliced sample (4 g) was extracted with 76 mL of 50% EtOH in a shaker at 45°C for 24 h. The extracts were then filtered (No. 4 filter paper; Whatman, Tokyo, Japan). The filtrates were evaporated on a rotary evaporator (RV 10 Digital, IKA, Staufen, Germany) to determine the dry weight, which was then diluted with 50% EtOH. TPC content was analyzed using a modification of the method described by Singleton et al. (1999). All experiments were performed in triplicates. Standard curves of gallic acid and naringin were constructed to calculate the TPC content and TFC, respectively.

Microbiological analysis of the dried blueberries was carried out by adding nine-fold the amount of peptone to 3.5 g of the sample and shaking at a constant rate for 2 min. The extracts were diluted using the 10-fold method, and 1 mL of each dilution was added to a Petri dish. Plate count agar and potato dextrose agar were poured and cultured at 35±1 and 25±1°C, respectively. The counts are expressed as log CFU/g.

Duncan’s multiple range test and Pearson’s correlation analysis were performed using SPSS (14.0, SPSS, Inc., Chicago, IL, USA) to identify significant differences between the results (p<0.05).

Hyperspectral data acquisition

Samples were placed on a white background plate and line-scanned using a Vis/NIR HSI system (IMEC Kapeldreef 75, 3001 Heverlee, Belgium) with a spectral range of 470- 900 nm and a resolution of 2,048×1,088. The number of spectral bands was 150, and the pixel size was 5.5 μm×5. The spectrometer was installed with six halogen lamps of 20 W each. The frame rate during spectral acquisition was 10, and the exposure time during the measurement was 10 s. All raw images were calibrated using white and dark images as follows:

${I}_{cal}=\frac{{I}_{raw}-{I}_{dark}}{{I}_{white}-{I}_{dark}}$
(2)

where Ical is a calibrated image, Iraw is a raw image, Idark is a dark image, and Iwhite is a white image. The Idark value is the spectrum obtained after closing the camera lens, and Iwhite is the spectrum of the white plate. To acquire sample spectra with removed background, the region of interest (ROI) function was applied in the ENVI software (version 4.8, Harris Co., Melbourne, FL, USA). The extraction of ROIs depended on human observation.

Spectral pre-processing

The extracted data should be used for analysis after removal of scattered light, unnecessary light sources in the surrounding environment, changes in the condition of the apparatus, and noise caused by the introduction of the sample. Bae et al. (2016) noted that it is possible to obtain more accurate spectra by performing various data preprocessing techniques such as normalization, smoothing, multiplicative scatter correction (MSC), standard normal variate (SNV) analysis, and the Savitzky-Golay (SG) method. Normalization makes the data distribution similar and corrects for the change of conditions. Smoothing replaced a certain range of values with the mean value in order to represent a spectrum with noise as a smooth curve. SNV and MSC pre-processing removed the influence of light scattering due to irregular shapes and sizes. SG-1 and SG-2 differentiation methods eliminated the movement of the baseline of the spectrum or emphasized the spectral characteristics.

Partial least-squared analysis (PLSA)

PLSR, a common chemometric method, is a popular multivariate calibration technique for quantitative analysis (Burnham et al., 1999). PLS-DA is a PLS method, which is a regression model analysis that identifies groups by setting virtual variables. The PLS regression method not only obtains a regression equation using the principal components of the spectral data, but also predicts new variables by considering the measured dependent variables. Because both independent and dependent variables are used, discrimination is performed more stably than analysis of principal components using only independent variables. The β-coefficient calculated in the PLS-DA model was used to predict the wavelength affecting the determination of the dryness of blueberries, and the PLS-DA was executed using Unscrambler (10.5, CAMO, Trondheim, Norway). As explained in section 2.2., in this study, the MC of blueberries was classified into five levels, and the five levels were confirmed by PLS-DA. Approximately 70% of the 670 blueberries were used to develop the calibration model, and 200 were applied to the developed model to confirm the discrimination accuracy.

PLSR was employed to predict the MC of dried blueberries. A regression model was derived using the HSI spectrum and measured MC of each blueberry. The PLSR model was calibrated using data from 70% of the 150 samples. The remaining 30% of the samples were used to validate the model. The accuracy of the developed model was estimated using the coefficients of calibration (R2c), full cross-validation (R2cv), root mean square error of calibration (RMSEC), and root mean square error of full cross-validation (RMSECV). The accuracy of the prediction model is represented by R2p. All model development procedures were conducted using Unscrambler X software (10.2, Camo, Oslo, Norway).

## Results and discussion

Physicochemical properties

The physicochemical data for the blueberries are presented in Table 1. From the analysis of the quality characteristics, MC was found to decrease significantly after the drying period (p<0.05), and TSS content did not show any significant change until the 6th day of storage, but displayed an increasing trend starting from the 9th day of storage. The pH was not affected by the drying time, whereas TPC and TFC showed a tendency to decrease significantly. Lightness (L*) decreased with time, whereas yellowness (b*) increased. As it dries, the blueberry becomes sticky and darker in color as the moisture and sugar inside the pulp move to the pericarp. Fruit firmness is one of the most critical quality indices of blueberries (Angeletti et al., 2010; Chen et al., 2015; Chen et al., 2017). Both the hardness and strength rapidly decreased until the 6th day of drying. This is because blueberries initially lose their hardness and strength as the texture of the flesh degrades. Softening is one of the main reasons for the high perishability and susceptible spoilage of blueberries during postharvest storage and influences fruit quality, storage life, and transport ability (Chen et al., 2015; Deng et al., 2014; Paniagua et al., 2013). Blueberry’s softening is primarily caused by the degradation of cell wall components, including water-soluble pectin, cellulose, and hemicellulose. During storage of blueberries, the water-soluble pectin content increases, and the contents of cellulose and hemicellulose decrease. (Liu et al., 2019). However, they became harder as the MC decreased rapidly from the 9th day of drying. This is consistent with the study of Kang et al. (2004), who reported that the hardness decreased and then increased again due to excessive evaporation of moisture in dried persimmons during the drying process.

Table 1. Physicochemical properties of blueberries according to drying period
Physicochemical property1) Drying period (day)
0 3 6 9 12
MC (%) 86.32±0.94e2) 82.11±1.67d 74.96±1.17c 64.80±2.08b 36.34±3.12a
TSS (Brix) 4.47±1.14ab 3.17±0.29ab 3.13±0.12ab 5.90±0.56b 9.10±1.82c
pH 4.38±0.23a 4.17±0.27a 4.21±0.07a 4.27±0.27a 4.28±0.24a
TPC (mg/g dw) 13.68±0.69d 11.88±0.33c 12.40±0.26c 10.63±0.09b 9.15±0.14a
TFC (mg/g dw) 7.52±0.19d 6.79±0.07c 6.82±0.24c 5.40±0.31b 3.75±0.10a
Color value3) L* 31.19±1.60d 30.05±1.46c 27.22±1.45b 27.16±1.21b 18.42±2.21a
a* 0.49±0.21b 0.36±0.13a 0.57±0.17b 0.48±0.22b 0.46±0.15ab
b* -4.42±1.04a -3.73±1.10b -2.69±0.97c -2.83±0.79c -1.76±0.44d
Texture Hardness 9,951.24±1,826.08c 3,822.49±756.14b 1,365.22±193.65a 1,700.13±436.72a 3,558.31±1,296.06b
Strength 3,035.26±361.57d 1,500.88±329.31c 193.82±45.92a 299.92±90.60a 787.28±337.88b
Bacteria ND4) ND ND ND 0.60±0.93
Mould and yeast 2.25±1.20 2.76±0.51 2.85±0.83 4.01±0.12 2.55±0.32

MC, moisture content; TSS, total soluble solids; TPC, total phenolic compound; TFC, total flavonoid content.

Mean±SD with different letters are significantly different at 5% level.

L*=lightness; a*=redness; b*=yellowness.

Not detect.

Microbiological analysis showed that bacteria were rarely produced throughout the experimental period, and the number of molds and yeasts was similar, except for the 9th day of drying. Jung et al. (2014) reported that total aerobic bacteria, yeast, and molds occurred at approximately 4 log CFU/g in the control group when stored for 12 days under refrigeration. This result is different from the results of the current study. In this study, it is considered that the bacterial occurrence was reduced because the MC and water activity both decreased as the blueberry dried.

Blueberries contain sucrose, glucose, and fructose, with sucrose as the main sugar (Jain et al., 2011; Rai et al., 2007). The TSS content of the blueberries varied between 55° and 65° Brix. Since blueberry fruit does not contain starch to support soluble sugar synthesis after harvest, the slight increase in TSS may be a consequence of cell wall degradation (Cordenunsi et al., 2003).

Classification results of PLS-DA full-wavelengths in accordance with various pre-processing

The classification results of the PLS-DA model pre- treated with various methods (smoothing, normalization, SNV, multiplicative scatter correction [MSC], SG-1 and SG-2) are shown in Table 2. In this study, the R2c and R2cv of the developed calibration model and cross-validation model using raw data without pre-treatment were higher than 0.91, and R2p determined the drying extent at a high accuracy rate of 0.81. Among the pre-processing methods used in this study, the accuracy of the predictive model was approximately 70% when applying the MSC and SNV preprocessing methods, which yielded the most inaccurate results. However, when the data were pre-processed using the smoothing moving average method, an R2p of 0.81 was shown. In the PLS-DA model for MC determination, smoothing pre-treatment of the spectrum was the most appropriate.

Table 2. Calibration and cross-validation values in accordance with spectra preprocessing using the partial least squares regression (PLSR) model
Pretreatment1) Calibration (n=470) Full cross-validation (n=470) Total Prediction (n=200)
Slope Offsetc RMSEC2) R2c Slope Offsetcv RMSECV R2cv R2
Raw 0.9147 0.2558 0.4129 0.9147 0.9134 0.2598 0.4192 0.9124 0.9136 80.5%
Normalisation 0.9140 0.2579 0.4146 0.9140 0.9130 0.2616 0.4207 0.9119 0.9130 75.5%
Smoothing 0.9142 0.2573 0.4142 0.9142 0.9129 0.2612 0.4202 0.9120 0.9131 81.0%
MSC 0.8839 0.3483 0.4819 0.8839 0.8826 0.3520 0.4862 0.8823 0.8831 69.5%
SNV 0.8845 0.3464 0.4805 0.8845 0.8833 0.3502 0.4848 0.8830 0.8838 70.0%
SG-1 0.9208 0.2375 0.3979 0.9208 0.9186 0.2438 0.3068 0.9176 0.9192 76.0%
SG-2 0.9082 0.2751 0.4283 0.9083 0.9043 0.2869 0.4386 0.9041 0.9062 72.0%

SNV, standard normal variate; MSC, multiple scatter correction; SG-1, Savitzky-Golay-1; SG-2, Savizky-Goloy-2.

RMSEC, root mean square error of calibration; RMSECV, root mean square error of cross-validation; R2c, the coefficient of determination in calibration; R2cv, the coefficient of determination in cross-validation;

The average value of R2c and R2cv.

Spectral interpretation

The original images obtained using the hyperspectral camera are shown in Fig. 1. The mean reflectance pre- treated smoothing of blueberries according to the drying period, as determined by HSI analysis, is shown in Fig. 2. Overall, the longer was the blueberries dried, the lower the reflectance. The β coefficient of the PLS-DA model can be used for meaningful wavelength selection (Cheng and Sun, 2015a). Fig. 2 shows the absolute value of the β coefficient of the smoothed spectra. There were five specific peaks at approximately 706, 790, 827, 868, and 894 nm. Using the β coefficient of the PLS-DA model, the informative wavelength was determined, and the quality characteristics affecting the absorption spectrum of the blueberry were measured. For example, as the primary structural components of the majority of food molecules contain functional groups such, as C-H (organic compounds), N[1]-H (proteins and amino acids) and O-H (water, carbohydrate and fat), these groups are closely related to the overtones and combination bands in the NIR spectral region (Pu et al., 2015).

Fig. 1. Hyperspectral image of blueberry acquired during drying process. BB-0, blueberries dried for 0 day (1.9±0.2 g); BB-3, blueberries dried for 3 days (1.5±0.2 g); BB-6, blueberries dried for 6 days (1.1±0.2 g); BB-9, blueberries dried for 9 days; BB-12, blueberries dried for 12 days.
Fig. 2. Smoothing pretreated hyperspectral mean spectra and beta coefficient of partial least-squares discriminant analysis model of blueberries with different degrees of drying. BB-0, blueberries dried for 0 day (1.9±0.2 g); BB-3, blueberries dried for 3 days (1.5±0.2 g); BB-6, blueberries dried for 6 days (1.1±0.2 g); BB-9, blueberries dried for 9 days; BB-12, blueberries dried for 12 days.

In the current study, the spectral region below approximately 690 nm did not differ between the groups (Fig. 1). The absorption bands at approximately 780 nm were attributed to the third overtone of the O-H stretching generated by water (Cheng and Sun, 2015b). He et al. (2013) also reported that the absorption band positioned at approximately 780 nm in the NIR region was due to the third overtone of O-H stretching.

Therefore, in Fig. 2, the peak around 790 nm is caused by the MC. A peak at approximately 850 nm is usually assigned to the stretching overtones of C-H and N-H, which are associated with proteins, lipids, and other organic compounds (Howard and Kjaergaard, 2006; Tarr and Zerbetto 1989). The TSS content represents the organic molecules, which contain C-H, O-H, C-O, and C-C bonds (Liu et al., 2010). Glucose has an informative region of 840-1,062 nm, where bands resulting from the first overtones of the C-H stretching modes of glucose exist (Workman and Weyer, 2007). Informative regions of 816-1,050 nm were obtained for fructose. These regions correspond to the regions for bands due to the combination of C-H stretching (Osborne et al., 1993). Sugar absorption bands, generally from O-H stretching, and the third and fourth overtones of C-H stretching, occur near strong water absorption regions (Giangiacomo, 2006). This suggests that 827, 868 and 894 nm, shown in Fig. 2, are the peaks caused by the TSS difference between the groups.

A substance substituted with a period (-OH), such as that of one of the hydrogens of the benzene ring (C6H6), is referred to as a phenol. Polyphenols are defined as having two or more hydroxyl groups. Examples include catechins, tannins, anthocyanins, and flavonoids. In this study, the TPC and TFC of blueberries were significantly altered during drying. Therefore, it is considered that this affects the hyperspectral image spectra.

Pearson’s correlation

Pearson’s correlation coefficient (r) is commonly used to determine the relationship between two variables (Cha and Kim, 2018). Pearson’s correlations among the physicochemical analysis results are shown in Table 3. Overall, MC and various quality characteristics showed high correlations. Those between MC and lightness, TPC, TFC, and TSS were 0.9 or more, and high correlations between hardness and strength (0.979) and TPC and TFC (0.977) were also determined. Lee et al. (2011) reported that differences in TSS might be related to differences in MC during the drying process, and Shin et al. (2015) observed that the difference in texture of dried blueberries was correlated with the difference in MC during the drying process, which is similar to the results of this study. Pearson’s correlation coefficients can induce quality characteristics that indirectly affect the β coefficient. Therefore, it is expected that the MC of blueberries can be predicted by the above correlation from the tendencies in attributes, such as lightness, TSS, TFC and TPC. The wavelength range at which the absolute value of the β coefficient was high was influenced by the MC, TSS, TFC and TPC content. It can be seen that these wavelength bands are highly correlated with each other.

Table 3. Pearson’s correlation between physicochemical properties
Physicochemical property MC TSS pH L* a* b* Hardness Strength TPC TFC Bacteria Mold and yeast
MC 1 -0.919 -0.015 0.984 -0.036 -0.903 0.385 0.494 0.931 0.982 - -0.126
TSS -0.919 1 0.359 -0.867 -0.04 0.672 -0.042 -0.168 -0.814 -0.912 - 0.043
pH -0.015 0.359 1 0.031 0.299 -0.285 0.767 0.633 0.23 0.039 - -0.275
L* 0.984 -0.867 0.031 1 -0.093 -0.919 0.384 0.493 0.882 0.938 - -0.002
a* -0.036 -0.04 0.299 -0.093 1 0.237 -0.116 -0.0262 0.204 0.092 - 0.059
b* -0.903 0.672 -0.285 -0.919 0.237 1 -0.700 -0.793 -0.868 -0.872 - 0.279
Hardness 0.385 -0.042 0.767 0.384 -0.116 -0.700 1 0.979 0.563 0.43 - -0.656
Strength 0.494 -0.168 0.633 0.493 -0.262 -0.793 0.979 1 0.615 0.519 - -0.654
TPC 0.931 -0.814 0.23 0.882 0.204 -0.868 0.563 0.615 1 0.977 - -0.359
TFC 0.982 -0.912 0.039 0.938 0.092 -0.872 0.43 0.519 0.977 1 - -0.254
Bacteria - - - - - - - - - - 1 -
Mold and yeast -0.126 0.043 -0.275 -0.002 0.059 0.279 -0.656 -0.654 -0.359 -0.254 - 1

MC, moisture contents; TSS, total soluble solids; TPC, total phenolic contents; TFC, total flavonoid contents.

Prediction of MC using the PLS model

The PLSR model was developed using HSI spectra and the respective measured MC values. In total, 105 of the 150 samples were used to create the calibration model. When the MC prediction model was established based on the remaining 30% of the samples that were not used in the model calibration, most of the preprocessing methods showed high accuracy. A prediction model for predicting MC was implemented using spectra with different pre- treatments. The measured MC and predicted moisture are shown in Fig. 3. When MSC and SNV pre-processing were applied, the R2 value was higher than 0.91, indicating a higher prediction accuracy than with the raw data. This result differs from that of the PLS-DA model. Finally, we developed an MC prediction model with high prediction accuracy. These results suggest that HSI techniques can be used to monitor and predict the MC of blueberries during drying. In the future, it is necessary to develop a machine learning system to implement a machine learning system in actural food production, distribution, and storage sites by conducting field demonstration research in which an algorithm using hyperspectral imaging anlaysis is applied to a machine.

Fig. 3. Predicting accuracy of blueberry moisture content of partial least-squares model using various preprocessed hyperspectral spectra. MC, moisture contents; SNV, standard normal variate; MSC, multiplicative scatter correction; SG, Savizky Golay.

## Conflict of interests

The authors declare no potential conflict of interest.

## References

1.

Angeletti P, Castagnasso H, Miceli E, Terminiello L, Concellon A, Chaves A, Vicente AR. Effect of preharvest calcium applications on postharvest quality, softening and cell wall degradation of two blueberry (Vacciniuum corymbosum) varieties. Postharvest Biol Tec, 58, 98-103 (2010)

2.

Bae H, Seo YW, Kim DY, Lohumi S, Park ES, Cho BK. Development of non-destructive sorting techniques for viability of watermelon seed using hyperspectral image processing. J Korean Soc Nondestruc Test, 36, 35-44 (2016)

3.

Burnham AJ, Macgregor JF, Viveros R. Latent variable multivariate regression modeling. Chemom Intell Lab Syst, 48, 167-180 (1999)

4.

Cha J, Kim J. Analysis of fine dust correlation between air quality and meteorological factors using SPSS. J Korea Inst Inf Commun Eng, 22, 722-727 (2018)

5.

Chen H, Cao S, Fang X, Mu H, Yang H, Wang X, Xu Q, Gao H. Changes in fruit firmness, cell wall composition and cell wall degrading enzymes in postharvest blueberries during storage. Sci Hortic, 188, 44-48 (2015)

6.

Chen Y, Hung YC, Chen M, Lin H. Effects of acidic electrolyzed oxidizing water on retarding cell wall degradation and delaying softening of blueberries during postharvest storage. LWT-Food Sci Technol, 84, 650-657 (2017)

7.

Cheng JH, Sun DW. Rapid and non-invasive detection of fish microbial spoilage by visible and near infrared hyperspectral imaging and multivariate analysis. LWT- Food Sci Technol, 62, 1060-1068 (2015b)

8.

Cheng JH, Sun DW. Rapid quantification analysis and visualization of Escherichia coli loads in grass crap fish flesh by hyperspectral imaging method. Food Bioproc Technol, 8, 951-959 (2015a)

9.

Cordenunsi BR, Nascimento JRO, Lajolo FM. Physico- chemical changes related to quality of five strawberry fruit cultivars during cool-storage. Food Chem, 83, 167-173 (2003)

10.

Deng J, Shi ZJ, Li XZ, Liu HM. Effects of cold storage and 1-methycyclopropene treatments on ripening and cell wall degrading in rebbiteye blueberry (Vaccinium ashei) fruit. Food Sci Technol Int, 75, 1-12 (2014)

11.

Fan SX, Li CY, Huang WQ, Chen LP. Detection of blueberry internal bruising over time using NIR hyperspectral reflectance imaging with optimum wavelengths. Postharvest Biol Tec, 134, 55-66 (2017)

12.

Giangiacomo R. Study of water-sugar interactions at increasing sugar concentration by NIR spectroscopy. Food Chem, 96, 371-379 (2006)

13.

Gowen AA, O’Donnell CP, Cullen PJ, Downey G, Frias JM. Hyperspectral imaging-an emerging process analytical tool for food quality and safety control. Trends Food Sci Tech, 18, 590-598 (2007)

14.

He HJ, Wu D, Sun DW. Non-destructive and rapid analysis of moisture distribution in farmed Atlantic salmon (Salmo salar) fillets using visible and near-infrared hyperspectral imaging. Innovative Food Sci Emerging Technol, 18, 237-245 (2013)

15.

Howard DL, Kjaergaard HG. Influence of intramolecular hydrogen bond strength on OH-stretching overtones. J Phys Chem A, 110, 10245-10250 (2006)

16.

Huang M, Wang Q, Zhang M, Zhu Q. Prediction of color and moisture content for vegetable soybean during drying using hyperspectral imaging technology. J Food Eng, 128, 24-30 (2014)

17.

Huang Y, Min S, Duan J, Wu L, Li Q. Identification of additive components in powdered milk by NIR imaging methods. Food Chem, 145, 278-283 (2014)

18.

Hwang SH, Ko SH. Quality characteristics of muffins containing domestic blueberry (V. corymbosum). J East Asian Soc Diet Life, 20, 727-734 (2010)

19.

Iqbal A, Sun DW, Allen P. An overview on principle, techniques and application of hyperspectral imaging with special reference to ham quality evaluation and control. Food Control, 46, 242-254 (2014)

20.

Jain SK, Verma RC, Murdia LK, Jain HK, Sharmak GP. Optimization of process parameters for osmotic dehydration of papaya cubes. J Food Sci Technol, 48, 211-217 (2011)

21.

Jung SH, Kang JH, Park SJ, Seong KH, Song KB. Quality changes in ‘Elliot’ blueberries and ‘Sulhyang’ strawberries packed with two different packaging materials during refrigerated storage. J Korean Soc Food Sci Nutr, 43, 901-908 (2014)

22.

Kang WW, Kim JK, Oh SL, Kim JH, Han JH, Yang JM, Choi JU. Physicochemical characteristics of Sangju traditional dried persimmons during drying process. J Korean Soc Food Sci Nutr, 33, 386-391 (2004)

23.

Lee HR, Jung BR, Park JY, Hwang IW, Kim SK, Choi JU, Lee SH, Chung SK. Antioxidant activity and total phenolic contents of grape juice products in the Korean market. Korean J Food Preserv, 15, 445-449 (2008)

24.

Lee JG, Lee BY. Effect of media composition on growth and rooting of highbush blueberry cuttings. J Korean Soc Hortic Sci, 25, 355-359 (2007)

25.

Lee SW, Moon HK, Lee WY, Kim JK. Physicochemical characteristics of cold-air dried persimmons and traditional dried persimmons. Korean J Food Preserv, 18, 481-487 (2011)

26.

Leiva-Valenzuela GA, Lu R, Aguilera JM. Assessment of internal quality of blueberries using hyperspectral transmittance and reflectance images with whole spectra or selected wavelengths. Innov Food Sci Emerg, 24, 2-13 (2014)

27.

Lin X, Xu JL, Sun DW. Investigation of moisture content uniformity of microwave-vacuum dried mushroom (Agaricus bisporus) by NIR hyperspectral imaging. LWT-Food Sci Technol, 109, 108-117 (2019)

28.

Liu B, Wang K, Shu X, Liang J, Fan X, Sun L. Changes in fruit firmness, quality traits and cell wall constituents of two highbush blueberries (Vaccinium corymbosum L.) during postharvest cold storage. Sci Hortic, 246, 557-562 (2019)

29.

Liu Y, Sun X, Ouyang A. Nondestructive measurement of soluble solid content of navel orange fruit by visible- NIR spectrometric technique with PLSR and PCA- BPNN. LWT-Food Sci Technol, 43, 602-607 (2010)

30.

Nguyen-Do-Trong N, Dusabumuremyi JC, Saeys W. Cross-polarized VNIR hyperspectral reflectance imaging for non-destructive quality evaluation of dried banana slices, drying process monitoring and control. J Food Eng, 238, 85-94 (2018)

31.

Osborne BG, Fearn T, Hindle PH. Practical NIR Spectroscopy: With Applications in Food and Beverage Analysis. Longman Scientific and Technical, Harlow, UK, p. 227 (1993)

32.

Paniagua AC, East AR, Heyes JA. Interaction of temperature control deficiencies and atmosphere conditions during blueberry storage on quality outcomes. Postharvest Biol Technol, 95, 50-59 (2014)

33.

Park HM, Yang SJ, Kang EJ, Lee DH, Kim DI, Hong JH. Quality characteristics and granule manufacture of mulberry and blueberry fruit extracts. Korean J Food Cook Sci, 28, 375-381 (2012)

34.

Park SJ, Choi YB, Ko JR, Rha YA, Lee HY. Effects of drying methods on the quality and physiological activities of blueberry (Vaccinium ashei). Culi Sci Hos Res, 20, 55-64 (2014)

35.

Pu H, Kamruzzaman M, Sun DW. Selection of feature wavelengths for developing multispectral imaging systems for quality, safety and authenticity of muscle foods-a review. Trends Food Sci Tech, 45, 86-104 (2015)

36.

Pu YY, Feng YZ, Sun DW. Recent progress of hyperspectral imaging on quality and safety inspection of fruits and vegetables: A review. Compr Rev Food Sci F, 14, 176-188 (2015)

37.

Pu YY, Sun DW. Vis-NIR hyperspectral imaging in visualizing moisture distribution of mango slices during microwave-vacuum drying. Food Chem, 188, 271-278 (2015)

38.

Qin JW, Chao KL, Kim MS, Lu RF, Burks TF. Hyperspectral and multispectral imaging for evaluating food safety and quality. J Food Eng, 118, 157-171 (2013)

39.

Rai S, Pal RK, Jayachandran KS. Optimization of process parameters for osmotic dehydration of pineapple slices. Indian J Hortic, 64, 304-308 (2007)

40.

Shin DS, Yoo YM, Kim HY, Han GJ. Determine the effects of drying temperature on the quality change and antioxidant activity characteristics of blueberry. Korean J Food Preserv, 22, 505-511 (2015)

41.

Singleton VL, Orthofer R, Lamuela-Raventos RM. Analysis of total phenols and other oxidation substrates and antioxidants by means of Folin-Ciocalteu reagent. Methods Enzymol, 299, 152-178 (1999)

42.

Sugiyama T, Sugiyama J, Tsuta M, Fujita K, Shibata M, Kokawa M, Araki T, Nabetani H, Sagara Y. NIR spectral imaging with discriminant analysis for detecting foreign materials among blueberries. J Food Eng, 101, 244-252 (2010)

43.

Tarr AW, Zerbetto F. Absolute intensities of CH-stretching overtone in chloroform and deuterochloroform. Chem Phys Lett, 154, 273-279 (1989)

44.

Wei Y, Wu F, Xu J, Sha J, Li X. Visual Detection of the moisture content of tea leaves with hyperspectral imaging technology. J Food Eng, 248, 89-96 (2019)

45.

Westwood MN. Temperate-zone Pomology: Physiology and Culture. 3th ed, Timber Press, Portland, OR, USA p 523 (1993)

46.

Workman J, Weyer L. Practical Guide to Interpretative Near-Infrared Spectroscopy. CRC Press, Boca Raton, FL, USA (2007)

47.

Xu R, Takeda F, Krewer G, Li C. Measure of mechanical impacts in commercial blueberry packing lines and potential damage to blueberry fruit. Postharvest Biol, Technol, 110, 103-113 (2015)

48.

Xu Y, Zhang H, Zhang C, Wu P, Li J, Xia Y, Fan S. Rapid prediction and visualization of moisture content in single cucumber (Cucumis sativus L.) seed using hyperspectral imaging technology. Infrared Phys Technol, 102, 103034 (2019)

49.

Yang C, Lee WS, Gader P. Hyperspectral band selection for detecting different blueberry fruit maturity stages. Comput Electron Agric, 109, 23-31 (2014)