Next Article in Journal
Effects of Processing Methods and Conditioning Temperatures on the Cassava Starch Digestibility and Growth Performance of Broilers
Next Article in Special Issue
Distinguishing Natural Infections of the Bovine Mammary Gland by Staphylococcus from Streptococcus spp. Using Quantitative Milk Proteomics
Previous Article in Journal
A Descriptive Methodology for Studying the Ontogeny of Object Play and Breed Differences in Dogs (Canis lupus familiaris)
Previous Article in Special Issue
Precision Livestock Farming: What Does It Contain and What Are the Perspectives?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Alternatives to Locomotion Scoring for Detecting Lameness in Dairy Cattle in Tanzania: Infrared Thermography

by
Chacha W. Werema
1,2,*,
Linda J. Laven
1,
Kristina R. Mueller
1 and
Richard A. Laven
1
1
School of Veterinary Science, Massey University, Private Bag 11 222, Palmerston North 4442, New Zealand
2
College of Veterinary Medicine and Biomedical Sciences, Sokoine University of Agriculture, Morogoro 67 115, Tanzania
*
Author to whom correspondence should be addressed.
Animals 2023, 13(8), 1372; https://doi.org/10.3390/ani13081372
Submission received: 10 March 2023 / Revised: 6 April 2023 / Accepted: 7 April 2023 / Published: 17 April 2023
(This article belongs to the Special Issue Cattle Health, Production, Population Medicine and Therapeutics)

Abstract

:

Simple Summary

Locomotion scoring requires skilled, trained observers to accurately detect lameness. Many studies have thus evaluated infrared thermography as an alternative lameness detection method as it does not require a skilled observer. However, there are few reports of the use of infrared thermography in cattle in tropical environments like Tanzania. This study, therefore, aimed to assess whether using an infrared camera to measure the foot skin temperature of hind limbs could potentially be used as an alternative on Tanzanian dairy farms. Three study farms were visited twice each during the afternoon milking on consecutive days. Locomotion scoring using a 4-point scale (0–3) was conducted on the first day as the cows exited the milking parlour after being milked. On the following day, the hind limbs of the cows were thermally imaged while they were standing in the milking parlour, using a forward-looking infrared camera. Mean foot skin temperature increase was associated with an increase in locomotion score; for example, the mean temperature was higher for cows with a locomotion score of 3 than those with a score of 2. Therefore, the present study confirmed that measuring foot skin temperature using an infrared camera has the potential to be employed for detecting lameness on Tanzanian dairy farms. However, improvements in accuracy and reductions in infrared camera costs are needed.

Abstract

Lameness detection is a significant challenge. Locomotion scoring (LS), the most widely used system for detecting lameness, has several limitations, including its subjective nature and the existence of multiple systems, each with its own advantages and disadvantages. Therefore, this study aimed to evaluate whether the foot skin temperature (FST) of hind limbs, as measured using infrared thermography (IRT), could potentially be used as an alternative on Tanzanian dairy farms. Each of the three study farms were visited twice during the afternoon milking on consecutive days, with a total of 170 cows assessed. DairyNZ LS (4-point scale (0–3)) was undertaken on the first day as the cows exited the milking parlour after being milked, while on the following day, the plantar aspect of the hind limbs of the cows was thermally imaged while they were standing in the milking parlour, using a handheld T650sc forward-looking infrared camera. Mean FST was higher for cows with a locomotion score of 1 than those with a score of 0; higher for cows with a locomotion score of 2 than those with a score of 1; and higher for cows with a locomotion score of 3 than those with a score of 2, with each one-unit locomotion score increase being associated with a 0.57 °C increase in mean temperature across all zones. The optimal cut-off point of 38.0 °C for mean temperature across all zones was identified using a receiver operator characteristic curve. This cut-off point had a sensitivity of 73.2% and a specificity of 86.0% for distinguishing cows with a locomotion score ≥ 2 (clinical lameness). The prevalence of clinical lameness across all three farms was 33%, which meant that only 72% of cows with a mean FST across all zones ≥ 38.0 °C had been identified as clinically lame using LS. This study confirmed that IRT has the potential to be used to detect lameness on Tanzanian dairy farms. However, before it can be widely used, improvements in accuracy, especially specificity, are needed, as are reductions in equipment (IR camera) costs.

1. Introduction

Early lameness detection accompanied by effective treatment is crucial to minimise the pain and discomfort associated with lameness [1,2,3], as well as to decrease the risk of irreversible claw damage [4].
Visual locomotion scoring (LS) is the most commonly used active lameness detection method on dairy farms [5]. A wide range of different systems have been used; Schlageter-Tello et al. [5] identified 25 different LS systems that had been published in the peer-reviewed literature. Across all systems, the most critical challenge of LS systems is their subjective nature, with both within- and between-observer variation being high, especially when training is limited [6,7,8,9]. Thus, training is particularly important when farmers and farm staff are undertaking lameness detection using LS.
An objective technique for identifying lame cows would thus be useful for on-farm detection of lameness, especially where farm staff training is difficult to achieve (e.g., because of lack of infrastructure or due to limited farmer knowledge). This is the situation in Tanzania, where active lameness detection using LS is extremely rare. In Tanzania, farm size is generally small and infrastructure limited, making training at the farm level expensive and challenging.
One potential alternative to visual LS is infrared thermography (IRT), a non-invasive technique that records body surface temperature and produces a pictographic representation of the scanned anatomical area [10]. Body extremities and surface temperature mainly depend on blood perfusion and tissue metabolism [11]. Thus, changes in blood flow can influence the amount of radiated heat and can thus be detected by IRT [10]. One reason for changes in tissue blood flow is the inflammatory response. Thus, IRT has the potential to detect systemic infectious lameness such as foot-and-mouth disease [12], and localised infectious lameness such as bovine digital dermatitis [13,14,15], as well as non-infectious claw horn disease [16,17,18].
However, hoof temperature can be affected by factors other than foot diseases, such as physiological status [19,20], environmental factors [21,22], and activity level [23]. Thus, the production system and management of the cows in a herd can significantly affect hoof temperature, potentially altering the accuracy of IRT as a means of lameness detection. Thus to determine the usefulness of IRT in a system, IRT needs to be tested in that system [24]. Unfortunately, no data are available for tropical pastoral systems such as those which are common in many regions of Tanzania. Therefore, the aim of this study was to investigate the association between the foot skin temperature (FST) of hind limbs at afternoon milking and locomotion scoring on three Tanzanian dairy farms where cows grazed tropical pasture between morning and afternoon milking.

2. Materials and Methods

2.1. Study Area and Animals

This study was undertaken in the Morogoro region of eastern Tanzania. This region has a sub-humid tropical climate with two wet seasons per year (long rains from March to June and short rains from October to December). Mean temperature ranges from 27 to 33.7 °C and 14.2 to 21.7 °C during the dry and wet seasons, respectively.
A convenient selection of three dairy farms was made for this study. All three herds had >50 cows, were milked twice daily, grazed natural pasture off-farm between morning and afternoon milking, and had a flat concrete surface outside the milking parlour which was suitable for locomotion scoring. Calving was all year round on the three study farms, and the herds were composed of mostly multiparous milking cows, as no regular replacement of cows was done on any of the farms. All three farms had only dairy breed cattle—a mixture of breeds and crossbreeds of European dairy cattle (Friesian, Ayrshire, Jersey, and their crossbreeds) and grazed their cattle on natural pastures (principally Hyparrhenia spp., Megathyrsus maximus, Cenchrus ciliaris, and Brachiaria spp.) for approximately 8 h after morning milking.
The distance between the grazing area and milking parlour on all farms was approximately 5 to 10 km with no constructed trackways. After the afternoon milking session, cows were allowed to graze around the farm before being taken to the free stall barns in the evening (6 pm), where they were given hay (principally, from Chloris gayana and Pennisetum purpureum) and homemade concentrates. All three farms had herringbone milking parlours. Farms 1 and 2 used milking machines (reverting to hand milking when there was no power), while farm 3 employed hand-milking only. During the study period, the milking herds had 52, 58, and 60 cows on farms 1, 2, and 3, respectively.
Routine hoof trimming and systematic locomotion scoring were not undertaken on any of the farms, and none of the farms had accurate lameness treatment records.

2.2. Study Visits

In March 2020, farms were visited during the afternoon milking on two consecutive days. Locomotion scoring was undertaken on the first day, and on the second day the plantar aspects of both hind feet of all cows were imaged using an infrared camera.

2.3. Locomotion Scoring

All locomotion scoring was undertaken by the first author (CWW) who was trained in the DairyNZ lameness score using a combination of video and supervised scoring of live animals (see Werema et al. [24] for further details). Prior to the study commencing in March 2020, CWW rewatched the training videos as a refresher [25,26].
On each farm, all milking cows were locomotion scored and their ear tag number recorded as they exited the milking parlour. The site chosen for scoring observations on each farm was a flat well-maintained concrete surface that was at least 25 m long, and which was cleaned after every milking session. This area allowed the observer to view at a distance so as not to disturb cow flow and locomotion.

2.4. Infrared Thermography

During the afternoon milking following the LS, IRT imaging was performed using a handheld T650sc Forward-looking Infrared camera (FLIR Systems, Oregon, United States) that at foot skin temperatures had an accuracy of ±0.4 °C. Emissivity was set at 0.95. Pictures were taken while the observer was standing at a distance of approximately 1 m from the cow, with the plantar aspects of both hind feet being imaged. Prior to imaging, no foot preparation was undertaken, except that the cows walked through a footbath containing only water on entering the collecting yard for the milking parlour approximately 30 min before milking.
Foot images were analysed using FLIR Tools software (FLIR Systems Oregon, United States) with estimates of the surface temperature obtained from seven zones on each hind foot as described by Werema et al. [24] (see Figure 1). The maximum temperature for each zone was used for analysis in line with previous infrared studies aimed at lameness detection in the cow [24,27,28]. Mean ambient temperature during this study was 29.4 °C.

2.5. Statistical Data Analyses

SPSS version 27 (IBM Corporation, Armonk, NY, USA) was used for all data analysis unless otherwise reported. Descriptive statistics were created for each zone temperature measure. The normality of foot temperature was visually assessed using Q-Q plots and histograms, followed by checking the skewness and kurtosis statistics. A generalised linear marginal repeated measures model was then used to evaluate the effect of the farm, foot, and zone within the foot on skin temperature. Farm and foot (right or left hind) were the independent variables, zone within foot the repeated variable and skin temperature the outcome variable. Interaction between farm and foot was tested and removed from the final model as it was non-significant (p > 0.05). Covariance structure was identified using Quasi-likelihood under the independence model criterion (QIC). Residuals were checked for normality using Q-Q plots and histograms. Post-hoc pairwise comparisons between zones were undertaken using the Šidák correction for multiple comparisons [29].
The association between locomotion score and foot temperature was tested using six temperature measures (see Table 1 for temperature definitions). Univariable analyses were first performed to assess the association between the outcome (foot temperature) and predictor (locomotion score) variables. This analysis identified significant heteroscedasticity when we compared foot temperature across locomotion scores. Therefore, we employed a generalised linear model with an identity link function and robust estimators [30] to analyse the association between foot skin temperature and locomotion score. In model fitting, each temperature definition was used as the outcome variable, with farm and locomotion scores as the predictor variables. Two-way interactions between farm and locomotion scores were added, and backwards selection was then used to remove any interactions where p > 0.05. Both main effects were kept in the final model irrespective of p-value.
Six receiver operator characteristic (ROC) curves were created, one for each temperature definition, with categorised locomotion score (lame (locomotion score ≥ 2) vs. not lame (locomotion score < 2)) to establish the sensitivity and specificity of IRT to predict locomotion score ≥ 2. Area under the curve (AUC) and coordinates of the curve (CC) were then used to evaluate a model’s predictive accuracy. Optimal temperature cut-off values were identified by maximising sensitivity plus specificity. The statistical package MedCalc Version 19.5.1 (MedCalc Software, Ostend, Belgium) was then employed to calculate positive and negative predictive values for those optimal cut-off points.

3. Results

All milking cows were locomotion scored on all three farms. The distribution of locomotion scores is summarised in Table 2. Across the three farms, 33% of cows were identified as being lame (locomotion score ≥ 2).

3.1. Effect of Foot and Foot Zone on Skin Temperature

On average the FST of the left foot was higher than the right foot (37.510 °C and 37.411 °C, respectively). However, it is unlikely that this is biologically significant as mean difference was only 0.099 °C (95% CI: 0.003–0.194). However, biologically significant differences between zones were identified. For example, the difference between the zone with the lowest mean temperature (zone 6) and the zone with the highest mean temperature (zone 4) was 1.57 °C (95% CI: 1.392–1.748). Results for all zones and their comparisons are summarised in Table 3 and Table 4, respectively. Mean temperatures were higher for zones on the lateral claw than their equivalent zones on the medial claw (see Table 3); these differences were 0.243 °C, 0.284 °C, and 0.163 °C for zones 1 vs. 5, 2 vs. 6, and 3 vs. 7, respectively.

3.1.1. Infrared Thermography versus Locomotion Scoring

The mean FST for different locomotion scores is summarised in Table 5. For all six FST measures (see Table 1 for definitions), the temperature was higher for cows with a locomotion score of 1 than those with a score of 0; higher for cows with a locomotion score of 2 than those with a score of 1; and higher for cows with a locomotion score of 3 than those with a score of 2. The comparisons of all six-foot skin temperature measures versus different locomotion scores are presented in Table 6. As an example, for mean temperature, the mean difference between cows with scores 0 and 1 was 0.91 °C (95% CI: 0.498–1.313); between scores 1 and 2 cows, it was 0.83 °C (95% CI: 0.431–1.232); and between scores 2 and 3 cows, it was 0.65 °C (95% CI: 0.026–1.271).

3.1.2. Association of Foot Temperatures and Locomotion Scores

A linear association was demonstrated between individual cow locomotion scores and foot skin temperature for all six temperature measures. Detailed data are presented in Table 7. As an example, for mean temperature, each one-unit locomotion score increase was associated with a 0.57 °C (95% CI: 0.46–0.70) rise in mean temperature.

3.1.3. A Receiver Operating Characteristic (ROC) Analysis

A receiver operating curve is illustrated in Figure 2. Optimal threshold values, area under the curve and calculated parameters for temperature measures are summarised in Table 8.

4. Discussion

The objective of the current study was to assess the usefulness of infrared thermography (IRT) as a tool for detecting lameness in dairy cattle partly grazed and housed in a tropical region (Tanzania) versus visual locomotion scoring (LS).

4.1. Suitability of IRT for Measuring Foot Skin Temperature during Milking

On all three farms, unlike in New Zealand [24], there was sufficient time during milking to collect thermal images from all cows. Hand milking increased the time available for image collection as it was slower than when the machine was used; but during hand milking, the milker obstructed the cow’s feet, and the cows were more restless, leading to frequent changes in foot posture. Thus, during hand milking, images were collected prior to a cow being milked.
One limitation of this study is that only hind feet were imaged, as cows were photographed during milking. Not all lameness is hind-feet related; although data from housed cows suggest that hind limb lameness accounts for >90% of cases [31], in pasture-based dairy cattle in New Zealand, only 67% of lameness was in the hind limb [32]. No such data are available from Tanzanian dairy cattle. However, it is likely that only recording the FST of hind limbs lowered the sensitivity of IRT for lameness detection in these cattle. Further research is required to identify how much the sensitivity was lowered.
Another potential issue with IRT under Tanzanian conditions is the imaging of dirty feet. Prior to milking on all three farms, the cows were walked through a footbath containing water. This was routine practice on all three farms and was undertaken with the intention of reducing dirt on feet, although it is unlikely that there was any impact on foot cleanliness [33]. Although this study was undertaken in the dry season, the feet of all cows were generally clean and were not washed prior to IRT. In the rainy season, the accumulation of mud on the feet during grazing might influence the accuracy of IRT, especially if the feet are not cleaned before imaging. However, Stokes et al. [27] reported that in the UK, maximum sensitivity plus specificity of IRT was achieved in feet that were not cleaned. Further longitudinal research is required to investigate whether this is also the case in Tanzania.

4.2. Effect of Claw and Zone on Skin Foot Temperature

In these cattle, the lateral claw had a higher mean temperature than the medial claw, with the highest difference of 0.28 °C being between zones 2 and 6 (see Table 3). This difference is somewhat larger than the difference of 0.1 °C we found under New Zealand conditions using the same camera and protocol for IRT [24], but it is smaller than the differences of 1.0 to 1.7 °C that have been reported in studies on housed cows [17,19]. The reason for this difference is unclear, but may be related to different protocols and different equipment.
The effect of the claw was consistent across zones, with the mean temperature for zones on the lateral claw being higher than their equivalents on the medial claw (see Table 3). Additionally, the order was consistent across claws, such as if on the lateral claw, zone 1 had a higher mean temperature than zone 2, the same applied to zones 5 and 6 on the medial claw. However, the order was not the same as that reported by Werema et al. [24]. In the current study, zones 1 and 5 (coronary band) had a higher mean temperature than zones 3 and 7 (below accessory digit), whereas Werema et al. [24] reported that zones 3 and 7 had a higher temperature than zones 1 and 5. The reason for this difference is unclear. It may be related to differences in hindlimb disease (e.g., increased infectious lameness in Tanzanian cattle). However, as we did not record hoof or limb lesions in either study, this has to remain a suggestion.

4.3. Infrared Thermography as a Predictor of Locomotion Score

In this study, mean foot skin temperature increased as locomotion scores increased (Table 5), consistent with previous studies of IRT and locomotion scores that used different protocols and devices [24,34,35]. Consistent with Werema et al. [24], the current study found that each of the locomotion scores (0, 1 and ≥2) had significantly different mean FST. The key difference between our two studies in regard to FST and locomotion score was that in Tanzanian cattle, the optimal cut-off point was 38.0 °C, whereas in New Zealand, our optimal cut-off was 34.5 °C [24]. However, both these cut-offs are much higher than the cut-offs used in studies on housed cattle in the Northern Hemisphere, e.g., 23.3 °C [35], 25.25 °C [21], 25.5 °C, [34], and 27.0 °C [27]. The reason for the large difference in optimal cut-off between our studies and those previous studies is unclear. However, it may be related to protocol differences (e.g., measuring different sites, using different equipment or cleaning feet before imaging), differences in the environment such as ambient temperature, or differences in the cause of lameness (e.g., infectious vs. non-infectious lameness). Further data on the factors driving these differences in the optimal cut-off points are required. However, our data suggest that optimal cut-offs are likely to be protocol and production system specific. Thus, cut-offs should not be transferred from one study in one system to another, even if similar protocols are used.
The specificity and sensitivity of IRT at our 38.0 cut-off point for identifying cows with a locomotion score of ≥2 were both moderate (86.0 and 73.2%, respectively), and lower than the specificity and sensitivity we reported in New Zealand of 92.4 and 80.0%, respectively [24]. However, as with the difference between the two studies in relation to optimal cut-off, the difference between our two studies is much less than the differences between previous studies, which have reported the specificity and sensitivity of IRT. For example, Lin et al. [35] reported a sensitivity of 78.5% and specificity of 39.2%, while the equivalent figures for Rodríguez et al. [34] were 46.7% and 89.7%, respectively; and for Stokes et al. [27], they were 80% and 73%, respectively. As with optimal cut-off, we need more data on the factors which are driving this variability in specificity and sensitivity, but the much smaller differences between our two studies suggest that differences in protocols may be driving much of the difference between studies.
For alternative methods of lameness detection, it is the specificity that is the most important. If a technique is easy to apply, then moderate sensitivity can be overcome by repeated measurement. However, even relatively high specificity (>90%) can result in farmer fatigue if a high proportion of identified cows are not actually lame [36]. This fatigue is determined by the positive predictive value (PPV), the proportion of the positive tests which are true positive. The PPV is determined by specificity and prevalence; thus, although the specificity was lower in our Tanzanian study than in New Zealand, the much higher lameness prevalence (33 vs. 14%) meant that PPV was higher in the Tanzanian study.
Nevertheless, a PPV of 72% at the optimal cut-off of 38.0 °C is not ideal as ~1/4 of the cattle identified will not have changes in gait and posture associated with clinical lameness. However, such issues may be tolerated in the small dairy herds in Tanzania, where the number of false positive cattle identified at any one time may be relatively small. Further research is required to better understand Tanzanian dairy farmers’ approach to lameness and whether this technological approach is likely to be acceptable on Tanzanian dairy farms.
One fundamental limitation of the current study is equipment cost. The FLIR camera is not likely to be affordable for most farmers, especially in lower-middle-income countries like Tanzania. However, the costs of infrared cameras are decreasing, and the development of smartphone apps which utilise IRT technology may make the technology affordable, even for Tanzanian farmers. Thus, our data suggest that if these changes continue, IRT may, in the future, be a useful objective method of detecting lameness in Tanzanian dairy cows.

5. Conclusions

The current study showed that, with a trained observer, locomotion scoring after milking can be a useful method of lameness detection under Tanzania conditions. However, the study farms were chosen based on having a suitable area for observation, so may not be typical of Tanzanian dairy farms. Further evaluation of LS on Tanzanian dairy farms is warranted, but the requirement for a trained observer may remain a key issue. The present study results demonstrated that FST measured by IRT was able to differentiate between cows with different locomotion scores. However, the sensitivity, and especially the specificity of IRT need to be improved before it can be recommended for lameness detection under Tanzanian conditions. Additionally, the development of cheap, accurate smartphone IRT apps is needed for such technology to become affordable even on large Tanzanian dairy farms.

Author Contributions

C.W.W. and the supervisory team conceptualised the research. C.W.W. collected data and did data analysis. R.A.L. validated data analysis. C.W.W. wrote the first draft of the manuscript, which was contributed to and approved by all authors. L.J.L., K.R.M. and R.A.L. supervised the project. R.A.L. was responsible for funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

Massey University Foundation: Massey University Postgraduate Student Research Fund.

Institutional Review Board Statement

The observation and image capture described in this study does not meet the definition of manipulation in the Tanzania Animal Welfare Act 2008. Therefore, ethical approval for animal manipulations was not required.

Informed Consent Statement

Farmers were clients of the Sokoine University of Agriculture (SUA)—Animal Teaching Hospital, and were interested in participating when informed about this project. We did not get specific written consent from the farmers as they were clients of the SUA—Animal Teaching Hospital, and were happy to have the researchers on their farms.

Data Availability Statement

Data are available at request from the corresponding author.

Acknowledgments

The authors would like to thank the farmers and staff involved. Similarly, C.W.W. acknowledges the financial support from Richard Laven and the School of Veterinary Science, Massey University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leach, K.A.; Tisdall, D.A.; Bell, N.J.; Main, D.C.J.; Green, L.E. The effects of early treatment for hindlimb lameness in dairy cows on four commercial UK farms. Vet. J. 2012, 193, 626–632. [Google Scholar] [CrossRef] [PubMed]
  2. Ghotoorlar, S.M.; Ghamsari, S.M.; Nowrouzian, I.; Ghidary, S.S. Lameness scoring system for dairy cows using force plates and artificial intelligence. Vet. Rec. 2012, 170, 126. [Google Scholar] [CrossRef] [PubMed]
  3. Pedersen, S.; Wilson, J. Early detection and prompt effective treatment of lameness in dairy cattle. Livestock 2021, 26, 115–121. [Google Scholar] [CrossRef]
  4. Newsome, R.; Green, M.J.; Bell, N.J.; Chagunda, M.G.G.; Mason, C.S.; Rutland, C.S.; Sturrock, C.J.; Whay, H.R.; Huxley, J.N. Linking bone development on the caudal aspect of the distal phalanx with lameness during life. J. Dairy Sci. 2016, 99, 4512–4525. [Google Scholar] [CrossRef] [PubMed]
  5. Schlageter-Tello, A.A.; Bokkers, E.A.; Koerkamp, P.W.G.; Van Hertem, T.; Viazzi, S.; Romanini, C.E.; Halachmi, I.; Bahr, C.; Berckmans, D.; Lokhorst, K. Manual and automatic locomotion scoring systems in dairy cows: A review. Prev. Vet. Med. 2014, 116, 12–25. [Google Scholar] [CrossRef]
  6. Ramanoon, S.Z.; Sadiq, M.B.; Shaik Mossadeq, W.; Mansor, R.; Syed-Hussain, S. The impact of lameness on dairy cattle welfare: Growing need for objective methods of detecting lame cows and assessment of associated pain. In Animal Welfare; IntechOpen: London, UK, 2018; pp. 51–72. [Google Scholar]
  7. Van Nuffel, A.; Zwertvaegher, I.; Van Weyenberg, S.; Pluym, L.; Pastell, M.; Thorup, V.M.; Sonck, B.; Saeys, W. Lameness detection in dairy cows: Part 1. How to distinguish between non-lame and lame cows based on differences in locomotion or behavior. Animals 2015, 5, 838–860. [Google Scholar] [CrossRef]
  8. Schlageter-Tello, A.A.; Bokkers, E.A.; Koerkamp, P.W.G.; Van Hertem, T.; Viazzi, S.; Romanini, C.E.; Halachmi, I.; Bahr, C.; Berckmans, D.; Lokhorst, K. Comparison of locomotion scoring for dairy cows by experienced and inexperienced raters using live or video observation methods. Anim. Welf. 2015, 24, 69–79. [Google Scholar] [CrossRef]
  9. Schlageter-Tello, A.; Van Hertem, T.; Bokkers, E.A.M.; Viazzi, S.; Bahr, C.; Lokhorst, K. Performance of human observers and an automatic 3-dimensional computer-vision-based locomotion scoring method to detect lameness and hoof lesions in dairy cows. J. Dairy Sci. 2018, 101, 6322–6335. [Google Scholar] [CrossRef]
  10. Eddy, A.L.; Van Hoogmoed, L.M.; Snyder, J.R. The role of thermography in the management of equine lameness. Vet. J. 2001, 162, 172–181. [Google Scholar] [CrossRef]
  11. Berry, R.J.; Kennedy, A.D.; Scott, S.L.; Kyle, B.L.; Schaefer, A.L. Daily variation in the udder surface temperature of dairy cows measured by infrared thermography: Potential for mastitis detection. Can. J. Anim. Sci. 2003, 83, 687–693. [Google Scholar] [CrossRef]
  12. Rainwater-Lovett, K.; Pacheco, J.M.; Packer, C.; Rodriguez, L.L. Detection of foot-and-mouth disease virus infected cattle using infrared thermography. Vet. J. 2009, 180, 317–324. [Google Scholar] [CrossRef]
  13. Alsaaod, M.; Syring, C.; Dietrich, J.; Doherr, M.G.; Gujan, T.; Steiner, A. A field trial of infrared thermography as a non-invasive diagnostic tool for early detection of digital dermatitis in dairy cows. Vet. J. 2014, 199, 281–285. [Google Scholar] [CrossRef]
  14. Harris-Bridge, G.; Young, L.; Handel, I.; Farish, M.; Mason, C.; Mitchell, M.A.; Haskell, M.J. The use of infrared thermography for detecting digital dermatitis in dairy cattle: What is the best measure of temperature and foot location to use? Vet. J. 2018, 237, 26–33. [Google Scholar] [CrossRef]
  15. Fabbri, G.; Fiore, E.; Piccione, G.; Giudice, E.; Gianesella, M.; Morgante, M.; Armato, L.; Bonato, O.; Giambelluca, S.; Arfuso, F. Detection of digital and interdigital dermatitis in Holstein Friesian dairy cows by means of infrared thermography. Large Anim. Rev. 2020, 26, 113–116. [Google Scholar]
  16. Alsaaod, M.; Büscher, W. Detection of hoof lesions using digital infrared thermography in dairy cows. J. Dairy Sci. 2012, 95, 735–742. [Google Scholar] [CrossRef]
  17. Wilhelm, K.; Wilhelm, J.; Fürll, M. Use of thermography to monitor sole haemorrhages and temperature distribution over the claws of dairy cattle. Vet. Rec. 2015, 176, 146. [Google Scholar] [CrossRef] [PubMed]
  18. Orman, A.; Endres, M.I. Use of thermal imaging for identification of foot lesions in dairy cattle. Acta Agric. Scand. Sect. A—Anim. Sci. 2016, 66, 1–7. [Google Scholar] [CrossRef]
  19. Nikkhah, A.; Plaizier, J.C.; Einarson, M.S.; Berry, R.J.; Scott, S.L.; Kennedy, A.D. Infrared thermography and visual examination of hooves of dairy cows in two stages of lactation. J. Dairy Sci. 2005, 88, 2749–2753. [Google Scholar] [CrossRef]
  20. Bobić, T.; Mijić, P.; Gantner, V.; Glavaš, H.; Gregić, M. The effects of parity and stage of lactation on hoof temperature of dairy cows using a thermovision camera. J. Cent. Eur. Agric. 2018, 19, 777–783. [Google Scholar] [CrossRef]
  21. Main, D.C.; Stokes, J.E.; Reader, J.D.; Whay, H.R. Detecting hoof lesions in dairy cattle using a hand-held thermometer. Vet. Rec. J. Br. Vet. Assoc. 2012, 171, 504. [Google Scholar] [CrossRef]
  22. Wang, F.-K.; Shih, J.-Y.; Juan, P.-H.; Su, Y.-C.; Wang, Y.-C. Non-invasive cattle body temperature measurement using infrared thermography and auxiliary sensors. Sensors 2021, 21, 2425. [Google Scholar] [CrossRef]
  23. Gloster, J.; Ebert, K.; Gubbins, S.; Bashiruddin, J.; Paton, D.J. Normal variation in thermal radiated temperature in cattle: Implications for foot-and-mouth disease detection. BMC Vet. Res. 2011, 7, 73. [Google Scholar] [CrossRef]
  24. Werema, C.W.; Laven, L.; Mueller, K.; Laven, R. Evaluating Alternatives to Locomotion Scoring for Lameness Detection in Pasture-Based Dairy Cows in New Zealand: Infra-Red Thermography. Animals 2021, 11, 3473. [Google Scholar] [CrossRef]
  25. DairyNZ. DairyNZ Lameness Scoring. Available online: https://www.dairynz.co.nz/animal/cow-health/lameness/lameness-scoring/ (accessed on 5 April 2023).
  26. AHDB. AHDB Dairy UK Mobility with AHDB Dairy. Available online: https://ahdb.org.uk/knowledge-library/mobilityscoring-for-dairy-cows (accessed on 5 April 2023).
  27. Stokes, J.E.; Leach, K.A.; Main, D.C.J.; Whay, H.R. An investigation into the use of infrared thermography (IRT) as a rapid diagnostic tool for foot lesions in dairy cattle. Vet. J. 2012, 193, 674–678. [Google Scholar] [CrossRef] [PubMed]
  28. Whay, H.; Bell, M.; Main, D. Validation of lame limb identification through thermal imaging. In Proceedings of the 13th International Symposium and 5th Conference on Lameness in Ruminants, Maribor, Slovenia, 11–15 February 2004. [Google Scholar]
  29. Abdi, H. Bonferroni and Šidák corrections for multiple comparisons. Encycl. Meas. Stat. 2007, 3, 103–107. [Google Scholar]
  30. Astivia, O.L.O.; Zumbo, B.D. Heteroskedasticity in Multiple Regression Analysis: What it is, How to Detect it and How to Solve it with Applications in R and SPSS. Pract. Assess. Res. Eval. 2019, 24, 1. [Google Scholar]
  31. Murray, R.; Downham, D.; Clarkson, M.; Faull, W.; Hughes, J.; Manson, F.; Merritt, J.; Russell, W.; Sutherst, J.; Ward, W. Epidemiology of lamenessin dairy cattle: Description and analysis of foot lesions. Vet. Rec. 1996, 138, 586–591. [Google Scholar] [CrossRef]
  32. Chesterton, R.; Lawrence, K.; Laven, R. A descriptive analysis of the foot lesions identified during veterinary treatment for lameness on dairy farms in north Taranaki. N. Z. Vet. J. 2008, 56, 130–138. [Google Scholar] [CrossRef]
  33. Fjeldaas, T.; Knappe-Poindecker, M.; Bøe, K.; Larssen, R. Water footbath, automatic flushing, and disinfection to improve the health of bovine feet. J. Dairy Sci. 2014, 97, 2835–2846. [Google Scholar] [CrossRef] [PubMed]
  34. Rodríguez, A.R.; Olivares, F.J.; Descouvieres, P.T.; Werner, M.P.; Tadich, N.A.; Bustamante, H.A. Thermographic assessment of hoof temperature in dairy cows with different mobility scores. Livest. Sci. 2016, 184, 92–96. [Google Scholar] [CrossRef]
  35. Lin, Y.C.; Mullan, S.; Main, D.C. Optimising lameness detection in dairy cattle by using handheld infrared thermometers. Vet. Med. Sci. 2018, 4, 218–226. [Google Scholar] [CrossRef] [PubMed]
  36. O’Leary, N.; Byrne, D.; O’Connor, A.; Shalloo, L. Invited review: Cattle lameness detection with accelerometers. J. Dairy Sci. 2020, 103, 3895–3911. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Infrared thermography image of the plantar aspect of the left hind foot overlaid to depict the seven zones for which estimates of surface skin temperature were measured. On the lateral claw; zone 1: coronary band (CB), zone 2: above the coronary band (ACB), zone 3: below accessory digit (BAD), zone 4: interdigital space (IDS). Zones 5 to 7 are the equivalent to zones 1 to 3 but on the medial claw.
Figure 1. Infrared thermography image of the plantar aspect of the left hind foot overlaid to depict the seven zones for which estimates of surface skin temperature were measured. On the lateral claw; zone 1: coronary band (CB), zone 2: above the coronary band (ACB), zone 3: below accessory digit (BAD), zone 4: interdigital space (IDS). Zones 5 to 7 are the equivalent to zones 1 to 3 but on the medial claw.
Animals 13 01372 g001
Figure 2. Receiver operating characteristic (ROC) curves were used to determine the optimal threshold values for the infrared thermography’s sensitivity and specificity, assuming locomotion scores ≥ 2 as locomotion score = 2 on three dairy farms in Morogoro (n = 170). Note: coronary band (CB), above the coronary band (ACB), below accessory digit (BAD).
Figure 2. Receiver operating characteristic (ROC) curves were used to determine the optimal threshold values for the infrared thermography’s sensitivity and specificity, assuming locomotion scores ≥ 2 as locomotion score = 2 on three dairy farms in Morogoro (n = 170). Note: coronary band (CB), above the coronary band (ACB), below accessory digit (BAD).
Animals 13 01372 g002
Table 1. Definitions of the infrared thermography estimates utilised for analysis adopted from [24].
Table 1. Definitions of the infrared thermography estimates utilised for analysis adopted from [24].
Foot/Zone *Description
Mean temperatureAverage temperature, across both feet (all 14 zones)
Hottest zoneHighest zone temperature, across both feet (all 14 zones)
Hottest zone 4The highest zone 4 temperature on either foot
Hottest coronary band (CB)The highest zone 1 or 5 temperature on either foot
Hottest above the coronary band zone (ACB)The highest zone 2 or 6 temperature on either foot
Hottest zone below the accessory digit (BAD)The highest zone 3 or 7 temperature on either foot
* For all analyses, the maximum temperature for each zone was used for the analysis.
Table 2. Distribution of locomotion score and their percentage in brackets on three farms in Morogoro, Tanzania (n = 170).
Table 2. Distribution of locomotion score and their percentage in brackets on three farms in Morogoro, Tanzania (n = 170).
Farm NameLocomotion Score
0123Total
111 (21.2)24 (46.1)12 (23.1)5 (9.6)52
215 (25.9)25 (43.1)15 (25.9)3 (5.1)58
313 (21.7)26 (43.3)14 (23.3)7 (1.7)60
Total39 (23.0)75 (44.1)41 (24.1)15 (8.8)170
Table 3. The mean temperature (degrees centigrade) and 95% confidence interval for zones 1–7 (see Table 1 for definitions) for both hindlimbs on three dairy farms in Morogoro (n = 170).
Table 3. The mean temperature (degrees centigrade) and 95% confidence interval for zones 1–7 (see Table 1 for definitions) for both hindlimbs on three dairy farms in Morogoro (n = 170).
ZoneMean95% Confidence Interval
LowerUpper
137.8237.7837.87
237.0236.9637.09
337.4637.4537.47
438.3138.2938.32
537.5837.5237.65
636.7436.7036.77
737.3037.2637.33
Note: zones 1–3 are on the lateral claw, and 5–7 are equivalents on the medial claw.
Table 4. Comparisons of mean foot skin temperature (degrees centigrade) and 95% confidence interval for zones 1–7 (see Table 1 for definitions) for both hindlimbs on three dairy farms in Morogoro (n = 170).
Table 4. Comparisons of mean foot skin temperature (degrees centigrade) and 95% confidence interval for zones 1–7 (see Table 1 for definitions) for both hindlimbs on three dairy farms in Morogoro (n = 170).
(I) Zone(J) ZoneMean Difference (I − J)95% Confidence Interval
Lower BoundUpper Bound
120.800.531.08
30.360.090.64
50.24−0.030.52
61.090.811.36
70.530.250.80
260.280.010.56
320.440.170.72
60.720.451.00
70.16−0.110.44
410.480.210.76
21.291.011.56
30.850.571.12
50.730.451.00
61.571.301.85
71.010.741.30
520.560.290.83
30.12−0.160.40
60.840.571.12
70.280.010.56
720.280.000.55
60.560.290.84
Table 5. The mean temperature (degrees centigrade) from plantar aspects across both feet (all 14 zones) and 95% confidence interval for different locomotion scores (0–3) DairyNZ on three farms in Morogoro (n = 170).
Table 5. The mean temperature (degrees centigrade) from plantar aspects across both feet (all 14 zones) and 95% confidence interval for different locomotion scores (0–3) DairyNZ on three farms in Morogoro (n = 170).
Locomotion ScoreMean95% Confidence Interval
LowerUpper
036.4536.1236.76
137.3537.1937.51
238.1938.0038.38
338.8138.6338.99
Table 6. Comparisons of mean foot skin temperature for different zones versus locomotion scores (0–3) DairyNZ on three dairy farms in Morogoro (n = 170).
Table 6. Comparisons of mean foot skin temperature for different zones versus locomotion scores (0–3) DairyNZ on three dairy farms in Morogoro (n = 170).
Temperature Measure(I) LS(J) LSMean Difference (I − J)95% Confidence Interval
Lower BoundUpper Bound
MT100.910.501.31
201.741.282.20
10.830.431.23
302.391.763.01
11.480.902.06
20.650.031.27
Hottest zone100.830.421.24
201.561.102.03
10.730.331.13
302.261.632.90
11.430.842.02
20.700.071.33
Hottest zone 4100.890.441.35
201.571.062.09
10.680.231.12
302.361.663.06
11.470.822.11
20.790.091.48
Hottest CB100.810.371.25
201.601.102.09
10.790.361.22
302.331.653.00
11.520.892.14
20.730.061.40
Hottest zone ACB100.820.381.26
201.681.182.19
10.860.431.30
302.261.582.94
11.440.802.07
20.57−0.101.25
Hottest zone BAD100.870.471.27
201.641.182.10
10.770.371.17
302.361.742.98
11.490.912.07
20.720.111.34
MT = mean temperature, CB = coronary band, ACB = above the coronary band, BAD = below the accessory digits.
Table 7. Temperature measures estimates and their 95% confidence intervals on three dairy farms in Morogoro (n = 170).
Table 7. Temperature measures estimates and their 95% confidence intervals on three dairy farms in Morogoro (n = 170).
Model Parameter 95% Confidence Interval
Lower BoundUpper Bound
Mean temperature (intercept)38.6638.3938.93
Locomotion score *0.570.460.70
Hottest zone (intercept)40.0639.8140.31
Locomotion score *0.570.460.70
Hottest zone 4 (intercept)39.8139.5640.07
Locomotion score *0.700.570.87
Hottest coronary band (intercept)39.4739.1839.75
Locomotion score *0.660.530.81
Hottest zone ACB (intercept)38.4738.1838.76
Locomotion score *0.670.540.83
Hottest zone BAD (intercept)39.1338.9039.37
Locomotion score *0.560.460.70
* Effect of increased locomotion score of 1 unit in DairyNZ lameness score. See Table 1 for the definition of temperature measurement. ACB = above the coronary band, BAD = below the accessory digits.
Table 8. Optimal cut-off points for skin foot temperature measurements (degrees centigrade) for determining lame cows (Dairy NZ lameness score ≥ 2) on three dairy farms in Morogoro (n = 170).
Table 8. Optimal cut-off points for skin foot temperature measurements (degrees centigrade) for determining lame cows (Dairy NZ lameness score ≥ 2) on three dairy farms in Morogoro (n = 170).
Temperature Measure 1Optimal Threshold (°C)AUC
(95% CI)
Specificity
(95% CI)
Sensitivity
(95% CI)
PPV
(95% CI) *
NPV
(95% CI) *
Mean temperature38.00.88 (0.83–0.93)86.0 (78.2–91.8)73.2 (59.7–84.2)71.9 (61.3–80.6)86.7 (80.8–91.0)
Hottest zone39.00.85 (0.79–0.91)80.7 (72.3–87.5)78.6 (65.6–88.4)66.7 (57.3–74.9)88.5 (82.2–92.7)
Hottest zone 438.90.84 (0.78–0.90)79.8 (71.3–86.8)71.4 (57.8–82.7)63.5 (53.8–72.2)85.1 (78.8–89.7)
Hottest CB38.60.85 (0.79–0.91)82.5 (74.2–88.9)75.0 (61.6–85.6)67.7 (57.8–76.3)87.0 (80.9–91.4)
Hottest zone ACB38.10.87 (0.81–0.92)86.8 (79.2–92.4)69.6 (55.9–81.2)72.2 (61.1–81.1)85.3 (79.6–89.7)
Hottest zone BAD38.70.86 (0.80–0.92)93.9 (87.8–97.5)64.3 (50.4–76.6)83.7 (71.0–91.5)84.3 (79.0–88.4)
1 See Table 1 for the definition of temperature measure; AUC, area under the curve; CI, confidence interval; PPV, positive predictive value and NPV, negative predictive value; CB, coronary band; ACB, above the coronary band; BAD, below accessory digit (BAD). *, PPV and NPV calculated at a prevalence of lameness of 32.94%.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Werema, C.W.; Laven, L.J.; Mueller, K.R.; Laven, R.A. Assessing Alternatives to Locomotion Scoring for Detecting Lameness in Dairy Cattle in Tanzania: Infrared Thermography. Animals 2023, 13, 1372. https://doi.org/10.3390/ani13081372

AMA Style

Werema CW, Laven LJ, Mueller KR, Laven RA. Assessing Alternatives to Locomotion Scoring for Detecting Lameness in Dairy Cattle in Tanzania: Infrared Thermography. Animals. 2023; 13(8):1372. https://doi.org/10.3390/ani13081372

Chicago/Turabian Style

Werema, Chacha W., Linda J. Laven, Kristina R. Mueller, and Richard A. Laven. 2023. "Assessing Alternatives to Locomotion Scoring for Detecting Lameness in Dairy Cattle in Tanzania: Infrared Thermography" Animals 13, no. 8: 1372. https://doi.org/10.3390/ani13081372

APA Style

Werema, C. W., Laven, L. J., Mueller, K. R., & Laven, R. A. (2023). Assessing Alternatives to Locomotion Scoring for Detecting Lameness in Dairy Cattle in Tanzania: Infrared Thermography. Animals, 13(8), 1372. https://doi.org/10.3390/ani13081372

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop