If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Lameness severely impairs cattle's locomotion, and it is among the most important threats to animal welfare, performance, and productivity in the modern dairy industry. However, insight into the pathological alterations of claw biomechanics leading to lameness and an understanding of the biomechanics behind development of claw lesions causing lameness are limited. Biplane high-speed fluoroscopic kinematography is a new approach for the analysis of skeletal motion. Biplane high-speed videos in combination with bone scans can be used for 3-dimensional (3D) animations of bones moving in 3D space. The gold standard, marker-based animation, requires implantation of radio-opaque markers into bones, which impairs the practicability for lameness research in live animals. Therefore, the purpose of this study was to evaluate the comparative accuracy of 2 noninvasive, markerless animation techniques (semi-automatic and manual) in 3D animation of the bovine distal limb. Tantalum markers were implanted into each of the distal, middle, and proximal phalanges of 5 isolated bovine distal forelimbs, and biplane high-speed x-ray videos of each limb were recorded to capture the simulation of one step. The limbs were scanned by computed tomography to create bone models of the 6 digital bones, and 3D animation of the bones' movements were subsequently reconstructed using the marker-based, the semi-automatic, and the manual animation techniques. Manual animation translational bias and precision varied from 0.63 ± 0.26 mm to 0.80 ± 0.49 mm, and rotational bias and precision ranged from 2.41 ± 1.43° to 6.75 ± 4.67°. Semi-automatic translational values for bias and precision ranged from 1.26 ± 1.28 mm to 2.75 ± 2.17 mm, and rotational values varied from 3.81 ± 2.78° to 11.7 ± 8.11°. In our study, we demonstrated the successful application of biplane high-speed fluoroscopic kinematography to gait analysis of bovine distal limb. Using the manual animation technique, kinematics can be measured with sub-millimeter accuracy without the need for invasive marker implantation.
Scientific analysis of the underlying causes, particularly a biomechanical understanding of claw–floor interactions, is key for successful reduction of claw tissue damage and the detection and prevention of lameness (
Röntgenanatomische und querschnittsanatomische Untersuchungen unter Berücksichtigung magnetresonanztomographischer Befunde an der Hintergliedmaße des Rindes. Faculty of Veterinary Medicine,
). Two x-ray generators in combination with 2 fluoroscopic image intensifiers capture high-speed video sequences (Figure 1). These sequences serve as templates for 3D bone animations of bone reconstructions and enable direct and precise measurements of bone elements in live animals during locomotion (
; Figure 2, top). It requires surgical implantation of at least 3 radio-opaque markers into each bone of interest to facilitate 3D animations. Although marker-based registration offers high accuracy (0.12 ± 0.08 mm;
), noninvasive markerless techniques ought to be applied to live animals for animal welfare reasons and to enable evaluation of regions that are surgically difficult to access. For noninvasive animation, 2 animation techniques are available: a semi-automatic (Figure 2, middle;
). Applied to the bovine distal limb, HFK may allow visualization and measurement of pathological alterations of claw biomechanics. These new insights may improve the prevention of lameness; for example, by investigating claw–floor interactions on various types of flooring. They may also enable the evaluation of influences of claw-trimming methods on claw biomechanics. Digital motion as well as the influence of ground conditions may be measured precisely.
Hence, the first aim of this study was to establish the application of biplane HFK to bovine distal limb and to provide a proof of concept in describing 3D bone movement reliably and precisely. The second aim was to evaluate 2 noninvasive markerless animation techniques (semi-automatic and manual technique) with regard to future applicability in live dairy cows. We hypothesize that (1) the accuracies of the two markerless animation techniques do not differ significantly; (2) the accuracy measurements of left and right phalanges do not differ significantly; and (3) both noninvasive, markerless animation techniques can match all 6 phalanges throughout simulated steps.
MATERIALS AND METHODS
The phalanges of 5 isolated distal limbs of German Holstein dairy cows were analyzed with 3 different animation techniques with regard to bias and precision (Figure 3).
Three left and 2 right bovine distal forelimbs, separated just below the metacarpal joint, were obtained from a local abattoir. Only limbs that were externally intact and showed no pathologies in computed tomography (CT) scans were included in this study. The digits were classified according to their proximity to the image intensifiers and cameras. The phalanges were defined as left (L; far from the camera) and right (R; closer to the camera). Claw length was measured according to
(mean left claw length: 116.02 mm; mean right claw length: 119.48 mm). Proximal and middle phalanges were surgically implanted with 3 spherical tantalum markers of 1-mm diameter (X-medics Scandinavia Smba, Frederiksberg, Denmark), and distal phalanges were implanted with 4. Markers were inserted through 2-mm venous catheters (Braun, Melsungen, Germany) into canals drilled through skin incisions into the bone or into the horn capsule with a drill of 2-mm diameter and held in place with cyanoacrylate glue (TEDI GmbH & Co. KG, Dortmund, Germany). The markers were evenly spread across each bone (Figure 4).
Computed Tomography and Bone Models
The CT scans of the 5 bovine limbs were performed on a multi-slice helical CT-scanner (Philips Healthcare, DA Best, the Netherlands) with a median in-plane resolution of 0.475 mm, a slice thickness of 1 mm, and an overlapping increment of 0.5 mm. Model preparation was done using a sharp bone filter (Filter Type and Convolution Kernel D). The distal, middle, and proximal phalanges were manually separated from each other for bone modeling using the image processing software MeVisLab version 2.3.1 (MeVis Medical Solutions AG, Bremen, Germany). The separation allowed generation of distinct 3D surface renderings of isolated bone segments using dedicated image analysis software (VTK 3.0, Kitware Inc., Clifton Park, NY). Further enhancement, as well as inspection of the 3D models was performed if necessary (ParaView 4.1.0, Kitware Inc.). Any trace of the markers was digitally removed within the CT data to avoid bias in the markerless animations.
Biplane High-Speed Fluoroscopy
Two x-ray tubes (Philips Medio 65 CP-H X-Ray Generator) and 2 fluoroscopes (Philips Typ BX 3i-2123) with image intensifiers and high-speed cameras (Optronis Cam Record CR600x2, Kehl, Germany) captured video sequences (Figure 1, Supplemental Video V1; https://doi.org/10.3168/jds.2017-12563) of the bones during the simulation of one step (heel strike, midstance, push off). Each limb was moved manually in the trial field of the biplane fluoroscopic setup in a step-like motion. The direction and walking velocity as they would occur in in vivo experiments were re-created. The system operated in continuous mode at 77 to 81 kV and 100 mA with an interbeam angle of 60°. Video sequences with an average length of 2.4 s were taken at 500 frames/s, 0.5 ms shutter, and a resolution of 1,024 × 1,280 pixels. The x-ray source-to-image distance was 1.4 m (Figure 1). This setup was selected with regard to applicability to live dairy cows, providing sufficient space for an in vivo walkway. Image distortion was corrected using a perforated steel sheet with a defined hole-diameter in XMALab version 1.2.17 (
). For manual and semi-automatic animation, bone contrast on each image was enhanced using Adobe Photoshop CS2 version 9.0 Smart Sharpen filter (Adobe Systems Inc., San Jose, CA).
Marker-Based 3D Animation
Marker-based registration was done with XMALab (Figure 2, top). The software measured marker trajectories representing motion of the implanted markers over time by tracking the markers in both x-ray videos. Together with marker coordinates extracted from CT data, translations and rotations could then be calculated for every bone (
; Figure 2, middle). After manual adjustment of the DRR to the x-ray images to obtain a good visual match for both biplane views, the software automatically completes the matching of the remaining video based on 3D bone shape and texture. The position and orientation of a bone is estimated by maximizing the match between the DRR and the biplane x-ray images. Using this technique, the 3D position and orientation of each phalanx were determined independently for all frames of each trial. After completion, the accordance of the DRR with every 10th x-ray image was visually verified. If a mismatch was detected, matching started again from the last frame of accordance. If the mismatch persisted, manual correction was repeated up to 3 times. If unsuccessful, the mismatch was carried through to evaluation. Semi-automatic technique data were used for 3D animations in Autodesk Maya as already described (Supplemental Video V3; https://doi.org/10.3168/jds.2017-12563).
). In this hierarchy, the middle phalanges were placed as the “lead” bone and the distal as well as the proximal phalanges attached. So, when registering the middle phalanges' congruence with its fluoroscopic image, the other 2 bone models followed, resulting in pre-alignment. Once the middle phalanx was positioned correctly, the proximal and the distal phalanges could be matched accurately. The intervals between matched frames were continuously reduced to every eighth frame of the video sequence.
Analysis of 3D Animation
To compare the 3 different tracking techniques in terms of 3D bone movement, coordinate systems were applied to the bones. Each bone's center point served as its origin, while the bones' axes of inertia formed the coordinate system's axes (
For the evaluation of the 2 markerless animation techniques, 5 bovine distal forelimbs were used. The data of corresponding bones (left proximal phalanges, right proximal phalanges, left middle phalanges, right middle phalanges, left distal phalanges, right distal phalanges) of all 5 legs were pooled and analyzed together. Only image sequences in which all 6 bones were visible at the same time and in which all 3 tracking techniques (both markerless techniques as well as the marker-based technique) accomplished a match were analyzed. For each bone, a different number of frames was evaluated (left proximal phalanges: 3,009; right proximal phalanges: 2,670; left middle phalanges: 2,816; right middle phalanges: 2,642; left distal phalanges: 1,977; right distal phalanges: 1,033).
For further statistical analysis, 2 right distal phalanges (3 and 5) were removed. During semi-automatic matching of the right distal phalanges of claw 3, the program stopped immediately after starting and achieved no match between the DRR and the x-ray image. For the right distal phalanges of claw 5, the semi-automatic technique worked successfully for one video. On the other camera, the DRR were arranged between the left and right distal phalanges (Figure 5).
To evaluate the accuracy of semi-automatic and manual animation, the translational and rotational data were analyzed in relation to marker-based registration. The translational accuracy was defined as the Euclidean distance between a bone's center point tracked with the semi-automatic and marker-based technique or the manual and marker-based technique in one frame. Euclidean differences were defined as
where NI was the 3D vector of the markerless tracking method and MB the vector of the marker-based tracking method. Accordingly, rotational accuracy was calculated using the same equation. Changes of the x-angle represented motion in the frontal plane, y-angle in the sagittal plane, and z-angle in the transverse plane. Accuracy of manual and semi-automatic techniques was quantified in terms of bias and precision (
). Bias was defined as the average Euclidean difference of one bone type (left proximal phalanges, right proximal phalanges, left middle phalanges, right middle phalanges, left distal phalanges, right distal phalanges) over the entire image sequences of all legs in the semi-automatic and marker-based or the manual and marker-based technique. Precision was the standard deviation of the bias. All values were assessed unfiltered and filtered at 20 Hz using a fifth-order, low-pass Butterworth filter. For both noninvasive animation techniques, the distribution and dimension of translational and rotational tracking errors over the different phalanges are illustrated in Figure 6, Figure 7 with violin plots including box plots with medians, 25th and 75th percentiles, and whiskers that extend to ±1.5 interquartile range. Additionally, mean values of standard deviations of intermarker distances per bone were evaluated as a measure of marker-tracking accuracy (
Statistical analysis of the data was accomplished with R version 3.2.4 (The R Foundation for Statistical Computing, Vienna, Austria). The Shapiro-Wilk normality test was used to determine whether values followed a normal distribution. Noninvasive animation technique values were compared with marker-based tracking using Wilcoxon's rank sum test for one sample, P < 0.01. Comparisons of manual animation values with semi-automatic technique values and left with right phalanges were performed using the Mann-Whitney U-test, P < 0.01.
We successfully applied HFK to the bovine distal limb. We performed 3D animations of the proximal, middle, and distal phalanges of all 5 distal limbs with marker-based and 2 markerless animation techniques (semi-automatic and manual; Figure 2). We had to exclude 2 right distal phalanges from statistical analysis. The semi-automatic technique achieved no match between the DRR of the right distal phalanx of claw 3 and the x-ray image. The semi-automatic technique's misalignment of the right distal phalanx of claw 5 resulted in a median deviation of its translational values of at least 12.4 mm compared with the other right distal phalanges (Figure 5), and it was therefore excluded.
As a measure of marker-tracking accuracy, we recorded standard deviations of intermarker distances of 0.04 mm (min 0.005 mm; max 0.14 mm) for all bones and limbs.
In translational and rotational analysis, low values indicate a high tracking accuracy and high values a low tracking accuracy. The Shapiro-Wilk normality test resulted in 0.00 for every group, showing that values did not follow a normal distribution.
The distribution and amount of unfiltered translational tracking errors of manual and semi-automatic techniques are shown in Figure 6 with violin plots in combination with box plots. Unfiltered translational bias and precision as well as medians and 95% confidence intervals (CI) for semi-automatic and manual animation are given in Table 1. Unfiltered manual animation bias and precision varied from 0.63 ± 0.26 mm (left distal phalanges) to 0.80 ± 0.49 mm (right distal phalanges). Manual animations' lowest medians occurred with 0.58 mm for right middle phalanges and the highest with 0.71 mm for right distal phalanges. Unfiltered semi-automatic technique bias and precision ranged from 1.26 ± 1.28 mm (right middle phalanges) to 2.75 ± 2.17 mm (left distal phalanges). Semi-automatic techniques' lowest medians arose with 0.76 mm for right middle phalanges and the highest with 2.09 mm for left distal phalanges. For both noninvasive animation techniques, medians were lower in comparison with bias.
Table 1Unfiltered translational bias and precision values with median and 95% CI for the different bone types of both markerless animation techniques
In a comparison of the lowest 95% CI value for the semi-automatic technique with the highest 95% CI value for the manual animation (mindiff) and vice versa (maxdiff), the magnitude of the difference between both markerless techniques ranged from 1.5 up to 4.1 in favor of manual animation. On average, manual animation showed medians that were 2.1 times lower than those of the semi-automatic technique. Manual animation showed continuously lower values for translation than the semi-automatic technique for all bones (P < 0.01). Noninvasive animation methods were different in translation from the marker-based method (P < 0.01).
In a comparison of the values of left and right phalanges, both markerless techniques showed differences in translation between each left and right phalanx (P < 0.01). For the semi-automatic animation technique, higher errors (median and 95% CI maximum, Table 1) occurred in left phalanges compared with corresponding right phalanges across all bones. The shape of the violin plots of the semi-automatic technique in Figure 6 also suggests a higher probability of higher errors in left phalanges. The manual method did not show this consistent distribution. Higher translational values were generated for right proximal and distal phalanges, whereas values for right middle phalanges were lower than the corresponding left side (Table 1 and Figure 6).
Figure 7 shows the distribution and dimensions of unfiltered rotational tracking errors of manual and semi-automatic animation. Table 2 shows unfiltered rotational bias and precision as well as medians and 95% CI for semi-automatic and manual animation. Unfiltered manual animation bias and precision ranged from 2.41 ± 1.43° (right proximal phalanges) to 6.75 ± 4.67° (left middle phalanges). Manual animation's lowest medians occurred with 2.09° for right proximal phalanges and the highest with 5.63° for the left middle phalanges. Unfiltered semi-automatic technique bias and precision varied from 3.81 ± 2.78° (right distal phalanges) to 11.7 ± 8.11° (left middle phalanges). Semi-automatic techniques' lowest medians were 3.26° for right distal phalanges and the highest were 10.25° for left middle phalanges. Medians were lower in comparison with bias for both noninvasive animation techniques.
Table 2Unfiltered rotational bias and precision values with median and 95% CI for the different bone types of both markerless animation techniques
In a comparison of semi-automatic techniques' lowest 95% CI value with manual animation's highest 95% CI value (mindiff) and vice versa (maxdiff), the accuracy between both markerless techniques ranged from 1.3 to 3.4 times in favor of manual animation. On average, manual animation showed 2.0 times higher rotational accuracy than semi-automatic animation. Manual animation showed consistently lower values for rotation than semi-automatic animation (P < 0.01). Rotational accuracy was lower than translational accuracy for all phalanges in both noninvasive animation techniques. Rotations in both markerless tracking techniques were different from the marker-based method (P < 0.01).
The analysis of corresponding bones' tracking accuracy revealed that manual animation showed a difference in rotational accuracy between corresponding middle phalanges (P < 0.01), but not for proximal phalanges (P = 0.87) or distal phalanges (P = 0.04). Semi-automatic technique showed differences in rotation between the left and right side for middle phalanges and distal phalanges (P < 0.01), but not for proximal phalanges (P = 0.03). In accordance with the results for translational values for the semi-automatic animation technique, higher rotational errors occurred in left phalanges compared with corresponding right phalanges across all bones. Manual techniques' lowest medians were also found in left phalanges compared with corresponding right phalanges across all bones (Table 2 and Figure 7).
In this study, we successfully applied biplane high-speed fluoroscopy to the bovine distal limb and evaluated 2 markerless animation techniques (semi-automatic and manual) for noninvasive kinematic measurements of bovine phalanges. To determine bias and precision as well as medians, we used marker-based tracking as the gold standard. As a measure of marker-tracking accuracy, standard deviations of intermarker distances were recorded. We observed a mean of standard deviations of unfiltered intermarker distances lower than the mean values of previous reports: 0.064 and 0.072 mm (filtered) (
). This finding is indicative of high marker-tracking accuracy and a rigid connection between markers and bones. It demonstrates that marker positions within bones did not change during limb motion. Marker-based tracking therefore provided a precise reference for the evaluation of both noninvasive animation techniques.
The manual simulation of the steps in this study was created to visually resemble in vivo steps in radiographic recordings. Care was taken to meet general radiographically visible criteria of a step (i.e., heel strike, midstance, push off) with velocity representative of a walking cow. In addition, superimposition of the phalanges was accurately re-created according to the expected situation in in vivo experiments. However, although vertical and accelerating–braking forces have been generated, they were not comparable to the in vivo situation. Because it was not technically possible to realistically simulate both appearance and forces of a bovine step at the same time, the parameter that was more likely to affect the results of this study (i.e., appearance) was prioritized to meet the required criteria. Therefore, this study provides a proof of concept for the application of noninvasive tracking techniques to the bovine distal limb. Future studies will have to determine their accuracy under real load-bearing conditions and when load bearing shifts during lameness.
For the evaluation of the 2 noninvasive animation methods, we used 5 bovine distal forelimbs (3 left and 2 right limbs). Limbs were not divided into left and right because no morphometrical differences exist between left and right phalanges (
. Besides, for the results of 3D bone model animation, superimposition of the 2 digits in the biplane x-ray videos is very important. Because of this superimposition, validation of HFK if the phalanges were farther or closer to the cameras was of major interest. Therefore, we decided to classify the digits as left (far from the camera) and right (closer to the camera). For statistical analysis we had to exclude 2 right distal phalanges (3 and 5). For the right distal phalanx of claw 3 the semi-automatic technique was not able to accomplish a match between the DRR and the x-ray image at all. For the right distal phalanx of claw 5 the semi-automatic technique arranged the DRR between left and right distal phalanges on one camera, whereas the other view worked satisfactorily. In oblique biplane recording of artiodactyle phalanges, one of the corresponding bones is the anterior and the other is the posterior in one view; in the second view, the perspectives are reversed. The misalignment of the right distal phalanx of claw 5 occurred in the view in which it was the posterior one and therefore more heavily influenced by superimposition, which may explain the semi-automatic technique's difficulties in matching the DRR with the x-ray image. The reason why only 2 of the 5 right distal phalanges were affected by misalignment may be due to the indistinct bone shape of these 2 right distal phalanges. Consequently, semi-automatic technique may have encountered problems with finding and differentiating bone contours. The removal of 2 of the right distal phalanges for statistical computing reduced the analyzed number of frames to 1,033, which was relatively small compared with the other phalanges, but the number was sufficient for statistical analysis.
For manual animation, the bone models were placed in a hierarchy with virtual joints (
). This attachment in an articular chain provided the benefit of transferring the rotations and translations of the “lead bone” (middle phalanges) to all bone models rather than rotoscoping each bone entirely independently. The resulting pre-alignment of the attached bone models saved time, but only if the lead bone was matched with high accuracy itself. We chose the middle phalanges as lead bones because of their position within the central beam. Although low values for middle phalanges translation were recorded, the high values for rotation (up to 2.7 times higher than the other phalanges) necessitate a reconsideration of the hierarchy used. These high rotational values for the middle phalanges may arise from its rotund shape and indistinct contours, which particularly occur at its distal parts due to the horn capsule. These high values may also be connected to the middle phalanges' position as the lead bone because it was the first bone model to be aligned manually to the biplane x-ray images. The bone models of the other phalanges followed according to the hierarchy and had to undergo less manipulation to be matched accurately. The time-saving benefit of a hierarchy is a great advantage, but for further studies we suggest placing the proximal phalanges at the top level. Because of their long shape and defined bone contour, the proximal phalanges' bone models could be aligned more accurately to the biplane x-ray videos, resulting in lower bias values for translation and rotation. With the proximal phalanx as the lead bone, its accurate alignment can be transferred to the middle and distal phalanges.
Although HFK is a complex laboratory-based method, it has the potential to enable new insights in claw biomechanics by direct and precise measurements of bone movement. In the future, HFK may also permit visualizing and measuring digital motion and allow for investigation of claw–floor interactions for various types of flooring or evaluation of the influences of claw-trimming methods. These new insights can help to improve the understanding of the development and prevention of lameness.
The size of the biplane x-ray system's field of view, which was defined by the overlapping region of the crossing x-ray beams (Figure 1), determined the extent to which skeletal motion could be measured. The intersecting volume of the 2 x-ray beams is approximately 6,000 cm3 (
). Therefore, we focused on heel strike, midstance, and push off of the limb. We chose a setup with an interbeam angle of 60° and a source-to-image distance of 1.4 m. This arrangement allowed us to keep a walkway of more than 1-m width clear of equipment for future in vivo application. However, an increased angle of 90° may improve the spatial resolution, and a different arrangement of x-ray tubes and cameras allowing for a clear distinction between the phalanges will probably have a positive effect on the accuracy of noninvasive animation techniques. The selected setup enables the direct applicability of this study to future studies on live animals, but future hardware developments will have to address this issue.
This study presents the first application of HFK to the bovine distal limb. The 2 markerless animation techniques (semi-automatic and manual) for noninvasive kinematic measurements of bovine phalanges were evaluated in relation to marker-based registration. Marker-based registration is the gold standard, with a methodological error of 0.12 ± 0.08 mm and 0.09 ± 0.08° (
). This error has to be added to the translational and rotational bias and precision of both the manual and semi-automatic measurements. The comparison of the reported results to previous studies is difficult because the previous investigations pertained to entirely different anatomical structures in other species. To our knowledge, only 1 reference exists for the accuracy of the manual technique validating markerless animation techniques in equines (
). We found translational accuracy to be in accordance with data of the equine study for proximal and middle phalanges, whereas the distal phalanges showed lower accuracy. Rotational accuracy median values were up to 1.6 times higher for all phalanges. For the semi-automatic technique, we observed translational and rotational bias values and precision of at least 0.45 mm and 2.82° to be higher than those previously published when using a shape-based tracking technique: 0.0 ± 0.56 mm and 0.0 ± 0.61° (
). The paired digits might explain the differences in translational and rotational bias and precision across semi-automatic and manual animation techniques. The superimposition and the resulting overlapping bone contours impaired semi-automatic technique's matching of the DRR to the x-ray images. Other factors may be irregular bone shape or texture that compromised the semi-automatic technique's performance of distinguishing between bone contours. Additionally, the horn capsule superimposes the lower part of the middle and the entire distal phalanges. In relation, impaired semi-automatic matching became particularly apparent in the left distal phalanges, resulting in high bias values.
The accuracies of both markerless animation techniques do differ (P < 0.01). Therefore, our first hypothesis stating that the accuracies of both markerless animation techniques do not differ significantly has to be rejected. Overall, manual animation showed higher translational and rotational tracking accuracy than semi-automatic animation (P < 0.01), which proves manual animation to be the preferable markerless tracking technique for bovine distal limbs. The differences in accuracy between both markerless techniques are directly visible in the created animations. For example, in the semi-automatic 3D animation considerable translational and rotational movements of the middle phalanges in relation to the distal and proximal phalanges are visible in Supplemental Video V3. In connection with the presented results, these movements are even more pronounced in the left bones, which are the ones farther from the camera. Because of the lower translational and rotational errors of the manual method, jerky movements of the individual bones in the manual animation (Supplemental Video V4) are less prominent compared with the ones in the semi-automatic animation. The animations in Supplemental Videos V1 to V4 (https://doi.org/10.3168/jds.2017-12563) all depict the same motion. Supplemental Video V1 shows the captured biplane high-speed fluoroscopic recordings, and V2 shows the marker-based animation of the bones in the biplane videos. Supplemental Video V3 depicts the semi-automatic animation of those bones and V4 the manual animation of the same bones. Supplemental Video V2 represents the gold standard. All movements that do not occur in V2 but are visible in V3 and V4 are therefore attributed to differences introduced by the noninvasive techniques.
We examined the tracking errors of corresponding bones to investigate if the accuracies of left and right phalanges differed. Both noninvasive animation techniques showed a difference in all corresponding bones for translational accuracy (P < 0.01). For rotational accuracy, we found only a difference for middle phalanges and semi-automatic technique's matching of distal phalanges (P < 0.01). Semi-automatic matching showed consistently better results for right phalanges, while manual animation values revealed no side preference. This result further supports use of the manual over the semi-automatic technique. Our second hypothesis stating that the accuracies of left and right phalanges do not differ significantly has to be partially rejected. The software accomplishes the matching between the DRR and the biplane x-ray images based on 3D bone shape and texture. The side preference of semi-automatic tracking may be explained with differences in image quality where one bone side possibly had distinctive contours, which increased detection by the software. The superimposition of bone structures and the resulting reduction in bone contour, which likely impaired semi-automatic matching, are the same for manual animation, but animator skill and anatomical knowledge could increase the accuracy of detail recognition. For both noninvasive animation techniques, bone contrast in each image was enhanced using an image sharpening filter. For further increase of animation accuracy, each image may be enhanced individually. The image contrast of the posterior and therefore the superimposed bone should be optimized for each camera independently from the other phalanges to improve their visibility.
The amount of time needed to create the animations is an important factor. Even when the bones were placed in a hierarchy with virtual joints, each bone had to be aligned individually for proper accordance of bone models with the x-ray image. Depending on the video sequence, manual animation of one distal limb took approximately 80 to 120 h. Semi-automatic animation took approximately the same amount of time because of low reliability of the software, numerous mismatches, and frequent necessity of verification of DRR accordance with the x-ray image. At the same time, the semi-automatic technique produced higher translational and rotational values but, as mentioned previously, was not capable of matching all 6 phalanges. Future software developments, such as improvement of bone contour detection for semi-automatic technique or enhancement of the images, may improve results. Nevertheless, our third hypothesis stated that both noninvasive markerless animation techniques are capable of matching all 6 phalanges throughout simulated steps has to be rejected for the semi-automatic technique. In contrast, our hypothesis is confirmed for manual animation. Manual animation achieved 3D animations of every phalanx, producing results with high accuracy, which offsets the effort and the expenditure of time. An appropriate bone hierarchy might increase manual animation's accuracy and may save time.
We successfully applied the method of biplane high-speed fluoroscopy to the bovine distal limb for the first time. We evaluated 2 noninvasive markerless animation techniques in a simulation model. Our results clearly show that bovine distal limb kinematics can be precisely measured with sub-millimeter accuracy without the need for invasive marker implantation. The ex vivo validation indicates that the manual animation method is the better choice. It is accurate to within 0.63 to 0.80 mm in translation and 2.41 to 6.75° in rotation. It provides a reference for the level of accuracy achievable for studies in vivo. It may enable new insights into claw biomechanics and claw–floor interactions, a key factor in understanding the pathogenesis of lameness as well as furthering its prevention. In the future, these animations may be applied to live animals to permit visualizing and evaluating movement of the distal phalanx inside the horn capsule. This technique represents a new, innovative approach for analyzing 3D bone and joint kinematics in 6 degrees of freedom for further locomotion research to maintain the well-being, performance, and productivity of dairy cows.
The project was supported by funds of the Landwirtschaftliche Rentenbank (Frankfurt am Main, Germany). The authors acknowledge Peter Boettcher (Department of Veterinary Medicine at the Freie Universitaet Berlin, Germany) and Dirk Brause, Juliane Munzel, Benjamin Oehme, and Nicole Röhrmann (Faculty of Veterinary Medicine of Leipzig University, Germany) for assisting with the collection of data and supporting figure preparation. The authors declare no competing or financial interests.
Röntgenanatomische und querschnittsanatomische Untersuchungen unter Berücksichtigung magnetresonanztomographischer Befunde an der Hintergliedmaße des Rindes. Faculty of Veterinary Medicine,