1
|
Ikebe S, Shimoto T, Higaki H, Banks SA. DRR acceleration using inexpensive GPUs for model-image registration based joint kinematic measurements. J Biomech 2023; 160:111824. [PMID: 37862924 DOI: 10.1016/j.jbiomech.2023.111824] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Revised: 08/30/2023] [Accepted: 10/03/2023] [Indexed: 10/22/2023]
Abstract
Model-image registration methods are commonly used in research to measure three-dimensional joint kinematics from single-plane and bi-plane x-ray images. These methods have the potential to be beneficial if used clinically, but current techniques are too slow or expensive to be clinically practical. One technical element of these methods for measuring natural bone motion is the use of digitally reconstructed radiographs (DRRs). DRRs can be very expensive to compute, or require expensive and fast computer hardware. In this technical development, a numerically efficient Siddon-Jacobs algorithm for computing DRRs was implemented on a consumer-grade graphics card using a programming language for parallel architectures. Compared to traditional voxel projection algorithms with a central-processing-unit-only implementation, the parallel computation implementation on the graphics card provided speedups of 650-1546 times faster rendering, while retaining equivalent performance for joint kinematics measurements. The use of consumer grade graphics hardware may contribute to making model-image registration measurements of joint kinematics practical for clinical use.
Collapse
Affiliation(s)
- Satoru Ikebe
- National Institute of Technology, Kitakyushu College, Kitakyushu, Fukuoka, Japan; University of Florida, Gainesville, FL, USA.
| | | | | | | |
Collapse
|
2
|
Lu HY, Lin CC, Shih KS, Lu TW, Kuo MY, Li SY, Hsu HC. Integration of statistical shape modeling and alternating interpolation-based model tracking technique for measuring knee kinematics in vivo using clinical interleaved bi-plane fluoroscopy. PeerJ 2023; 11:e15371. [PMID: 37334125 PMCID: PMC10276557 DOI: 10.7717/peerj.15371] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 04/18/2023] [Indexed: 06/20/2023] Open
Abstract
Background A 2D fluoroscopy/3D model-based registration with statistical shape modeling (SSM)-reconstructed subject-specific bone models will help reduce radiation exposure for 3D kinematic measurements of the knee using clinical alternating bi-plane fluoroscopy systems. The current study aimed to develop such an approach and evaluate in vivo its accuracy and identify the effects of the accuracy of SSM models on the kinematic measurements. Methods An alternating interpolation-based model tracking (AIMT) approach with SSM-reconstructed subject-specific bone models was used for measuring 3D knee kinematics from dynamic alternating bi-plane fluoroscopy images. A two-phase optimization scheme was used to reconstruct subject-specific knee models from a CT-based SSM database of 60 knees using one, two, or three pairs of fluoroscopy images. Using the CT-reconstructed model as a benchmark, the performance of the AIMT with SSM-reconstructed models in measuring bone and joint kinematics during dynamic activity was evaluated in terms of mean target registration errors (mmTRE) for registered bone poses and the mean absolute differences (MAD) for each motion component of the joint poses. Results The mmTRE of the femur and tibia for one image pair were significantly greater than those for two and three image pairs without significant differences between two and three image pairs. The MAD was 1.16 to 1.22° for rotations and 1.18 to 1.22 mm for translations using one image pair. The corresponding values for two and three image pairs were 0.75 to 0.89° and 0.75 to 0.79 mm; and 0.57 to 0.79° and 0.6 to 0.69 mm, respectively. The MAD values for one image pair were significantly greater than those for two and three image pairs without significant differences between two and three image pairs. Conclusions An AIMT approach with SSM-reconstructed models was developed, enabling the registration of interleaved fluoroscopy images and SSM-reconstructed models from more than one asynchronous fluoroscopy image pair. This new approach had sub-millimeter and sub-degree measurement accuracy when using more than one image pair, comparable to the accuracy of CT-based methods. This approach will be helpful for future kinematic measurements of the knee with reduced radiation exposure using 3D fluoroscopy with clinically alternating bi-plane fluoroscopy systems.
Collapse
Affiliation(s)
- Hsuan-Yu Lu
- Department of Biomedical Engineering, National Taiwan University, Taipei, Taiwan, R.O.C.
| | - Cheng-Chung Lin
- Department of Electrical Engineering, Fu-Jen Catholic University, New Taipei, Taiwan, R.O.C.
| | - Kao-Shang Shih
- Department of Orthopedics, Shin Kong Wu Ho-Su Memorial Hospital, Taipei, Taiwan, R.O.C.
| | - Tung-Wu Lu
- Department of Biomedical Engineering, National Taiwan University, Taipei, Taiwan, R.O.C.
- Department of Orthopaedic Surgery, School of Medicine, National Taiwan University, Taipei, Taiwan, R.O.C.
| | - Mei-Ying Kuo
- Department of Physical Therapy, China Medical University, Taichung, Taiwan, R.O.C.
| | - Song-Ying Li
- Department of Biomedical Engineering, National Taiwan University, Taipei, Taiwan, R.O.C.
| | - Horng-Chaung Hsu
- Department of Orthopaedic Surgery, China Medical University, Taichung, Taiwan, R.O.C.
| |
Collapse
|
3
|
Luan Y, Wang H, Fang C, Zhang M, Li J, Zhang N, Liu B, Su J, Cheng CK. Mechanical alignment tolerance of a cruciate-retaining knee prosthesis under gait loading—A finite element analysis. Front Bioeng Biotechnol 2023; 11:1148914. [PMID: 37064225 PMCID: PMC10098169 DOI: 10.3389/fbioe.2023.1148914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Accepted: 03/20/2023] [Indexed: 04/03/2023] Open
Abstract
Component alignment is one of the most crucial factors affecting total knee arthroplasty’s clinical outcome and survival. This study aimed to investigate how coronal, sagittal, and transverse malalignment affects the mechanical behavior of the tibial insert and to determine a suitable alignment tolerance on the coronal, sagittal, and transverse planes. A finite element model of a cruciate-retaining knee prosthesis was assembled with different joint alignments (−10°, −7°, −5°, −3°, 0°, 3°, 5°, 7°, 10°) to assess the effect of malalignment under gait loading. The results showed that varus or valgus, extension, internal rotation, and excessive external rotation malalignments increased the maximum Von Mises stress and contact pressure on the tibial insert. The mechanical alignment tolerance of the studied prosthesis on the coronal, sagittal, and transverse planes was 3° varus to 3° valgus, 0°–10° flexion, and 0°–5° external rotation, respectively. This study suggests that each prosthesis should include a tolerance range for the joint alignment angle on the three planes, which may be used during surgical planning.
Collapse
Affiliation(s)
- Yichao Luan
- Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education, Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Huizhi Wang
- School of Biomedical Engineering, Shanghai Jiao Tong University, Engineering Research Center of Digital Medicine, Ministry of Education, Shanghai, China
| | - Chaohua Fang
- School of Biomedical Engineering, Shanghai Jiao Tong University, Engineering Research Center of Digital Medicine, Ministry of Education, Shanghai, China
- Department of Sport Medicine, Ningbo No 6 Hospital, Ningbo, China
| | - Min Zhang
- Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education, Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Junwei Li
- Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education, Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Ningze Zhang
- Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education, Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Bolun Liu
- Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education, Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Jian Su
- Key Laboratory of Biomechanics and Mechanobiology, Ministry of Education, Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
| | - Cheng-Kung Cheng
- School of Biomedical Engineering, Shanghai Jiao Tong University, Engineering Research Center of Digital Medicine, Ministry of Education, Shanghai, China
- *Correspondence: Cheng-Kung Cheng,
| |
Collapse
|
4
|
Tang S, Zheng L, Luo Y, Wu R, Tian Q, Wang L. Validation of an MRI Technique for the 6-DOF Knee Kinematics Measurement. Front Bioeng Biotechnol 2022; 10:904012. [PMID: 36601392 PMCID: PMC9806800 DOI: 10.3389/fbioe.2022.904012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Accepted: 06/09/2022] [Indexed: 01/07/2023] Open
Abstract
Background: For total knee arthroplasty (TKA), the optimal rotational position of the femoral component is felt to be critically important. The current knee joint kinematics measurement technology is unable to identify the exact rotation axis of the knee joint, the main reasons being low measurement accuracy and insufficient three-dimensional data (2D-3D image matching technology). In order to improve the effect of TKA surgery, we proposed a knee joint kinematics measurement method, based on the MRI technology, and verified its measurement accuracy. We then employed this method to identify the personalized optimal rotation axis of the knee joint for TKA patients. Purposes: The purpose of the study was 1) to propose a method for measuring knee joint kinematics and verify its accuracy and 2) to propose a method for determining the optimal rotation axis of knee joint for TKA surgery, based on accurate kinematic measurement results. Materials and Methods: The experiment was divided into two parts: in vitro and in vivo. The purpose of the in vitro experiment was to verify the measurement accuracy of our method. We fixed two aquarium stones (approximately 10 cm * 10 cm * 10 cm in size, close to the size of the distal femur and proximal tibia) firmly on the fixed and moving arms of the goniometer/vernier caliper with glue and immersed the aquarium stones in the water to capture MRI images. The MRI images were then processed with MATLAB software, and the relative motion of the two aquarium stones was measured. The measurement accuracy of our method was verified via the scale reading of the moving arm on the goniometer/vernier caliper. In vivo, 36 healthy elderly participants (22 females, 14 males) were recruited from the local community; our method was then employed to measure the relative motion of the tibia and femur and to observe the rollback and screw home motion of the medial/lateral condyle of the femur, which was identified as specific kinematic features of the knee joint. Results: In vitro, all measurements were accurate to <1 mm and <1°. In vivo, all knee measurements showed rollback motion (the rollback distance of the medial femoral condyle was 18.1 ± 3.7 mm and that of the lateral condyle was 31.1 ± 7.3 mm) and screw home motion. Conclusion: In the application scenario of knee joint kinematics measurement, our method has an accuracy of <1° of rotation angle and <1 mm of translation for all reference points, and it can be employed to identify the most stable axis of the knee joint. Significance: Using our method to accumulate data on the knee rotation axis of more subjects to establish an average rotation axis of a given population may help in knee prosthesis design and reduce the patient dissatisfaction rate. Individually measuring the patient's rotation axis before TKA surgery and adjusting the prosthesis installation in TKA may further reduce the patient dissatisfaction rate, and automatic computer measurement may be realized in the future, but it is still time-consuming for now.
Collapse
Affiliation(s)
- Shixiong Tang
- Department of Radiology, The Second Xiangya Hospital, Central South University, Changsha, China,Clinical Research Center for Medical Imaging in Hunan Province, Changsha, China
| | - Liwen Zheng
- Department of Rehabilitation, The Second Xiangya Hospital, Central South University, Changsha, China,*Correspondence: Liwen Zheng,
| | - Yongheng Luo
- Department of Radiology, The Second Xiangya Hospital, Central South University, Changsha, China
| | - Ren Wu
- Department of Orthopedics, The Second Xiangya Hospital, Central South University, Changsha, China
| | - Qunyan Tian
- Department of Orthopedics, The Second Xiangya Hospital, Central South University, Changsha, China
| | - Lei Wang
- School of Information Science and Engineering, Central South University, Changsha, China
| |
Collapse
|
5
|
Marsilio L, Faglia A, Rossi M, Mainardi L, Manzotti A, Cerveri P. CEL-Unet: a novel CNN architecture for 3D Segmentation of Knee Bones affected by Severe Osteoarthritis for PSI-Based Surgical Planning. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2022; 2022:5039-5042. [PMID: 36085733 DOI: 10.1109/embc48229.2022.9871953] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Unet architectures are promising deep learning networks exploited to perform the automatic segmentation of bone CT images, in line with their ability to deal with pathological deformations and size-varying anatomies. However, bone degeneration, like the development of irregular osteophytes as well as mineral density alterations might interfere with this automated process and demand extensive manual refinement. The aim of this work is to implement an innovative Unet variant, the CEL-Unet, to improve the femur and tibia segmentation outcomes in osteoarthritic knee joints. In this network the decoding path is split into a region and contour-aware branch to increase the prediction reliability in such pathological conditions. The comparison between the segmentation results achieved with a standard Unet and its novel variant (CEL-Unet) was performed as follows: the Unet was trained with 5 different loss functions: Dice Loss, Focal Loss, Exponential Logarithmic Loss, Double Cross Entropy Loss and Distanced Cross Entropy loss. The CEL-Unet was instead trained with two loss functions, one for each of the network outputs, namely Mask and Edge, yielding the so-called Combined Edge Loss (CEL) function. A set of 259 knee CT scans was used to train the model and test segmentation performance. The CEL-Unet outperformed all other Unet-based models, reaching the highest Jaccard values of about 0.97 and 0.96 on femur and tibia, respectively. Clinical Relevance- With the increasing rate of Total Knee Arthoplasty deep learning-based methods can achieve fast accurate and automatic 3D segmentation of the knee joint bones to enhance new costumized pre-operative planning.
Collapse
|