Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 20 October 2021
Sec. Biomechanics

Three-Dimensional Subject-Specific Knee Shape Reconstruction with Asynchronous Fluoroscopy Images Using Statistical Shape Modeling

Hsuan-Yu Lu&#x;Hsuan-Yu Lu1Kao-Shang Shih,&#x;Kao-Shang Shih2,3Cheng-Chung LinCheng-Chung Lin4Tung-Wu Lu,
Tung-Wu Lu1,5*Song-Ying LiSong-Ying Li1Hsin-Wen KuoHsin-Wen Kuo1Horng-Chaung HsuHorng-Chaung Hsu6
  • 1Department of Biomedical Engineering, National Taiwan University, Taipei, Taiwan
  • 2Department of Orthopedics, Shin Kong Wu Ho-Su Memorial Hospital, Taipei, Taiwan
  • 3School of Medicine, Fu Jen Catholic University, Taipei, Taiwan
  • 4Department of Electrical Engineering, Fu Jen Catholic University, Taipei, Taiwan
  • 5Department of Orthopaedic Surgery, School of Medicine, National Taiwan University, Taipei, Taiwan
  • 6Department of Orthopaedic Surgery, China Medical University, Taipei, Taiwan

Background and objectives: Statistical shape modeling (SSM) based on computerized tomography (CT) datasets has enabled reasonably accurate reconstructions of subject-specific 3D bone morphology from one or two synchronous radiographs for clinical applications. Increasing the number of radiographic images may increase the reconstruction accuracy, but errors related to the temporal and spatial asynchronization of clinical alternating bi-plane fluoroscopy may also increase. The current study aimed to develop a new approach for subject-specific 3D knee shape reconstruction from multiple asynchronous fluoroscopy images from 2, 4, and 6 X-ray detector views using a CT-based SSM model; and to determine the optimum number of planar images for best accuracy via computer simulations and in vivo experiments.

Methods: A CT-based SSM model of the knee was established from 60 training models in a healthy young Chinese male population. A new two-phase optimization approach for 3D subject-specific model reconstruction from multiple asynchronous clinical fluoroscopy images using the SSM was developed, and its performance was evaluated via computer simulation and in vivo experiments using one, two and three image pairs from an alternating bi-plane fluoroscope.

Results: The computer simulation showed that subject-specific 3D shape reconstruction using three image pairs had the best accuracy with RMSE of 0.52 ± 0.09 and 0.63 ± 0.085 mm for the femur and tibia, respectively. The corresponding values for the in vivo study were 0.64 ± 0.084 and 0.69 ± 0.069 mm, respectively, which was significantly better than those using one image pair (0.81 ± 0.126 and 0.83 ± 0.108 mm). No significant differences existed between using two and three image pairs.

Conclusion: A new two-phase optimization approach was developed for SSM-based 3D subject-specific knee model reconstructions using more than one asynchronous fluoroscopy image pair from widely available alternating bi-plane fluoroscopy systems in clinical settings. A CT-based SSM model of the knee was also developed for a healthy young Chinese male population. The new approach was found to have high mode reconstruction accuracy, and those for both two and three image pairs were much better than for a single image pair. Thus, two image pairs may be used when considering computational costs and radiation dosage. The new approach will be useful for generating patient-specific knee models for clinical applications using multiple asynchronous images from alternating bi-plane fluoroscopy widely available in clinical settings. The current SSM model will serve as a basis for further inclusion of training models with a wider range of sizes and morphological features for broader applications.

Introduction

The knee joint plays an essential role in the normal function of the lower extremities, providing stability and mobility necessary for locomotion while bearing the body’s weight. Knowledge of the biomechanical interactions of the bones and their surrounding force-bearing tissues is thus crucial for a better understanding of the functional changes of the knee under normal, pathological and post-surgical conditions (Wilson et al., 1998; Eck et al., 2002; Kleipool and Blankevoort, 2010). To this end, it is critical to be able to describe accurately individual-specific morphology of bones. These descriptions play important roles in various clinical applications such as fluoroscopy-based kinematic measurement (Lu et al., 2008), pre-surgical planning (Krekel et al., 2006), customized finite element analysis (Fernandez et al., 2004), joint implant design (Cheng et al., 1999; Harrysson et al., 2007) and computer-aided orthopaedic surgery (Gamage et al., 2009).

An accurate description of the shapes of bones is particularly important in bone model-based tracking techniques for in vivo measurement of three-dimensional (3D) joint kinematics for evaluating functional changes owing to disorders and/or treatments (Georgoulis et al., 2003; McDonald et al., 2014; Kobayashi et al., 2015), and for deriving soft tissue deformations and forces (Li et al., 2004; Haraguchi et al., 2009; Wang et al., 2013; Ding and Khan, 2019). Among the existing methods, fluoroscopy-to-CT or MRI registration techniques have been shown to be effective and less-invasive in measuring the 3D kinematics of various joints during weight-bearing functional activities (Li et al., 2008; Anderst et al., 2009; Tsai et al., 2013; Lin et al., 2014a). These techniques use subject-specific, CT- or MRI-based bone models and determine their 3D pose by searching for the final pose of the bone model whose digitally reconstructed radiograph or projections on the image plane best match the fluoroscopy image(s) (Tsai et al., 2010). Techniques using single-plane fluoroscopy have been applied to 3-D kinematic measurements of normal (Lu et al., 2008; Yamaguchi et al., 2009), pathological (Dennis et al., 2005; Hamai et al., 2009; Kobayashi et al., 2015) and replaced joints (Dennis et al., 2003; Liu et al., 2007; Yamaguchi et al., 2011), but the measured translations normal to the image plane are substantially less accurate than those of the other components (Fregly et al., 2005; Lin et al., 2014b). Bi-plane approaches address the issue by simultaneously incorporating an additional X-ray image plane, achieving much higher accuracy (Li et al., 2008; Anderst et al., 2009; Brainerd et al., 2010; Tsai et al., 2013; Lin et al., 2014b; Kapron et al., 2014; Ito et al., 2015). Generally, CT-based models are more accurate than MRI-based models but the radiation exposure is of concern. Therefore, reconstruction of three-dimensional patient-specific bone models of accuracy comparable to that of CT-based models from planar radiographs or fluoroscopy images will be useful for reducing the radiation exposure.

Statistical shape modeling (SSM) techniques have been used in the development of fully automated bone segmentation methods (Lamecker et al., 2004; Josephson et al., 2005; Fripp et al., 2006), parametric descriptions of the bony geometry (Seber, 2009) and reconstruction of subject-specific bone models using 2D/3D registration (Baka et al., 2011; Valenti et al., 2016a; Valenti et al., 2016b; van IJsseldijk et al., 2016; Yu et al., 2016; Cerveri et al., 2017; Smoger et al., 2017; Yu et al., 2017; Zheng and Yu, 2017; Fotsin et al., 2019; Cerveri et al., 2020; Wu and Mahfouz, 2021). With the SSM, a bone model is described as the mean shape superimposed by a linear combination of the principal components of shape variations of the training dataset. Reconstruction of a bone model thus involves determining the coefficients of the linear combination that best match subject-specific information such as from one or more radiographs (Lamecker et al., 2006; Baka et al., 2011; Sarkalkan et al., 2014; Karade and Ravi, 2015). So far, existing SSM models are built based mainly on Caucasian populations (Baka et al., 2011; Baka et al., 2012; Karade and Ravi, 2015; Tsai et al., 2015). However, morphological differences have been noted between Chinese and Caucasian population (Cheng et al., 1999; Mahfouz et al., 2012). For example, Cheng et al showed significant differences in the ratios of anteroposterior and mediolateral dimensions between the resected surfaces of the tibial plateau in a Chinese patient group and the tibial component of a total knee replacement designed based on the Caucasian population (Cheng et al., 1999). They suggested that between-population differences in the knee morphology directly impact the design and implantation of total knee replacements. Another study by Mahfouz et al. (2012) also identified significant differences in the mean dimensions of the three-dimensional morphology of the distal femur and proximal tibia among different ethnic populations. These previous results show that knee shape variations among multiple ethnic groups are greater than those of a single ethnic group. Therefore, an SSM dataset aiming to address more than one ethnic group would need a sample size greater than that for a single ethnic group. Considering the challenge for a dataset big enough to cover the shape variations in multiple ethnic groups and the gap in the SSM for the Chinese population, there is a need to establish a CT-based SSM for a Chinese group, which can be further expanded to include a wider range of subjects for basic research and clinical applications.

The accuracy of the model’s reconstruction using SSM can be affected by the number of radiographic images and the algorithm used. Theoretically, the more the images used, the more accurate the reconstructed bone model. Alternating bi-plane fluoroscopy systems are widely available in hospitals, providing a single pair of asynchronous fluoroscopy images. More image pairs can be obtained but movement of either the bones or the clinical fluoroscopy systems cannot be avoided, leading to temporally and spatially asynchronous images. Previous studies have mainly used a single fluoroscopy image (Fleute and Lavallée, 1999; Lamecker et al., 2006; Zheng et al., 2006; Wu and Mahfouz, 2021) or a single pair of fluoroscopy images (Laporte et al., 2003; Tang and Ellis, 2005; Baka et al., 2011; Zhu and Li, 2011; Baka et al., 2012) for 3D model reconstruction. More recent computer simulation studies used multiple image pairs for model reconstruction but considered only synchronous images (Valenti et al., 2016a; Valenti et al., 2016b; Smoger et al., 2017), so the methods proposed in these studies cannot be applied to alternating bi-plane fluoroscopy systems in clinical settings. Moreover, the accuracies reported were not representative of those in the scenarios of clinical applications. Up to the present, no study was found to evaluate the effects of image number on the accuracy of SSM-reconstructed personalized bone models of the knee to determine systematically the optimum number of planar images for personalized SSM-model reconstruction via an in vivo experimental setup.

The current study aimed to develop a new approach for 3D subject-specific knee shape reconstruction from multiple asynchronous fluoroscopy images from 2, 4, and 6 X-ray detector views using a CT-based SSM model of the knee; and to determine the optimum number of planar images for the new approach by systematically evaluating the effects of image number on reconstruction accuracy via computer simulation and in vivo data. It was hoped that the new approach would be used not only with the current SSM model but also with other existing SSM models for 3D subject-specific knee shape reconstruction from asynchronous fluoroscopy images, and that the current CT-based SSM model could form a basis for further inclusion of a wider range of subjects for basic research and clinical applications.

Materials and Methods

Statistical Shape Modeling of the Knee

The general procedure of the SSM of the knee included 1) obtaining a set of CT-derived training shape models, 2) choosing a reference model with a predefined surface mesh; 3) establishing shape (mesh) correspondence between individual training models by transforming the reference model to individual training ones; and 4) determining the mean model and primary modes of shape variations using Principal Component Analysis (PCA) (Figure 1).

FIGURE 1
www.frontiersin.org

FIGURE 1. The general procedure of statistical shape modeling (SSM) for the human knee: 1) a set of CT-derived training shape models, 2) choosing a reference model; 3) establishing shape (mesh) correspondence and alignment between individual training models; 4) determining the mean shape model and shape variations; and 5) determining primary modes of shape variations using Principal Component Analysis (PCA).

Training Shape Models of the Knee

The training shape models of the knee were reconstructed from the CT data of the distal femur and the proximal tibia from 60 healthy Chinese males (age: 22.89 ± 2 years; body height: 172.64 ± 5.2 cm; body mass: 66.35 ± 10.6 kg) who gave informed written consent as approved by the Institutional Review Board. The CT datasets were acquired for a total length of about 420 mm with a voxel size ranging from 0.441 mm × 0.441 mm × 0.625–0.822 mm × 0.822 mm × 0.625 mm (PQ-5000, Picker International, United States). The regions of the femur and tibia were segmented semi-automatically and reconstructed to obtain subject-specific training shape models using ITK-SNAP 3.6.0 (University of Pennsylvania, United States).

Shape Correspondence

Shape correspondence between individual training shape models was established by applying a self-organizing network with the Growing and Adaptive Meshes (GAMEs) algorithm (Ferrarini et al., 2007) to a randomly-selected knee model. A reference model was then selected as the model that was closest to the mean shape of the individual models. The reference model was spatially aligned to each of the training shape models via the Iterative Closet Point (ICP) method (Besl and McKay, 1992), and subsequently deformed non-rigidly to match fully with the shape model using the Coherent Point Drift (CPD) method (Myronenko and Song, 2010), yielding training shape models with corresponding meshes. The CPD method has been proven to have robust and accurate performance with respect to noise, outliers and missing points (Myronenko and Song, 2010).

Shape Alignment and Shape Variation

All the training models with corresponding meshes were best-aligned using Generalized Procrustes Analysis (GPA), minimizing the surface distances between all the models (Goodall, 1991; Cootes and Taylor, 2004). These best-aligned models were then analyzed using Principal Component Analysis (PCA) (Wold et al., 1987) to give a set of eigenvalues λi and eigenvectors ϕi, with which the form of the SSM was given as follows.

ms=m¯+i=1cbiϕi(1)

where ms is a training shape model; m¯ is the mean shape of all training models; ϕi indicates the principal modes of the shape variation following the descending order of the corresponding eigenvalues λi; bi are shape parameters bounded within an interval of [-3λi ,3λi]] (Cootes and Taylor, 2004; Sarkalkan et al., 2014); and c is the number of principal modes used for generating ms such that the ratio of the accumulated variance of eigenvalues to the total variance reaches a given level such as 0.9–0.98 (Heimann and Meinzer, 2009).

Reconstruction of Personalized Statistical Shape Model

To generate personalized bone shape model using the trained SSM with asynchronous 2D images of the bone, a two-phase optimization scheme was developed as described below (Figure 2).

FIGURE 2
www.frontiersin.org

FIGURE 2. Two-phase optimization scheme for reconstructing a personalized statistical shape model of the knee. In the current study, m = 10 and n = 10 were taken.

Digitally Reconstructed Radiograph

The two-phase optimization scheme utilized a 2D fluoroscopy to 3D model registration technique. This required the generation of digitally reconstructed radiographs (DRR) of the volumetric model of the bone onto the image planes of the fluoroscopy system. Generally, the planar detector of the fluoroscopy system was modeled as a perspective projection system, and the parameters for the positions of the X-ray sources relative to the image plane and possible image distortions were obtained via a validated calibration procedure (Tsai, 1987; Baltzopoulos, 1995; Kaptein et al., 2011). Under the perspective projection model, DRR of the bone was then produced using the ray-tracing method with trilinear interpolation (Otake et al., 2011), casting rays from a point source of x-ray to each pixel position of an image plane through the volumetric bone model. Each of these rays went through a number of voxels of the volume, the linear attenuation coefficient values of which were then integrated along the ray and projected onto the imaging plane to obtain a DRR image resembling a radiograph (Siddon, 1985; Penney et al., 2001; De Bruin et al., 2008).

Phase 1: Pose and Shape Optimization

The first phase of the new method involved the search for the optimum set of design variables to maximize the similarity between the DRRs of the volumetric bone model and the multiple asynchronous 2D fluoroscopy images (Figure 2). The design variables were the pose parameters (i.e., six degrees-of-freedom of the bone) and the coefficients for the first m principal modes of the shape variation, i.e., m shape parameters bi, i=1,m. The initial guesses of the pose parameters of the mean shape models of the femur and tibia were manually assigned separately by an operator. For generating the DRRs of the volumetric bone model, the shape model defined by the shape parameters was voxelized to simulate the CT-based bone model by filling the interior of the shape model with voxels by taking intersections of the interior of the shape model, and a set of parallel virtual 2-D slices in the transverse plane with a pixel size of 1 × 1 mm and slice thickness of 1 mm (Patil and Ravi, 2005). The resulting virtual voxels interior to the shape model were assigned a value of 700 to simulate the Hounsfield unit (HU) value of bone, while voxels outside the contours were assigned −1,000 to simulate air.

For a given set of design variables, the DRRs of the voxelized bone model in space were generated and compared with each corresponding fluoroscopy image according to the a similarity measure called Weighted Edge-Matching Score (WEMS, fw) that emphasized the alignment of longer edges between the DRRs and fluoroscopy images (Tsai et al., 2010). For the ith fluoroscopy image, the edges of the image (Ef,i) were first detected using the Canny operator (Canny, 1986), and then dilated with a given band to give Bf,i. Similarly, the edges of the DRR were also detected to give EDRR,i. The separated edges in  EDRR,i were given weightings depending on their lengths and stored in a weighting image WDRR,i. The WEMS values to be minimized were defined as follows.

fw(x,y,z,α,β,γ,b1,,bm)i=(x,y)Bf,i(x,y)WDRR,i(x,y)(x,y)WDRR,i(x,y)(2)

For p fluoroscopic images, the combined similarity measure to be minimized is defined as follows.

fw,p= i=1pfw(x,y,z,α,β,γ,b1,,bm)i(3)

The initial shape parameters b1,,bm were set to be zero because the mean training model and initial pose parameters were manually given via a graphic user interface. The resulting optimization problem was solved using a genetic algorithm (Goldberg and Holland, 1988).

Phase 2: Sequential Refined Search

In the second phase, the shape model and its pose obtained in Phase 1 were refined further in two steps. At the first step, the shape model was further refined by including n additional shape parameters, i.e., bm+1,,bm+n, for a better match with p fluoroscopy images taking the pose and shape parameters obtained in Phase 1 as fixed parameters. The similarity between the model-projected DRRs and the corresponding fluoroscopy images for this refining search was defined using a metric called the gradient correlation (GC, fGC) (Penney et al., 1998). The sum of the fGC of p fluoroscopy images was taken as the cost function to be minimized, as follows.

fshape= i=1pfGC(bm+1, ,bm+n)i(4)

The resulting optimization problem was solved using a genetic algorithm (Goldberg and Holland, 1988), giving the final shape of the bone model described by b1,,bm+n. At the second step with the final bone model, the pose parameters were further refined such that the model-projected DRRs best matched the p fluoroscopy images, minimizing the sum of the fGC of p fluoroscopy images as follows.

fpose= i=1pfGC(x,y,z,α,β,γ)i(5)

The resulting optimization problem was also solved using a genetic algorithm (Goldberg and Holland, 1988), giving the final pose of the bone model.

Evaluation of the SSM

The performance of the proposed SSM method was evaluated using computer simulation and in vivo studies as described below.

Computer Simulation

Computer simulations using a leave-one-out cross-validation scheme were performed to evaluate the proposed SSM and personalized reconstruction procedure under ideal conditions using synchronized images. Six out of the 60 subjects from the training shape model were randomly chosen one at a time, and the corresponding CT-based model was used to simulate the personalized reconstruction procedure, but was ruled out in building the SSM. For each selected subject, multiple X-ray image pairs simulating those produced from the bi-plane fluoroscopy system were generated by perspective projection of the CT-based volumetric knee model onto the image planes of the simulated fluoroscope. The CT-based volumetric knee model was positioned at the isocenter of the simulated fluoroscopy system, and three combinations of asynchronous X-ray image pairs (DRR pairs) were produced for evaluating the personalized reconstruction process: 1) one image pair (two orthogonal images) simulating the bi-plane fluoroscopy system; 2) two image pairs (four images produced by rotating the simulated bi-plane fluoroscope around the isocenter by 45°); 3) three image pairs (six images produced by rotating the simulated bi-plane fluoroscope around the isocenter by 30° and 60°) (Figure 3). The 3D personalized model of the knee was reconstructed from each of the three combinations of asynchronous image pairs using the proposed reconstruction procedure.

FIGURE 3
www.frontiersin.org

FIGURE 3. Illustration of the views of simulated X-ray images for three image pair conditions. The first image pair was obtained from the starting reference detector position. The second and third image pairs were obtained by rotating the bi-plane fluoroscope around the isocenter from the reference detector position by 30° and 60°, respectively.

In vivo Study

Ten healthy young male volunteers (age: 23 ± 2.36 years; height: 174.6 ± 4.27 cm; body mass: 63.6 ± 9.92 kg) without any neuromusculoskeletal disease or any surgical history of the lower limbs participated in the current study. The subjects were fully informed of the experimental protocol as approved by the Research Ethics Committee of China Medical University and Hospital (approval number: CMUH107-REC2-078), and gave their written consent prior to the experiment. All the subjects were scanned by CT (PQ-5000, Picker International, United States) to reconstruct the volumetric model of the knee with a voxel size of 0.709 × 0.709 × 0.625 mm. A bi-plane fluoroscopy system (Allura XPER FD 20/20, Philips Medical Systems, Netherlands) was used to acquire the dynamic X-ray images at 512 × 512 resolution during the experiment. Prior to data acquisition, intrinsic and extrinsic parameters for the bi-plane projection models of the fluoroscopy system were determined following well-established calibration procedures (Kaptein et al., 2011). Four lead markers were attached on the distal thigh and proximal shank, and one lead marker was attached on the patella. The subject was then asked to stand on a rotating plate with the knee located at the isocenter of the bi-plane imaging system. Using the plate, the lower limb was rotated sequentially by 0°, 30°, 45°, and 60° about the vertical axis while the bi-plane X-ray images were acquired. The 3D coordinates of the lead markers were determined using radiostereometric analysis for each pair of fluoroscopy images (Karrholm, 1989). Transformations between the X-ray image pairs were obtained by co-registering the known coordinates of the lead markers, from which the three combinations of image pairs were obtained as in the computer simulation study.

Error Metrics

To evaluate the performance of the proposed SSM and subject-specific model reconstruction procedure using different numbers of image pairs, the shape differences between the reconstructed model and the corresponding gold-standard CT-based model were quantified using a metric based on point-to-surface distances as follows (Cignoni et al., 1998):

e(p,S)= minpϵSd(p,p)(6)

For each point p on the surface of the SSM-generated personalized bone model, its shortest Euclidean distance e to a point p on the surface S of the CT-derived bone model was calculated, and the root-mean-squared values of e over the surface points (RMSE) was then obtained as a measure of the performance of the proposed method.

Statistical Analysis

Means and standard deviations of the RMSE were obtained for both computer simulation and in vivo study. For computer simulations, univariate analysis with a polynomial linear test was performed to determine the trend of the reconstruction accuracy variables with increasing image pairs. For the in vivo study, paired t-tests were used to compare the accuracy between image pair conditions, and between the results of Phases 1 and 2 of the two-phase optimization method. All statistical analyses were performed using SPSS (SPSS Inc., Chicago, United States) at a significance level set at 0.05.

Results

Principal Component Analysis

From the principal component analysis of the training models, it was found that the first principal component of the femur and tibia accounted for 70.7 and 58.9% of the total shape variance, respectively, while the accumulated variance of the first ten, first twenty and first thirty consecutive principal components were 85.32, 91.13, and 94.91% for the femur respectively, and 84.14, 90.45, and 94.86% for the tibia, respectively. For the femur, the first principal component was related to the deformation of the overall shape of the femur, while the second one was related to the deformation of the medial and lateral condyles, and the third one to the medial and lateral aspects of the femoral shaft (Figure 4). For the tibia, the first principal component would deform the overall shape of the tibia while the second one would deform the medial and lateral condyles and the radius of the shaft, and the third one was related to the variance in the shapes of the shaft of the tibia (Figure 4).

FIGURE 4
www.frontiersin.org

FIGURE 4. The 3D shape variation of SSMs models of (A) the femur and (B) the tibia with colors representing the difference caused by each principal component on the reference model. The first few principal components account for most of the total shape variance.

Computer Simulation

The RMSE of the reconstructed knee models using one image pair for the femur and tibia were 0.62 ± 0.075 and 0.72 ± 0.076 mm, respectively. The corresponding values for two and three image pairs were 0.57 ± 0.088 and 0.67 ± 0.074 mm, and 0.52 ± 0.09 and 0.63 ± 0.085 mm, respectively, (Table 1). The RMSE of the reconstructed knee models decreased linearly as the number of fluoroscopy image pairs increased.

TABLE 1
www.frontiersin.org

TABLE 1. The means (standard deviations) of the RMSE (unit: mm) for the femur and tibia shape reconstruction using 1, 2, and 3 image pairs by computer simulation. p values of a univariate analysis with a polynomial linear test are also given for the trend of reconstruction accuracy with increasing image pairs.

In vivo Study

The RMSE of the reconstructed knee models using one image pair for the femur and tibia were 0.81 ± 0.126 and 0.83 ± 0.108 mm, respectively, Figure 5). The corresponding values for two and three image pairs were 0.68 ± 0.088 and 0.73 ± 0.04 mm, and 0.64 ± 0.084 and 0.69 ± 0.069 mm, respectively, (Figure 5). The RMSE of both the femur and tibia for one image pair were significantly greater than those for two and three image pairs, while no significant differences existed between two and three image pairs (Figure 5). The RMSE of the reconstructed knee models were significantly reduced after Phase 2 optimization for all image pair conditions for the tibia and for 1, 2, and image pairs for the femur when compared to those with only Phase 1 optimization (Table 2). Comparisons of these results with those of previous studies are also shown in Table 3.

FIGURE 5
www.frontiersin.org

FIGURE 5. The means (standard deviations) of the RMSE (unit: mm) for the in vivo femur (A) and tibia (B) shape reconstruction using 1, 2, and 3 image pairs. The asterisks indicate significant differences (p < 0.05).

TABLE 2
www.frontiersin.org

TABLE 2. Means (standard deviations) of the RMSE (unit: mm) for the in vivo femur and tibia shape reconstruction after Phase 1 and Phase 2 using the two-phase optimization method.

TABLE 3
www.frontiersin.org

TABLE 3. Comparisons of the bone type, number of reconstructed models, number of fluoroscopy images, methods of modeling (number of PCA if using SSM), absolute mean errors, RMSE, image synchronization, and the experiment type for 3D knee shape reconstruction between published studies and the current study.

Discussion

The current study aimed to develop a new two-phase optimization approach for SSM-based 3D subject-specific knee model reconstructions using more than one asynchronous fluoroscopy image pair from widely available alternating bi-plane fluoroscopy systems in clinical settings. The approach was implemented and evaluated using a CT-based SSM model of the knee for a healthy young Chinese male population. Both computer simulation and in vivo evaluations shoCw that the new two-phase optimization approach with SSM was capable of reconstructing subject-specific knee models with high accuracy, and two or three image pairs achieved a much better accuracy than using a single image pair for both the femur and tibia. It was also found the two-phase optimization was indeed producing more accurate results than single optimization phases. The performance of the SSM model of the knee was affected by several factors, primarily the population and number of training models, shape correspondence, shape alignment, and the selection of the principal components (Baka et al., 2011; Sarkalkan et al., 2014; Tsai et al., 2015). The number of training models in the current study was higher than that of most previous studies (from 22 to 43 training models) (Baka et al., 2011; Zhu and Li, 2011; Karade and Ravi, 2015), covering greater variability of the geometric features of the joint. The current SSM of the knee joint based on a healthy young Chinese male population has been shown to produce results with high reconstruction accuracy for other subjects in the same population. This may be expected to be better than using Caucasian-based SSM because ethnic differences in knee morphology have been observed (Baka et al., 2011; Baka et al., 2012; Karade and Ravi, 2015; Tsai et al., 2015). However, further studies will be needed to confirm whether the current Chinese-based SSM would enable subject-specific reconstruction of a Caucasian knee at the same accuracy. Further inclusion of models with a wider range of sizes and morphological features will also be needed for broader applications. In the current approach, the entire pre-processing for mesh correspondence between training models was fully automated and free from any manual interventions, avoiding the possible variability in node distributions resulting from manual digitization (Zhu and Li, 2011; Tsai et al., 2015), and thus contributing to the final accuracy.

The shape modes used for representing subject-specific models were chosen considering the accuracy and computational efficiency. The first few principal components with higher eigenvalues contributed to the major shape variances, but the number of shape modes chosen may vary among SSM models depending on the population or ethnic group from which the training datasets were obtained. For example, the femur was determined with the first 20, while the tibia was determined with the first 30 by Tsai et al. (2015); Baka et al. adopted the first 30 PCA to retain 95% of the variance (Baka et al., 2011); while the first nine and first seven modes were chosen for the proximal and distal femur to account for accumulated shape variances of about 85.9 and 86.1%, respectively, by Tang et al. (Tang and Ellis, 2005). These previous SSM models were all based on Caucasian groups or a Caucasian population. In the current SSM model for the Chinese population, accumulated shape variances of about 90% were achieved with the first 20 PCA for both the femur and tibia. While the more the number of shape parameters used, the higher the accuracy a training model can be described, the choice of the number of shape parameters can have direct impact on the accuracy and efficiency of the reconstruction of a subject-specific 3D knee model from planar images.

The new two-phase optimization approach has several important features to achieve high accuracy (sub-millimeter in RMSE) and efficiency in 3D subject-specific reconstruction of the knee. Firstly, the 2-phase optimization scheme used model-generated DRRs for the 3D/2D image registration in the optimization procedure, contributing to the observed high accuracy in reconstructing subject-specific knee models (Table 1; Figure 5). This is in contrast to previous methods using mainly the contours of the model projected onto the fluoroscopy imaging plane (Laporte et al., 2003; Baka et al., 2011; Zhu and Li, 2011; Baka et al., 2014; Karade and Ravi, 2015). While the assumed homogeneous density of the bone model did not reproduce the real CT radiodensity information, the resulting attenuation of bony contours and structural overlapping on the DRRs (e.g., bilateral condyles) helped improve the structural similarities with the fluoroscopy images. Secondly, by taking the shape and bone pose parameters as design variables and the summation of the similarity measures of all the image pairs as objective function, the new 2-phase optimization approach was successful in tackling the problems of temporal and spatial asynchronization of the bones (changes in bone poses) involved in the imaging at different time instances in the reconstruction of subject-specific bone models. Using two or three pairs of fluoroscopy images from a clinical alternating bi-plane fluoroscopy system, very high accuracy with an RMSE of less than 0.74 mm was achieved for both femur and tibia. This is better than most previously reported SSM-based methods using a single synchronized image pair for the femur with RMSE ranging from 1.33 to 1.68 mm (Fleute and Lavallée, 1999; Laporte et al., 2003; Baka et al., 2011; Karade and Ravi, 2015) (Table 3). Baka et al. (2014) used multiple frames of images of the knee during particular tasks from two fixed perspective views of a synchronous bi-plane fluoroscopy system. The shape customization process would thereby involve simultaneously searching for the optimal shape parameters and the pose parameters from a great number of image frames. In contrast, the current approach acquired fluoroscopy image pairs from 4 different views (0°, 30°, 45°, and 60°) with respect to the subject’s knee by rotating the subject with a custom-made rotating plate so the x-ray image pairs from various perspectives could capture more skeletal features (Figure 3). Also, the knee images were collected in the standing posture instead of during motion tasks as in Baka et al. for images with less motion blur. The transformations among the fluoroscopic views were determined during an experimental calibration procedure, so the shape customization process involved fewer unknown parameters (i.e., shape parameters and pose parameters in one instant). All these features were considered beneficial to the shape reconstruction process. Thirdly, by taking only the first 10 shape parameters with an accumulated variance of 85% in Phase 1, the new approach enabled a relatively fast search of a first estimate of the shape and pose of the bone. This was then followed by the refined search with the second 10 shape parameters with an accumulated variance of 90% in Phase 2 for a significantly increased accuracy (Table 2). With these new features, the 2-phase optimization approach was able to produce results with sub-millimeter accuracy in the 3D knee shape reconstruction as compared to previous studies (Table 3).

The current computer simulations showed that more synchronized image pairs improve reconstruction performance. Similar outcomes were also found in previous computer simulation studies (Tang and Ellis, 2005). However, under real life conditions, the addition of an extra pair of images did not increase the reconstruction accuracy as shown by the current in vivo study (Tables 1), presumably owing to the errors arising from the temporal and spatial asynchronization of the images. The current results suggest that the reconstruction accuracy could benefit from more images as long as the temporal and spatial asynchronization of these images are taken into account in the reconstruction procedure. Considering both reconstruction quality and computing efficiency, two image pairs would be the best choice for subject-specific model reconstruction when using clinically obtained asynchronous images.

The current study established the first SSM model of the knee for a healthy young Chinese male population. With the proposed 2-phase optimization approach, the reconstruction of subject-specific knee model from two or more pairs of temporally and spatially asynchronous fluoroscopy images was made possible, and was shown to produce highly accurate results. Further inclusion of other types of knee models, such as females and older people and those with disease or deformities, in the training database would be helpful for expanding the current database for future basic research and clinical applications such as studies of the features of specific diseases or deformities of the joint. Further inclusion of Caucasian knee models in the training database or as test models would also enable quantitative comparisons of the performance of the new approach in subject-specific model reconstructions for Chinese and Caucasian subjects using different datasets. With the accuracy of the SSM model and the 2-phase reconstruction method, the proposed approach will also be useful for studying image-based knee kinematics during functional activities, as well as for clinical applications using asynchronous fluoroscopy systems.

Conclusion

A new two-phase optimization approach was developed for SSM-based 3D subject-specific knee model reconstructions using more than one asynchronous fluoroscopy image pair from widely available alternating bi-plane fluoroscopy systems in clinical settings. A CT-based SSM model of the knee was also developed for a healthy young Chinese male population. Both computer simulation and in vivo evaluations show that the new optimization approach was capable of reconstructing subject-specific knee models with high accuracy, and two or three image pairs achieved a much better accuracy than using a single image pair for both the femur and tibia. Considering computational costs, two image pairs may be preferred over three image pairs. The new approach will be useful for generating patient-specific knee models from SSM models for clinical applications using multiple asynchronous images from alternating bi-plane fluoroscopy widely available in clinical settings. The current SSM model will serve as a basis for further inclusion of training models with a wider range of sizes and morphological features for broader applications.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Requests to access the datasets should be directed to T-WL, twlu@ntu.edu.tw.

Ethics Statement

The studies involving human participants were reviewed and approved by the Research Ethics Committee China Medical University and Hospital. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

Conceptualization: T-WL and K-SS; Methodology: T-WL, K-SS, and C-CL; Data curation: H-YL and S-YL; Formal analysis and investigation: H-YL and S-YL and H-WK; Writing—original draft preparation: T-WL and H-YL; Writing—review and editing: C-CL and K-SS; Funding acquisition: T-MW and H-CH; Resources: K-SS and C-CL; Supervision: T-WL All authors have read and agreed to the published version of the article.

Funding

The authors are grateful for the financial support from the Ministry of Science and Technology, Taiwan, R.O.C. (MOST 107-2221-E-002-057 and MOST 109-2221-E-039-010).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Anderst, W., Zauel, R., Bishop, J., Demps, E., and Tashman, S. (2009). Validation of Three-Dimensional Model-Based Tibio-Femoral Tracking During Running. Med. Eng. Phys. 31 (1), 10–16. doi:10.1016/j.medengphy.2008.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Baka, N., de Bruijne, M., van Walsum, T., Kaptein, B. L., Giphart, J. E., Schaap, M., et al. (2012). Statistical Shape Model-Based Femur Kinematics from Biplane Fluoroscopy. IEEE Trans. Med. Imaging 31 (8), 1573–1583. doi:10.1109/TMI.2012.2195783

PubMed Abstract | CrossRef Full Text | Google Scholar

Baka, N., Kaptein, B. L., de Bruijne, M., van Walsum, T., Giphart, J. E., Niessen, W. J., et al. (2011). 2D-3D Shape Reconstruction of the Distal Femur from Stereo X-ray Imaging Using Statistical Shape Models. Med. Image Anal. 15 (6), 840–850. doi:10.1016/j.media.2011.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Baka, N., Kaptein, B. L., Giphart, J. E., Staring, M., Bruijne, M. d., Lelieveldt, B. P. F., et al. (2014). Evaluation of Automated Statistical Shape Model Based Knee Kinematics from Biplane Fluoroscopy. J. Biomech. 47 (1), 122–129. doi:10.1016/j.jbiomech.2013.09.022

CrossRef Full Text | Google Scholar

Baltzopoulos, V. (1995). A Videofluoroscopy Method for Optical Distortion Correction and Measurement of Knee-Joint Kinematics. Clin. Biomech. 10 (2), 85–92. doi:10.1016/0268-0033(95)92044-m

PubMed Abstract | CrossRef Full Text | Google Scholar

Besl, P. J., and McKay, N. D. (1992). “Method for Registration of 3-D Shapes,” in Sensor Fusion IV: Control Paradigms and Data Structures (Boston, MA, United States: International Society for Optics and Photonics), 586–606.

Google Scholar

Brainerd, E. L., Baier, D. B., Gatesy, S. M., Hedrick, T. L., Metzger, K. A., Gilbert, S. L., et al. (2010). X-Ray Reconstruction of Moving Morphology (XROMM): Precision, Accuracy and Applications in Comparative Biomechanics Research. J. Exp. Zool. 9999A (5), a–n. doi:10.1002/jez.589

CrossRef Full Text | Google Scholar

Canny, J. (1986). A Computational Approach to Edge Detection. IEEE Trans. Pattern Anal. Mach. Intell. PAMI-8 (6), 679–698. doi:10.1109/tpami.1986.4767851

CrossRef Full Text | Google Scholar

Cerveri, P., Belfatto, A., and Manzotti, A. (2020). Predicting Knee Joint Instability Using a Tibio-Femoral Statistical Shape Model. Front. Bioeng. Biotechnol. 8, 253. doi:10.3389/fbioe.2020.00253

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerveri, P., Sacco, C., Olgiati, G., Manzotti, A., and Baroni, G. (2017). 2D/3D Reconstruction of the Distal Femur Using Statistical Shape Models Addressing Personalized Surgical Instruments in Knee Arthroplasty: A Feasibility Analysis. Int. J. Med. Robotics Comput. Assist. Surg. 13 (4), e1823. doi:10.1002/rcs.1823

CrossRef Full Text | Google Scholar

Cheng, C.-K., Lung, C.-Y., Lee, Y.-M., and Huang, C.-H. (1999). A New Approach of Designing the Tibial Baseplate of Total Knee Prostheses. Clin. Biomech. 14 (2), 112–117. doi:10.1016/s0268-0033(98)00054-0

CrossRef Full Text | Google Scholar

Cignoni, P., Rocchini, C., and Scopigno, R. (1998). “Metro: Measuring Error on Simplified Surfaces,” in Computer Graphics Forum (Pisa, Toscana, Italy: Wiley Online Library), 167–174.

CrossRef Full Text | Google Scholar

Cootes, T. F., and Taylor, C. J. (2004). Statistical Models of Appearance for Computer Vision. University of Manchester. Technical report.

Google Scholar

De Bruin, P. W., Kaptein, B. L., Stoel, B. C., Reiber, J. H. C., Rozing, P. M., and Valstar, E. R. (2008). Image-Based RSA: Roentgen Stereophotogrammetric Analysis Based on 2D-3D Image Registration. J. Biomech. 41 (1), 155–164. doi:10.1016/j.jbiomech.2007.07.002

CrossRef Full Text | Google Scholar

Dennis, D. A., Komistek, R. D., Mahfouz, M. R., Haas, B. D., and Stiehl, J. B. (2003). Conventry Award Paper: Multicenter Determination of In Vivo Kinematics after Total Knee Arthroplasty. Clin. Orthopaedics Relat. Res. 416, 37–57. doi:10.1097/01.blo.0000092986.12414.b5

CrossRef Full Text | Google Scholar

Dennis, D. A., Mahfouz, M. R., Komistek, R. D., and Hoff, W. (2005). In Vivo Determination of Normal and Anterior Cruciate Ligament-Deficient Knee Kinematics. J. Biomech. 38 (2), 241–253. doi:10.1016/j.jbiomech.2004.02.042

CrossRef Full Text | Google Scholar

Ding, B. T. K., and Khan, S. A. (2019). The Judet Quadricepsplasty for Elderly Traumatic Knee Extension Contracture: A Case Report and Review of the Literature. Biomedicine (Taipei) 9 (3), 21. doi:10.1051/bmdcn/2019090321

PubMed Abstract | CrossRef Full Text | Google Scholar

Eck, J. C., Humphreys, S. C., Lim, T. H., Jeong, S. T., Kim, J. G., Hodges, S. D., et al. (2002). Biomechanical Study on the Effect of Cervical Spine Fusion on Adjacent-Level Intradiscal Pressure and Segmental Motion. Spine (Phila Pa 1976) 27 (22), 2431–2434. doi:10.1097/00007632-200211150-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandez, J. W., Mithraratne, P., Thrupp, S. F., Tawhai, M. H., and Hunter, P. J. (2004). Anatomically Based Geometric Modelling of the Musculo-Skeletal System and Other Organs. Biomech. Model. mechanobiology 2 (3), 139–155. doi:10.1007/s10237-003-0036-1

CrossRef Full Text | Google Scholar

Ferrarini, L., Olofsen, H., Palm, W., Vanbuchem, M., Reiber, J., and Admiraalbehloul, F. (2007). GAMEs: Growing and Adaptive Meshes for Fully Automatic Shape Modeling and Analysis. Med. image Anal. 11 (3), 302–314. doi:10.1016/j.media.2007.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Filippi, S., Motyl, B., and Bandera, C. (2008). Analysis of Existing Methods for 3D Modelling of Femurs Starting from Two Orthogonal Images and Development of a Script for a Commercial Software Package. Comput. Methods Programs Biomed. 89 (1), 76–82. doi:10.1016/j.cmpb.2007.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Fleute, M., and Lavallée, S. (1999). “Nonrigid 3-D/2-D Registration of Images Using Statistical Models,” in International Conference on Medical Image Computing and Computer-Assisted Intervention, Cambridge (Springer), 138–147.

CrossRef Full Text | Google Scholar

Fotsin, T. J. T., Vázquez, C., Cresson, T., and De Guise, J. (2019). “Shape, Pose and Density Statistical Model for 3D Reconstruction of Articulated Structures from X-ray Images,” in 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Berlin (IEEE), 2748–2751.

Google Scholar

Fregly, B. J., Rahman, H. A., and Banks, S. A. (2005). Theoretical Accuracy of Model-Based Shape Matching for Measuring Natural Knee Kinematics with Single-Plane Fluoroscopy. J. Biomechanical Eng. 127 (4), 692–699. doi:10.1115/1.1933949

CrossRef Full Text | Google Scholar

Fripp, J., Warfield, S. K., Crozier, S., and Ourselin, S. (2006). “Automatic Segmentation of the Knee Bones Using 3D Active Shape Models,” in 18th International Conference on Pattern Recognition (ICPR'06), Hong Kong (IEEE), 167–170.

CrossRef Full Text | Google Scholar

Gamage, P., Xie, S. Q., Delmas, P., and Xu, P. (2009). “3D Reconstruction of Patient Specific Bone Models from 2D Radiographs for Image Guided Orthopedic Surgery,” in 2009 Digital Image Computing: Techniques and Applications, Melbourne (IEEE), 212–216.

CrossRef Full Text | Google Scholar

Georgoulis, A. D., Papadonikolakis, A., Papageorgiou, C. D., Mitsou, A., and Stergiou, N. (2003). Three-Dimensional Tibiofemoral Kinematics of the Anterior Cruciate Ligament-Deficient and Reconstructed Knee during Walking. Am. J. Sports Med. 31 (1), 75–79. doi:10.1177/03635465030310012401

CrossRef Full Text | Google Scholar

Goldberg, D. E., and Holland, J. H. (1988). Genetic Algorithms and Machine Learning. Machine Learn. 3, 95–99. doi:10.1023/A:1022602019183

CrossRef Full Text | Google Scholar

Goodall, C. (1991). Procrustes Methods in the Statistical Analysis of Shape. J. R. Stat. Soc. Ser. B (Methodological) 53 (2), 285–321. doi:10.1111/j.2517-6161.1991.tb01825.x

CrossRef Full Text | Google Scholar

Hamai, S., Moro-oka, T.-A., Miura, H., Shimoto, T., Higaki, H., Fregly, B. J., et al. (2009). Knee Kinematics in Medial Osteoarthritis during In Vivo Weight-Bearing Activities. J. Orthop. Res. 27 (12), 1555–1561. doi:10.1002/jor.20928

CrossRef Full Text | Google Scholar

Haraguchi, N., Armiger, R. S., Myerson, M. S., Campbell, J. T., and Chao, E. Y. S. (2009). Prediction of Three-Dimensional Contact Stress and Ligament Tension in the Ankle during Stance Determined from Computational Modeling. Foot Ankle Int. 30 (2), 177–185. doi:10.3113/fai-2009-017710.3113/fai.2009.0177

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrysson, O. L., Hosni, Y. A., and Nayfeh, J. F. (2007). Custom-Designed Orthopedic Implants Evaluated Using Finite Element Analysis of Patient-specific Computed Tomography Data: Femoral-Component Case Study. BMC Musculoskelet. Disord. 8 (1), 91. doi:10.1186/1471-2474-8-91

PubMed Abstract | CrossRef Full Text | Google Scholar

Heimann, T., and Meinzer, H.-P. (2009). Statistical Shape Models for 3D Medical Image Segmentation: A Review. Med. image Anal. 13 (4), 543–563. doi:10.1016/j.media.2009.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ito, K., Hosoda, K., Shimizu, M., Ikemoto, S., Kume, S., Nagura, T., et al. (2015). Direct Assessment of 3D Foot Bone Kinematics Using Biplanar X-Ray Fluoroscopy and an Automatic Model Registration Method. J. Foot Ankle Res. 8 (1), 1–10. doi:10.1186/s13047-015-0079-4

CrossRef Full Text | Google Scholar

Josephson, K., Ericsson, A., and Karlsson, J. (2005). “Segmentation of Medical Images Using Three-Dimensional Active Shape Models,” in Scandinavian Conference on Image Analysis, Joensuu (Springer), 719–728.

CrossRef Full Text | Google Scholar

Kapron, A. L., Aoki, S. K., Peters, C. L., Maas, S. A., Bey, M. J., Zauel, R., et al. (2014). Accuracy and Feasibility of Dual Fluoroscopy and Model-Based Tracking to Quantify In Vivo Hip Kinematics during Clinical Exams. J. Appl. Biomech. 30 (3), 461–470. doi:10.1123/jab.2013-0112

CrossRef Full Text | Google Scholar

Kaptein, B. L., Shelburne, K. B., Torry, M. R., and Erik Giphart, J. (2011). A Comparison of Calibration Methods for Stereo Fluoroscopic Imaging Systems. J. Biomech. 44 (13), 2511–2515. doi:10.1016/j.jbiomech.2011.07.001

CrossRef Full Text | Google Scholar

Karade, V., and Ravi, B. (2015). 3D Femur Model Reconstruction from Biplane X-Ray Images: A Novel Method Based on Laplacian Surface Deformation. Int. J. CARS 10 (4), 473–485. doi:10.1007/s11548-014-1097-6

CrossRef Full Text | Google Scholar

Kärrholm, J. (1989). Roentgen Stereophotogrammetry: Review of Orthopedic Applications. Acta Orthopaedica Scand. 60 (4), 491–503. doi:10.3109/17453678909149328

CrossRef Full Text | Google Scholar

Kleipool, R. P., and Blankevoort, L. (2010). The Relation between Geometry and Function of the Ankle Joint Complex: A Biomechanical Review. Knee Surg. Sports Traumatol. Arthrosc. 18 (5), 618–627. doi:10.1007/s00167-010-1088-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, T., Suzuki, E., Yamazaki, N., Suzukawa, M., Akaike, A., Shimizu, K., et al. (2015). In Vivo talocrural Joint Contact Mechanics with Functional Ankle Instability. Foot & Ankle Specialist 8 (6), 445–453. doi:10.1177/1938640015585967

PubMed Abstract | CrossRef Full Text | Google Scholar

Krekel, P. R., Botha, C. P., Valstar, E. R., de Bruin, P. W., Rozing, P. M., and Post, F. H. (2006). “Interactive Simulation and Comparative Visualisation of the Bone-Determined Range of Motion of the Human Shoulder,” in Simulation und Visualisierung 2006 (SimVis 2006), Aachen, 275–288.

Google Scholar

Lamecker, H., Seebass, M., Hege, H.-C., and Deuflhard, P. (2004). “A 3D Statistical Shape Model of the Pelvic Bone for Segmentation,” in Medical Imaging 2004: Image Processing (California, CA, United States: International Society for Optics and Photonics), 1341–1351.

CrossRef Full Text | Google Scholar

Lamecker, H., Wenckebach, T. H., and Hege, H.-C. (2006). “Atlas-Based 3D-Shape Reconstruction from X-ray Images,” in 18th International Conference on Pattern Recognition (ICPR'06) (IEEE), 371–374.

CrossRef Full Text | Google Scholar

Laporte, S., Skalli, W., de Guise, J. A., Lavaste, F., and Mitton, D. (2003). A Biplanar Reconstruction Method Based on 2D and 3D Contours: Application to the Distal Femur. Comput. Methods Biomech. Biomed. Eng. 6 (1), 1–6. doi:10.1080/1025584031000065956

CrossRef Full Text | Google Scholar

Li, G., DeFrate, L. E., Sun, H., and Gill, T. J. (2004). In Vivo Elongation of the Anterior Cruciate Ligament and Posterior Cruciate Ligament during Knee Flexion. Am. J. Sports Med. 32 (6), 1415–1420. doi:10.1177/0363546503262175

CrossRef Full Text | Google Scholar

Li, G., Van de Velde, S. K., and Bingham, J. T. (2008). Validation of a Non-Invasive Fluoroscopic Imaging Technique for the Measurement of Dynamic Knee Joint Motion. J. Biomech. 41 (7), 1616–1622. doi:10.1016/j.jbiomech.2008.01.034

CrossRef Full Text | Google Scholar

Lin, C.-C., Lu, T.-W., Wang, T.-M., Hsu, C.-Y., Hsu, S.-J., and Shih, T.-F. (2014a). In Vivo Three-Dimensional Intervertebral Kinematics of the Subaxial Cervical Spine during Seated Axial Rotation and Lateral Bending via a Fluoroscopy-To-CT Registration Approach. J. Biomech. 47 (13), 3310–3317. doi:10.1016/j.jbiomech.2014.08.014

CrossRef Full Text | Google Scholar

Lin, C.-C., Lu, T.-W., Wang, T.-M., Hsu, C.-Y., and Shih, T.-F. (2014b). Comparisons of Surface vs. Volumetric Model-Based Registration Methods Using Single-Plane vs. Bi-Plane Fluoroscopy in Measuring Spinal Kinematics. Med. Eng. Phys. 36 (2), 267–274. doi:10.1016/j.medengphy.2013.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Cheng, J., Komistek, R. D., Mahfouz, M. R., and Sharma, A. (2007). In Vivo Evaluation of Dynamic Characteristics of the Normal, Fused, and Disc Replacement Cervical Spines. Spine 32 (23), 2578–2584. doi:10.1097/brs.0b013e318158cdf8

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, T.-W., Tsai, T.-Y., Kuo, M.-Y., Hsu, H.-C., and Chen, H.-L. (2008). In Vivo Three-Dimensional Kinematics of the Normal Knee during Active Extension Under Unloaded and Loaded Conditions Using Single-Plane Fluoroscopy. Med. Eng. Phys. 30 (8), 1004–1012. doi:10.1016/j.medengphy.2008.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahfouz, M., Abdel Fatah, E. E., Bowers, L. S., and Scuderi, G. (2012). Three-Dimensional Morphology of the Knee Reveals Ethnic Differences. Clin. Orthop. Relat. Res. 470 (1), 172–185. doi:10.1007/s11999-011-2089-2

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonald, C. P., Chang, V., McDonald, M., Ramo, N., Bey, M. J., and Bartol, S. (2014). Three-Dimensional Motion Analysis of the Cervical Spine for Comparison of Anterior Cervical Decompression and Fusion versus Artificial Disc Replacement in 17 Patients. J. Neurosurg. Spine 20 (3), 245–255. doi:10.3171/2013.11.SPINE13392

CrossRef Full Text | Google Scholar

Myronenko, A., and Xubo Song, X. (2010). Point Set Registration: Coherent Point Drift. IEEE Trans. Pattern Anal. Mach. Intell. 32 (12), 2262–2275. doi:10.1109/TPAMI.2010.46

PubMed Abstract | CrossRef Full Text | Google Scholar

Otake, Y., Armand, M., Armiger, R. S., Kutzer, M. D., Basafa, E., Kazanzides, P., et al. (2011). Intraoperative Image-Based Multiview 2D/3D Registration for Image-Guided Orthopaedic Surgery: Incorporation of Fiducial-Based C-Arm Tracking and GPU-Acceleration. IEEE Trans. Med. Imaging 31 (4), 948–962. doi:10.1109/TMI.2011.2176555

PubMed Abstract | CrossRef Full Text | Google Scholar

Patil, S., and Ravi, B. (2005). “Voxel-Based Representation, Display and Thickness Analysis of Intricate Shapes,” in Ninth International Conference on Computer Aided Design and Computer Graphics (CAD-CG'05), Hong Kong (IEEE), 6.

Google Scholar

Penney, G. P., Batchelor, P. G., Hill, D. L. G., Hawkes, D. J., and Weese, J. (2001). Validation of a Two- to Three-Dimensional Registration Algorithm for Aligning Preoperative CT Images and Intraoperative Fluoroscopy Images. Med. Phys. 28 (6), 1024–1032. doi:10.1118/1.1373400

PubMed Abstract | CrossRef Full Text | Google Scholar

Penney, G. P., Weese, J., Little, J. A., Desmedt, P., Hill, D. L. G., and hawkes, D. J. (1998). A Comparison of Similarity Measures for Use in 2-D-3-D Medical Image Registration. IEEE Trans. Med. Imaging 17 (4), 586–595. doi:10.1109/42.730403

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarkalkan, N., Weinans, H., and Zadpoor, A. A. (2014). Statistical Shape and Appearance Models of Bones. Bone 60, 129–140. doi:10.1016/j.bone.2013.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Seber, G. A. (2009). Multivariate Observations. Auckland, New Zealand: John Wiley & Sons.

Google Scholar

Siddon, R. L. (1985). Fast Calculation of the Exact Radiological Path for a Three-Dimensional CT Array. Med. Phys. 12 (2), 252–255. doi:10.1118/1.595715

PubMed Abstract | CrossRef Full Text | Google Scholar

Smoger, L. M., Shelburne, K. B., Cyr, A. J., Rullkoetter, P. J., and Laz, P. J. (2017). Statistical Shape Modeling Predicts Patellar Bone Geometry to Enable Stereo-Radiographic Kinematic Tracking. J. Biomech. 58, 187–194. doi:10.1016/j.jbiomech.2017.05.009

CrossRef Full Text | Google Scholar

Tang, T. S., and Ellis, R. E. (2005). “2D/3D Deformable Registration Using a Hybrid Atlas,” in International Conference on Medical Image Computing and Computer-Assisted Intervention, Palm Springs (Springer), 223–230.

CrossRef Full Text | Google Scholar

Tsai, R. (1987). A Versatile Camera Calibration Technique for High-Accuracy 3D Machine Vision Metrology Using Off-The-Shelf TV Cameras and Lenses. IEEE J. Robot. Automat. 3 (4), 323–344. doi:10.1109/jra.1987.1087109

CrossRef Full Text | Google Scholar

Tsai, T.-Y., Li, J.-S., Wang, S., Li, P., Kwon, Y.-M., and Li, G. (2015). Principal Component Analysis in Construction of 3D Human Knee Joint Models Using a Statistical Shape Model Method. Comput. Methods Biomech. Biomed. Eng. 18 (7), 721–729. doi:10.1080/10255842.2013.843676

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsai, T.-Y., Li, J.-S., Wang, S., Lin, H., Malchau, H., Li, G., et al. (2013). A Novel Dual Fluoroscopic Imaging Method for Determination of THA Kinematics: In-Vitro and In-Vivo Study. J. Biomech. 46 (7), 1300–1304. doi:10.1016/j.jbiomech.2013.02.010

CrossRef Full Text | Google Scholar

Tsai, T.-Y., Lu, T.-W., Chen, C.-M., Kuo, M.-Y., and Hsu, H.-C. (2010). A Volumetric Model-Based 2D to 3D Registration Method for Measuring Kinematics of Natural Knees with Single-Plane Fluoroscopy. Med. Phys. 37 (3), 1273–1284. doi:10.1118/1.3301596

PubMed Abstract | CrossRef Full Text | Google Scholar

Valenti, M., De Momi, E., Yu, W., Ferrigno, G., Akbari Shandiz, M., Anglin, C., et al. (2016a). Fluoroscopy-Based Tracking of Femoral Kinematics with Statistical Shape Models. Int. J. CARS 11 (5), 757–765. doi:10.1007/s11548-015-1299-6

CrossRef Full Text | Google Scholar

Valenti, M., Ferrigno, G., Martina, D., Yu, W., Zheng, G., Shandiz, M. A., et al. (2016b). Gaussian Mixture Models Based 2D-3D Registration of Bone Shapes for Orthopedic Surgery Planning. Med. Biol. Eng. Comput. 54 (11), 1727–1740. doi:10.1007/s11517-016-1460-6

PubMed Abstract | CrossRef Full Text | Google Scholar

van IJsseldijk, E. A., Valstar, E. R., Stoel, B. C., Nelissen, R. G. H. H., Baka, N., van’t Klooster, R., et al. (2016). Three Dimensional Measurement of Minimum Joint Space Width in the Knee from Stereo Radiographs Using Statistical Shape Models. Bone Jt. Res. 5 (8), 320–327. doi:10.1302/2046-3758.58.2000626

CrossRef Full Text | Google Scholar

Wang, S., Park, W. M., Gadikota, H. R., Miao, J., Kim, Y. H., Wood, K. B., et al. (2013). A Combined Numerical and Experimental Technique for Estimation of the Forces and Moments in the Lumbar Intervertebral Disc. Comput. Methods Biomech. Biomed. Eng. 16 (12), 1278–1286. doi:10.1080/10255842.2012.668537

CrossRef Full Text | Google Scholar

Wilson, D. R., Feikes, J. D., and O’Connor, J. J. (1998). Ligaments and Articular Contact Guide Passive Knee Flexion. J. Biomech. 31 (12), 1127–1136. doi:10.1016/S0021-9290(98)00119-5

CrossRef Full Text | Google Scholar

Wold, S., Esbensen, K., and Geladi, P. (1987). Principal Component Analysis. Chemometrics Intell. Lab. Syst. 2 (1-3), 37–52. doi:10.1016/0169-7439(87)80084-9

CrossRef Full Text | Google Scholar

Wu, J., and Mahfouz, M. R. (2021). Reconstruction of Knee Anatomy from Single-Plane Fluoroscopic X-ray Based on a Nonlinear Statistical Shape Model. J. Med. Imaging 8 (1), 016001. doi:10.1117/1.jmi.8.1.016001

CrossRef Full Text | Google Scholar

Yamaguchi, S., Sasho, T., Kato, H., Kuroyanagi, Y., and Banks, S. A. (2009). Ankle and Subtalar Kinematics during Dorsiflexion-Plantarflexion Activities. Foot Ankle Int. 30 (4), 361–366. doi:10.3113/fai.2009.0361

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaguchi, S., Tanaka, Y., Kosugi, S., Takakura, Y., Sasho, T., and Banks, S. A. (2011). In Vivo kinematics of Two-Component Total Ankle Arthroplasty during Non-Weightbearing and Weightbearing Dorsiflexion/Plantarflexion. J. Biomech. 44 (6), 995–1000. doi:10.1016/j.jbiomech.2011.02.078

CrossRef Full Text | Google Scholar

Yu, W., Chu, C., Tannast, M., and Zheng, G. (2016). Fully Automatic Reconstruction of Personalized 3D Volumes of the Proximal Femur from 2D X-ray Images. Int. J. CARS 11 (9), 1673–1685. doi:10.1007/s11548-016-1400-9

CrossRef Full Text | Google Scholar

Yu, W., Tannast, M., and Zheng, G. (2017). Non-Rigid Free-Form 2D-3D Registration Using a B-Spline-Based Statistical Deformation Model. Pattern recognition 63, 689–699. doi:10.1016/j.patcog.2016.09.036

CrossRef Full Text | Google Scholar

Zheng, G., Rajamani, K. T., Zhang, X., Dong, X., Styner, M., and Nolte, L.-P. (2006). “Kernel Regularized Bone Surface Reconstruction from Partial Data Using Statistical Shape Model,” in 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference, Shanghai (IEEE), 6579–6582. doi:10.1109/iembs.2005.1616008

CrossRef Full Text | Google Scholar

Zheng, G., and Yu, W. (2017). “Statistical Shape and Deformation Models Based 2D-3D Reconstruction,” in Statistical Shape and Deformation Analysis (Elsevier), 329–349. doi:10.1016/b978-0-12-810493-4.00015-8

CrossRef Full Text | Google Scholar

Zhu, Z., and Li, G. (2011). Construction of 3D Human Distal Femoral Surface Models Using a 3D Statistical Deformable Model. J. Biomech. 44 (13), 2362–2368. doi:10.1016/j.jbiomech.2011.07.006

CrossRef Full Text | Google Scholar

Keywords: statistical shape model, subject-specific, knee joint, digitally reconstructed radiographs, two-phase optimization

Citation: Lu H-Y, Shih K-S, Lin C-C, Lu T-W, Li S-Y, Kuo H-W and Hsu H-C (2021) Three-Dimensional Subject-Specific Knee Shape Reconstruction with Asynchronous Fluoroscopy Images Using Statistical Shape Modeling. Front. Bioeng. Biotechnol. 9:736420. doi: 10.3389/fbioe.2021.736420

Received: 05 July 2021; Accepted: 28 September 2021;
Published: 20 October 2021.

Edited by:

Bernardo Innocenti, Université libre de Bruxelles, Belgium

Reviewed by:

Yun Peng, NuVasive, United States
Fabio Galbusera, Galeazzi Orthopedic Institute (IRCCS), Italy

Copyright © 2021 Lu, Shih, Lin, Lu, Li, Kuo and Hsu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Tung-Wu Lu, twlu@ntu.edu.tw

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.