Differences in trabecular bone texture between knees with and without radiographic osteoarthritis detected by fractal methods

Subchondral bone: An emerging target for the treatment of articular surface lesions of the knee

Purpose

When dealing with the health status of the knee articular surface, the entire osteochondral unit has gained increasing attention, and in particular the subchondral bone, which plays a key role in the integrity of the osteochondral unit. The aim of this article was to discuss the current evidence on the role of the subchondral bone.

Methods

Experts from different geographical regions were involved in performing a review on highly discussed topics about the subchondral bone, ranging from its etiopathogenetic role in joint degeneration processes to its prognostic role in chondral and osteochondral defects, up to treatment strategies to address both the subchondral bone and the articular surface.

Discussion

Subchondral bone has a central role both from an aetiologic point of view and as a diagnostic tool, and its status was found to be relevant also as a prognostic factor in the follow-up of chondral treatment. Finally, the recognition of its importance in the natural history of these lesions led to consider subchondral bone as a treatment target, with the development of osteochondral scaffolds and procedures to specifically address osteochondral lesions.

Conclusion

Subchondral bone plays a central role in articular surface lesions from different points of view. Several aspects still need to be understood, but a growing interest in subchondral bone is to be expected in the upcoming future towards the optimization of joint preservation strategies.

Level of Evidence

Level V, expert opinion.

Evaluation of Bone Mineral Changes in Panoramic Radiographs of Hypothyroid and Hyperthyroid Patients Using Fractal Dimension Analysis

Objective Hyperthyroidism and hypothyroidism are endocrinopathies that cause a decrease in bone mineral density. The aim of this study is to investigate possible bone changes in the mandible caused by hyperthyroidism and hypothyroidism using fractal analysis (FA) on panoramic radiographs. Material and Methods Panoramic radiographs of a total of 180 patients, including 120 patient groups (60 hyperthyroid, 60 hypothyroid) and 60 healthy control groups, were used. Five regions of interests (ROI) were determined from panoramic radiographs and FA was performed. ROI1: geometric midpoint of mandibular notch and mandibular foramen, ROI2: geometric midpoint of mandibular angle, ROI3: anterior of mental foramen, ROI4: basal cortical area from distal mental foramen to distal root of first molar, ROI5: geometric center of mandibular foramen and mandibular ramus. Results While a significant difference was observed between the patient and control groups regarding ROI1 and ROI2 (p < 0.05); there was no significant difference between the groups in relation to ROI3, ROI4, and ROI5. All FA values were lower in the hyperthyroid group than in the hypothyroid group. Conclusion Fractal analysis proves to be an effective method for early detection of bone mass changes. In the present study, it was concluded that while the mandibular cortical bone was intact, trabecular rich regions were affected by osteoporosis caused by thyroid hormones. Necessary precautions should be taken against the risk of osteoporosis in patients with thyroid hormone disorders.

Trabecular bone texture analysis of conventional radiographs in the assessment of knee osteoarthritis: review and viewpoint

BACKGROUND

Trabecular bone texture analysis (TBTA) has been identified as an imaging biomarker that provides information on trabecular bone changes due to knee osteoarthritis (KOA). Consequently, it is important to conduct a comprehensive review that would permit a better understanding of this unfamiliar image analysis technique in the area of KOA research. We examined how TBTA, conducted on knee radiographs, is associated to (i) KOA incidence and progression, (ii) total knee arthroplasty, and (iii) KOA treatment responses. The primary aims of this study are twofold: to provide (i) a narrative review of the studies conducted on radiographic KOA using TBTA, and (ii) a viewpoint on future research priorities.

METHOD

Literature searches were performed in the PubMed electronic database. Studies published between June 1991 and March 2020 and related to traditional and fractal image analysis of trabecular bone texture (TBT) on knee radiographs were identified.

RESULTS

The search resulted in 219 papers. After title and abstract scanning, 39 studies were found eligible and then classified in accordance to six criteria: cross-sectional evaluation of osteoarthritis and non-osteoarthritis knees, understanding of bone microarchitecture, prediction of KOA progression, KOA incidence, and total knee arthroplasty and association with treatment response. Numerous studies have reported the relevance of TBTA as a potential bioimaging marker in the prediction of KOA incidence and progression. However, only a few studies have focused on the association of TBTA with both OA treatment responses and the prediction of knee joint replacement.

CONCLUSION

Clear evidence of biological plausibility for TBTA in KOA is already established. The review confirms the consistent association between TBT and important KOA endpoints such as KOA radiographic incidence and progression. TBTA could provide markers for enrichment of clinical trials enhancing the screening of KOA progressors. Major advances were made towards a fully automated assessment of KOA.

Investigation of mandibular fractal dimension on digital panoramic radiographs in bruxist individuals

OBJECTIVE

This study aimed to evaluate changes in mandibular trabecular bone structure in bruxism using fractal analysis on digital panoramic radiographs obtained with automatic dosing.

STUDY DESIGN

In this prospective study, fractal analysis was performed on radiographs of 126 bruxists and 126 non-bruxists. Eight paired mandibular regions of interest were selected: the bilateral condylar and gonial regions, and the bilateral dentate regions between the apical areas of the first molar and second premolar and between the first premolar and canine. Fractal dimensions (FDs) were calculated at each site.

RESULTS

Mean FD values in the bilateral gonial regions of the bruxists were significantly lower than those of controls (P ≤ .049). In both groups, FD values of the right dentate region anterior to the mental foramen were significantly lower than those on the left side (P ≤ .042). Females exhibited significantly lower FD values in both condylar regions in both groups (P ≤ .039) and in the right dentate regions in the controls (P ≤ .022). Correlations between age and FD in all regions were positive but nonsignificant in both groups (P > .05).

CONCLUSIONS

FD values of mandibular trabecular bone are affected by bruxism in the gonial region and by laterality and sex differences in the condylar and dentate regions.

Evaluation of the mandibular trabecular bone in patients with bruxism using fractal analysis

OBJECTIVE

The aims of this study were (1) to investigate the effect of bruxism on the fractal dimension (FD) of the mandibular trabecular bone through digital panoramic radiographs, and (2) to evaluate the effectiveness of fractal analysis as a diagnostic test for bruxism.

METHODS

One hundred and six bruxer and 106 non-bruxer patients were included in the study. Three bilateral regions of interest (ROI) were selected: ROI-1, the mandibular condyle; ROI-2, the mandibular angle; ROI-3, the-area between the apical regions of the mandibular second premolar and the first molar teeth. FD values for the bruxer and non-bruxer groups were compared for each ROI.

RESULTS

Only the FD measurements for the right mandibular condyle (ROI-1) showed a statistically significant difference (p = 0.041) between the bruxer and non-bruxer individuals. FD values measured in the bruxers (1.40 ± 0.09) were lower than in the non-bruxers (1.42 ± 0.08).

CONCLUSION

Fractal analysis may be a useful method for discerning trabecular differences in the condylar areas of bruxer individuals. In future studies, the unilateral mastication habits, the characteristics of dental wear, and the occlusal bite forces of individuals should be documented.

Subchondral bone deterioration in femoral heads in patients with osteoarthritis secondary to hip dysplasia: A case-control study

Objectives

Residual hip dysplasia is the most common underlying condition leading to secondary osteoarthritis (OA) of the hip. Subchondral bone alterations in OA secondary to hip dysplasia (HD-OA) are poorly investigated. The aim of the present study was to analyse the microarchitecture, bone remodelling and pathological alterations of subchondral bone in femoral heads from patients with HD-OA.

Methods

Subchondral bone specimens were extracted from both weight-bearing and non–weight-bearing regions of femoral heads from 20 patients with HD-OA and 20 patients with osteoporotic femoral neck fracture, during hip replacement surgery. Micro-CT and histological examination were performed to assess the microarchitecture and histopathological changes.

Results

The weight-bearing subchondral bone showed significantly more sclerotic microarchitecture and higher bone remodelling level in HD-OA as compared with osteoporosis. In the non–weight-bearing region, the two diseases shared similar microarchitectural characteristics, but higher bone remodelling level was detected in HD-OA. Distinct regional differences were observed in HD-OA, whereas the two regions exhibited similar characteristics in osteoporosis. In addition, HD-OA displayed more serious pathological alterations, including subchondral bone cyst, metaplastic cartilaginous tissue, bone marrow oedema and fibrous tissue, especially in the weight-bearing region.

Conclusions

Osteoarthritic deteriorations of subchondral bone induced by hip dysplasia spread throughout the whole joint, but exhibit region-dependent variations, with the weight-bearing region more seriously affected. Biomechanical stress might exert a pivotal impact on subchondral bone homeostasis in hip dysplasia.

The translational potential of this article

The histomorphometric findings in the project indicate an early intervention for the development of hip dysplasia in clinic.

Bone Density and Texture from Minimally Post-Processed Knee Radiographs in Subjects with Knee Osteoarthritis

Plain radiography is the most common modality to assess the stage of osteoarthritis. Our aims were to assess the relationship of radiography-based bone density and texture between radiographs with minimal and clinical post-processing, and to compare the differences in bone characteristics between controls and subjects with knee osteoarthritis or medial tibial bone marrow lesions (BMLs). Tibial bone density and texture was evaluated from radiographs with both minimal and clinical post-processing in 109 subjects with and without osteoarthritis. Bone texture was evaluated using fractal signature analysis. Significant correlations (p < 0.001) were found in all regions (between 0.94 and 0.97) for calibrated bone density between radiographs with minimal and clinical post-processing. Correlations varied between 0.51 and 0.97 (p < 0.001) for FD_Ver texture parameter and between − 0.10 and 0.97 for FD_Hor. Bone density and texture were different (p < 0.05) between controls and subjects with osteoarthritis or BMLs mainly in medial tibial regions. When classifying healthy and osteoarthritic subjects using a machine learning-based elastic net model with bone characteristics, area under the receiver operating characteristics (ROCAUC) curve was 0.77. For classifying controls and subjects with BMLs, ROCAUC was 0.85. In conclusion, differences in bone density and texture can be assessed from knee radiographs when using minimal post-processing.

Comparison of bone texture between normal individuals and patients with Kashin-Beck disease from plain radiographs in knee

To compare tibial bone texture between Kashin-Beck disease (KBD) patients and normal individuals from plain radiographs using an advanced image analysis. Plain knee radiographs were obtained from KBD patients (n = 49) and age-matched healthy controls (n = 98). KBD were graded with diagnostic criteria WS/T 207-2010. The textural values related to bone structure from medial and lateral tibial subchondral and trabecular bones were evaluated using entropy of Laplacian-based image (E_Lap), entropy of local binary patterns (E_LBP), homogeneity indices (HI) of local angles (HI_Mean, HI_Perp and HI_Paral), and fractal dimensions from horizontal (FD_Hor) and vertical (FD_Ver) structures. KBD patients were shorter in height and lighter in weight, and their tibial width was wider than controls. Anatomical angle of KBD patients showed more genu valgus. Total KBD patients and subgroups had higher E_Lap, HI_Mean, HI_Perp and HI_Paral in detected tibial subchondral and trabecular bones than controls, except E_Lap in lateral subchondral bone. E_LBP, FD_Hor and FD_Ver from the detected tibial bone in KBD patients and subgroups were lower than controls, except FD_Ver in lateral trabecular bone. Our results indicate that micro-scale in bone texture in KBD-affected knees can be quantitatively examined from plain radiographs using an advanced image analysis.

Cancellous bone and theropod dinosaur locomotion. Part I-an examination of cancellous bone architecture in the hindlimb bones of theropods

This paper is the first of a three-part series that investigates the architecture of cancellous ('spongy') bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and has previously been used to infer locomotor biomechanics in extinct tetrapod vertebrates, especially primates. Despite great promise, cancellous bone architecture has remained little utilized for investigating locomotion in many other extinct vertebrate groups, such as dinosaurs. Documentation and quantification of architectural patterns across a whole bone, and across multiple bones, can provide much information on cancellous bone architectural patterns and variation across species. Additionally, this also lends itself to analysis of the musculoskeletal biomechanical factors involved in a direct, mechanistic fashion. On this premise, computed tomographic and image analysis techniques were used to describe and analyse the three-dimensional architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs for the first time. A comprehensive survey across many extant and extinct species is produced, identifying several patterns of similarity and contrast between groups. For instance, more stemward non-avian theropods (e.g. ceratosaurs and tyrannosaurids) exhibit cancellous bone architectures more comparable to that present in humans, whereas species more closely related to birds (e.g. paravians) exhibit architectural patterns bearing greater similarity to those of extant birds. Many of the observed patterns may be linked to particular aspects of locomotor biomechanics, such as the degree of hip or knee flexion during stance and gait. A further important observation is the abundance of markedly oblique trabeculae in the diaphyses of the femur and tibia of birds, which in large species produces spiralling patterns along the endosteal surface. Not only do these observations provide new insight into theropod anatomy and behaviour, they also provide the foundation for mechanistic testing of locomotor hypotheses via musculoskeletal biomechanical modelling.

Subchondral bone in osteoarthritis: association between MRI texture analysis and histomorphometry

OBJECTIVE

Magnetic resonance imaging (MRI) texture analysis is a method of analyzing subchondral bone alterations in osteoarthritis (OA). The objective of this study was to evaluate the association between MR texture analysis and ground-truth subchondral bone histomorphometry at the tibial plateau.

DESIGN

The local research ethics committee approved the study. All subjects provided written, informed consent. This was a cross-sectional study carried out at our institution between February and August 2014. Ten participants aged 57-84 with knee OA scheduled for total knee arthroplasty (TKA) underwent pre-operative MRI of the symptomatic knee at 3T using a high spatial-resolution coronal T1 weighted sequence. Tibial plateau explants obtained at the time of TKA underwent histological preparation to allow calculation of bone volume fraction (BV.TV), trabecular thickness (Tb.Th), trabecular separation (Tb.Sp) and trabecular number (Tb.N). Texture analysis was performed on the tibial subchondral bone of MRI images matched to the histological sections. Regression models were created to assess the association of texture analysis features with BV.TV, Tb.Th, Tb.Sp and Tb.N.

RESULTS

MRI texture features were significantly associated with BV.TV (R² = 0.76), Tb.Th (R² = 0.47), Tb.Sp (R² = 0.75) and Tb.N (R² = 0.60, all P < 0.001). Simple gray-value histogram based texture features demonstrated the highest standardized regression coefficients for each model.

CONCLUSION

MRI texture analysis features were significantly associated with ground-truth subchondral bone histomorphometry at the tibial plateau.

The Use of Tomosynthesis in the Global Study of Knee Subchondral Insufficiency Fractures

RATIONALE AND OBJECTIVES

Subchondral insufficiency fractures (SIF), previously termed spontaneous osteonecrosis of the knee, are marked by a sudden onset of severe pain. Other than the size of the lesion, prediction for progression to joint replacement is difficult. The objective was to determine if quantitative analysis of bone texture using digital tomosynthesis imaging would be useful in predicting more rapid progression to joint replacement.

MATERIALS AND METHODS

Tomosynthesis studies of 30 knees with documented SIF were quantified by fractal, mean intercept length (MIL), and line fraction deviation analyses. Fractal dimension, lacunarity, MIL, and line fraction deviation variables measured from these analyses were then correlated to short interval progression to joint replacement surgery.

RESULTS

Higher odds for joint replacement were related to higher values of the standard deviation of slope lacunarity and to morphometric measures (eg, MIL).

CONCLUSIONS

Using digital tomosynthesis images for bone texture assessment may help distinguish condylar bone response in SIF, potentially acting as a clinically relevant predictive tool. In the future, contrasting SIF to the more gradual long-term process of osteoarthritis, there may be a better understanding of the different mechanisms for the two conditions.

Subchondral tibial bone texture analysis predicts knee osteoarthritis progression: data from the Osteoarthritis Initiative: Tibial bone texture & knee OA progression

OBJECTIVES

To examine whether trabecular bone texture (TBT) parameters assessed on computed radiographs could predict knee osteoarthritis (OA) progression.

METHODS

This study was performed using data from the Osteoarthritis Initiative (OAI). 1647 knees in 1124 patients had bilateral fixed flexion radiographs acquired 48 months apart. Images were semi-automatically segmented to extract a patchwork of regions of interest (ROI). A fractal texture analysis was performed using different methods. OA progression was defined as an increase in the joint space narrowing (JSN) over 48 months. The predictive ability of TBT was evaluated using logistic regression and receiver operating characteristic (ROC) curve. An optimization method for features selection was used to reduce the size of models and assess the impact of each ROI.

RESULTS

Fractal dimensions (FD's) were predictive of the JSN progression for each method tested with an area under the ROC curve (AUC) up to 0.71. Baseline JSN grade was not correlated with TBT parameters (R < 0.21) but had the same predictive capacity (AUC 0.71). The most predictive model included the clinical covariates (age, gender, body mass index (BMI)), JSN and TBT parameters (AUC 0.77). From a statistical point of view we found higher differences in TBT parameters computed in medial ROI between progressors and non-progressors. However, the integration of TBT results from the whole patchwork including the lateral ROIs in the model provided the best predictive model.

CONCLUSIONS

Our findings indicate that TBT parameters assessed in different locations in the joint provided a good predictive ability to detect knee OA progression.

Advanced Knee Structure Analysis (AKSA): a comparison of bone mineral density and trabecular texture measurements using computed tomography and high-resolution peripheral quantitative computed tomography of human knee cadavers

Background

A change of loading conditions in the knee causes changes in the subchondral bone and may be a cause of osteoarthritis (OA). However, quantification of trabecular architecture in vivo is difficult due to the limiting spatial resolution of the imaging equipment; one approach is the use of texture parameters. In previous studies, we have used digital models to simulate changes of subchondral bone architecture under OA progression. One major result was that, using computed tomography (CT) images, subchondral bone mineral density (BMD) in combination with anisotropy and global homogeneity could characterize this progression.

The primary goal of this study was a comparison of BMD, entropy, anisotropy, variogram slope, and local and global inhomogeneity measurements between high-resolution peripheral quantitative CT (HR-pQCT) and CT using human cadaveric knees. The secondary goal was the verification of the spatial resolution dependence of texture parameters observed in the earlier simulations, two important prerequisites for the interpretation of in vivo measurements in OA patients.

Method

The applicability of texture analysis to characterize bone architecture in clinical CT examinations was investigated and compared to results obtained from HR-pQCT. Fifty-seven human knee cadavers (OA status unknown) were examined with both imaging modalities. Three-dimensional (3D) segmentation and registration processes, together with automatic positioning of 3D analysis volumes of interest (VOIs), ensured the measurement of BMD and texture parameters at the same anatomical locations in CT and HR-pQCT datasets.

Results

According to the calculation of dice ratios (>0.978), the accuracy of VOI locations between methods was excellent. Entropy, anisotropy, and global inhomogeneity showed significant and high linear correlation between both methods (0.68 < R² < 1.00). The resolution dependence of these parameters simulated earlier was confirmed by the in vitro measurements.

Conclusion

The high correlation of HR-pQCT- and CT-based measurements of entropy, global inhomogeneity, and anisotropy suggests interchangeability between devices regarding the quantification of texture. The agreement of the experimentally determined resolution dependence of global inhomogeneity and anisotropy with earlier simulations is an important milestone towards their use to quantify subchondral bone structure. However, an in vivo study is still required to establish their clinical relevance.

Automated selection of bone texture regions on hand radiographs: Data from the Osteoarthritis Initiative

Manual selection of finger trabecular bone texture regions on hand X-ray images is time-consuming, tedious, and observer-dependent. Therefore, we developed an automated method for the region selection. The method selects square trabecular bone regions of interest above and below the second to fifth distal and proximal interphalangeal joints. Two regions are selected per joint (16 regions per hand). The method consists of four integral parts: (1) segmentation of a radiograph into hand and background, (2) identification of finger regions, (3) localization of center points of heads of distal phalanges and the distal interphalangeal, proximal interphalangeal, and metacarpophalangeal joints, and (4) placement of the regions of interest under and above the distal and proximal interphalangeal joints. A gold standard was constructed from regions selected by two observers on 40 hand X-ray images taken from Osteoarthritis Initiative cohort. Datasets of 520 images were generated from the 40 images to study the effects of hand and finger positioning. The accuracy in regions selection and the agreement in calculating five directional fractal parameters were evaluated against the gold standard. The accuracy, agreement, and effects of hand and finger positioning were measured using similarity index (0 for no overlap and 1 for entire overlap) and interclass correlation coefficient as appropriate. A high accuracy in selecting regions (similarity index ≥ 0.79) and a good agreement in fractal parameters (interclass correlation coefficient ≥ 0.58) were achieved. Hand and finger positioning did not affect considerably the region selection (similarity index ≥ 0.70). These results indicate that the method developed selects bone regions on hand X-ray images with accuracy sufficient for fractal analyses of bone texture.

Knee joint subchondral bone structure alterations in active athletes: a cross-sectional case-control study

OBJECTIVE

It has been shown that trabecular bone structure parameters extracted from radiographs known as fractal signature analysis (FSA) are able to predict structural outcomes such as radiographic osteoarthritis (OA) progression. Little is known about their involvement in early disease or about differences between subjects exposed to increased joint loading such as young active athletes compared to non-athletes. Aim was to compare horizontal and vertical dimensions of bone texture considering athlete status, gender, previous anterior cruciate ligament (ACL) surgery and age.

DESIGN

Included were 685 patients of which 135 consecutive athletes (82% soccer players) 18-36 years old and 550 non-athletes controls in the same age range had knee radiography for assessment of subacute or chronic knee complaints. Regions of interest (ROI) were placed in the subchondral medial and lateral tibial plateaus. Fractal signatures were calculated in the horizontal and vertical dimensions. Curve fitting algorithms were applied taking into account all four risk factors in the same model adjusting for each other.

RESULTS

For the horizontal dimensions significant differences were observed for gender (estimate (E) 0.098 (95% confidence interval(CI)) (-0.009, 0.008), P < .0001), previous ACL surgery (E -0.031, 95% CI (-0.043, -0.019), P < .0001) and highest age group (E -0.039, 95% CI (-0.048, -0.029), P < .0001). For vertical dimensions, significant differences were shown for athletes (E -0.012, 95% CI (-0.020, -0.004), P < .0001), gender (E 0.056, 95% CI (0.049, 0.062), P < .0001), and age range from 28 to 32 years (E -0.028, 95% CI (-0.037, -0.019), P < .0001).

CONCLUSIONS

Trabecular bone structure differs between athletes and non-athletes, in regard to previous ACL surgery, gender and higher age.

Correlation of Subchondral Bone Density and Structure from Plain Radiographs with Micro Computed Tomography Ex Vivo

Osteoarthritis causes changes in the subchondral bone structure and composition. Plain radiography is a cheap, fast, and widely available imaging method. Bone tissue can be well seen from plain radiograph, which however is only a 2D projection of the actual 3D structure. Therefore, the aim was to investigate the relationship between bone density- and structure-related parameters from 2D plain radiograph and 3D bone parameters assessed from micro computed tomography (µCT) ex vivo. Right tibiae from eleven cadavers without any diagnosed joint disease were imaged using radiography and with µCT. Bone density- and structure-related parameters were calculated from four different locations from the radiographs of proximal tibia and compared with the volumetric bone microarchitecture from the corresponding regions. Bone density from the plain radiograph was significantly related with the bone volume fraction (r = 0.86; n = 44; p < 0.01). Mean homogeneity index for orientation of local binary patterns (HI_angle,mean) and fractal dimension of vertical structures (FD_Ver) were related (p < 0.01) with connectivity density (HI_angle,mean: r = −0.73, FD_Ver: r = 0.69) and trabecular separation (HI_angle,mean: r = 0.73, FD_Ver: r = −0.70) when all ROIs were pooled together (n = 44). Bone density and structure in tibia from standard clinically available 2D radiographs are significantly correlated with true 3D microstructure of bone.

Multi-dimensional reliability assessment of fractal signature analysis in an outpatient sports medicine population

The aim of this study has been to test reproducibility of fractal signature analysis (FSA) in a young, active patient population taking into account several parameters including intra- and inter-reader placement of regions of interest (ROIs) as well as various aspects of projection geometry. In total, 685 patients were included (135 athletes and 550 non-athletes, 18-36 years old). Regions of interest (ROI) were situated beneath the medial tibial plateau. The reproducibility of texture parameters was evaluated using intraclass correlation coefficients (ICC). Multi-dimensional assessment included: (1) anterior-posterior (A.P.) vs. posterior-anterior (P.A.) (Lyon-Schuss technique) views on 102 knees; (2) unilateral (single knee) vs. bilateral (both knees) acquisition on 27 knees (acquisition technique otherwise identical; same A.P. or P.A. view); (3) repetition of the same image acquisition on 46 knees (same A.P. or P.A. view, and same unitlateral or bilateral acquisition); and (4) intra- and inter-reader reliability with repeated placement of the ROIs in the subchondral bone area on 99 randomly chosen knees. ICC values on the reproducibility of texture parameters for A.P. vs. P.A. image acquisitions for horizontal and vertical dimensions combined were 0.72 (95% confidence interval (CI) 0.70-0.74) ranging from 0.47 to 0.81 for the different dimensions. For unilateral vs. bilateral image acquisitions, the ICCs were 0.79 (95% CI 0.76-0.82) ranging from 0.55 to 0.88. For the repetition of the identical view, the ICCs were 0.82 (95% CI 0.80-0.84) ranging from 0.67 to 0.85. Intra-reader reliability was 0.93 (95% CI 0.92-0.94) and inter-observer reliability was 0.96 (95% CI 0.88-0.99). A decrease in reliability was observed with increasing voxel sizes. Our study confirms excellent intra- and inter-reader reliability for FSA, however, results seem to be affected by acquisition technique, which has not been previously recognized.

Trabecular orientation in the human femur and tibia and the relationship with lower-limb alignment for patients with osteoarthritis of the knee

Wolff׳s Law suggests that the orientation of trabeculae in human bone changes in response to altered loading patterns. The aim of this study was to investigate trabecular orientation in both the femur and tibia and to compare this with the mechanical axis of the leg. The study involved analysis of radiographs from patients with osteoarthritis of the knee (n=91). For each patient, the trabecular orientation in both the distal femur and proximal tibia was measured from a standard anteroposterior radiograph of the knee and the mechanical axis of the leg was calculated from a long leg view taken while weight bearing. There was a significant correlation between the mechanical axis and the trabecular orientation in each of the regions considered in the femur (r=-0.41, -0.30, 0.52, and 0.23) and tibia (r=-0.27 and 0.31). Multiple regression analysis, with mechanical axis as the dependent variable, produced an R(2) of 0.62. Greater trabecular anisotropy (i.e. greater alignment) was observed in the medial femur and tibia compared to the lateral side (p<0.01). The results give an insight into the trabecular changes that may take place during development of osteoarthritis and following surgery. In particular, we propose that the orientation of the trabeculae in both the distal femur and proximal tibia will reflect the angle of mechanical loading through the knee.

Joint unloading implant modifies subchondral bone trabecular structure in medial knee osteoarthritis: 2-year outcomes of a pilot study using fractal signature analysis

Background

Knee osteoarthritis (OA) is largely attributable to chronic excessive and aberrant joint loading. The purpose of this pilot study was to quantify radiographic changes in subchondral bone after treatment with a minimally invasive joint unloading implant (KineSpring^® Knee Implant System).

Methods

Nine patients with unilateral medial knee OA resistant to nonsurgical therapy were treated with the KineSpring System and followed for 2 years. Main outcomes included Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) pain, function, and stiffness subscores and independent core laboratory determinations of joint space width and fractal signature of the tibial cortex.

Results

WOMAC scores, on average, improved by 92% for pain, 91% for function, and 79% for stiffness over the 2-year follow-up period. Joint space width in the medial compartment of the treated knee significantly increased from 0.9 mm at baseline to 3.1 mm at 2 years; joint space width in the medial compartment of the untreated knee was unchanged. Fractal signatures of the vertically oriented trabeculae in the medial compartment decreased by 2.8% in the treated knee and increased by 2.1% in the untreated knee over 2 years. No statistically significant fractal signature changes were observed in the horizontally oriented trabeculae in the medial compartment or in the horizontal or vertical trabeculae of the lateral compartment in the treated knee.

Conclusion

Preliminary evidence suggests that the KineSpring System may modify knee OA disease progression by increasing joint space width and improving subchondral bone trabecular integrity, thereby reducing pain and improving joint function.

Characterization of knee osteoarthritis-related changes in trabecular bone using texture parameters at various levels of spatial resolution-a simulation study

Articular cartilage and subchondral bone are the key tissues in osteoarthritis (OA). The role of the cancellous bone increasingly attracts attention in OA research. Because of its fast adaptation to changes in the loading distribution across joints, its quantification is expected to improve the diagnosis and monitoring of OA. In this study, we simulated OA progression-related changes of trabecular structure in a series of digital bone models and then characterized the potential of texture parameters and bone mineral density (BMD) as surrogate measures to quantify trabecular bone structure. Five texture parameters were studied: entropy, global and local inhomogeneity, anisotropy and variogram slope. Their dependence on OA relevant structural changes was investigated for three spatial resolutions typically used in micro computed tomography (CT; 10 μm), high-resolution peripheral quantitative CT (HR-pQCT) (90 μm) and clinical whole-body CT equipment (250 μm). At all resolutions, OA-related changes in trabecular bone architecture can be quantified using a specific (resolution dependent) combination of three texture parameters. BMD alone is inadequate for this purpose but if available reduces the required texture parameter combination to anisotropy and global inhomogeneity. The results are summarized in a comprehensive analysis guide for the detection of structural changes in OA knees. In conclusion, texture parameters can be used to characterize trabecular bone architecture even at spatial resolutions below the dimensions of a single trabecula and are essential for a detailed classification of relevant OA changes that cannot be achieved with a measurement of BMD alone.

Quantification of differences in bone texture from plain radiographs in knees with and without osteoarthritis

OBJECTIVE

To quantify differences in bone texture between subjects with different stages of knee osteoarthritis (OA) and age- and gender-matched controls from plain radiographs using advanced image analysis methods.

DESIGN

Altogether 203 knees were imaged using constant X-ray parameters and graded according to Kellgren-Lawrence (KL) grading scale (KL0: n = 110, KL1: n = 28, KL2: n = 27, KL3: n = 31, KL4: n = 7). Bone density-related and structure-related parameters were calculated from medial and lateral tibial subchondral bone plate and trabecular bone and from femur. Density-related parameters were derived from grayscale values and structure-related parameters from Laplacian- and local binary patterns (LBP)-based images.

RESULTS

Reproducibilities of structure-related parameters were better than bone density-related parameters. Bone density-related parameters were significantly (P < 0.05) higher in KL2-4 groups than in control group (KL0) in medial tibial subchondral bone plate and trabecular bone. LBP-based structure parameters differed significantly between KL0 and KL2-4 groups in medial subchondral bone plate, between KL0 and KL1-4 groups in medial and lateral trabecular bone, and between KL0 and KL1-4/KL2-4 in medial and lateral femur. Laplacian-based parameters differed significantly between KL0 and KL2-4 groups in medial side regions-of-interest (ROIs).

CONCLUSIONS

Our results indicate that the changes in bone texture in knee OA can be quantitatively evaluated from plain radiographs using advanced image analysis. Based on the results, increased bone density can be directly estimated if the X-ray imaging conditions are constant between patients. However, structural analysis of bone was more reproducible than direct evaluation of grayscale values, and is therefore better suited for quantitative analysis when imaging conditions are variable.

Trabecular bone texture detected by plain radiography is associated with an increased risk of knee replacement in patients with osteoarthritis: a 6 year prospective follow up study

OBJECTIVE

To examine the association between trabecular bone texture and knee joint replacement (KJR) measured using a variance orientation transform (VOT) method.

METHODS

The association of trabecular bone texture and KJR was examined prospectively over 6 years in 123 subjects with symptomatic knee osteoarthritis (OA): data regarding KJR was available for 114 (93%). At baseline, weight-bearing anteroposterior tibio-femoral radiographs were acquired. Trabecular bone texture regions were selected from the medial and lateral subchondral tibia. The VOT method was applied to each region and five fractal bone texture parameters, i.e., mean fractal dimension (FDMEAN), fractal dimensions in the horizontal (FDH) and vertical (FDV) directions, and along the roughest part of trabecular bone (FD(Sta)), and texture aspect ratio (Str) were calculated. The association between groups with increasing baseline fractal parameters (defined using tertiles) with risk of JR was examined using logistic regression.

RESULTS

28 (25%) participants' study knees underwent KJR over 6 years. Participants with KJR had lower medial FD(MEAN) and FD(H) parameters (P = 0.02 for difference). With increasing FD(MEAN), adjusted for age, gender, body mass index (BMI), osteophyte grade, joint space narrowing (JSN) grade and WOMAC pain score, the odds of KJR was reduced (P = 0.04 for trend).

CONCLUSION

This study suggests that the texture of medial tibial trabecular bone measured from plain radiographs is related to the risk of KJR: with increasing FD(MEAN) (the overall measure of bone texture roughness) the risk of KJR was reduced, independent of other clinical predictors of joint replacement. Tibial trabecular bone texture may be a useful marker of disease progression and a target of therapy in OA.

Subchondral bone trabecular integrity predicts and changes concurrently with radiographic and magnetic resonance imaging-determined knee osteoarthritis progression

OBJECTIVE

To evaluate subchondral bone trabecular integrity (BTI) on radiographs as a predictor of knee osteoarthritis (OA) progression.

METHODS

Longitudinal (baseline, 12-month, and 24-month) knee radiographs were available for 60 female subjects with knee OA. OA progression was defined by 12- and 24-month changes in radiographic medial compartment minimal joint space width (JSW) and medial joint space area (JSA), and by medial tibial and femoral cartilage volume on magnetic resonance imaging. BTI of the medial tibial plateau was analyzed by fractal signature analysis using commercially available software. Receiver operating characteristic (ROC) curves for BTI were used to predict a 5% change in OA progression parameters.

RESULTS

Individual terms (linear and quadratic) of baseline BTI of vertical trabeculae predicted knee OA progression based on 12- and 24-month changes in JSA (P < 0.01 for 24 months), 24-month change in tibial (P < 0.05), but not femoral, cartilage volume, and 24-month change in JSW (P = 0.05). ROC curves using both terms of baseline BTI predicted a 5% change in the following OA progression parameters over 24 months with high accuracy, as reflected by the area under the curve measures: JSW 81%, JSA 85%, tibial cartilage volume 75%, and femoral cartilage volume 85%. Change in BTI was also significantly associated (P < 0.05) with concurrent change in JSA over 12 and 24 months and with change in tibial cartilage volume over 24 months.

CONCLUSION

BTI predicts structural OA progression as determined by radiographic and MRI outcomes. BTI may therefore be worthy of study as an outcome measure for OA studies and clinical trials.

A dissimilarity-based multiple classifier system for trabecular bone texture in detection and prediction of progression of knee osteoarthritis

There is a growing need for classification systems that can accurately detect and predict knee osteoarthritis (OA) from plain radiographs. For this purpose, a system based on a support vector machine (SVM) classifier and distances measured between trabecular bone (TB) texture images was developed and tested in previous work. Unlike other systems, it allows an image classification without the calculation and selection of numerous texture features, and it is invariant to a range of imaging conditions encountered in a routine X-ray screening of knees. Although the system exhibited 85.4% classification accuracy in OA detection, which was higher than those obtained from other systems, its performance could be further improved. To achieve this, a dissimilarity-based multiple classifier (DMC) system is developed in this study. The system measures distances between TB texture images and generates a diverse ensemble of classifiers using prototype selection, bootstrapping of training set and heterogeneous classifiers. A measure of competence is used to select accurate (i.e. better-than-random) classifiers from the ensemble, which are then combined through the majority voting rule. To evaluate the newly developed system in OA detection (prediction of OA progression), TB texture images selected on standardised radiographs of healthy and OA (non-progressive and progressive OA) knees were used. The results obtained showed that the DMC system has higher classification accuracies for the detection (90.51% with 87.65% specificity and 93.33% sensitivity) and prediction (80% with 82.00% specificity and 77.97% sensitivity) than other systems, indicating its potential as a decision-support tool for the assessment of radiographic knee OA.

Assessment of bone mineral density and radiographic texture analysis at the tibial subchondral bone

Microstructural changes of subchondral bone constitute one of the figures characterising osteoarthritis on a structural level. Subchondral bone mineral density may reflect the complex relationship between bone and cartilage submitted to movement and loading. In this review, the authors discussed the interest of tibial subchondral bone mineral density assessment in the perspective of its diagnostic, etiopathogenic and prognostic value in osteoarthritis. In addition, the sources of variability linked to the measurement of tibial subchondral bone mineral density are precised. Trabecular bone structure characterisation by radiographic texture analyses may also represent a new promising tool to evaluate the microarchitectural changes that occur with initiation and progression of osteoarthritis. In this paper, the authors also highlighted the interest of different radiographic texture analyses and their clinical relevance in the field of osteoarthritis.

Knee malalignment is associated with an increased risk for incident and enlarging bone marrow lesions in the more loaded compartments: the MOST study

OBJECTIVE

To examine the relationship of knee malalignment with occurrence of incident and enlarging bone marrow lesions (BMLs) and regression of BMLs.

METHODS

Subjects from the Multicenter Osteoarthritis Study aged 50-79 years with or at high risk of knee osteoarthritis were studied. Full-limb radiographs were taken at baseline and hip-knee-ankle mechanical axis was measured. Baseline and 30-month magnetic resonance imaging (MRI) of knees (n = 1782) were semiquantitatively assessed for BMLs. Outcome was defined as a change in BML score in femoral/tibial condyle in medial/lateral compartments. Medial compartment in varus alignment and lateral compartment in valgus alignment were combined to form 'more loaded' compartment, while lateral compartment in valgus and medial compartment in varus were combined to form 'less loaded' compartment. Relative risk (RR) of BML score increase or decrease in relation to malalignment was estimated using a log linear regression model with the Poisson assumption, adjusting for age, gender, body mass index, physical activity scale for the elderly, race and clinic site. Further, results were stratified by ipsilateral meniscal and cartilage status at baseline.

RESULTS

Baseline varus alignment was associated with higher risk of BML score increase from baseline to follow-up in the medial compartment [adjusted RRs (95%CI): 1.5 (1.2-1.9)] and valgus alignment in the lateral compartment [1.4 (1.0-2.1)]. Increase in BML score was more likely in the more loaded compartments [1.7 (1.4-2.0)] in malaligned knees. Regardless of ipsilateral cartilage or meniscus status, adjusted RR for BML score increase was higher in the more loaded compartments of malaligned knees than those with neutral alignment. Decrease in BML score was less likely in the more loaded compartments in malaligned knees [0.8 (0.7-1.0)].

CONCLUSION

Knee malalignment is associated with increased risk of incident and enlarging BMLs in the more loaded compartments of the tibiofemoral joint.

Prediction of progression of radiographic knee osteoarthritis using tibial trabecular bone texture

OBJECTIVE

To develop a system for predicting the progression of radiographic knee osteoarthritis (OA) using tibial trabecular bone texture.

METHODS

We studied 203 knees with (n = 68) or without (n = 135) radiographic tibiofemoral OA in 105 subjects (90 men and 15 women with a mean age of 54 years) in whom 2 sets of knee radiographs were obtained 4 years apart. We determined medial and lateral compartment tibial trabecular bone texture using an automated region selection method. Three texture parameters were calculated: roughness, degree of anisotropy, and direction of anisotropy based on a signature dissimilarity measure method. We evaluated tibiofemoral OA progression using a radiographic semiquantitative outcome: an increase in the medial joint space narrowing (JSN) grade. We examined the predictive ability of trabecular bone texture in knees with and those without preexisting radiographic OA, with adjustment for age, sex, and body mass index, using logistic regression (generalized estimating equations) and receiver operating characteristic curves.

RESULTS

The prediction of increased medial JSN in knees with or without preexisting radiographic OA was the most accurate for medial trabecular bone texture; the area under the curve (AUC) was 0.77 and 0.75, respectively. For lateral trabecular bone texture, the AUC was 0.71 in knees with preexisting OA and 0.72 in knees without preexisting OA.

CONCLUSION

We have developed a system, based on analyzing tibial trabecular bone texture, which yields good prediction of loss of tibiofemoral joint space. The predictive ability of the system needs to be further validated.

Risk factors for medial meniscal pathology on knee MRI in older US adults: a multicentre prospective cohort study

OBJECTIVES

Meniscal pathology in which the aetiology is often unclear is a frequent finding on knee MRI. This study investigates potential risk factors for medial meniscal lesions or extrusion in middle-aged and elderly persons.

METHODS

Prospective cohort study using population-based subjects from Birmingham, Alabama and Iowa City, Iowa, USA (the Multicenter Osteoarthritis Study). 644 men and women aged 50-79 years with or at high risk of knee osteoarthritis (Kellgren and Lawrence grade 0-2) but with normal medial meniscal status at baseline were studied. Paired baseline and 30-month 1.0 T knee MRI were scored for meniscal lesions and extrusion (pathology) and the following systemic, knee-specific and compartment-specific potential risk factors were evaluated: age, sex, body mass index, bony enlargement of finger joints, knee trauma, leg-length inequality and knee alignment.

RESULTS

Of 791 knees, 77 (9.7%) had medial meniscal pathology at 30 months follow-up. 61 of the 77 (81%) had no report of trauma during follow-up. Including all potential risk factors in the multivariable model, the adjusted OR for medial meniscal pathology was 4.14 (95% CI 2.06 to 8.31) for knee trauma during follow-up, 1.64 (1.00 to 2.70) for five or more bony enlargements of finger joints (vs ≤ 4) and 2.00 (1.18 to 3.40) for varus alignment (vs not varus) at baseline examination. Obesity was a risk factor for the development of meniscal extrusion, OR 3.04 (1.04 to 8.93) but not for meniscal lesions, OR 1.15 (0.52 to 2.54).

CONCLUSIONS

Apart from knee trauma, possible generalised osteoarthritis, expressed as multiple bony enlargements of finger joints, varus alignment and obesity are risk factors for medial meniscal pathology.

A computer analysis method for correlating knee X-rays with continuous indicators

PURPOSE

To develop an image analysis method that can automatically find correlations between a set of plain radiographs and continuous clinical or physiological indicators.

METHODS

Knee X-rays taken from the Baltimore Longitudinal Study of Aging are used in this study. The computer analysis method is based on the WND-CHARM image feature set filtered by using the Pearson correlation of each feature with the continuous variable, and the estimated value is determined by a weighted nearest neighbor interpolation.

RESULTS

Experimental results using 300 radiographs show that the proposed method can correlate knee X-rays with physiological indicators such as sex, age, height, weight, and BMI. For instance, the Pearson correlation between the X-ray images and the height and weight were 0.59 and 0.62, respectively.

CONCLUSIONS

Using computer analysis, X-ray images can be correlated to continuous physiological variables that might not have a direct and straightforward link to the visual content of the radiograph. This approach of radiology image analysis can be used in population studies for detecting biomarkers and also in genome-wide association studies for studying the link between genes and anatomy.

Trabecular bone texture detected by plain radiography and variance orientation transform method is different between knees with and without cartilage defects

The objective of this work is to evaluate differences in trabecular bone (TB) texture between subjects with and without tibiofemoral cartilage defects using a variance orientation transform (VOT) method. A case-control study was performed in subjects without radiographic knee osteoarthritis (OA) (K&L grade <2) matched on sex, BMI, age, knee compartment, and meniscectomy where cases (n = 28) had cartilage defects (grade ≥2) and controls (n = 28) had no cartilage defects (grade <2). Cartilage defects were assessed from MRI using validated methods. The VOT was applied to TB regions selected on medial and lateral compartments in knee X-rays and fractal signatures (FS) in the horizontal (FS(H) ) and vertical (FS(V) ) directions, and along the roughest part of TB (FS(Sta) ) and texture aspect ratio signatures (StrS), at different trabecular image sizes (0.30-0.70 mm) were calculated. Compared with controls, FS(V) for cases were higher (p < 0.011) at image sizes 0.30-0.40 mm and 0.45-0.55 mm in the medial compartment. In the lateral compartment, FS(H) and FS(Sta) for cases were higher (p < 0.028) than those for controls at 0.30-0.40 mm and 0.45-0.55 mm, while FS(V) was higher (p < 0.02) at 0.30-0.40 mm. TB texture roughness was greater in subjects with cartilage defects than in subjects without, suggesting thinning and fenestration of TB occur early in OA and that the VOT identifies changes in TB in knees with early cartilage damage. No differences in StrS (p > 0.05) were found.

Dependence of trabecular structure on bone quantity: a comparison between osteoarthritic and non-pathological bone

BACKGROUND

The mechanical characterization of trabecular bone is related to its structure. In order to describe the trabecular structure and to study the mechanical behavior of the trabecular tissue, several parameters are presented in the literature. Some studies suggest a possible dependence of the structure on bone volume fraction; this dependence could bias the validity of previous studies. The problem increases its complexity when pathological bone such as osteoarthritic tissue is studied, where the organization of the trabecular structure could be different if compared to the non-pathological tissue. The primary aim of this study was to evaluate the dependence between trabecular structure and bone volume fraction. The secondary aim was to compare osteoarthritic and non-pathological bone considering the correlation between structure and bone volume fraction.

METHODS

Sixty trabecular bone specimens were extracted from femoral heads of two groups of 30 Caucasian donors; an osteoarthritic group and a non-pathological group. Several structural parameters, such as bone volume fraction, direct trabecular thickness, fabric tensor eigenvalues and their normalizations, were calculated from micro-CT analysis. A statistical analysis was carried out to identify the dependences between structural parameters and bone volume fraction. The comparison between osteoarthritic bone and non-pathological bone was also performed.

FINDINGS

Only the normalized eigenvalues of the fabric tensor were not correlated to bone volume fraction (R<0.5). The first and second normalized eigenvalues were significantly different between osteoarthritic bone and non-pathological bone (respectively P<0.05 and P<0.001).

INTERPRETATION

In conclusion, orientation and anisotropy of the trabecular structure do not depend on bone volume fraction. Moreover, differences in the first and second normalized fabric tensor eigenvalues suggest in the osteoarthritic group a structure more oriented along the main trabecular direction.

Meniscal pathology on MRI increases the risk for both incident and enlarging subchondral bone marrow lesions of the knee: the MOST Study

OBJECTIVES

To investigate the association between meniscal pathology and incident or enlarging bone marrow lesions (BML) in knee osteoarthritis.

METHODS

The authors studied subjects from the Multicenter Osteoarthritis Study aged 50-79 years either with knee osteoarthritis or at high risk of the disease. Baseline and 30-months magnetic resonance images of knees (n=1344) were scored for subchondral BML. Outcome was defined as an increase in BML score in either the tibial or femoral condyle in medial and lateral compartments, respectively. The authors defined meniscal pathology at baseline as the presence of either meniscal lesions or meniscal extrusion. The risk of an increase in BML score in relation to meniscal status in the same compartment was estimated using a log linear regression model adjusted for age, sex, body mass index, physical activity level and mechanical axis. In secondary analyses the investigators stratified by ipsilateral tibiofemoral cartilage status at baseline and compartments with pre-existing BML.

RESULTS

The adjusted relative risk of incident or enlarging BML ranged from 1.8; 95% CI 1.3 to 2.3 for mild medial meniscal pathology to 5.0; 95% CI 3.2 to 7.7 for major lateral meniscal pathology (using no meniscal pathology in the same compartment as reference). Stratification by cartilage or BML status at baseline had essentially no effect on these estimates.

CONCLUSIONS

Knee compartments with meniscal pathology have a substantially increased risk of incident or enlarging subchondral BML over 30 months. Higher relative risks were seen in those with more severe and with lateral meniscal pathology.

Progression analysis and stage discovery in continuous physiological processes using image computing

We propose an image computing-based method for quantitative analysis of continuous physiological processes that can be sensed by medical imaging and demonstrate its application to the analysis of morphological alterations of the bone structure, which correlate with the progression of osteoarthritis (OA). The purpose of the analysis is to quantitatively estimate OA progression in a fashion that can assist in understanding the pathophysiology of the disease. Ultimately, the texture analysis will be able to provide an alternative OA scoring method, which can potentially reflect the progression of the disease in a more direct fashion compared to the existing clinically utilized classification schemes based on radiology. This method can be useful not just for studying the nature of OA, but also for developing and testing the effect of drugs and treatments. While in this paper we demonstrate the application of the method to osteoarthritis, its generality makes it suitable for the analysis of other progressive clinical conditions that can be diagnosed and prognosed by using medical imaging.

A signature dissimilarity measure for trabecular bone texture in knee radiographs

PURPOSE

The purpose of this study is to develop a dissimilarity measure for the classification of trabecular bone (TB) texture in knee radiographs. Problems associated with the traditional extraction and selection of texture features and with the invariance to imaging conditions such as image size, anisotropy, noise, blur, exposure, magnification, and projection angle were addressed.

METHODS

In the method developed, called a signature dissimilarity measure (SDM), a sum of earth mover's distances calculated for roughness and orientation signatures is used to quantify dissimilarities between textures. Scale-space theory was used to ensure scale and rotation invariance. The effects of image size, anisotropy, noise, and blur on the SDM developed were studied using computer generated fractal texture images. The invariance of the measure to image exposure, magnification, and projection angle was studied using x-ray images of human tibia head. For the studies, Mann-Whitney tests with significance level of 0.01 were used. A comparison study between the performances of a SDM based classification system and other two systems in the classification of Brodatz textures and the detection of knee osteoarthritis (OA) were conducted. The other systems are based on weighted neighbor distance using compound hierarchy of algorithms representing morphology (WND-CHARM) and local binary patterns (LBP).

RESULTS

Results obtained indicate that the SDM developed is invariant to image exposure (2.5-30 mA s), magnification (x1.00 - x1.35), noise associated with film graininess and quantum mottle (< 25%), blur generated by a sharp film screen, and image size (> 64 x 64 pixels). However, the measure is sensitive to changes in projection angle (> 5 degrees), image anisotropy (> 30 degrees), and blur generated by a regular film screen. For the classification of Brodatz textures, the SDM based system produced comparable results to the LBP system. For the detection of knee OA, the SDM based system achieved 78.8% classification accuracy and outperformed the WND-CHARM system (64.2%).

CONCLUSIONS

The SDM is well suited for the classification of TB texture images in knee OA detection and may be useful for the texture classification of medical images in general.

Differences in trabecular bone texture between knees with and without radiographic osteoarthritis detected by directional fractal signature method

OBJECTIVE

To evaluate differences in tibial trabecular bone (TB) texture between subjects with and without radiographic knee osteoarthritis (OA) using a variance orientation transform (VOT) method.

DESIGN

Subjects with knee OA (Kellgren & Lawrence grade > or =2) and controls without OA (both n=26, seven women) were matched by sex, age, body mass index and compartment. The VOT method was applied to TB X-ray images and fractal signature and dimension in horizontal (FS(H), FD(H)) and vertical (FS(V), FD(V)) directions and along the roughest part of TB (FS(Sta), FD(Sta)), texture aspect ratio (Str) and signature (StrS), and mean FD (FD(MEAN)) were calculated. The VOT method was compared against an augmented Hurst orientation transform (HOT) method using paired t tests, intraclass correlation coefficients (ICCs) and coefficients of variation (CVs%). Longitudinal sensitivity to OA bone changes was not assessed.

RESULTS

For the reproducibility of texture parameters, ICCs were >0.75 and CVs% were <8.2% for both methods. Compared with controls, FD(MEAN), FD(H), FD(V) and FD(Sta) for OA knees were lower (P<0.001), while Str was higher in both medial (P=0.03) and lateral (P=0.02) compartments. FS(H), FS(Sta) were lower for OA knees than for controls at sizes 0.3-0.7 mm (P<0.001) in both compartments. In lateral compartment, FS(V) for OA knees was lower than for controls at sizes 0.3-0.5 mm (P<0.001) and 0.55-0.70 mm (P<0.02), while in medial compartment at sizes 0.3-0.7 mm (P<0.001). Compared with controls, StrS for OA knees was higher at sizes 0.3, 0.55-0.70 mm in medial (P<0.03) and lateral (P<0.04) compartments.

CONCLUSIONS

The VOT method is comparable to HOT method in the reproducibility of texture parameters and the ability to discriminate between non-OA and OA TB textures. However, unlike the HOT method, it quantifies texture roughness along the roughest part of the tibial bone, texture anisotropy at individual trabecular sizes and it works over a larger range of trabecular sizes. The VOT method may be a valuable tool for studying OA changes in TB.

Trabecular morphometry by fractal signature analysis is a novel marker of osteoarthritis progression

OBJECTIVE

To evaluate the effectiveness of using subchondral bone texture observed on a radiograph taken at baseline to predict progression of knee osteoarthritis (OA) over a 3-year period.

METHODS

A total of 138 participants in the Prediction of Osteoarthritis Progression study were evaluated at baseline and after 3 years. Fractal signature analysis (FSA) of the medial subchondral tibial plateau was performed on fixed flexion radiographs of 248 nonreplaced knees, using a commercially available software tool. OA progression was defined as a change in joint space narrowing (JSN) or osteophyte formation of 1 grade according to a standardized knee atlas. Statistical analysis of fractal signatures was performed using a new model based on correlating the overall shape of a fractal dimension curve with radius.

RESULTS

Fractal signature of the medial tibial plateau at baseline was predictive of medial knee JSN progression (area under the curve [AUC] 0.75, of a receiver operating characteristic curve) but was not predictive of osteophyte formation or progression of JSN in the lateral compartment. Traditional covariates (age, sex, body mass index, knee pain), general bone mineral content, and joint space width at baseline were no more effective than random variables for predicting OA progression (AUC 0.52-0.58). The predictive model with maximum effectiveness combined fractal signature at baseline, knee alignment, traditional covariates, and bone mineral content (AUC 0.79).

CONCLUSION

We identified a prognostic marker of OA that is readily extracted from a plain radiograph using FSA. Although the method needs to be validated in a second cohort, our results indicate that the global shape approach to analyzing these data is a potentially efficient means of identifying individuals at risk of knee OA progression.

Early detection of radiographic knee osteoarthritis using computer-aided analysis

OBJECTIVE

To determine whether computer-based analysis can detect features predictive of osteoarthritis (OA) development in radiographically normal knees.

METHOD

A systematic computer-aided image analysis method weighted neighbor distances using a compound hierarchy of algorithms representing morphology (WND-CHARM) was used to analyze pairs of weight-bearing knee X-rays. Initial X-rays were all scored as normal Kellgren-Lawrence (KL) grade 0, and on follow-up approximately 20 years later either developed OA (defined as KL grade=2) or remained normal.

RESULTS

The computer-aided method predicted whether a knee would change from KL grade 0 to grade 3 with 72% accuracy (P<0.00001), and to grade 2 with 62% accuracy (P<0.01). Although a large part of the predictive signal comes from the image tiles that contained the joint, the region adjacent to the tibial spines provided the strongest predictive signal.

CONCLUSION

Radiographic features detectable using a computer-aided image analysis method can predict the future development of radiographic knee OA.

Directional fractal signature analysis of trabecular bone: evaluation of different methods to detect early osteoarthritis in knee radiographs

There is ongoing research directed towards the development of cheap and reliable decision support systems for the detection and prediction of osteoarthritis (OA) in knee joints. Fractal analysis of trabecular bone texture X-ray images is one of the most promising approaches. It is cheap and non-invasive. However, difficulties arise when the fractal signature methods are used to quantify bone roughness and anisotropy on individual scales. This is because the fractal methods are able to quantify bone texture only in the vertical and horizontal directions, and previous studies showed that OA bone changes can occur in any direction. To address these difficulties, three directional fractal signature methods were developed in this study, i.e. a fractal signature Hurst orientation transform (FSHOT) method, a variance orientation transform (VOT) method, and a blanket with rotating-grid (BRG) method. These methods were tested and the best performing method was selected. Unlike other methods, the newly developed techniques are able to calculate fractal dimensions (FDs) on individual scales (i.e. fractal signature) in all possible directions. The accuracy of the methods developed in measuring texture roughness and anisotropy on individual scales was evaluated. The effects of imaging conditions such as image noise, blur, exposure, magnification, and projection angle and the effects of translation of the bone region of interest on texture parameters were also evaluated. Computer-generated fractal surface images with known FDs and X-ray images obtained for a human tibia head were used. Results obtained show that the VOT method performs better than the FSHOT and BRG methods.

Biometric identification using knee X-rays

Identification of people often makes use of unique features of the face, fingerprints and retina. Beyond this, a similar identifying process can be applied to internal parts of the body that are not visible to the unaided eye. Here we show that knee X-rays can be used for the identification of individual persons. The image analysis method is based on the wnd-charm algorithm, which has been found effective for the diagnosis of clinical conditions of knee joints. Experimental results show that the rank-10 identification accuracy using a dataset of 425 individuals is ~56%, and the rank-1 accuracy is ~34%. The dataset contained knee X-rays taken several years apart from each other, showing that the identifiable features correspond to specific persons, rather than the present clinical condition of the joint.

Automated selection of trabecular bone regions in knee radiographs

Osteoarthritic (OA) changes in knee joints can be assessed by analyzing the structure of trabecular bone (TB) in the tibia. This analysis is performed on TB regions selected manually by a human operator on x-ray images. Manual selection is time-consuming, tedious, and expensive. Even if a radiologist expert or highly trained person is available to select regions, high inter- and intraobserver variabilities are still possible. A fully automated image segmentation method was, therefore, developed to select the bone regions for numerical analyses of changes in bone structures. The newly developed method consists of image preprocessing, delineation of cortical bone plates (active shape model), and location of regions of interest (ROI). The method was trained on an independent set of 40 x-ray images. Automatically selected regions were compared to the "gold standard" that contains ROIs selected manually by a radiologist expert on 132 x-ray images. All images were acquired from subjects locked in a standardized standing position using a radiography rig. The size of each ROI is 12.8 x 12.8 mm. The automated method results showed a good agreement with the gold standard [similarity index (SI) = 0.83 (medial) and 0.81 (lateral) and the offset =[-1.78, 1.27]x[-0.65,0.26] mm (medial) and [-2.15, 1.59]x[-0.58, 0.52] mm (lateral)]. Bland and Altman plots were constructed for fractal signatures, and changes of fractal dimensions (FD) to region offsets calculated between the gold standard and automatically selected regions were calculated. The plots showed a random scatter and the 95% confidence intervals were (-0.006, 0.008) and (-0.001, 0.011). The changes of FDs to region offsets were less than 0.035. Previous studies showed that differences in FDs between non-OA and OA bone regions were greater than 0.05. ROIs were also selected by a second radiologist and then evaluated. Results indicated that the newly developed method could replace a human operator and produces bone regions with an accuracy that is sufficient for fractal analyses of bone texture.