- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
THR Simulator – the software for generating radiographs of THR prosthesis
© Wu et al; licensee BioMed Central Ltd. 2009
Received: 24 March 2008
Accepted: 16 January 2009
Published: 16 January 2009
Measuring the orientation of acetabular cup after total hip arthroplasty is important for prognosis. The verification of these measurement methods will be easier and more feasible if we can synthesize prosthesis radiographs in each simulated condition. One reported method used an expensive mechanical device with an indeterminable precision. We thus develop a program, THR Simulator, to directly synthesize digital radiographs of prostheses for further analysis.
Under Windows platform and using Borland C++ Builder programming tool, we developed the THR Simulator. We first built a mathematical model of acetabulum and femoral head. The data of the real dimension of prosthesis was adopted to generate the radiograph of hip prosthesis. Then with the ray tracing algorithm, we calculated the thickness each X-ray beam passed, and then transformed to grey scale by mapping function which was derived by fitting the exponential function from the phantom image. Finally we could generate a simulated radiograph for further analysis.
Using THR Simulator, the users can incorporate many parameters together for radiograph synthesis. These parameters include thickness, film size, tube distance, film distance, anteversion, abduction, upper wear, medial wear, and posterior wear. These parameters are adequate for any radiographic measurement research. This THR Simulator has been used in two studies, and the errors are within 2° for anteversion and 0.2 mm for wearing measurement.
We design a program, THR Simulator that can synthesize prosthesis radiographs. Such a program can be applied in future studies for further analysis and validation of measurement of various parameters of pelvis after total hip arthroplasty.
Measuring the orientation of acetabulum cup and the wearing of insert on plain radiograph of patients who underwent total hip arthroplasty is important for prognosis. Verifying the orientation measurement [1–6] and wearing [7–14] methods are both important, which may require a simulator to mimic every situation for such an analysis. Mechanical simulator has once been reported in a study to measure the wearing of acetabular insert . Although such a mechanical device is straightforward, there are disadvantages including expensive price, undetermined precision, as well as requiring image processing from radiograph to digital form. Every processing step may cause error and interfere with the final precision.
Many reported methods used Fourier transformation to fasten the process in generating the radiographs from computed tomography data [15, 16]. However, Fourier transformation may decrease the precision, which is the first priority in the measurement analysis on plain radiographs. On the other hand, ray tracing, which is popular in computer game, may be suitable for this transformation. Unfortunately, current built libraries only provide reflection images instead of transparent images that are needed in our analysis. Therefore we have to build up our whole software program before practical application in the plain radiographs.
Another problem is physics, i.e., once X-rays beam pass through the prosthesis, they then generate the image on the radiogram film. The grey scale on the radiogram film is determined by the amount of the X-ray passed that is dependent on the thickness of the metal in the pathway. Such condition follows Beer-Lambert law.
Penetration = e-kbc
k: molar absorbability
b: path length
The parameter k is different among various metals and radiation energy (kv in X-ray). In real X-ray machine, the distribution of kv follows the rule of normal distribution, which is different among various X-ray machines. After calculating the amount of X-rays passed, we need another formula to transfer them to grey scale.
Our goal is to build a simulated total hip prosthesis. Virtually, femoral head equals to a ball.
x 2+y 2+z 2<r f 2
(x, y, z) is the point of the simulated three-dimensional Cartesian coordinate system. r f is the radius of femoral head.
In our program, we make the ball move to simulate wearing of insert.
(x-d x )2+(y-d y )2+(z-d z )2<r f 2
d x , d y , d z are femoral head movement in three directions.
Virtually, acetabulum is composed of two balls and one plane.
x 2+y 2+z 2<r ao 2
x 2+y 2+z 2>r io 2
r ao means radius of acetabulum's outer shell, r io means radius of acetabulum's inner shell, (a, b, c) means the normal vector of the acetabulum which can be derived from inclination and anteversion of acetabulum. Liaw et al. derived the following formula for this process.
(a, b, c) = (sinφ × cosθ, -cosφ × cosθ, sinθ)
Vector (a, b, c) means the normal vector of the acetabulum, φ means the inclination of acetabulum, θ means the anteversion of acetabulum, positive θ means anteversion, and negative θ means retroversion.
Theoretically, the X-ray source is set at (0,0,-d t ). d t means tube distance (the X-ray tube to the acetabulum center). The points at film are (x f , y f , d f ). (x f , y f ) means point at film. d f means distance from film to the acetabulum center.
(x, y, z) = (t x f , t y f , t(d f +d t )-d t ) 0<t <1
The ray-tracing algorithm means calculating every ray from X-ray source to film. We used formula (7), simulating every X-ray from source to film. We then determined the total length the X-ray beam passed through femoral head by calculating the length between the two extreme solutions of formulas (3) and (7). Finally, we came out at the total length the X-ray beam passed through acetabulum by calculating length between the two extreme solutions of formulas (4), (5), (6), and (7). We use analytical mathematics for these calculations. The detailed process is illustrated in Appendix section.
In summary, we first built a mathematical model of acetabulum and femoral head, formulas (3) to (7). The real dimension data was adopted to generate the proper prosthesis figure. Then we calculated the total thickness of metal the X-ray beams passed by the ray-tracing algorithm, and then transformed these digital data to grey scale by mapping function. Finally, the gray scale generated for the prosthesis was shifted according to the aforementioned method. We could generate various radiographs according to the different parameters used.
The functional parameters in the THR Simulator include the following:
Thickness: refers to the thickness of the acetabulum shell.
Film size: refers to the X-ray film size it simulates.
Tube distance: refers to the distance from X-ray tube to the acetabulum center.
Film distance: refers to the distance from X-ray film to the acetabulum center.
Anteversion: refers to the version of the acetabular cup. The user can choose either version (anatomical, radiological, or operational) to simulate. Negative value means retroversion.
Abduction: refers to the abduction (or inclination) of the acetabular cup.
Upper movement: refers to the femoral head moving upward.
Medial movement: refers to the femoral head moving medially.
Posterior movement: refers to the femoral head moving posteriorly.
Picture size: refers to the size of picture file to simulate.
Results and discussion
Under Windows platform and using Borland C++ Builder programming tool, we developed the THR Simulator shown in Figure 1. In the THR Simulator, the user can incorporate many parameters before generating a simulated radiograph, i.e., with different parameters adopted, we can generate different radiographs. These can be used for further analysis.
Previously reported methods for measurement of acetabular cup orientation include mechanical simulator used for verifying insert wearing measuring methods , as well as other methods adopting Fourier transformation algorithm [15, 16]. These methods have their own disadvantages. The mechanical simulator can be used to generate radiographs directly but has some inherent problems. The simulation process takes much time, about 50 seconds with Pentium III 500 MHz notebook. On the other hand, other methods using Fourier transformation algorithm can do it in real time [15, 16] at the expense of lower precision. Thus we abandoned Fourier transformation algorithm and developed our own precise algorithm.
Parameters of total hip prostheses (U2, United Orthopedic Corporation, Hsinchu, Taiwan).
Acetabular shell diameter(mm)
Acetabular insert thickness (mm)
Femoral head diameter (mm)
Acetabular shell thickness (mm)
To make simulated radiographs more real, users can superimpose the synthesized radiographs onto real radiographs. We do not routinely recommend doing so because this action may make users misread patient's position and thus confuse the standardization process.
We use upper, medial, and posterior movements to indicate wearing of the three directions. Because femoral heads may not locate in the center of the acetabulum, the user should adjust these three vectors to fit every situation.
Our program can accommodate a three-dimensional wear vector by incorporating movements in the inferior-superior (upper), medial-lateral, and anterior-posterior directions. Inferior-superior and medial-lateral wears change the location of the femoral head relative to the cup while anterior-posterior wear changes the apparent size of the femoral head. However, the change in size is negligible.
Calculating wearing volume is another interesting issue. Kosak et al. has published a mathematical model to calculate it with the three wearing directions.
We designed new software THR Simulator that can generate radiographs after total hip arthroplasty. The strength is its accuracy and precision. The limitation is that it can not synthesize the details of the prosthesis and surrounding bone. We hope it can be used in future studies about measurements of geometrical parameters of pelvis after total hip arthroplasty.
Availability and requirements
Mathematical detail of calculating metal thickness
Equation (7) shows the line from X-ray source to the film.
We first change Formula (3) to Equation (3e).
(x-d x )2+(y-d y )2+(z-d z )2 = r f 2
Then we find the solution of t between Equations (7) and (3e). If there is no solution or only one solution, the thickness is zero, otherwise it means that the X-ray beam pass through it. From t we can calculate the real point using Formula (3), and then we can calculate the distance between the two points, which is the metal thickness the X-ray passed.
Similarly, we change Formula (4), (5), and (6) to Equations (4e), (5e), and (6e).
x 2+y 2+z 2 = r ao 2
x 2+y 2+z 2 = r io 2
ax+by+cz = 0
Then we find the solution of t between Equations (7) and (4e). If there is no solution or only one solution, the thickness is zero, otherwise it means that the X-ray beam pass through it. We keep the solutions of t in the solution set. Similarly, we find the solution of t from Equations (7) and (5e). Then we keep the solutions of t in the solution set if there are solutions.
We solve t from Equation (7) and Formula (6) and get the range of t, i.e. t>rx or t<rx, and rx is the solution of t from Equations (7) and (6e).
We exclude the solutions in the solution set outside the range of t. If the number of solutions in the solution set is odd, we append rx into the set.
Now we sort the solutions in the set, and pair the solutions with their neighborhood. We apply these paired t solutions to Equation (7) and get pairs of points.
Finally we add the distance of these pair points to we previously calculated distance intra-head. The sum is the total thickness the X-ray beam passed through the metal.
We thank United Orthopedic Corporation, Hsinchu, Taiwan for providing us technical data of U2 total hip arthroplasty.
This study was supported by the grant of NSC96-2320-B-087-001, Taiwan, ROC.
- Liaw CK, et al: A new tool for measuring cup orientation in total hip arthroplasties from plain radiographs. Clin Orthop Relat Res. 2006, 451: 134-9. 10.1097/01.blo.0000223988.41776.fa.View ArticlePubMedGoogle Scholar
- Fabeck L, et al: A method to measure acetabular cup anteversion after total hip replacement. Acta Orthop Belg. 1999, 65 (4): 485-91.PubMedGoogle Scholar
- Ackland MK, Bourne WB, Uhthoff HK: Anteversion of the acetabular cup. Measurement of angle after total hip replacement. J Bone Joint Surg Br. 1986, 68 (3): 409-13.PubMedGoogle Scholar
- Pradhan R: Planar anteversion of the acetabular cup as determined from plain anteroposterior radiographs. J Bone Joint Surg Br. 1999, 81 (3): 431-5. 10.1302/0301-620X.81B3.9067.View ArticlePubMedGoogle Scholar
- Visser JD, Konings JG: A new method for measuring angles after total hip arthroplasty. A study of the acetabular cup and femoral component. J Bone Joint Surg Br. 1981, 63B (4): 556-9.PubMedGoogle Scholar
- Widmer KH: A simplified method to determine acetabular cup anteversion from plain radiographs. J Arthroplasty. 2004, 19 (3): 387-90. 10.1016/j.arth.2003.10.016.View ArticlePubMedGoogle Scholar
- Ebramzadeh E, et al: Accuracy of measurement of polyethylene wear with use of radiographs of total hip replacements. J Bone Joint Surg Am. 2003, 85-A (12): 2378-84.PubMedGoogle Scholar
- Livermore J, Ilstrup D, Morrey B: Effect of femoral head size on wear of the polyethylene acetabular component. J Bone Joint Surg Am. 1990, 72 (4): 518-28.PubMedGoogle Scholar
- Hardinge K, et al: Measurement of hip prostheses using image analysis. The maxima hip technique. J Bone Joint Surg Br. 1991, 73 (5): 724-8.PubMedGoogle Scholar
- Shaver SM, et al: Digital edge-detection measurement of polyethylene wear after total hip arthroplasty. J Bone Joint Surg Am. 1997, 79 (5): 690-700.PubMedGoogle Scholar
- Martell JM, Berdia S: Determination of polyethylene wear in total hip replacements with use of digital radiographs. J Bone Joint Surg Am. 1997, 79 (11): 1635-41.PubMedGoogle Scholar
- Collier MB, et al: Evaluation of contemporary software methods used to quantify polyethylene wear after total hip arthroplasty. J Bone Joint Surg Am. 2003, 85-A (12): 2410-8.PubMedGoogle Scholar
- Devane PA: Measurement of polyethylene wear in metal-backed acetabular cups. I. Three-dimensional technique. Clin Orthop Relat Res. 1995, 303-16. 319Google Scholar
- Geerdink CH, et al: The determination of linear and angular penetration of the femoral head into the acetabular component as an assessment of wear in total hip replacement: A COMPARISON OF FOUR COMPUTER-ASSISTED METHODS. J Bone Joint Surg Br. 2008, 90 (7): 839-46. 10.1302/0301-620X.90B7.20305.View ArticlePubMedGoogle Scholar
- Calhoun PS, et al: Three-dimensional volume rendering of spiral CT data: theory and method. Radiographics. 1999, 19 (3): 745-64.View ArticlePubMedGoogle Scholar
- Muniyandi M, et al: Real-time PC based X-ray simulation for interventional radiology training. Stud Health Technol Inform. 2003, 94: 233-9.PubMedGoogle Scholar
- Liaw CK, et al: A simple mathematical standardized measurement of acetabulum anteversion after total hip arthroplasty. Computational and Mathematical Methods in Medicine. 2008, 9 (2): 105-119. 10.1080/17486700701865265.View ArticleGoogle Scholar
- Liaw CK, et al: Automatic digital PE wear measurement. AAOS annual meeting. Washington. 2005Google Scholar
- Kosak R, et al: Polyethylene wear in total hip prostheses: the influence of direction of linear wear on volumetric wear determined from radiographic data. Skeletal Radiol. 2003, 32 (12): 679-86. 10.1007/s00256-003-0685-2.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2474/10/8/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.