Next Article in Journal
Coordinated Control Strategy of Electro-Hydraulic Composite Braking Torque for the Distributed Electric Vehicles
Next Article in Special Issue
An Improved Data-Driven Calibration Method with High Efficiency for a 6-DOF Hybrid Robot
Previous Article in Journal
Machine Learning in CNC Machining: Best Practices
Previous Article in Special Issue
Elastostatic Stiffness Modeling and Performance Evaluation of a 2UPR–2PRU Redundantly Actuated Parallel Manipulator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A General Pose Recognition Method and Its Accuracy Analysis for 6-Axis External Fixation Mechanism Using Image Markers

Key Laboratory of Mechanism Theory and Equipment Design of Ministry of Education, Tianjin University, Tianjin 300350, China
*
Author to whom correspondence should be addressed.
Machines 2022, 10(12), 1234; https://doi.org/10.3390/machines10121234
Submission received: 14 November 2022 / Revised: 6 December 2022 / Accepted: 14 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Development and Applications of Parallel Robots)

Abstract

:
The 6-axis external fixation mechanism with Gough-Stewart configuration has been widely applied to the correction of long bone deformities in orthopedics. Pose recognition of the mechanism is essential for trajectory planning of bone correction, but is usually implemented by the surgeons’ experience, resulting in a relatively low level of correction accuracy. This paper proposes a pose recognition method based on novel image markers, and implements accuracy analysis. Firstly, a pose description of the mechanism is established with several freely installed markers, and the layout of the markers is also parametrically described. Then, a pose recognition method is presented by identifying the orientation and position parameters using the markers. The recognition method is general in that it encompasses all possible marker layouts, and the recognition accuracy is investigated by analyzing variations in the marker layout. On this basis, layout principles for markers that achieve a desired recognition accuracy are established, and an error compensation strategy for precision improvement is provided. Finally, experiments were conducted. The results show that volume errors of pose recognition were 0.368 ± 0.130 mm and 0.151 ± 0.045°, and the correction accuracy of the fracture model after taking compensation was 0.214 ± 0.573 mm and −0.031 ± 0.161°, validating the feasibility and accuracy of the proposed methods.

1. Introduction

In recent years, a new type of bone external fixation device in orthopedics, known as the 6-axis external fixation mechanism [1,2,3,4], has attracted extensive attention from the medical engineering field. This device can adjust fractured bone in six degrees of freedom while reducing the workload of surgeons. The 6-axis external fixation mechanism is based on the Gough-Stewart configuration [5,6], and is divided into manual fixators (Figure 1a–c) and robots (Figure 1d). It is usually composed of two fixation rings and six support struts, which respectively correspond to the moving platform, base platform, and SPS limbs (S represents the spherical joint, and P represents the prismatic driving joint). Each fixation ring connects with a bone segment through several percutaneous bone pins. By adjusting the lengths of the six support struts, the relative pose of the bone segments is adjusted to realize the deformity correction. The 6-axis external fixation mechanism has been utilized successfully in a variety of orthopedic treatments, such as minimally invasive surgery for fractures [7], osteotomy and correction for bone deformities [8], and distraction osteogenesis for bone defects [9].
Due to the individuality of each patient, the clinical application of the 6-axis external fixation mechanism requires developing an appropriate correction plan and determining the strut lengths for correction. From the robotics perspective, these two problems are trajectory planning in the operational space and kinematics solving in the joint space, respectively. The effectiveness of deformity correction is highly dependent on the accuracy of the mechanism’s trajectory, which is obtained through three stages: pose recognition by fluoroscopy measurement, trajectory planning by software, and motion execution by the mechanism. Pose recognition is fundamental because it determines the kinematic state of the mechanism and bone, building the foundation for the subsequent stages.
Most commercial fixators rely on manual measurement to recognize the pose. They require surgeons to draw lines and measure a certain number of image parameters on two orthogonal X-ray images [2,10]. Since it is difficult to achieve the orthogonality of image shooting and the geometrical requirements of fixator mounting under clinical conditions, the accuracy of the correction is quite low [11,12]. Moreover, measurement errors and the perspective effect in X-ray images cannot be avoided. The literature [9,13,14,15,16] shows an average correction accuracy of 4.0–7.3 mm in translation and 2.1–9.7° in angulation clinically, which is much lower than the accuracy of the products. Liu et al. utilized a drawing method to optimize the measurement of installation parameters [17] in order to increase the measurement accuracy of X-ray image parameters, whereas Ahrend et al. [18] used additional structural parts to achieve the orthogonal shooting of X-ray images. Their attempts have yielded some success, but pose recognition remains the bottleneck of the effectiveness of external fixation technology.
With the development of computed tomography (CT) and three-dimensional (3D) reconstruction, two types of 3D image guiding technology for orthopedics have emerged. One of them is surgical navigation [19,20,21], which rigidly connects each object to be tracked (including bone segments, robot manipulator and surgical instruments) with an optical frame tool. Intraoperative fluoroscopic images of these objects and 3D models are required to guide the surgeon’s operation. This type of technology is relatively mature, however, it requires expensive navigation equipment and the completion of all procedures during surgery. It is better suited for treatments requiring high real-time responsiveness, such as nail placement or joint shaping.
For the bone deformity corrections in question, another type of technology, the feature point positioning method, is better suited. The bone pose and the external fixator structure are recognized based on postoperative CT images, and the fixator can be adjusted at any time during hospitalization to achieve correction. Simpson et al. [22] extracted the coordinates of several feature points on the fixation rings, adjusted the 3D bone models to virtually reduce the fracture, and then solved the corresponding transformation matrix to obtain the strut lengths. Their model experiments achieved a correction accuracy of −2.5–3.0 mm and 0.8–4.4°. In contrast to the conventional measuring method, their method provided a more intuitive way of selecting graphical points, avoiding cumbersome measurements and reducing fixator mounting requirements. However, subjective errors still exist during point-picking, resulting in inadequate correction. Using metal marker balls on both sides of each support strut (twelve in total), Tang et al. [23] designed a 6-axis external fixation robot. The marker balls served as landmarks for the CT image, whose positions were determined by 3D reconstruction and point cloud fitting. Virtual fracture reduction was also utilized in their method to obtain strut lengths, however, the length of each strut was determined by the marker balls on both sides. In a later model experimental study [24], the accuracy of deformity correction was 0.8–2.5 mm in translation, 0.5–6.2° in lateral angulation, and 1.5–5.0° in axial rotation. Utilizing the point cloud data of marker balls improved the precision of position parameter identification. However, there were still some issues, including the redundant number of markers, the exclusiveness of the marker structure, and operational errors caused by segmenting the 3D reconstruction model.
Compared with conventional measurement methods, using image markers can simplify the operation and improve the accuracy of pose recognition. Our preliminary study found that the influence patterns of error factors and the universality of the recognition method need to be investigated. The existing recognition methods [22,23,24] adopt numerical solutions which are sensitive to error, and they require extra processes to deal with multiple solutions or no solution. For the accuracy problem of external fixation mechanism, related research [25,26,27] focuses mostly on the mechanism’s positioning precision. Their experiments indicate the mechanical system precision of the external fixation robots within 1 mm and 1°. For the bottleneck of pose recognition accuracy, however, it is essential to study the identification errors thoroughly and provide appropriate improvement strategies. Therefore, this paper presents a general method for pose recognition of the 6-axis external fixation mechanism using image markers, which is applicable for multiple existing products. By analyzing the effect of marker layout variations, pose recognition accuracy is acquired. The marker layout principles and error compensation strategy are further discussed in order to improve the accuracy for the deformity correction.
The rest of this paper is organized as follows. Section 2 addresses pose description of the 6-axis external fixation mechanism and parametric layout description of the novel designed marker. In Section 3, a general method for pose recognition of the mechanism using markers is proposed. In Section 4, the principles of the marker layout are established to guarantee pose recognition accuracy, using error modeling and analysis for pose recognition. An error compensation strategy is developed and model experiments are carried out in Section 5, before Section 6 presents the conclusions.

2. Descriptions of the Mechanism and Markers

Since each fixation ring connects rigidly to the bone segments, the relative movement of bone segments can be considered equal to that of the fixation rings under clinical precision. Therefore, the kinematic state of bones and rings can be uniformly described with the mechanism’s pose, avoiding measuring image parameters from the irregular geometry of bones. Image markers and their layout on the mechanism are also described.

2.1. Image Marker Design

As shown in Figure 2, an image marker was designed to locate the 6-axis external fixation mechanism (denoted in the following as the 6-SPS mechanism, for convenience) in the CT image space. The image marker consists of a marker ball made of aluminum alloy, and a connection rod made of photosensitive resin. The connection rod is manufactured by rapid prototyping, and its bottom buckle is suitable for the connection holes of the fixation ring in various products. After the marker is installed, the structure of the connection rod ensures the distance between marker ball’s center and the installing plane of the fixation ring as h mk . Utilizing the difference of X-ray attenuation between the connecting rod and the marker ball, the region of the marker ball in the CT image can be segmented automatically using threshold division technology. Consequently, the 3D reconstruction model of the marker balls is obtained independently and without manual labor.
To identify the position vector m of a marker ball in the CT image space, perform a 3D reconstruction of the marker ball and extract the surface point cloud data, then, fit the data into a sphere by a data fitting algorithm (for example, the least square method). The center of the sphere determines the position vector m .

2.2. Pose Description of the Mechanism

The two fixation rings of a 6-SPS mechanism are clinically named the proximal ring and the distal ring, following anatomical direction. Install N image markers on each ring for pose recognition; 2 N markers are used in total. The N markers on the same ring are inserted on the same side of the ring plane, and their placement is left to the surgeon’s choice. Because the identified marker position m contains errors, the marker layout will impact the accuracy of the pose recognition, which will be discussed in Section 4.
Assume 2 N = 6 . Figure 3a illustrates the proximal ring, the distal ring, the support strut l   ( l = 1 , 2 , , 6 ) , the bone segments, and the image markers. The proximal ring, the distal ring, and the l th support strut are designated as the moving platform, the base platform, and the l th limb, respectively. Denote the reference frame of the image space as the image frame { W } . Establish a frame { A } fixed to the moving platform with its origin O A at the center of the upper plane of the platform, z A -axis perpendicular upward to the platform, x A -axis toward the geometric front of the platform, and y A -axis determined by the right-hand rule. Similarly, a frame { B } is established for the base platform. The points A l and B l are the rotation centers of the two spherical joints on the l th limb, which connects to the moving platform and the base platform, respectively.
The pose of the 6-SPS mechanism is described by the relative pose of frame { A } with respect to frame { B } , that is, the mechanism’s operational space pose. The poses of the moving platform and the base platform in the CT image space need to be determined first. Number the markers on the moving platform from 1 to N , beginning from the x A -axis and increasing around the positive direction of z A -axis. The markers on the base platform are numbered from N + 1 to 2 N similarly. The point M i is the center of the marker ball i   ( i = 1 , 2 , , 2 N ) , whose position vector is the marker position m i identified previously, and will be used to determine the pose X A = ( p A T θ A T ) T and X B = ( p B T θ B T ) T of frame { A } and frame { B } , respectively. The position vectors p A , p A and orientation vectors θ A , θ B of pose X A and X B consist of parameters.
p A = ( x A y A z A ) T ,   p B = ( x B y B z B ) T ,   θ A = ( α A β A γ A ) T ,   θ B = ( α B β B γ B ) T
Wherein, orientation parameters α , β and γ are the three sequential rotation angles of the X-Y-Z fixed angles rotation, which characterize the orientation of the frame { A } or { B } by the orientation matrix
R = [ c β c γ s α s β c γ c α s γ c α s β c γ + s α s γ c β s γ s α s β s γ + c α c γ c α s β s γ s α c γ s β s α c β c α c β ]
where s and c are the sine and cosine operators, respectively. According to the general body direction of patients in CT scans, each orientation parameter is within the range ( 90 ° , 90 ° ) . Thus they can be determined by the entries R 21 at row 2 column 1, R 31 at row 3 column 1, and R 32 at row 3 column 2 of the orientation matrix R as
β = arcsin R 31 ,   α = arcsin ( R 32 / cos β ) ,   γ = arcsin ( R 21 / cos β )
Once the pose parameters of the two platforms are identified, the mechanism’s pose can be calculated as
p E = R B 1 ( p A p B ) ,   R E = R B 1 R A
here R B 1 is the inverse matrix of the orientation matrix R B . The position vector p E = ( x E y E z E ) T and orientation vector θ E = ( α E β E γ E ) T derived from R E describe the pose of the mechanism’s operational space. Denote the pose of the 6-SPS mechanism as X E = ( p E T θ E T ) T , establishing a pose description of the mechanism.

2.3. Layout Description of Markers

Using the identified position m i of an image marker, one point M i on the mechanism can be located in the CT image space. In order to recognize the pose of the mechanism, the layout information of the markers is required. The marker layout is parametrically described here, including the marker’s number N , and the installation parameters t i , φ i and h i of each marker i .
Take the marker layout on frame { A } as an illustration. As shown in Figure 3, denote the projection of point M i onto the x A y A -plane as point M t , i , the distance between points O A and M t , i as radius t i , the rotation angle from the x A -axis to the vector O A M t , i around the z A -axis as sweep angle φ i , and the signed distance from the x A y A -plane to point M i as height h i . The height h i is determined by the marker’s dimension h mk and the platform’s thickness h plt , that h i = h mk for the marker installed on the plane or h i = h mk h plt for beneath the plane. The difference between sweep angles φ i + 1 and φ i is defined as the included angle Δ φ i = φ i + 1 φ i . Using the installation parameters, the position vector of point M i can be expressed as
r i A = ( t i cos φ i t i sin φ i h i ) T
which describes the layout of each marker.

3. Pose Recognition of the Mechanism

This section proposes a pose recognition method by identifying the position and orientation parameters of the 6-SPS mechanism. Based on the marker layout parameters, this recognition method adopts an analytical approach to solving the issue of numerical solutions. The pose recognition of each platform is presented by taking frame { A } as an example, and then the pose of the mechanism’s operational space is determined.

3.1. Establishing Position Relationships of Markers

Based on the vector chain formed by the origin O A and position vector r i A in image frame { W } (Figure 3a), the closed-loop vector equation regarding the marker i is
p A + R A r i A m i = 0
If N = 2 markers are used, two points on the platform are settled, with one rotational freedom around the axis passing through these points. In other words, the solution of Equation (6), which includes two vector equations, is uncertain. Therefore, a minimum of three markers should be used for the pose recognition of a frame. Using N 3 markers, the pose X A of frame { A } can be identified by solving Equation (6). In order to implement analytical solutions, marker groups and geometric conditions are utilized.

3.2. Setting Up Marker Groups

Set up N groups of markers, such that three markers j 1 = j , j 2 = j + floor ( N / 3 ) , and j 3 = j + floor ( 2 N / 3 ) are included in the group j   ( j = 1 , 2 , , N ) . Herein, floor ( ) is the round-down function, and the addition computation of the marker numbers j 2 and j 3 follows an N-based overflow. Taking N = 6 markers as an example, the group j = 3 includes markers numbered j 1 = 3 , j 2 = 5 and j 3 = 1 . Determine the install parameters t j 1 , φ j 1 , …, r j 3 . Each group can be used to identify one set of pose parameters.

3.3. Identifying Pose Parameters of Each Group

Consider the parameter identification of group j . Denote the unit vectors of x A -axis and z A -axis as u ^ A ( j ) and w ^ A ( j ) , respectively. The unit vector v ^ A ( j ) of y A -axis and the orientation matrix R A ( j ) of frame { A } satisfy
R A ( j ) = [ u ^ A ( j ) v ^ A ( j ) w ^ A ( j ) ] = [ u ^ A ( j ) w ^ A ( j ) × u ^ A ( j ) w ^ A ( j ) ]
Based on the identified positions m j 1 , m j 2 and m j 3 of marker j 1 , j 2 and j 3 , use their relative positions
m Δ j 1 = m j 1 m j 2 ,   m Δ j 2 = m j 2 m j 3
to determine the unit vector of z A -axis as
w ^ A ( j ) = ( m Δ j 1 × m Δ j 2 ) / m Δ j 1 × m Δ j 2
where     represents the 2-norm of a vector. Taking the frame { A } as a reference, the position vector t j 1 A of the projection point M t , j 1 is invariant, and its unitization yields
t ^ j 1 A = ( cos φ j 1 sin φ j 1 0 ) T
In the image frame { W } , the description of vector t ^ j 1 A becomes R A t ^ j 1 A . The unitized relative position vector
m ^ Δ j 1 = m Δ j 1 / m Δ j 1
and vector R A t ^ j 1 A satisfy an axis-angle rotation relationship (Figure 3b): the vector R A t ^ j 1 A is the result of vector m ^ Δ j 1 rotating in the positive direction about z A -axis at angle λ j 1 . Expressing the axis-angle rotation with a transformation matrix R w ^ A ( j ) , λ j 1 yields
R A t ^ j 1 A = R w ^ A ( j ) , λ j 1 m ^ Δ j 1
Utilizing the Rodrigues formula, it is feasible to calculate the transformation matrix as
R w ^ A ( j ) , λ j 1 = I + [ w ^ A ( j ) × ] sin λ j 1 + [ w ^ A ( j ) × ] 2 ( 1 cos λ j 1 )
where I represents 3 × 3 identity matrix, and [ w ^ A ( j ) × ] is the skew-symmetric matrix of vector w ^ A ( j ) . The angle λ j 1 of the axis-angle rotation can be calculated with the included angle Δ φ j 1 = φ j 2 φ j 1 as
λ j 1 = arcsin ( ( t j 2 / M t , j 1 M t , j 2 ) sin Δ φ j 1 ) ,   M t , j 1 M t , j 2 = t j 1 2 + t j 2 2 2 t j 1 t j 2 cos Δ φ j 1
Substituting Equations (7) and (10) into Equation (12) leads to a linear equation system about the vector u ^ A ( j ) :
( I cos φ j 1 + [ w ^ A ( j ) × ] ) u ^ A ( j ) = R w ^ A ( j ) , λ j 1 m ^ Δ j 1
Denote the coefficient matrix of Equation (15) as G A ( j ) . If cos φ j 1 0 , the inverse matrix G A ( j ) 1 of coefficient matrix G A ( j ) exists, and the solution of u ^ A ( j ) is
u ^ A ( j ) = G A ( j ) 1 R w ^ A ( j ) , λ j 1 m ^ Δ j 1
Elseif cos φ j 1 = 0 , the vector t ^ j 1 A becomes ( 0 1 0 ) T , and Equation (12) turns into
v ^ A ( j ) = R w ^ A ( j ) , λ j 1 m ^ Δ j 1 .
Now the rest vector v ^ A ( j ) or u ^ A ( j ) can be calculated by the right-hand rule, using the vectors u ^ A ( j ) (from Equation (16)) and w ^ A ( j ) , or using v ^ A ( j ) (from Equation (17)) and w ^ A ( j ) . Determine the matrix R A ( j ) by Equation (7) with vectors u ^ A ( j ) , v ^ A ( j ) and w ^ A ( j ) , identify the orientation parameters θ A ( j ) = ( α A ( j ) β A ( j ) γ A ( j ) ) T of group j .
For identification of the position parameters, substituting the previously obtained orientation matrix R A ( j ) , the identified position m j 1 and the installation position r j 1 A of marker j 1 into Equation (6) leads to
p A ( j ) = ( x A ( j ) y A ( j ) z A ( j ) ) T = m j 1 R A ( j ) r j 1 A

3.4. Recognizing the Mechanism’s Pose

The above step has identified N groups of position and orientation parameters. Take the mean value of the identified position parameter p A ( j ) and orientation parameter θ A ( j ) to determine the recognized pose of frame { A } :
X A = ( p A T θ A T ) T = 1 N j = 1 N ( p A ( j ) T θ A ( j ) T ) T
For the pose X B of frame { B } , utilize the markers i = N + 1 , N + 2 , 2 N on the base platform and identify the position and orientation parameters similarly. After identifying the pose of moving platform X A and that of base platform X B in the image space, finish the pose recognition of the mechanism’s operational space X E by Equation (4).

4. Marker Layout Principles

The main error factors affecting the accuracy of pose recognition using image markers are CT scanning resolution limit, reconstruction point cloud distortion, image marker dimension error, and installation position error. In the clinical environment, it is difficult to reduce these error factors directly. This section establishes an error mapping model between the error factors and the pose recognition errors, and discusses the appropriate range of the marker layout.

4.1. Error Modeling of the Pose Recognition

The error model is about the pose recognition of a frame ( { A } or { B } ), and the subscript is omitted for simplicity. For the N markers utilized by the pose recognition of a frame, use image error parameters Δ m i = ( δ x m , i δ y m , i δ z m , i ) T to represent the systematic error of CT scan and 3D reconstruction, and installation error parameters Δ η i = ( δ t i δ φ i δ h i ) T to represent the dimensional error of the marker’s manufacture and installation. These error parameters are relatively small compared to the structural dimensions and installation parameters of the markers. Therefore, the linear approach is used to establish the differential error model of pose recognition. Denote the 6 N error parameters as a vector
( Δ m T Δ η T ) T = ( Δ m 1 T Δ m N T Δ η 1 T Δ η N T ) T
With regard to the identified pose parameters X = ( p T θ T ) T , calculating Jacobian matrix produces the error mapping model between the error parameters ( Δ m T Δ η T ) T and the pose recognition errors Δ X = ( Δ p T Δ θ T ) T :
Δ X = [ p m 1 p m N p η 1 p η N θ m 1 θ m N θ η 1 θ η N ] ( Δ m 1 Δ m N Δ η 1 Δ η N ) = J ( Δ m Δ η )
Each submatrix of Jacobian matrix J in Equation (21) represents the derivative matrix of the corresponding identified pose parameters with respect to error parameters, for instance,
p m 1 = [ x x m , 1 x y m , 1 x z m , 1 y x m , 1 y y m , 1 y z m , 1 z x m , 1 z y m , 1 z z m , 1 ] ,   θ η N = [ α t N α φ N α h N β t N β φ N β h N γ t N γ φ N γ h N ]

4.2. Analyzing the Effect of Marker Layout Variations

By considering the Jacobian matrix J as a function of the marker’s number and installation parameters, the trend of pose recognition errors versus the variation of marker layout can be found. According to the Monte Carlo method, assume that the error parameters conform to a certain probability distribution, and consider the pose recognition as a random sampling. Hence, the accuracy of pose recognition using a particular marker layout can be expressed as a mathematical expectation, which is derived from the recognition errors of multiple sampling.
Based on the error sensitivity analysis method [28,29], the Jacobian matrix J suggests that the installation parameter variations of radius t i and height h i have little effect on the error mapping. Therefore, the effect of the marker’s number N and sweep angle φ i is investigated here. A single-factor analysis approach is developed for the effect of marker layout variations. Based on computer simulation, the approach includes the following steps (Figure 4):
(1)
Configure the probability distribution to generate error parameters for each random sampling. A total of G t simulation groups were set univariately about the marker’s number N and sweep angle φ i . For the groups regarding the variation of number, take N = 3 , 4 , 5 , . For the groups regarding the variation of sweep angle, adopt combinations of discrete angles φ i = ψ , 2 ψ , 3 ψ , based on an interval ψ , and ensure that all included angles satisfy the condition Δ φ 1 + Δ φ 2 + + Δ φ N = 360 ° .
(2)
Generate K s samples of random error parameters ( Δ m T Δ η T ) k T for the simulation group g   ( g = 1 , 2 , , G t ) , and calculate the Jacobian matrix J g for the particular marker layout of this group (given parameters N , t i , φ i and h i ). The corresponding pose recognition errors of K s samples are
Δ X k = J g ( Δ m T Δ η T ) k T   ( k = 1 , 2 , , K s )
(1)
Take the absolute value of each sample, and calculate their expectation as
Δ X ¯ g = 1 K s k = 1 K s | Δ X k |
The expectation Δ X ¯ g represents the pose recognition accuracy about the marker layout adopted by group g .
(4)
Repeat steps 2–3 to assess the pose recognition accuracy of all simulation groups. Analyze the trend of pose recognition accuracy versus the variation of layout parameters, and further determine the marker layout principles.

4.3. Determining the Marker Layout Principles

According to our previous tests and relevant literature [30,31,32], assume that the image error parameters δ x m , i , δ y m , i and δ z m , i follow a normal distribution of N ( 0 , 0.3 2 ) . Take the marker’s installation deviations as position errors along the three principal axes of frame { A } or { B } . Regarding the manufacturing tolerance, the distribution of installation deviations is determined as ε x , i , ε y , i , ε z , i N ( 0 , ( 0.4 / 6 ) 2 ) based on the three-sigma rule. Then, the installation error parameters are calculated as
{ δ t i = ε x , i cos φ i + ε y , i sin φ i δ φ i = ( ε x , i sin φ i + ε y , i cos φ i ) ( 180 ° / π t i ) δ h i = ε z , i

4.3.1. The Number Principle

The simulation groups about the marker’s number take N = 3 , 4 , , 10 . The installation radius t i = 100 mm and height h i = 20 mm are fixed, and the sweep angle φ i is set following the markers evenly distributed along the platform’s circumferential direction. For each simulation group, K s = 2000 samples are generated, and the accuracy of pose recognition is evaluated. The results are shown in Figure 5.
In the same group, the recognition error of position parameter z and orientation parameter γ are slightly smaller than those of other position and orientation parameters. This agrees with the findings of the error sensitivity analysis, in that the sensitivity coefficients mapping error parameters ( Δ m T Δ η T ) k T to pose recognition errors Δ z and Δ γ are relatively low. Among the different groups, increasing the number of markers reduces the recognition errors partly, but not by order of magnitude. The pose recognition errors using N = 3 markers are
Δ X ¯ ( 3 mks ) = ( 0.165 mm 0.164 mm 0.137 mm 0.103 ° 0.106 ° 0.079 ° ) T
while increasing the marker number to N = 6 reduces the pose recognition errors by 27% in position and 25% in orientation. Considering the practicability and accuracy comprehensively, establish the principle of marker number as N = 3 .

4.3.2. The Sweep Angle Principle

Based on the preceding established number principle N = 3 , the different combinations of sweep angles φ i   ( i = 1 , 2 , 3 ) are taken for the simulation groups about the marker’s sweep angle. The remaining installation parameters are t i = 100 mm and h i = 20 mm . According to the rotational symmetry of the platform, the first marker can be placed at the x -axis of the platform, that is φ 1 = 0 . The sweep angles φ 2 and φ 3 are discretized with an interval ψ = 5 ° in a possible range, making the included angles satisfying
Δ φ 1 + Δ φ 2 + Δ φ 3 = 360 °
The accuracy of each group’s pose recognition was then evaluated by taking K s = 2000 samples. For clarity, the pose recognition accuracy was expressed as volume errors of position Δ p v and orientation Δ θ v , and the results are shown in Figure 6.
Δ p v = Δ x 2 + Δ y 2 + Δ z 2 ,   Δ θ v = Δ α 2 + Δ β 2 + Δ γ 2
The pose recognition errors decrease with the values of included angles Δ φ 1 and Δ φ 2 converging and their sum Δ φ 1 + Δ φ 2 approaching 240°, showing that evenly distributing markers on a platform reduces recognition errors. To determine the principle of the marker’s sweep angle, the included angles Δ φ 1 , Δ φ 2 and Δ φ 3 of each simulation group were arranged ascendingly intra-group, and then the groups were sorted according to the minimum-median-maximum priority of the included angle (Figure 7).
The scatters between two adjacent scales on the horizontal axis of Figure 7 are the groups whose minimum included angle equals the smaller scale value. The scatters from left to right represent the groups sorted with the minimum-median-maximum priority. Accoridng to the vertical axis value of the scatters, when the minimum and median included angles were both greater than 45°, the position and orientation recognition errors were less than 0.9 mm and 0.6°. If the minimum included angle reaches 15°, a difference between the maximum and median included angles of less than 150° and permitted recognition errors of less than 1 mm and 0.8°.
Based on the preceding analysis, the following are the layout principles for markers: use N = 3 markers for each platform, the minimum included angle Δ φ min 45 ° and median included angle Δ φ med 45 ° ; for limited installation conditions, Δ φ min 15 ° and the difference between maximum and median included angles Δ φ max Δ φ med 150 ° .

5. Error Compensation Strategy and Experiments

By installing markers following the proposed marker layout principles, the pose recognition errors can be kept within a small range. Compensating the recognized pose errors can further improve the accuracy of the correction results. Although current clinical practice lacks instruments for directly measuring pose errors, the commercial fixators ensure a certain degree of precision in the dimensions of their structure. In addition, kinematic calibrations can be performed on external fixation robots to reduce positioning errors to 0.05% of their operational space’s scale [33,34]. Therefore, recognition errors can be compensated using the dimension parameters of the 6-SPS mechanism. Model experiments of fracture reduction are also carried out in this section.

5.1. Compensation for Pose Recognition Errors

The error compensation strategy implements compensation in the mechanism’s operational space. Denote the mechanism’s recognized pose as X r . Within the precision of clinical practice, the mechanism’s theoretical pose X t can be considered equivalent to the mechanism’s actual pose and be used for compensation. X t is calculated through the forward kinematics model of the mechanism with nominal dimensions and limb lengths. The nominal dimensions are replaced with their identified values if the mechanism has been calibrated. Referring to Figure 3a, utilize the points A l and B l to establish the closed-loop vector equation of limb l :
p t + R t a l A b l B q t , l s ^ t , l = 0 ( l = 1 , 2 , , 6 )
where p t and R t are the position vector and orientation matrix corresponding to the theoretical pose X t , which are to be solved. a l A and b l B are the position vectors of point A l in frame { A } and point B l in frame { B } , respectively, which are determined by the mechanism’s dimensions. q t , l s ^ t , l is the product of length q t , l and unit axis vector s ^ t , l of the limb l . The limb length q t , l can be determined from the actual position of prismatic joint P (for example, the reading from motor encoders), while the vector s ^ t , l can be eliminated by taking norm of Equation (29). Writing the normed equations in a functional form we have
f l ( X t ) = p t + R t a l A b l B q t , l = 0
Combine Equation (30) from limb 1 to limb 6 to establish the mechanism’s forward kinematics model, and the solution X t can be obtained through numerical methods.
The recognition errors are calculated by the theoretical pose X t as
Δ X r = X r X t
For the target pose X a ( t ) at time t of the planned trajectory, use the recognition errors Δ X r to determine the compensated pose X c ( t ) as
X c ( t ) = X a ( t ) Δ X r
The compensated pose X c ( t ) is finally applied for the mechanism to achieve deformity correction. Regarding the execution of the mechanism, establish the inverse kinematics model based on Equation (30) to determine the corresponding joint space trajectory:
q c , l ( t ) = p c ( t ) + R c ( t ) a l A b l B
where the position vector p c ( t ) and orientation matrix R c ( t ) are derived from the compensated pose X c ( t ) .

5.2. Model Experiments of Fracture Reduction

To evaluate the methods proposed, we conducted model experiments of fracture reduction, as shown in Figure 8. The experiments used a self-developed external fixation mechanism. Two aluminum cylinders simulated a fracture’s proximal and distal bone segments. Each adjacent end of the cylinders featured a notch and a circumferential scale for correction planning and accuracy assessment. The fracture model with the mechanism was assembled, and the experiments were carried out as follows. The complete procedure of the experiment is presented in the Supplementary Video Material.
(1)
Adjust the lengths of limbs 1–6 arbitrarily to construct a random pose of the mechanism, simulating the initial fracture state. Theoretical pose X t is calculated according to the strut lengths q t , l provided by the motor encoders. Meanwhile, the relative initial pose X md 0 of the cylinders is determined by measuring the point pairs from circumferential scales.
(2)
Install the markers following the proposed marker layout principles, and then perform a CT scan and 3D reconstruction of the entire model. The 3D model of the marker balls was obtained by automatic segmentation based on the feature of markers described in Section 2. The 3D models of cylinders and mechanism platforms were used for correction planning in visual, thus they could be obtained easily through approximate manual segmentation. The positions m i of markers were then identified, and the mechanism’s recognized pose X r was determined by the proposed pose recognition method. Subsequently, the accuracy of pose recognition was analyzed by calculating the recognition errors Δ X r . In addition, a visual inspection was performed by registering standard 3D models to the reconstruction models of the platforms, using the recognized poses of frames { A } and { B } .
(3)
Using self-developed software, the 3D models of bone segments undergoing reduction motion were manipulated to design the fracture correction plan. The notch on the cylinder served as a reference point for the correction target. The correction motion was determined and the trajectory with target pose X a ( t ) was generated. After error compensation, the mechanism executed the compensated trajectory X c ( t ) and reduced the fracture model. The effectiveness of fracture reduction was evaluated by measuring the relative final pose X md 1 of the cylinders.
A caliper with 0.02 mm resolution was utilized for measuring. A uCT760 instrument (United Imaging, Shanghai, China) and MIMICS 20.0 software (Materialise, Leuven, Belgium) were utilized for CT scan and 3D reconstruction. Pose recognition, trajectory planning, and motion control were accomplished in self-developed software. The mechanism has been calibrated to an accuracy of 0.0574 mm and 0.0294°.
Ten reduction experiments on the fracture model have been conducted. The statistical results of pose recognition errors Δ X r are shown in Table 1, indicating that the proposed pose recognition method is highly accurate. The volume errors were 0.368 ± 0.130 mm in position and 0.151 ± 0.045°in orientation (mean ± standard deviation), close to the simulation results of Section 4, which validates the effectiveness of the proposed marker layout principles.
The relative poses of the fracture model before reduction X md 0 and after reduction X md 1 are shown in Table 2. After reduction, the average of relative pose was 0.214 ± 0.573 mm in position and −0.031 ± 0.161°in orientation, with a maximum deviation of 1.24 mm in z and −0.26°in β . As a result of the cylinders’ relative axial movement being obstructed by the contact of their end faces, the minimum deviation in z was 0. The mean and standard deviation of relative orientation γ after reduction was close to that of α and β , indicating that the proposed pose recognition method can correct the rotational deformity around the bone axis with higher accuracy than the X-ray measurement methods.
According to the model experiments, the proposed method for pose recognition has a high level of accuracy and, with the addition of error compensation, can achieve good reduction results. Moreover, preliminary clinical trials [35] based on the findings of this paper have been conducted. Twenty-one patients with tibial fractures were treated using 6-axis external fixators. The reduction results, as indicated by radiographic indices, were 0–1.72 mm in translation and 0–1.25° in angulation. Our method has significant advantages over existing clinical reports in terms of correction accuracy, and merits further application and promotion.

6. Conclusions

Towards the application of 6-axis external fixation mechanism in orthopedics, this paper proposes a general method for pose recognition using image markers, analyzes its accuracy, and puts forward accuracy improvements. Conclusions are as follows:
(1)
Measuring the pose of 6-axis external fixation mechanism in CT image space served as the foundation of deformity correction planning. The position and orientation parameters were utilized to describe the mechanism’s pose. Image markers were designed and implemented to eliminate subjective measurement errors of pose recognition, and their layout on the mechanism is parametrically described.
(2)
Utilizing CT scan and 3D reconstruction, an analytical method was developed for recognizing the mechanism’s pose based on 2 N 6 markers. The proposed method encompasses all possible marker layouts that can be implemented in practice, thereby expanding its applicability. In addition, the proposed method has more stable parameter identification compared to numerical methods.
(3)
The geometric error model of pose recognition was established. The effect of marker layout variations on the pose recognition errors were investigated. Based on the Monte Carlo method, the probability distribution of error parameters was set, and the single-factor analysis of layout parameters was carried out. The principles of marker layout were established to guide clinical application.
(4)
Ten groups of fracture model reduction experiments were conducted. A self-developed 6-axis external fixation mechanism was utilized to execute deformity correction. The results showed that the maximum errors of pose recognition were 0.40 mm in position and 0.22 ° in orientation, and the average accuracy of correction was 0.214 ± 0.573 mm and −0.031 ± 0.161° after compensation. It was demonstrated that the pose recognition method and accuracy improvements could achieve precise and safe correction of bone deformities.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/machines10121234/s1.

Author Contributions

Supervision, project administration and funding acquisition, Y.S.; methodology, software, validation and writing—original draft preparation, S.L.; data curation, visualization, and writing—review and editing, B.L.; data curation, conceptualization, investigation and resources, T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Tianjin Science and Technology Planning Project (No. 20201193).

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author Yimin Song ([email protected]) upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Taylor, J.C. Perioperative Planning for Two- and Three-Plane Deformities. Foot Ankle Clin. 2008, 13, 69–121. [Google Scholar] [CrossRef] [PubMed]
  2. Ferreira, N.; Birkholtz, F. Radiographic analysis of hexapod external fixators: Fundamental differences between the Taylor Spatial Frame and TrueLok-Hex. J. Med. Eng. Technol. 2015, 39, 173–176. [Google Scholar] [CrossRef] [PubMed]
  3. Solomin, L.N.; Paley, D.; Shchepkina, E.; Vilensky, V.A.; Skomoroshko, P.V. A comparative study of the correction of femoral deformity between the Ilizarov apparatus and Ortho-SUV frame. Int. Orthop. 2013, 38, 865–872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Gigi, R.; Mor, J.; Lidor, I.; Ovadia, D.; Segev, E. Auto Strut: A novel smart robotic system for external fixation device for bone deformity correction, a preliminary experience. J. Child. Orthop. 2021, 15, 130–136. [Google Scholar] [CrossRef] [PubMed]
  5. Gough, V.E. Universal tyre test machine. In Proceedings of the FISITA 9th International Automobile Technical Congress, London, UK, 30 April–5 May 1962; pp. 117–137. [Google Scholar]
  6. Stewart, D. A platform with six degrees of freedom. Proc. Inst. Mech. Eng. 1965, 180, 371–386. [Google Scholar] [CrossRef]
  7. Henderson, D.J.; Barron, E.; Hadland, Y.; Sharma, H.K. Functional Outcomes After Tibial Shaft Fractures Treated Using the Taylor Spatial Frame. J. Orthop. Trauma 2015, 29, e54–e59. [Google Scholar] [CrossRef]
  8. Floerkemeier, T.; Stukenborg-Colsman, C.; Windhagen, H.; Waizy, H. Correction of Severe Foot Deformities Using the Taylor Spatial Frame. Foot Ankle Int. 2011, 32, 176–182. [Google Scholar] [CrossRef]
  9. Marangoz, S.; Feldman, D.S.; Sala, D.A.; Hyman, J.E.; Vitale, M.G. Femoral Deformity Correction in Children and Young Adults Using Taylor Spatial Frame. Clin. Orthop. Relat. Res. 2008, 466, 3018–3024. [Google Scholar] [CrossRef] [Green Version]
  10. Sökücü, S.; Demir, B.; Lapçin, O.; Yavuz, U.; Kabukçuoğlu, Y.S. Perioperative versus postoperative measurement of Taylor Spatial Frame mounting parameters. Acta Orthop. Traumatol. Turc. 2014, 48, 491–494. [Google Scholar] [CrossRef]
  11. Ahrend, M.-D.; Baumgartner, H.; Ihle, C.; Histing, T.; Schröter, S.; Finger, F. Influence of axial limb rotation on radiographic lower limb alignment: A systematic review. Arch. Orthop. Trauma Surg. 2021, 142, 3349–3366. [Google Scholar] [CrossRef]
  12. Jamali, A.A.; Meehan, J.P.; Moroski, N.M.; Anderson, M.J.; Lamba, R.; Parise, C. Do small changes in rotation affect measurements of lower extremity limb alignment? J. Orthop. Surg. Res. 2017, 12, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Arvesen, J.E.; Watson, J.T.; Israel, H. Effectiveness of Treatment for Distal Tibial Nonunions With Associated Complex Deformities Using a Hexapod External Fixator. J. Orthop. Trauma 2017, 31, e43–e48. [Google Scholar] [CrossRef] [PubMed]
  14. Eren, I.; Eralp, L.; Kocaoglu, M. Comparative clinical study on deformity correction accuracy of different external fixators. Int. Orthop. 2013, 37, 2247–2252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Naqui, S.Z.H.; Thiryayi, W.; Foster, A.; Tselentakis, G.; Evans, M.; Day, J.B. Correction of Simple and Complex Pediatric Deformities Using the Taylor-Spatial Frame. J. Pediatr. Orthop. 2008, 28, 640–647. [Google Scholar] [CrossRef]
  16. Manner, H.M.; Huebl, M.; Radler, C.; Ganger, R.; Petje, G.; Grill, F. Accuracy of complex lower-limb deformity correction with external fixation: A comparison of the Taylor Spatial Frame with the Ilizarov ring fixator. J. Child. Orthop. 2007, 1, 55–61. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, Y.; Yushan, M.; Liu, Z.; Liu, J.; Ma, C.; Yusufu, A. Application of elliptic registration and three-dimensional reconstruction in the postoperative measurement of Taylor spatial frame parameters. Injury 2020, 51, 2975–2980. [Google Scholar] [CrossRef]
  18. Ahrend, M.-D.; Finger, F.; Grünwald, L.; Keller, G.; Baumgartner, H. Improving the accuracy of patient positioning for long-leg radiographs using a Taylor Spatial Frame mounted rotation rod. Arch. Orthop. Trauma. Surg. 2020, 141, 55–61. [Google Scholar] [CrossRef]
  19. Dagnino, G.; Georgilas, I.; Tarassoli, P.; Atkins, R.; Dogramadzi, S. Vision-based real-time position control of a semi-automated system for robot-assisted joint fracture surgery. Int. J. Comput. Assist. Radiol. Surg. 2015, 11, 437–455. [Google Scholar] [CrossRef] [Green Version]
  20. Li, C.; Wang, T.; Hu, L.; Zhang, L.; Du, H.; Zhao, L.; Wang, L.; Tang, P. A visual servo-based teleoperation robot system for closed diaphyseal fracture reduction. Proc. Inst. Mech. Eng. Part H J. Eng. Med. 2015, 229, 629–637. [Google Scholar] [CrossRef]
  21. Dagnino, G.; Georgilas, I.; Morad, S.; Gibbons, P.; Tarassoli, P.; Atkins, R.; Dogramadzi, S. Intra-operative fiducial-based CT/fluoroscope image registration framework for image-guided robot-assisted joint fracture surgery. Int. J. Comput. Assist. Radiol. Surg. 2017, 12, 1383–1397. [Google Scholar] [CrossRef]
  22. Simpson, A.L.; Ma, B.; Slagel, B.; Borschneck, D.P.; Ellis, R.E. Computer-assisted distraction osteogenesis by Ilizarov’s method. Int. J. Med. Robot. Comput. Assist. Surg. 2008, 4, 310–320. [Google Scholar] [CrossRef] [PubMed]
  23. Tang, P.; Hu, L.; Du, H.; Gong, M.; Zhang, L. Novel 3D hexapod computer-assisted orthopaedic surgery system for closed diaphyseal fracture reduction. Int. J. Med. Robot. Comput. Assist. Surg. 2011, 8, 17–24. [Google Scholar] [CrossRef] [PubMed]
  24. Hu, L.; Zhang, J.; Li, C.; Wang, Y.; Yang, Y.; Tang, P.; Fang, L.; Zhang, L.; Du, H.; Wang, L. A femur fracture reduction method based on anatomy of the contralateral side. Comput. Biol. Med. 2013, 43, 840–846. [Google Scholar] [CrossRef] [PubMed]
  25. Li, C.; Wang, T.; Hu, L.; Zhang, L.; Du, H.; Wang, L.; Luan, S.; Tang, S.L.A.P. Accuracy Analysis of a Robot System for Closed Diaphyseal Fracture Reduction. Int. J. Adv. Robot. Syst. 2014, 11, 169. [Google Scholar] [CrossRef]
  26. Faschingbauer, M.; Heuer, H.J.D.; Seide, K.; Wendlandt, R.; Münch, M.; Jürgens, C.; Kirchner, R. Accuracy of a hexapod parallel robot kinematics based external fixator. Int. J. Med. Robot. Comput. Assist. Surg. 2015, 11, 424–435. [Google Scholar] [CrossRef]
  27. Liu, M.; Li, J.; Sun, H.; Guo, X.; Xuan, B.; Ma, L.; Xu, Y.; Ma, T.; Ding, Q.; An, B. Study on the Modeling and Compensation Method of Pose Error Analysis for the Fracture Reduction Robot. Micromachines 2022, 13, 1186. [Google Scholar] [CrossRef]
  28. Sun, T.; Zhai, Y.; Song, Y.; Zhang, J. Kinematic calibration of a 3-DoF rotational parallel manipulator using laser tracker. Robot. Comput. Manuf. 2016, 41, 78–91. [Google Scholar] [CrossRef]
  29. Song, Y.; Zhang, J.; Lian, B.; Sun, T. Kinematic calibration of a 5-DoF parallel kinematic machine. Precis. Eng. 2016, 45, 242–261. [Google Scholar] [CrossRef]
  30. Liang, X.; Lambrichts, I.; Sun, Y.; Denis, K.; Hassan, B.; Li, L.; Pauwels, R.; Jacobs, R. A comparative evaluation of Cone Beam Computed Tomography (CBCT) and Multi-Slice CT (MSCT). Part II: On 3D model accuracy. Eur. J. Radiol. 2010, 75, 270–274. [Google Scholar] [CrossRef]
  31. Prevrhal, S.; Fox, J.C.; Shepherd, J.A.; Genant, H.K. Accuracy of CT-based thickness measurement of thin structures: Modeling of limited spatial resolution in all three dimensions. Med. Phys. 2002, 30, 1–8. [Google Scholar] [CrossRef]
  32. Van Eijnatten, M.; van Dijk, R.; Dobbe, J.; Streekstra, G.; Koivisto, J.; Wolff, J. CT image segmentation methods for bone used in medical additive manufacturing. Med. Eng. Phys. 2018, 51, 6–16. [Google Scholar] [CrossRef] [PubMed]
  33. He, Z.; Lian, B.; Li, Q.; Zhang, Y.; Song, Y.; Yang, Y.; Sun, T. An Error Identification and Compensation Method of a 6-DoF Parallel Kinematic Machine. IEEE Access 2020, 8, 119038–119047. [Google Scholar] [CrossRef]
  34. He, Z.; Song, Y.; Lian, B.; Sun, T. Kinematic Calibration of a 6-DoF Parallel Manipulator with Random and Less Measurements. IEEE Trans. Instrum. Meas. 2022. [Google Scholar] [CrossRef]
  35. Liu, Q.; Liu, Y.; Li, H.; Fu, X.; Zhang, X.; Liu, S.; Zhang, J.; Zhang, T. Marker- three dimensional measurement versus traditional radiographic measurement in the treatment of tibial fracture using Taylor spatial frame. BMC Musculoskelet. Disord. 2022, 23, 155. [Google Scholar] [CrossRef]
Figure 1. Commercial external fixation devices including (a) Taylor Spatial Frame, (b) TL-HEX, (c) Ortho-SUV, and (d) Auto Strut.
Figure 1. Commercial external fixation devices including (a) Taylor Spatial Frame, (b) TL-HEX, (c) Ortho-SUV, and (d) Auto Strut.
Machines 10 01234 g001
Figure 2. (a) Structure of the image marker. (b) Image marker installed on a fixation ring. (c) CT image region segmentation. (d) 3D reconstruction result.
Figure 2. (a) Structure of the image marker. (b) Image marker installed on a fixation ring. (c) CT image region segmentation. (d) 3D reconstruction result.
Machines 10 01234 g002
Figure 3. (a) Schematic of the mechanism and image markers. (b) Projection of the marker centers on x A y A -plane.
Figure 3. (a) Schematic of the mechanism and image markers. (b) Projection of the marker centers on x A y A -plane.
Machines 10 01234 g003
Figure 4. Flow chart of the single-factor analysis.
Figure 4. Flow chart of the single-factor analysis.
Machines 10 01234 g004
Figure 5. Recognition errors in (a) position and (b) orientation using different numbers of markers.
Figure 5. Recognition errors in (a) position and (b) orientation using different numbers of markers.
Machines 10 01234 g005
Figure 6. Recognition errors in (a) position and (b) orientation by using different combinations of included angles.
Figure 6. Recognition errors in (a) position and (b) orientation by using different combinations of included angles.
Machines 10 01234 g006
Figure 7. Recognitions errors of sorted marker groups, with position error in blue color and orientation error in pink color.
Figure 7. Recognitions errors of sorted marker groups, with position error in blue color and orientation error in pink color.
Machines 10 01234 g007
Figure 8. (a) Initial state of the 6-SPS mechanism fixing a fracture model. (b) 3D reconstruction models and registered standard models. (c) Axially aligned state and (d) reduction complete state of the fracture model.
Figure 8. (a) Initial state of the 6-SPS mechanism fixing a fracture model. (b) 3D reconstruction models and registered standard models. (c) Axially aligned state and (d) reduction complete state of the fracture model.
Machines 10 01234 g008
Table 1. Pose recognition errors of the model experiment.
Table 1. Pose recognition errors of the model experiment.
x (mm)y (mm)z (mm)α (°)β (°)γ (°)
Range−0.32–0.36−0.39–0.40−0.23–0.33−0.10–0.22−0.15–0.21−0.12–0.09
Mean0.062−0.0360.0280.0190.024−0.015
Standard deviation0.2270.2740.1850.1020.1070.064
Table 2. Relative pose of the model before and after reduction.
Table 2. Relative pose of the model before and after reduction.
x (mm)y (mm)z (mm)α (°)β (°)γ (°)
Range (before)−12.87–9.78−13.19–16.806.18–20.75−7.41–6.76−8.32–4.54−2.43–3.55
Range (after)−1.04–0.65−0.79–1.050–1.24−0.23–0.18−0.26–0.24−0.21–0.21
Mean (after)−0.1480.2630.529−0.073−0.0570.037
S.D. (after)0.6180.4960.4120.1780.1580.139
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, S.; Song, Y.; Lian, B.; Sun, T. A General Pose Recognition Method and Its Accuracy Analysis for 6-Axis External Fixation Mechanism Using Image Markers. Machines 2022, 10, 1234. https://doi.org/10.3390/machines10121234

AMA Style

Liu S, Song Y, Lian B, Sun T. A General Pose Recognition Method and Its Accuracy Analysis for 6-Axis External Fixation Mechanism Using Image Markers. Machines. 2022; 10(12):1234. https://doi.org/10.3390/machines10121234

Chicago/Turabian Style

Liu, Sida, Yimin Song, Binbin Lian, and Tao Sun. 2022. "A General Pose Recognition Method and Its Accuracy Analysis for 6-Axis External Fixation Mechanism Using Image Markers" Machines 10, no. 12: 1234. https://doi.org/10.3390/machines10121234

APA Style

Liu, S., Song, Y., Lian, B., & Sun, T. (2022). A General Pose Recognition Method and Its Accuracy Analysis for 6-Axis External Fixation Mechanism Using Image Markers. Machines, 10(12), 1234. https://doi.org/10.3390/machines10121234

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