Journal of Dairy Science
Volume 90, Issue 7 , Pages 3118-3125, July 2007

Estimation of Storage Time of Yogurt with Artificial Neural Network Modeling

Suleyman Demirel University, Food Engineering Department, 32200, Isparta, Turkey

Received 13 September 2006; accepted 23 February 2007.

Article Outline

Abstract 

Changes in the physical, chemical, and microbiological structure of yogurt determine the storage and shelf life of the product. In this study, microbial counts and pH values of yogurt during storage were determined at d 1, 7, and 14. Simultaneously, image processing of yogurt was digitized by using a machine vision system (MVS) to determine color changes during storage, and the obtained data were modeled with an artificial neural network (ANN) for prediction of shelf life of set-type whole-fat and low-fat yogurts. The ANN models were developed using back-propagation networks with a single hidden layer and sigmoid activation functions. The input variables of the network were pH; total aerobic, yeast, mold, and coliform counts; and color analysis values measured by the machine vision system. The output variable was the storage time of the yogurt. The modeling results showed that there was excellent agreement between the experimental data and predicted values, with a high determination coefficient (R2 = 0.9996) showing that the developed model was able to analyze nonlinear multivariant data with very good performance, fewer parameters, and shorter calculation time. The model might be an alternative method to control the expiration date of yogurt shown in labeling and provide consumers with a safer food supply.

Key words: artificial neural network, yogurt, machine vision, prediction of shelf life

 

Back to Article Outline

Introduction 

Yogurt is a semisolid fermented milk product consumed in most parts of the world. Changes in the physical, chemical, and microbiological structure of yogurt determine the storage and shelf life of the product. Alteration of these properties causes color, aroma, and texture deterioration of yogurt, which are considered important quality criteria by consumers.

Artificial neural networks (ANN) are new information processing techniques offering solutions to problems that have not been clearly formulated. Much of the excitement surrounding neural networks is their unique ability to learn by experience. In the past few years, neural networks have shown increased power over many other statistical methods when solving nonlinear prediction problems (Bochereau et al., 1992). Artificial neural networks are general nonlinear models based on a simplification of human brain function; more than other modeling strategies, they have the capability to internally self-adapt and encapsulate complex nonlinear relationships between input and output variables without the need for an a priori rigid model structure (Hornik et al., 1989; Thibault and Grandjean, 1992). As in nature, the network function is determined largely by the connections between elements. A neural network can be trained to perform a particular function by adjusting the values of the connections (weights) between the elements. Commonly, neural networks are adjusted or trained so that a particular input leads to a specific target output.

There have been many applications of ANN reported for the interpretation of images in the agri-food industry (Bishop, 1994). Fuzzy logic, genetic algorithms, and similar regression analyses have been used for prediction of shelf life of food, but generally the best results have been obtained by using ANN (Doganis et al., 2006). Artificial neural networks have been used, successfully, as a modeling tool in several food processing applications such as sensory analysis and quality control (color analysis, textural evaluation, human preferences, and so on), classifications, microbiology, and drying applications (Ni and Gunasekaran, 1998; Edwards and Cobb, 1999; Farkas et al., 2000; Hussian et al., 2002).

Increasing demand for fast, reliable, and objective techniques to determine shelf life of food products with high quality and safety standards encourages scientists to use image analysis techniques for food quality evaluation (Park and Chen, 2000; Kavdir and Guyer, 2002; Brosnan and Sun, 2004; Sun and Du, 2004). Ni and Gunasekaran (1998) and Brosnan and Sun (2004) reviewed the use of computer vision technology for food analysis. In image analysis for food products, color is an influential attribute of visual information and a powerful descriptor of measurement. Research shows that color has been successful in classifying a variety of food products (Brosnan and Sun, 2004; Sun, 2004). Besides the use of machine vision systems for inspection and grading of fruits and vegetables (Brosnan and Sun, 2004), it has been used routinely in the quality assessment of meat, cheese, and pizza (McDonald and Chen, 1990; Gerrard et al., 1996; Jamieson, 2002). Analysis of the characteristics of Cheddar and Mozzarella cheeses with machine vision during cooking gave promising results, suggesting that the method provided an objective and easy approach for analyzing the functional properties of cheese (Wang and Sun, 2002). Ni and Gunasekaran (1998) developed an image-processing algorithm to recognize individual cheese shreds and automatically measure the shred length.

The objective of this study was to determine color changes during storage of set-type whole-fat and low-fat yogurts using a machine vision system and to develop a prediction model using ANN for prediction of shelf life of yogurt.

Back to Article Outline

Materials and Methods 

Raw Material 

Commercial yogurt samples were obtained from a local market in Isparta (southwestern Turkey). The samples included plain, set-type whole-fat (3% fat) and low-fat (1.5% fat) yogurt and were stored in the refrigerator (4 to 7°C) for 14 d. Analyses were performed during storage on d 1, 7, and 14.

Chemical Analyses 

Measurements of pH values of the yogurt samples were performed at 17 to 20°C using a pH meter (Inolab WTW, Measurement System, Fort Myers, FL). Fat content was measured by the Gerber method (British Standards Institution, 1955).

Microbiological Analyses 

Total aerobic counts (TAC) were determined on plate count agar (Merck, Darmstadt, Germany) incubated at 32°C for 3 d. Coliform, yeast, and mold counts of yogurts were determined using the pour-plate method. Coliform counts (CC) were estimated using violet red bile agar (Merck) incubated at 37°C for 48h. Yeast and mold counts (YMC) were carried out on potato dextrose agar (Merck) incubated at 25°C for 5 d.

Machine Vision System 

Hardware 

The machine vision system described by Demir et al. (2002) and Balaban et al. (1997) was used as an alternative for color measurements of yogurt samples. Images were digitized using pixels with information on the levels of 3 primary colors (RGB: red, green, and blue) for an image. This system can give a complete description and amount of all the colors present in the samples with different sizes, shapes, surface textures, and colors. The hardware consisted of a light box, digital color video camera, and computer. The light box was used to standardize illumination and comprised an aluminum framework covered by 100% white safety glazing fabrics. The light box had top and bottom lighting with 2 fluorescent lights (illuminant D50, simulating noon summer sun; General Electric, Cleveland, OH) inside the box. Three yogurt samples were placed in the light box with the top lighting, and the door of the light box was closed. Yogurt samples were observed on d 1, 7, and 14 with a charge-coupled device (CCD) color video camera (Sony DFW-V500, Sony Corp., Tokyo, Japan) with 64-bit color, 525 lines, and a horizontal resolution of >460 lines. During image acquisition a standardized color tile was placed next to the yogurt samples in the light box. Images of samples were captured and saved for further color analysis (Luzuriaga et al., 1997). A schematic presentation of the machine vision system is shown in Figure 1.

  • View full-size image.
  • Figure 1. 

    Machine vision system used for color analysis of yogurt. MVS = machine vision system; CAS = color analysis software; CCD charge-coupled device camera; ANN = artificial neural network.

Software 

A color analysis program called Color Expert, running under the Microsoft Windows XP environment, was used to extract color information from a color image. The user interface program was written in Microsoft Visual Basic professional (version 6.0, Microsoft Inc., Redmond, WA), and also uses the Pinnacle Graphics Server (version 5.0, Bits Per Second Ltd., Brighton, UK), and the Formula One spreadsheet from Visual Tools (version 4.1, Visual Components, Lenexa, KS; Luzuriaga et al., 1997). Color data from images of yogurt samples were analyzed by looking at every pixel. The system included a database that stored reference values to correctly evaluate yogurt quality, and saved all results from the analysis in one spreadsheet.

ANN Modeling 

The network was adjusted based on a comparison of the output and the target, until the network output matched the target. Typically many such input/target output pairs are used to train a network (Demuth and Beale, 2001). Batch training of a network proceeds by making weight and bias changes based on an entire set (batch) of input vectors. Incremental training changes the weight and biases of a network as needed after presentation of each individual input vector. Incremental training is sometimes referred to as online or adaptive training.

In this study, Levenberg-Marquardt back-propagation (BP) network training was applied until the evaluation standard was reached (Demuth and Beale, 2001). These networks are most widely used and are considered the workhorses of ANN (Rumelhart et al., 1986). A BP network is a multilayer perceptron consisting of 1) an input layer with nodes representing input variables, 2) an output layer with nodes representing the dependent variables (i.e., what is being modeled), and 3) one or more hidden layers containing nodes to help capture the nonlinearity in the data. Using supervised learning (with the error correction learning rule), these networks can learn the mapping from one data space to another using examples. The term “back-propagation” refers to the way the error computed at the output side is propagated backward from the output layer, to the hidden layer, and finally to the input layer (Hassoun, 1995).

A neural network model consists of an input layer, an output layer, and one or more hidden layers. In this study, the ANN was designed with only one hidden layer (Figure 2). One hidden layer is usually sufficient to approximate any continuous nonlinear function, although more complex networks must be used in special applications (Bucinski et al., 2004).

  • View full-size image.
  • Figure 2. 

    Detail of the function of a neuron. The inputs (y) into a neuron are multiplied by their corresponding connection weights (w) and summed. This sum is then transformed through the sigmoid function to produce a single output that may be passed on to other neurons.

The function of the neuron is schematically shown in Figure 2. The inputs (yi) into a neuron are multiplied by their corresponding connection weights (wij) and summed:

[1]

This sum is then transformed using the sigmoid function to produce a single output, yj, which may be passed on to other neurons:

[2]

In the same way, the yk (prediction values) were calculated taking the corresponding output values of the previous hidden layer neurons (yj) as inputs and multiplying them by the connection weights wjk.

In the training step, the connection parameters, w, were adjusted such that, for the given input data set, the ANN-predicted output data set matched with the real output data set. At the beginning of the training phase, the network connection weights were random values between 0 and 1. For the given set of inputs to the network, the response of each neuron in the output layer (yk) was then calculated and compared with the corresponding real output response (rk). Then, the prediction error associated with the output response (Ek) was computed according to

[3]

The weights were adjusted to reduce prediction errors through a BP algorithm where Ek was back-distributed to the previous layers across the network. The optimization of the connection weights was performed by minimizing the error according to

[4]
where is the initial connection weight and μ is the learning coefficient. This coefficient controls the degree by which the connection weights were modified during the training phase. In this study, model selection was completed for the training set by fixing the learning coefficient and momentum to 0.001 and 0.1, respectively. Then, the trained networks were used to run a set of test data.

Solving Eq. [4] leads to

[5]
for the connection weights corresponding to the output layer or for the connection weights corresponding to the hidden layer:
[6]

Logarithmic sigmoid transfer function was used as the activation function for hidden and output layers. The development of the ANN model involved 2 basic steps, training/learning and testing/validation. During the training stage, the network used the inductive-learning principle to learn from a set of examples called the training set. Input variables of the ANN model were included microbial counts (TAC, YMC, CC), pH values, and obtained color data of yogurt during storage from the machine vision system, which was the color index value performed by the color analysis program looking at each pixel of the input yogurt images.

Out of 30 cases in the prediction time data set, 22 cases (∼73%) were used for training/testing and 8 cases (∼27%) were used for validation, randomly (Figure 3).

Normalized values in the range of 0 to 1 were used as input and output values. Normalization is a transformation performed on a single data input to distribute the data evenly and scale it into an acceptable range for the network. In normalizing data, the goal is to ensure that the statistical distribution of values for each net input and output is roughly uniform. In addition, the values should be scaled to match the range of the input neurons. Therefore, each input should be normalized in addition to the transformations performed on network inputs.

Proximity of the correlation coefficient to 1 and closeness of the mean square error to 0 shows the efficiency of the shelf life-prediction ANN models. The input layer had 16 neurons. The hidden layer nodes were chosen according to ANN performance. Training was continued until 10,000 epochs had been executed. For a given set of inputs to the network, the response to each neuron in the output layer was calculated and compared with the corresponding desired output response. Different ANN models were developed by using 20 different neurons in 1 hidden layer. The optimum model result was obtained in a hidden layer with 5 neurons; this model is shown in Figure 4.

  • View full-size image.
  • Figure 4. 

    Selected neural network structure: 16 inputs were used for the artificial neural network (ANN): 12 inputs corresponding to color blocks resulting from machine vision system analysis; 1 input corresponding to pH values; 3 inputs corresponding to total aerobic count (TAC), yeast and mold counts (YMC), and coliform counts (CC) resulting from microbial analysis. For the ANN used to predict the shelf life, 1 output was used and the output value was the “shelf-life”.

Statistical Analysis 

The degree of association between pH and microbial counts upon storage was determined by using PROC GLM of SAS (SAS Institute Inc., Cary, NC). Separation of the means (P<0.05) was accomplished by Duncan's multiple range test. All experiments and analyses were replicated 3 times. All analyses within experimental replications were performed in duplicate.

Back to Article Outline

Results and Discussion 

The extent of growth of certain classes of microorganisms and changes in sensory parameters are widely used to mark end of shelf life of foods (Muir, 1996; Lewis and Dale, 2000). In this study, pH, microbial counts (TAC, YMC, CC), and color data were used to mark shelf life (d 1 through 14) of set-type whole- and low-fat yogurts. Then, the ANN model was developed with the obtained data to predict the shelf life of yogurt.

Changes in pH and microbial counts during storage are shown in Table 1. Although the pH of low-fat yogurt did not change significantly (P>0.05), the pH of whole-fat yogurt decreased from 4.20 to 3.76 during storage (P<0.05). Total aerobic plate counts in low-fat yogurt samples increased by 2 log cycles from 3.43 to 5.12 log cfu/g after 14 d of storage. Total plate counts were lower for whole-fat yogurt samples, possibly due to the decrease in pH during storage. Yeasts and molds counts of low-fat yogurt increased during storage with counts ranging within 2 log cycles (from 3.45 to 5.63 log cfu/g; P<0.05), whereas YMC became undetectable for whole-fat yogurts after d 7. Flavor defects in yogurt, quarg, and (initially) yeast-free labneh have been reported to develop when counts of yeasts and molds reach levels of 105 cfu/g (Tamime and Robinson, 1999; Kadamany et al., 2003). The coliform counts of low-fat yogurt were high on d 1 (2.58 log cfu/g). Coliforms were present at a level of 1.26 log cfu/g in whole-fat yogurt and decreased progressively during storage, with growth becoming undetectable on the d 7 of storage, presumably due to the inhibitory effect of increased acid production (Tamime and Robinson, 1999). The presence of coliforms in the samples is indicative of postpasteurization contamination at one or more stages during processing.

Table 1. The pH measurements and microbial counts of low-fat and whole-fat yogurt cultures during storage for 1, 7, and 14 d1
Storage dayYogurt samplespHTotal aerobic counts, log (cfu/g)Yeast and mold counts, log (cfu/g)Coliform count, log (cfu/g)
1Low fat4.41±0.01a3.43±0.07a3.45±0.45a2.58±2.11a
Whole fat4.20±0.16b1.72±0.16b1.07±1.19b1.26±1.01b
7Low fat4.16±0.01c5.04±0.19c2.96±1.47c1.06±0.24c
Whole fat3.86±0.03b1.95±0.99dND2ND
14Low fat4.29±0.00c5.12±2.55a5.63±0.11dND
Whole fat3.76±0.01bNDNDND

a–dValues with different superscripts within the same column were statistically significant (P<0.05).

1Means±standard error of the means.

2Not detectable.

The color data measured by the machine vision system were reported as a percentage of the total area of the surface by a given color block (data not shown) as described in Luzuriaga et al. (1997). There was an increase in the percentage of the total area of color blocks having color of pale-greenish-yellow, grayish-yellow, light grayish-green, and yellowish-green, out of the 4,096 color blocks from d 1 to 14. The color block with grayish-greenish yellow color was dominant in both whole- and low-fat yogurts at the end of storage. The presence of these colors is associated with microbial spoilage of the food product. The color analysis data were parallel to pH and microbial count data.

The development of the ANN model involved 2 basic steps, training/learning and testing/validation. During the training stage the network used the inductive-learning principle to learn from a set of data called the training set. Out of 30 cases in the prediction time data set, 22 cases (∼73%) were used for training/testing and 8 cases (∼27%) were used for validation, randomly (Figure 3). Sixteen inputs were used for ANN: 12 inputs corresponding to color blocks resulting from machine vision analysis; 1 input corresponding to pH values; and 3 inputs corresponding to TAC, YMC, and CC resulting from microbial analysis. For the ANN used to predict the shelf life, 1 output was used and the output value was the “shelf-life” (Figure 4).

Determining the optimal ANN topology consisted of selecting the number of neurons that gave a minimum final error in a minimal number of iterations during training of the ANN. Moreover, the final use of the ANN must be considered to select the adequate ANN topology. It is necessary to take into account that light topologies (fewer neurons) involve smaller learning data sets than heavy topologies (more neurons). The learning step is therefore easier and faster and the final ANN calculation time is also shorter. Twenty configurations of ANN topology were conducted (Table 2). As seen from Table 2, configurations 14, 3, and 5 resulted in correlation coefficients (R2) of 0.992, 0.997, and 0.999, respectively. Because the correlation coefficient of configuration 5 was closest to 1, it was selected as the optimal ANN topology to give the best prediction.

Table 2. Correlation coefficient (R2) and mean square error (MSE) values of neurons in the topology optimization
Neurons in hidden layerMSER2
12.586080.8928
23.267110.6633
32.836370.997
43.886880.8994
51.410490.999
61.194080.935
74.587210.8912
86.548160.9214
97.596670.6281
104.555570.8185
111.138050.8086
121.306130.9739
131.733340.8952
141.011990.992
151.635290.6819
163.750310.4629
171.096030.8422
181.635290.7877
191.831380.1869
201.152060.9569

The performance of the ANN model for the training and validation data sets is presented in Figures 5 and 6, respectively. These figures show a comparison of model predictions with the experimental values of shelf life. The results show that the training procedure for shelf-life prediction was very successful and that a perfect match was obtained between the measured and the calculated values. The R2 values of the ANN model were 1 for the training set and 0.9996 for the testing/validation set.

  • View full-size image.
  • Figure 5. 

    Correlation of experimental vs. neural network values of time with training data set using the optimal network, with 1 hidden layer, 5 neurons per hidden layer, and a training data set of 22 cases.

  • View full-size image.
  • Figure 6. 

    Correlation of experimental vs. neural network values of time with validation data set using the optimal network, with 1 hidden layer, 5 neurons per hidden layer, and a validation data set of 8 cases.

The performance of ANN models was clearly superior for training the data sets and was reasonably good for testing the data sets. The ANN predictions for testing data sets can be considered to be good predictions, because the ANN predicted values for unseen data sets, whereas regression equations of the statistical models in McClure et al. (1994) and Zaika et al. (1994) were developed for the complete data sets (i.e., all data sets were used for developing the statistical models). Unless the number of data points is much larger than the number of parameters in the polynomial regression equations, there is no assurance that predictions for unseen data would be accurate (Specht, 1991; Barayani and Roberts, 1995). Bratchell et al. (1990) examined how the quantity of data sets can affect model predictions and reported that the response surfaces produced by the model developed from partial data sets were quite different from those produced by the model developed from complete data sets. Therefore, the predictions of statistical models may not be accurate for unseen data.

The simplified algebraic equations derived from the ANN model for the prediction of shelf life with pH, microbial counts, and color analysis are presented in Appendix A. With a priori knowledge of a yogurt's pH, microbial counts, and color analysis, control of its expiration date and prediction of shelf life can be computed using these equations compiled in a simple subroutine.

Back to Article Outline

Conclusions 

Combining neural networks with objective readings can result in powerful quality predictions. The modeling results showed that there was excellent agreement between the experimental data and predicted values, with a high determination coefficient (R2 = 0.9996) showing that the developed model was able to analyze nonlinear multivariant data with very good performance, fewer parameters, and shorter calculation time. Although ANN has been used in many applications in the agri-food industry, this is the first study modeling the prediction of shelf life of yogurt using ANN. The use of ANN provides an inexpensive and easy technique for evaluation of yogurt quality parameters.

The conceptual ANN model provides a database and an alternative generic framework for the modeling of shelf life of yogurt. This model has potential to be used as an alternative method to control the expiration date, estimate the shelf life of yogurt, and ensure the safety of the product.

Back to Article Outline

Acknowledgments 

We thank Necla Demir for granting permission to use her equipment and facility and we acknowledge her technical guidance. This research was supported by Süleyman Demirel University, Research Foundation (Project No: 0947-YL-04).

Back to Article Outline

Appendix 

Equations defining E1 to E6 and F1 to F5 are the summation and activation functions of each neuron, respectively, where c1 to c12 represent color values obtained from the machine vision system.

[A1]
[A2]
[A3]
[A4]
[A5]
[A6]
[A7]
[A8]
[A9]
[A10]
[A11]
[A12]

Back to Article Outline

References 

  1. Balaban MO, Luzuriaga DA, Yeralan S, Hasan R. Automated Quality Evaluation of Shrimp. In:  Jowitt R editors. Engineering and Food In ICEF 7. Sheffield, UK: Sheffield Academic Press; 1997;p. A105–A108
  2. Barayani J, Roberts TA. Mathematics of predictive food microbiology. Int. J. Food Microbiol. 1995;26:199–218
  3. Bishop MC. Neural network and their applications. Rev. Sci. Instrum. 1994;64:1803–1831
  4. Bratchell N, McClure PJ, Kelly TM, Roberts TA. Predicting microbial growth: Graphical methods for comparing models. Int. J. Food Microbiol. 1990;11:279–288
  5. British Standards Institution. Gerber Method for the Determination of Fat in Milk and Milk Products. London, UK: British Standards Institution; 1955;Standard No. 696, part 2
  6. Brosnan T, Sun DW. Improving quality inspection of food products by computer vision––A review. J. Food Eng. 2004;61:3–16
  7. Bochereau L, Bourgine P, Palagos B. A Method for prediction by combining data analysis and neural networks: Application to prediction of apple quality using near infra-red spectra. J. Agric. Eng. Res. 1992;51:207–216
  8. Bucinski A, Zielinski H, Kozlowska H. Artificial neural networks for prediction of antioxidant capacity of cruciferous sprouts. Trends Food Sci. Technol. 2004;15:161–169
  9. Demir N, Balaban MO, Brecht JK, Sims CA, Gulay MS. Objective quality assessment of modified atmosphere stored zucchini (Cucurbita pepo L., cv. Elite) slices using electronic nose. In: Proc. IFT Annual Meeting. Anaheim, CA. 2002;p. 181
  10. Demuth H, Beale M. Neural Network toolbox. User Guide. Version 4. Natick, MA: The MathWorks Inc.; 2001;
  11. Doganis P, Alexandridis A, Patrinos P. Time series sales forecasting for short shelf-life food products based on artificial neural networks and evolutionary computing. J. Food Eng. 2006;75:196–204
  12. Edwards EJ, Cobb AH. The effect of prior storage on the potential of potato tubers (Solanum tuberosum L.) to accumulate glycoalkaloids and chlorophylls during light exposure, including artificial neural network modelling. J. Sci. Food Agric. 1999;79:1289–1297
  13. Farkas I, Remenyi P, Biro A. A neural network topology for modelling grain drying. Comp. Electron. Agric. 2000;26:147–158
  14. Gerrard DE, Gao X, Tan J. Beef marbling and color score determination by image processing. J. Food Sci. 1996;61:145–148
  15. Hassoun MH. Fundamentals of Artificial Neural Networks. Cambridge, MA: MIT Press; 1995;
  16. Hornik K, Stinchcombe M, White H. Multilayer feedforward networks are universal approximators. Neural Netw. 1989;2:359–366
  17. Hussian MA, Shaur M, Ng NW. Predictions of pores formation (porosity) in foods during drying: Generic models by the use of hybrid neural network. J. Food Eng. 2002;51:239–248
  18. Jamieson V. Physics raises food standards. Physics World. 2002;1:21–22
  19. Kadamany EA, Khattar M, Haddad T, Toufeili I. Estimation of shelf-life of concentrated yogurt by monitoring selected microbiological and physicochemical changes during storage. Lebensm. Wiss. Technol. 2003;36:407–414
  20. Kavdir I, Guyer DE. Apple sorting using artificial neural networks and spectral imaging. Trans. ASAE. 2002;45:1995–2005
  21. Lewis M, Dale RH. Chilled yogurt and other dairy desserts. In:  Man D,  Jones A editor. Shelf Life Evaluation of Foods. Gaithersburg, MD: Aspen Publishers; 2000;p. 89–109
  22. Luzuriaga D, Balaban M, Yeralan S. Analysis of visual quality attributes of white shrimp by machine vision. J. Food Sci. 1997;61(1):113–118
  23. Ni H, Gunasekaran S. Food quality prediction with neural networks. Food Technol. 1998;5:60–65
  24. McClure PJ, de Blackburn CW, Cole MB, Curtis PS, Jones JE, Legan JD, et al. Modelling the growth, survival and death of microorganisms in foods: The UK food micromodel approach. Int. J. Food. Microbiol. 1994;23:265–275
  25. McDonald T, Chen YR. Separating connected muscle tissues in images of beef carcass rib eyes. Trans. ASAE. 1990;33:2059–2065
  26. Muir DD. The shelf life of dairy products: 2. Raw milk and fresh products. J. Soc. Dairy Technol. 1996;49:44–48
  27. Park B, Chen YR. Real-time dual-wavelength image processing for poultry safety inspection. J. Food Process Eng. 2000;23:329–351
  28. Rumelhart DE, Hinton GE, Williams RJ. Learning internal representation by error propagation. In:  Rumelhart DE,  McClelland JL editor. Chapter 8 in Parallel Distributed Processing: Exploration in the Microstructure of Cognition. Vol. 1. Cambridge, MA: MIT Press; 1986;
  29. Specht DF. A general regression neural network. IEEE Trans. Neural Netw. 1991;2:568–576
  30. Sun D-W. Computer vision—An objective, rapid and non-contact quality evaluation tool for the food industry. J. Food Eng. 2004;61:1–2
  31. Sun D-W, Du CJ. Segmentation of complex food images by stick growing and merging algorithm. J. Food Eng. 2004;61:17–26
  32. Tamime AY, Robinson RK. Yogurt science and technology. 2nd ed.. Boca Raton, FL: CRC Pres; 1999;
  33. Thibault J, Grandjean BPA. Process control using feedforward neural networks. J. Syst. Eng. 1992;2:198–212
  34. Wang HH, Sun DW. Assessment of cheese browning affected by baking conditions using computer vision. J. Food Eng. 2002;56:339–345
  35. Zaika LL, Moulden E, Weimer L, Phillips JG, Buchanan RL. Model for the combined effects of temperature, initial pH, sodium chloride and sodium nitrite concentrations on anaerobic growth of Shigella flexneri. Int. J. Food Microbiol. 1994;23:345–358

PII: S0022-0302(07)71759-9

doi:10.3168/jds.2006-591

Journal of Dairy Science
Volume 90, Issue 7 , Pages 3118-3125, July 2007