An integrated system for 3D hip joint reconstruction from 2D X-rays: a preliminary validation study

Comparison of 3D Bone Position Estimation Using QR Code and Metal Bead Markers

To improve the accuracy of a 3D bone position estimation system that displays 3D images in response to changes in the position of fluoroscopic images, modified markers using quick response (QR) codes were developed. The aims of this study were to assess the accuracy of the estimated bone position on 3D images with reference to QR code markers on fluoroscopic images and to compare its accuracy with metal bead markers. Bone positions were estimated from reference points on a fluoroscopic image compared with those on a 3D image. The positional relationships of QR code and metal bead markers on the fluoroscopic image were compared with those on the 3D image in order to establish whether a 3D image may be drawn by tracking positional changes in radius models. Differences were investigated by comparing the distance between markers on the fluoroscopic image and that on the 3D image, which was projected on the monitor. The error ratio, which was defined as the difference in the measurement between the fluoroscopic and 3D images divided by the fluoroscopic measurement, was compared between QR code and metal bead markers. Error ratios for the QR code markers were 5.0 ± 2.0%, 6.4 ± 7.6%, and 1.0 ± 0.8% in the anterior–posterior view, ulnar side lateral view, and posterior–anterior view, respectively. Error ratios for the metal bead markers were 1.3 ± 1.7%, 13.8 ± 14.5%, and 4.7 ± 5.7% in the anterior–posterior view, ulnar side lateral view, and posterior–anterior view, respectively. The error ratio for the metal bead markers was smaller in the initial position (p < 0.01). However, the error ratios for the QR code markers were smaller in the lateral position and the posterior–anterior position (p < 0.05). In QR code marker tracking, tracking was successful even with discontinuous images. The accuracy of a 3D bone position estimation was increased by using the QR code marker system. QR code marker tracking facilitates real-time comparisons of dynamic changes in preoperative 3D and intraoperative fluoroscopic images.

Hybrid curvature-geometrical detection of landmarks for the automatic analysis of the reduction of supracondylar fractures of the femur

BACKGROUND AND OBJECTIVE

The analysis of the features of certain tissues is required by many procedures of modern medicine, allowing the development of more efficient treatments. The recognition of landmarks allows the planning of orthopedic and trauma surgical procedures, such as the design of prostheses or the treatment of fractures. Formerly, their detection has been carried out by hand, making the workflow inaccurate and tedious. In this paper we propose an automatic algorithm for the detection of landmarks of human femurs and an analysis of the quality of the reduction of supracondylar fractures.

METHODS

The detection of anatomical landmarks follows a knowledge-based approach, consisting of a hybrid strategy: curvature and spatial decomposition. Prior training is unrequired. The analysis of the reduction quality is performed by a side-to-side comparison between healthy and fractured sides. The pre-clinical validation of the technique consists of a two-stage study: Initially, we tested our algorithm with 14 healthy femurs, comparing the output with ground truth values. Then, a total of 140 virtual fractures was processed to assess the validity of our analysis of the quality of reduction. A two-sample t test and correlation coefficients between metrics and the degree of reduction have been employed to determine the reliability of the algorithm.

RESULTS

The average detection error of landmarks was maintained below 1.7 mm and 2^∘ (p< 0.01) for points and axes, respectively. Regarding the contralateral analysis, the resulting P-values reveal the possibility to determine whether a supracondylar fracture is properly reduced or not with a 95% of confidence. Furthermore, the correlation is high between the metrics and the quality of the reduction.

CONCLUSIONS

This research concludes that our technique allows to classify supracondylar fracture reductions of the femur by only analyzing the detected anatomical landmarks. A initial training set is not required as input of our algorithm.

An Experimental Study of a 3D Bone Position Estimation System Based on Fluoroscopic Images

To compare a 3D preoperative planning image and fluoroscopic image, a 3D bone position estimation system that displays 3D images in response to changes in the position of fluoroscopic images was developed. The objective of the present study was to evaluate the accuracy of the estimated position of 3D bone images with reference to fluoroscopic images. Bone positions were estimated from reference points on a fluoroscopic image compared with those on a 3D image. The four reference markers positional relationships on the fluoroscopic image were compared with those on the 3D image to evaluate whether a 3D image may be drawn by tracking positional changes in the radius model. Intra-class correlations coefficients for reference marker distances between the fluoroscopic image and 3D image were 0.98–0.99. Average differences between measured values on the fluoroscopic image and 3D bone image for each marker corresponding to the direction of the bone model were 1.1 ± 0.7 mm, 2.4 ± 1.8 mm, 1.4 ± 0.8 mm, and 2.0 ± 1.6 mm in the anterior-posterior view, ulnar side lateral view, posterior-anterior view, and radial side lateral view, respectively. Marker positions were more accurate in the anterior-posterior and posterior-anterior views than in the radial and ulnar side lateral views. This system helps in real-time comparison of dynamic changes in preoperative 3D and intraoperative fluoroscopy images.

Rapid X-Ray-Based 3-D Finite Element Modeling of Medial Knee Joint Cartilage Biomechanics During Walking

Finite element (FE) modeling is becoming an increasingly popular method for analyzing knee joint mechanics and biomechanical mechanisms leading to osteoarthritis (OA). The most common and widely available imaging method for knee OA diagnostics is planar X-ray imaging, while more sophisticated imaging methods, e.g., magnetic resonance imaging (MRI) and computed tomography (CT), are seldom used. Hence, the capability to produce accurate biomechanical knee joint models directly from X-ray imaging would bring FE modeling closer to clinical use. Here, we extend our atlas-based framework by generating FE knee models from X-ray images (N = 28). Based on measured anatomical landmarks from X-ray and MRI, knee joint templates were selected from the atlas library. The cartilage stresses and strains of the X-ray-based model were then compared with the MRI-based model during the stance phase of the gait. The biomechanical responses were statistically not different between MRI- vs. X-ray-based models when the template obtained from X-ray imaging was the same as the MRI template. However, if this was not the case, the peak values of biomechanical responses were statistically different between X-ray and MRI models. The developed X-ray-based framework may pave the way for a clinically feasible approach for knee joint FE modeling.

Advances in FAI Imaging: a Focused Review

PURPOSE OF REVIEW

Femoroacetabular impingement (FAI) is one of the main causes of hip pain in young adults and poses clinical challenges which have placed it at the forefront of imaging and orthopedics. Diagnostic hip imaging has dramatically changed in the past years, with the arrival of new imaging techniques and the development of magnetic resonance imaging (MRI). This article reviews the current state-of-the-art clinical routine of individuals with suspected FAI, limitations, and future directions that show promise in the field of musculoskeletal research and are likely to reshape hip imaging in the coming years.

RECENT FINDINGS

The largely unknown natural disease course, especially in hips with FAI syndrome and those with asymptomatic abnormal morphologies, continues to be a problem as far as diagnosis, treatment, and prognosis are concerned. There has been a paradigm shift in recent years from bone and soft tissue morphological analysis towards the tentative development of quantitative approaches, biochemical cartilage evaluation, dynamic assessment techniques and, finally, integration of artificial intelligence (AI)/deep learning systems. Imaging, AI, and hip preserving care will continue to evolve with new problems and greater challenges. The increasing number of analytic parameters describing the hip joint, as well as new sophisticated MRI and imaging analysis, have carried practitioners beyond simplistic classifications. Reliable evidence-based guidelines, beyond differentiation into pure instability or impingement, are paramount to refine the diagnostic algorithm and define treatment indications and prognosis. Nevertheless, the boundaries of morphological, functional, and AI-aided hip assessment are gradually being pushed to new frontiers as the role of musculoskeletal imaging is rapidly evolving.

Automatic detection of landmarks for the analysis of a reduction of supracondylar fractures of the humerus

An accurate identification of bone features is required by modern orthopedics to improve patient recovery. The analysis of landmarks enables the planning of a fracture reduction surgery, designing prostheses or fixation devices, and showing deformities accurately. The recognition of these features was previously performed manually. However, this long and tedious process provided insufficient accuracy. In this paper, we propose a geometrically-based algorithm that automatically detects the most significant landmarks of a humerus. By employing contralateral images of the upper limb, a side-to-side study of the landmarks is also conducted to analyze the goodness of supracondylar fracture reductions. We conclude that a reduction can be classified by only considering the detected landmarks. In addition, our technique does not require a prior training, thus becoming a reliable alternative to treat this kind of fractures.

Multimodal molecular 3D imaging for the tumoral volumetric distribution assessment of folate-based biosensors

The aim of this study was to characterize the in vivo volumetric distribution of three folate-based biosensors by different imaging modalities (X-ray, fluorescence, Cerenkov luminescence, and radioisotopic imaging) through the development of a tridimensional image reconstruction algorithm. The preclinical and multimodal Xtreme imaging system, with a Multimodal Animal Rotation System (MARS), was used to acquire bidimensional images, which were processed to obtain the tridimensional reconstruction. Images of mice at different times (biosensor distribution) were simultaneously obtained from the four imaging modalities. The filtered back projection and inverse Radon transformation were used as main image-processing techniques. The algorithm developed in Matlab was able to calculate the volumetric profiles of ^99mTc-Folate-Bombesin (radioisotopic image), ¹⁷⁷Lu-Folate-Bombesin (Cerenkov image), and FolateRSense™ 680 (fluorescence image) in tumors and kidneys of mice, and no significant differences were detected in the volumetric quantifications among measurement techniques. The imaging tridimensional reconstruction algorithm can be easily extrapolated to different 2D acquisition-type images. This characteristic flexibility of the algorithm developed in this study is a remarkable advantage in comparison to similar reconstruction methods.

3D–2D registration in mobile radiographs: algorithm development and preliminary clinical evaluation

An image-based 3D-2D registration method is presented using radiographs acquired in the uncalibrated, unconstrained geometry of mobile radiography. The approach extends a previous method for six degree-of-freedom (DOF) registration in C-arm fluoroscopy (namely 'LevelCheck') to solve the 9-DOF estimate of geometry in which the position of the source and detector are unconstrained. The method was implemented using a gradient correlation similarity metric and stochastic derivative-free optimization on a GPU. Development and evaluation were conducted in three steps. First, simulation studies were performed that involved a CT scan of an anthropomorphic body phantom and 1000 randomly generated digitally reconstructed radiographs in posterior-anterior and lateral views. A median projection distance error (PDE) of 0.007 mm was achieved with 9-DOF registration compared to 0.767 mm for 6-DOF. Second, cadaver studies were conducted using mobile radiographs acquired in three anatomical regions (thorax, abdomen and pelvis) and three levels of source-detector distance (~800, ~1000 and ~1200 mm). The 9-DOF method achieved a median PDE of 0.49 mm (compared to 2.53 mm for the 6-DOF method) and demonstrated robustness in the unconstrained imaging geometry. Finally, a retrospective clinical study was conducted with intraoperative radiographs of the spine exhibiting real anatomical deformation and image content mismatch (e.g. interventional devices in the radiograph that were not in the CT), demonstrating a PDE = 1.1 mm for the 9-DOF approach. Average computation time was 48.5 s, involving 687 701 function evaluations on average, compared to 18.2 s for the 6-DOF method. Despite the greater computational load, the 9-DOF method may offer a valuable tool for target localization (e.g. decision support in level counting) as well as safety and quality assurance checks at the conclusion of a procedure (e.g. overlay of planning data on the radiograph for verification of the surgical product) in a manner consistent with natural surgical workflow.

Cup Implant Planning Based on 2-D/3-D Radiographic Pelvis Reconstruction-First Clinical Results

GOAL

In the following, we will present a newly developed X-ray calibration phantom and its integration for 2-D/3-D pelvis reconstruction and subsequent automatic cup planning. Two different planning strategies were applied and evaluated with clinical data.

METHODS

Two different cup planning methods were investigated: The first planning strategy is based on a combined pelvis and cup statistical atlas. Thereby, the pelvis part of the combined atlas is matched to the reconstructed pelvis model, resulting in an optimized cup planning. The second planning strategy analyzes the morphology of the reconstructed pelvis model to determine the best fitting cup implant.

RESULTS

The first planning strategy was compared to 3-D CT-based planning. Digitally reconstructed radiographs of THA patients with differently severe pathologies were used to evaluate the accuracy of predicting the cup size and position. Within a discrepancy of one cup size, the size was correctly identified in 100% of the cases for Crowe type I datasets and in 77.8% of the cases for Crowe type II, III, and IV datasets. The second planning strategy was analyzed with respect to the eventually implanted cup size. In seven patients, the estimated cup diameter was correct within one cup size, while the estimation for the remaining five patients differed by two cup sizes.

CONCLUSION

While both planning strategies showed the same prediction rate with a discrepancy of one cup size (87.5%), the prediction of the exact cup size was increased for the statistical atlas-based strategy (56%) in contrast to the anatomically driven approach (37.5%).

SIGNIFICANCE

The proposed approach demonstrated the clinical validity of using 2-D/3-D reconstruction technique for cup planning.

Comparisons of preoperative three-dimensional planning and surgical reconstruction in primary cementless total hip arthroplasty

Reconstruction of important parameters such as femoral offset and torsion is inaccurate, when templating is based on plain x-rays. We evaluate intraoperative reproducibility of pre-operative CT-based 3D-templating in a consecutive series of 50 patients undergoing primary cementless THA through an anterior approach. Pre-operative planning was compared to a postoperative CT scan by image fusion. The implant size was correctly predicted in 100% of the stems, 94% of the cups and 88% of the heads (length). The difference between the planned and the postoperative leg length was 0.3 + 2.3 mm. Values for overall offset, femoral anteversion, cup inclination and anteversion were 1.4 mm ± 3.1, 0.6° ± 3.3°, -0.4° ± 5° and 6.9° ± 11.4°, respectively. This planning allows accurate implant size prediction. Stem position and cup inclination are accurately reproducible.

X-ray image calibration and its application to clinical orthopedics

X-ray imaging is one of the most commonly used medical imaging modality. Albeit X-ray radiographs provide important clinical information for diagnosis, planning and post-operative follow-up, the challenging interpretation due to its 2D projection characteristics and the unknown magnification factor constrain the full benefit of X-ray imaging. In order to overcome these drawbacks, we proposed here an easy-to-use X-ray calibration object and developed an optimization method to robustly find correspondences between the 3D fiducials of the calibration object and their 2D projections. In this work we present all the details of this outlined concept. Moreover, we demonstrate the potential of using such a method to precisely extract information from calibrated X-ray radiographs for two different orthopedic applications: post-operative acetabular cup implant orientation measurement and 3D vertebral body displacement measurement during preoperative traction tests. In the first application, we have achieved a clinically acceptable accuracy of below 1° for both anteversion and inclination angles, where in the second application an average displacement of 8.06±3.71 mm was measured. The results of both applications indicate the importance of using X-ray calibration in the clinical routine.

The virtual skeleton database: an open access repository for biomedical research and collaboration

BACKGROUND

Statistical shape models are widely used in biomedical research. They are routinely implemented for automatic image segmentation or object identification in medical images. In these fields, however, the acquisition of the large training datasets, required to develop these models, is usually a time-consuming process. Even after this effort, the collections of datasets are often lost or mishandled resulting in replication of work.

OBJECTIVE

To solve these problems, the Virtual Skeleton Database (VSD) is proposed as a centralized storage system where the data necessary to build statistical shape models can be stored and shared.

METHODS

The VSD provides an online repository system tailored to the needs of the medical research community. The processing of the most common image file types, a statistical shape model framework, and an ontology-based search provide the generic tools to store, exchange, and retrieve digital medical datasets. The hosted data are accessible to the community, and collaborative research catalyzes their productivity.

RESULTS

To illustrate the need for an online repository for medical research, three exemplary projects of the VSD are presented: (1) an international collaboration to achieve improvement in cochlear surgery and implant optimization, (2) a population-based analysis of femoral fracture risk between genders, and (3) an online application developed for the evaluation and comparison of the segmentation of brain tumors.

CONCLUSIONS

The VSD is a novel system for scientific collaboration for the medical image community with a data-centric concept and semantically driven search option for anatomical structures. The repository has been proven to be a useful tool for collaborative model building, as a resource for biomechanical population studies, or to enhance segmentation algorithms.