Digital-image analysis to predict weight and yields of boneless subprimal beef cuts

Análise de imagem digital para a previsão de pesos e rendimentos de cortes de carne bovina

Abstracts

Por várias décadas, a avaliação de carcaça bovina em sistemas de tipificação ou em pesquisas tem dependido de escores subjetivos e medidas obtidas manualmente, mas ultimamente tem havido um crescente interesse por novas tecnologias capazes de aumentar a acurácia das estimativas. Este trabalho teve como objetivo desenvolver equações para a previsão de pesos e rendimentos de cortes bovinos, através da análise de imagem digital (VIA) de uma seção do contrafilé da 12ª costela. As equações de previsão do peso dos cortes do traseiro especial (CUTS) apresentaram coeficientes de determinação (CD) de 0,84 e de 0,87 – 0,88, quando as variáveis independentes usadas eram os parâmetros VIA e o peso da meia carcaça (HC) ou o peso total do traseiro especial (TP), respectivamente. As equações de previsão do rendimento dos cortes do traseiro especial (CUTS%) representaram de 37,1 a 46,8% e de 21,3 a 30,6% da variação total, quando a principal variável independente utilizada na equação era HC ou TP, respectivamente. Nas equações de previsão dos pesos individuais dos cortes do traseiro especial o CD variou de 0,40 – 0,72 e de 0,43 – 0,74, usando as variáveis HC ou TP, respectivamente. O sistema de análise de imagem digital utilizado apresentou boa repetibilidade, podendo ser considerado um procedimento confiável para a estimativa do peso em cortes do traseiro especial e de alguns dos seus cortes individuais.

indústria da carne; carcaça bovina; análise de imagem; área de olho de lombo; espessura de gordura

NOTE

Digital-image analysis to predict weight and yields of boneless subprimal beef cuts

Análise de imagem digital para a previsão de pesos e rendimentos de cortes de carne bovina

Gustavo Adolfo Teira^I; Eduardo Tinois^II; Roberto de Alencar Lotufo^III; Pedro Eduardo de Felício^IV

^IDepto. de Industria de Alimentos, Universidad Nacional de Entre Ríos, 3200 - Concordia, ER, Argentina

^IIUNICAMP/CEB, C.P. 6040 - 13084-971 - Campinas, SP - Brasil

^IIIUNICAMP/FEEC - Depto.de Engenharia de Computação e Automação Industrial , C.P. 6101 - 13083-970 - Campinas, SP - Brasil

^IVUNICAMP/FEA - Depto. de Tecnologia de Alimentos, C.P. 6121 - 13081-970 - Campinas, SP - Brasil

^{Correspondence}

ABSTRACT

For several decades, beef carcass evaluation for grading or research purposes has relied upon subjective visual scores, and manually taken measurements, but in recent times there has been a growing interest in new technologies capable of improving accuracy of estimates. Equations to predict weight and yield of beef pistol subprimal cuts were developed in this work using digital image analysis (VIA) of the 12

^th rib steak. Equations to predict total pistol subprimal cuts weight (CUTS) had coefficients of determination (CD) of 0.84, or 0.87 to 0.88, when the independent variables were the VIA parameters and the half carcass weight (HC) or the total pistol weight (TP), respectively. The predicted values for the total seven subprimal cuts, as a percentage of half carcass weight (CUTS%), presented CD values ranging from 0.37 to 0.47, or 0.21 to 0.31, using HC or TP as a principal independent variable. Likewise, the equation for weight of the individual subprimal cuts had CD values ranging from 0.40 to 0.72, or 0.43 to 0.74 using HC or TP, respectively. In this research, the developed VIA procedure has demonstrated good repeatability and accuracy to estimate the total pistol subprimal weights, and some individual subprimal weights.

Key words: beef industry, beef carcass, image analysis, rib eye area, fat thickness

RESUMO

Por várias décadas, a avaliação de carcaça bovina em sistemas de tipificação ou em pesquisas tem dependido de escores subjetivos e medidas obtidas manualmente, mas ultimamente tem havido um crescente interesse por novas tecnologias capazes de aumentar a acurácia das estimativas. Este trabalho teve como objetivo desenvolver equações para a previsão de pesos e rendimentos de cortes bovinos, através da análise de imagem digital (VIA) de uma seção do contrafilé da 12^a costela. As equações de previsão do peso dos cortes do traseiro especial (CUTS) apresentaram coeficientes de determinação (CD) de 0,84 e de 0,87 – 0,88, quando as variáveis independentes usadas eram os parâmetros VIA e o peso da meia carcaça (HC) ou o peso total do traseiro especial (TP), respectivamente. As equações de previsão do rendimento dos cortes do traseiro especial (CUTS%) representaram de 37,1 a 46,8% e de 21,3 a 30,6% da variação total, quando a principal variável independente utilizada na equação era HC ou TP, respectivamente. Nas equações de previsão dos pesos individuais dos cortes do traseiro especial o CD variou de 0,40 – 0,72 e de 0,43 – 0,74, usando as variáveis HC ou TP, respectivamente. O sistema de análise de imagem digital utilizado apresentou boa repetibilidade, podendo ser considerado um procedimento confiável para a estimativa do peso em cortes do traseiro especial e de alguns dos seus cortes individuais.

Palavras-chave: indústria da carne, carcaça bovina, análise de imagem, área de olho de lombo, espessura de gordura

INTRODUCTION

The value of beef cuts should reflect both consumer demand and reliability, and marketing standards which emphasize quality, uniformity and consistency of the meat. Therefore, an evaluation procedure for predicting weights and yields of carcasses and beef retail cuts becomes of great importance for the beef industry (Cross & Belk, 1994).

Several studies have described different methods that have great potential for industrial applications: ultrasound, video-image analysis, nuclear magnetic resonance, optical probes and bioelectrical impedance. These technologies have been used successfully in laboratories. However their utilization in the beef industry is restricted. Only fast, consistent, accurate, non-destructive and reasonable priced methods are of interest to the industry (Swatland, 1995).

Digital image analysis has been identified as the most efficient procedure to evaluate yield of beef carcasses (Cross et al., 1983). It has also been utilized for other purposes: determination of color and fat thickness (Monin, 1998), marbling scores in beef (Albrecht et al., 1996) and pork (Scholz et al., 1996) and water retention capacity in beef (Irie et al., 1996). Li et al. (1999) reported that tenderness of beef could be predicted on the basis of color, marbling and texture measures obtained by image analysis.

The use of this technique has been reported to develop an objective system for carcass classification which has been used as a standard of payment to European Union beef producers based mainly on the determination of conformation and finish traits (Borggaard et al., 1996). Other studies have reported the use of video image analysis of steaks removed from the 12^th rib for prediction of weight and yield of subprimal cuts (Cross et al., 1983; Wassenberg et al., 1986; Shackelford et al., 1998).

No information on the use of such systems for the determination of traits of carcasses and cuts produced from Zebu steers, predominantly slaughtered in Brazil, is available. In addition, it is necessary to follow current trends concerning international beef commercialization. The objectives of this work were to predict weight and yield of whole and individual pistol cuts by regression equations using variables obtained through image analysis of the 12^th rib steak, and to test the accuracy and repeatability of the procedure.

MATERIAL AND METHODS

Selection of carcasses - Carcasses produced from 51 Nelore crossbred steers, slaughtered at Bertin beef processing plant (Lins, SP-Brazil) were selected according to various carcass weight, maturity and fatness to obtain a trait variability as large as possible.

Cuts utilized - Carcasses were chilled for 24 hours (1 – 2°C) and broken between the 5^th and the 6^th rib into hindquarter (HQ) and forequarter (FQ). The former was separated from flank and plate (FP) using a saw and knife, 50 mm to the side of the Longissimus muscle to obtain the pistol (TP). This was further boned to obtain the following cuts: top round, bottom round, eye round, knuckle, rump, cube roll & strip loin, and tenderloin side strap off. Connective tissue, subcutaneous fat and trimmings were separated following the standard procedure for cuts exported to Germany by the beef plant.

Sample preparation - A 2.5 cm steak sample was removed from the 12^th rib before boning each cube roll & strip loin, and the rib eye area of the steak was delineated for later measurement with a compensating polar planimeter. The sample weight was not considered in the total weight of the cut.

Image analysis - Images were captured using a Sony^Ò digital camera, MVC-FD71. Bone-in steak samples were placed flat on a non-glare black surface and illuminated with two tungsten, 150 W halogen bulbs. A reference white square paper was positioned on the black surface, next to the object such that the same distance and focus were kept when images were captured.

Digital images, saved in a 3.5" floppy disk, were processed using Matlab Mathematical Morphology Toolbox (Barrera et al., 1998). With these tools, a program for segmentation of the original digital image and extraction of selected parameters was developed. Figure 1a: input image (a known square is included in the image to calibrate the pixel dimensions).

To compute the total rib steak area (Figure 1b, TRSA, cm²) the total rib was extracted from the background using a standard thresholding technique. The image illustrates only the contour of the rib and the square. The rib eye area (Figure 1c, REA, cm²) was detected by an automatic procedure based on the so-called watershed transform (Beucher & Meyer, 1993; Lotufo & Falcao, 2000). First the input image was pre-processed by an edge enhancement procedure. Then the best closed contour was found based on a prior knowledge of the approximate position of the rib eye. The watershed transform was responsible to find the best closed contour on the edge enhanced image. The rib eye area percentage (REA%) resulted from REA/TRSA x 100.

Fat thickness (FT, mm) is shown in Figure 1d. The black line indicates where the thickness was measured. This line was automatically drawn perpendicular to the rib eye contour. It was traced from a point located on a straight line segment determined by the ends of the rib eye at ¼ of its length from the left.

A rectangle was automatically generated from the points of the rib eye ends (Figure 1e). It shows the location and subsequent extraction of fat area, determined by subcutaneous fat in the lower part of rib eye area (Figure 1f, FA, cm²). This parameter was measured from an image generated using a thresholding technique to detect the fat area. Dark pixels were not considered and white pixels (fat) were counted. The fat area percentage (FA%) resulted from FA/TRSA x 100.

Statistical analysis - VIA parameters and measured side carcass traits were statistically analyzed to select the best linear regression equations for weight and yield prediction. These equations were assessed using the coefficient of determination (CD), statistical Cp (Mallows, 1973) and mean square error (MSE) with Stepwise/Maxr procedure, SAS Institute (1985). A second image was captured to test repeatability of the VIA procedure following the Bland & Altman (1986) methodology.

RESULTS AND DISCUSSION

Table 1 shows mean values, standard deviation, maximum and minimum values obtained for the 51 half carcasses used in the present study. The great variability shown by the dependent and independent variables were expected on account of the applied selected methodology.

The pistol subprimal cut weights (CUTS), their weight as a percentage of the half carcass weight (CUTS%) and the weight of individual pistol subprimal cuts did not presented correlations (P > 0.01) with the VIA parameters. These parameters however in some cases were correlated (P < 0.01): TRSA and REA (r = 0.62), REA and REA% (r = 0.82), FT and FA (r = 0.59), FT and FA% (r = 0.59) and FA and FA% (r = 0.98). Rib eye area, measured by both VIA (REA) and manually with a planimeter (AOLPLAN) presented a significant correlation (P < 0.01) of 0.70. The graphic representation of the observed values of CUTS and CUTS% variables against independent variables used did not show non-linear relationships.

Weight prediction of the seven subprimal cuts - To obtain the prediction equation of the CUTS variable, the two best individual indicators, HC or TP and all VIA parameters were used (

Tables 2 and

3). When the HC variable was used, the best equations accounted for 83.5 to 84.0% of the observed variation of the subprimal cut weights. When HC was replaced by TP, the coefficient of determination of the best equations varied from 0.87 to 0.88.

In order to select the most adequate models of regression equations, belonging to the same subset of variables, a Cp graph (Mallow statistics) was plotted against p (number of parameters included in the equation). If there was no bias, Cp should be equal to p. In practice, models of Cps near to, but preferably less than p, should be chosen (Pike, 1986).

This kind of graphical representations, subsets corresponding to equations 1 to 6 (Table 2) and 7 to 12 (Table 3), indicates that the inclusion of VIA parameters results in obtaining the most appropriate models.

Comparison of predictions in relation to carcass weight - The simultaneous graphic representation of both observed and predicted values of the CUTS variable, using the different regression equations, versus chilled carcass weight (kg) ranges (A:230-250; B:250-270; C:270-290; D:290-310; E:310-330) conveniently selected, indicated the results of the different models under such conditions.

These equations including the HC variable tended to overestimate values for heavier carcasses (E range) and equations including the TP variable tended to underestimate values for lighter carcasses (A range). Equation n° 8 (Table 3) was considered the most adequate for any weight range carcass. These observations were validated through the Tukey test of residue mean comparison (observed values minus predicted values) for each equation and weight range.

Prediction of total yield of seven pistol subprimal cuts - Prediction of yield accounted for 37.1 to 46.0% or 21.3 to 30.6% of the variation when the major independent variable was HC or TP, respectively (Tables 4 and 5).

Prediction of weight of individual pistol cuts - It is important to predict the weight of individual cuts in order to determine the true value of each carcass. That is why VIA parameters were used to develop different prediction equations. The results obtained (Table 6 and 7) presented coefficients of determination that varied from 0.40 to 0.72 or from 0.43 to 0.74 when the major independent variable was HC or TP, respectively. These prediction equations were considered satisfactory for some cuts (knuckle, rump) but they presented lower CDs than those reported by Shackelford et al. (1998). German export standards used at the Bertin beef plant may account for these lower values since the removal of all visible fat is required, connective tissue and the trimming of the cuts which may vary from one employee to another.

Repeatability of the digital image analysis procedure - Twenty five digital image sample duplicates were utilized to test repeatability of the procedure. VIA parameters determined from duplicates were used in regression equation n^o6 to predict the CUTS variable. The new estimated values were compared to those predicted by the same equation using VIA parameters which were obtained with the 25 original images. The repeatability of the system remained within the established limits (± 1.96 S.D).

CONCLUSIONS

Weight prediction equations with a high degree of accuracy for pistol cuts were obtained using the digital image analysis procedure, equations n^o6 and n^o12 being especially recommended. When grouped according to chilled carcass weight range, prediction presented a tendency to over or underestimate extreme values (ranges A and E). However equation n° 8 made acceptable predictions for all considered weight ranges. Yield predictions of pistol subprimal cuts had a relatively low degree of accuracy. Prediction of individual pistol subprimal cut weights (kg) was good in some cases (knuckle, rump) but all values were below the expected. Current trends towards standardization of beef cuts especially for export require further work to be carried out along these lines to develop a "prototype" with features suitable for its implementation in industrial continuous lines.

ACKNOWLEDGEMENTS

To CAPES/SPU-Argentina for the scholarship, of the first author. Fundación Antorchas, Argentina, for partially funding this project. Engº Pedro Bertin Beloto, R&D manager of Bertin Beef Processing Plant, Lins, S.P.; Prof. Dr. Armando Infante, UNICAMP/IMECC, and Department of English from the School of Food Science at the University of Entre Ríos, Concordia, ER, Argentina.

Received January 09, 2002

• ALBRECHT, E.; WEGNER, J.; ENDER, K. A new technique for objective evaluation of marbling in beef. Fleischwirtschaft, v.4, p.11-16, 1996.
• BARRERA, J.; BANON, G.J.; LOTUFO, R.A.; HIRATA JR., R. MMach: a mathematical morphology toolbox for the KHOROS system. Journal of Electronic Imaging, v.7, p.174-210, 1998.
• BEUCHER, S; MEYER, F. The morphological approach to segmentation: the watershed transformation. In: DOUGHERTY, E.R. (Ed.) Mathematical morphology in image processing. New York: Marcel Dekker, 1993. cap.12, p.433-481.
• BLAND, J.M.; ALTMAN, D.G. Statistical methods for assessing agreement between two methods of clinical measurement. The Lancet, v.1, p.307-310, 1986.
• BORGGAARD, C.; MADSEN, N.T.; THODBERG, H.H. In line image analysis in the slaughter industry, illustrated by beef carcass classification. Meat Science, v.43, p.s151-s163, 1996.
• CROSS, H.R.; BELK, K.E. Objective measurements of carcass and meat quality. Meat Science, v.36, p.191-202, 1994.
• CROSS, H.R.; GILLILAND, D.A.; DURLAND, P.R.; SEIDEMAN, S. Beef carcass evaluation by use of a video image analysis system. Journal of Animal Science, v.7, p.908-917, 1983.
• IRIE M.; IZUMO, A.; MOHRI, S. Rapid method for determining water-holding capacity in meat using video image analysis and simple formulae. Meat Science, v.42, p.95-102, 1996.
• LI, J.; TAN, J.; MARTZ, F.A.; HEYMANN, H. Image texture features as indicators of beef tenderness. Meat Science, v.53, p.17-22, 1999.
• LOTUFO, R.; FALCAO, A. The ordered queue and the optimality of the watershed approaches. In: GOUTSIAS, J.; VINCENT, L.; BLOOMBERG, D. (Ed.) Morphology and its application to image and signal processing Dordrecht: Kluwer Academic, 2000. p.341-350. (Computational Imaging and Vision, 12)
• MALLOWS, C.L. Some comments on Cp. Technometrics, v.15, p.661-675, 1973.
• MONIN, G. Recent methods for predicting quality of whole meat. Meat Science, v.49, p.s231-s243, 1998.
• PIKE, D.J. A practical approach to regression. In: PIGGOT, J.R. (Ed.) Statistical procedures in food research. London: Elsevier Applied Science, 1986. cap.3, p.61-100.
• SAS INSTITUTE. User's guide: Statistic. 5.ed. Cary: Statistical Analysis System Institute, 1985. 956p.
• SCHOLZ, A.; PAULKE, T.; EGER, H. Determining the degree of marbling in the pig. Use of computer-supported video picture analysis. Fleischwirtschaft International, v.1, p.9-11, 1996.
• SHACKELFORD, S.D.; WHEELER, T.L.; KOOHMARAIE, M. Coupling of image analysis and tenderness classification to simultaneously evaluate carcass cutability, longissimus area, subprimal cut weights, and tenderness of beef. Journal of Animal Science, v.76, p.2631-2640, 1998.
• SWATLAND, H.J. Objective assessment of meat yield and quality. Trends in Food Science and Technology, v.6, p.117-120, 1995.
• WASSENBERG, R.L.; ALLEN, D.M.; KEMP, K.E. Video image analysis prediction of total kilograms and percent primal lean and fat yield of beef carcasses. Journal of Animal Science, v.62, p.1609-1616, 1986.

Publication Dates

• Publication in this collection
  15 May 2003
• Date of issue
  2003

History

• Received
  09 Jan 2002