Computer-based methods for measuring joint space and estimating erosion volume in the finger and wrist joints of patients with rheumatoid arthritis

Application of ImageJ software for quantification of Hand Joint Space Narrowing in Patients with Rheumatoid Arthritis

BACKGROUND

ImageJ software is used to quantify the joint space width (JSW) of hand and wrist in patients with rheumatoid arthritis (RA) as well as in the healthy control group.

METHOD

Forty-one RA patients and 31 healthy controls are included in this study. All of 72 participants underwent digital radiography of the bilateral hand and wrist, then all the images were opened by ImageJ software to measure the width of wrist and hand joint space (total 2160 joints). Joint space narrowing (JSN) was defined if the width was less than the mean - 2SD of the control group.

RESULT

The mean JSW of all sites of wrist and hand joints of RA patients was significantly reduced as compared to those in the control group (p<0.001). There were 37/41 (90.24%) RA patients who had JSN in at least one joint in hand or wrist. In total, 70.89% of joints on the right and 68.46% of joints on the left wrist and hand had JSN.

CONCLUSION

ImageJ software was simple and convenient , which helps rheumatologists quantify the width of joint space for diagnosis and follow-up in RA patients.

Computer-aided diagnosis based on hand thermal, RGB images, and grip force using artificial intelligence as screening tool for rheumatoid arthritis in women

Rheumatoid arthritis (RA) is an autoimmune disorder that typically affects people between 23 and 60 years old causing chronic synovial inflammation, symmetrical polyarthritis, destruction of large and small joints, and chronic disability. Clinical diagnosis of RA is stablished by current ACR-EULAR criteria, and it is crucial for starting conventional therapy in order to minimize damage progression. The 2010 ACR-EULAR criteria include the presence of swollen joints, elevated levels of rheumatoid factor or anti-citrullinated protein antibodies (ACPA), elevated acute phase reactant, and duration of symptoms. In this paper, a computer-aided system for helping in the RA diagnosis, based on quantitative and easy-to-acquire variables, is presented. The participants in this study were all female, grouped into two classes: class I, patients diagnosed with RA (n = 100), and class II corresponding to controls without RA (n = 100). The novel approach is constituted by the acquisition of thermal and RGB images, recording their hand grip strength or gripping force. The weight, height, and age were also obtained from all participants. The color layout descriptors (CLD) were obtained from each image for having a compact representation. After, a wrapper forward selection method in a range of classification algorithms included in WEKA was performed. In the feature selection process, variables such as hand images, grip force, and age were found relevant, whereas weight and height did not provide important information to the classification. Our system obtains an AUC ROC curve greater than 0.94 for both thermal and RGB images using the RandomForest classifier. Thirty-eight subjects were considered for an external test in order to evaluate and validate the model implementation. In this test, an accuracy of 94.7% was obtained using RGB images; the confusion matrix revealed our system provides a correct diagnosis for all participants and failed in only two of them (5.3%). Graphical abstract.

Detection of Fine Radiographic Progression in Finger Joint Space Narrowing Beyond Human Eyes: Phantom Experiment and Clinical Study with Rheumatoid Arthritis Patients

The visual assessment of joint space narrowing (JSN) on radiographs of rheumatoid arthritis (RA) patients such as the Genant-modified Sharp score (GSS) is widely accepted but limited by its subjectivity and insufficient sensitivity. We developed a software application which can assess JSN quantitatively using a temporal subtraction technique for radiographs, in which the chronological change in JSN between two radiographs was defined as the joint space difference index (JSDI). The aim of this study is to prove the superiority of the software in terms of detecting fine radiographic progression in finger JSN over human observers. A micrometer measurement apparatus that can adjust arbitrary joint space width (JSW) in a phantom joint was developed to define true JSW. We compared the smallest detectable changes in JSW between the JSDI and visual assessment using phantom images. In a clinical study, 222 finger joints without interval score change on GSS in 15 RA patients were examined. We compared the JSDI between joints with and without synovial vascularity (SV) on power Doppler ultrasonography during the follow-up period. True JSW difference was correlated with JSDI for JSW differences ranging from 0.10 to 1.00 mm at increments of 0.10 mm (R² = 0.986 and P < 0.001). Rheumatologists were difficult to detect JSW difference of 0.30 mm or less. The JSDI of finger joints with SV was significantly higher than those without SV (P = 0.030). The software can detect fine differences in JSW that are visually unrecognizable.

The shape language in application to the diagnosis of cervical vertebrae pathology

In this paper the possibility of classification of X-ray images of the cervical vertebrae is studied. The images should be classified into one of the following classes-the images of healthy vertebrae and the images of vertebrae with syndesmophytes. The vertebra contours, described unambiguously by using the generalized shape language, are the basis of the analysis. As a result, the contour is represented as a chain of sinquads that determine switches. The found switches are the characteristic points of the analyzed contour. In these points additional features of the contour are determined. On the basis of these features two aforementioned classes of images are defined as fuzzy sets. Such an approach allows us to create a hierarchical algorithm of classification based on the syntactic and fuzzy description of the contour.

Computer-Based Radiographic Quantification of Joint Space Narrowing Progression Using Sequential Hand Radiographs: Validation Study in Rheumatoid Arthritis Patients from Multiple Institutions

We have developed a refined computer-based method to detect joint space narrowing (JSN) progression with the joint space narrowing progression index (JSNPI) by superimposing sequential hand radiographs. The purpose of this study is to assess the validity of a computer-based method using images obtained from multiple institutions in rheumatoid arthritis (RA) patients. Sequential hand radiographs of 42 patients (37 females and 5 males) with RA from two institutions were analyzed by a computer-based method and visual scoring systems as a standard of reference. The JSNPI above the smallest detectable difference (SDD) defined JSN progression on the joint level. The sensitivity and specificity of the computer-based method for JSN progression was calculated using the SDD and a receiver operating characteristic (ROC) curve. Out of 314 metacarpophalangeal joints, 34 joints progressed based on the SDD, while 11 joints widened. Twenty-one joints progressed in the computer-based method, 11 joints in the scoring systems, and 13 joints in both methods. Based on the SDD, we found lower sensitivity and higher specificity with 54.2 and 92.8%, respectively. At the most discriminant cutoff point according to the ROC curve, the sensitivity and specificity was 70.8 and 81.7%, respectively. The proposed computer-based method provides quantitative measurement of JSN progression using sequential hand radiographs and may be a useful tool in follow-up assessment of joint damage in RA patients.

Automatic Quantification of Radiographic Wrist Joint Space Width of Patients With Rheumatoid Arthritis

Objective: Wrist joint space narrowing is a main radiographic outcome of rheumatoid arthritis (RA). Yet, automatic radiographic wrist joint space width (JSW) quantification for RA patients has not been widely investigated. The aim of this paper is to present an automatic method to quantify the JSW of three wrist joints that are least affected by bone overlapping and are frequently involved in RA. These joints are located around the scaphoid bone, viz. the multangular-navicular, capitate-navicular-lunate, and radiocarpal joints. Methods: The joint space around the scaphoid bone is detected by using consecutive searches of separate path segments, where each segment location aids in constraining the subsequent one. For joint margin delineation, first the boundary not affected by X-ray projection is extracted, followed by a backtrace process to obtain the actual joint margin. The accuracy of the quantified JSW is evaluated by comparison with the manually obtained ground truth. Results: Two of the 50 radiographs used for evaluation of the method did not yield a correct path through all three wrist joints. The delineated joint margins of the remaining 48 radiographs were used for JSW quantification. It was found that 90% of the joints had a JSW deviating less than 20% from the mean JSW of manual indications, with the mean JSW error less than 10%. Conclusion: The proposed method is able to automatically quantify the JSW of radiographic wrist joints reliably.

SIGNIFICANCE

The proposed method may aid clinical researchers to study the progression of wrist joint damage in RA studies.Objective: Wrist joint space narrowing is a main radiographic outcome of rheumatoid arthritis (RA). Yet, automatic radiographic wrist joint space width (JSW) quantification for RA patients has not been widely investigated. The aim of this paper is to present an automatic method to quantify the JSW of three wrist joints that are least affected by bone overlapping and are frequently involved in RA. These joints are located around the scaphoid bone, viz. the multangular-navicular, capitate-navicular-lunate, and radiocarpal joints. Methods: The joint space around the scaphoid bone is detected by using consecutive searches of separate path segments, where each segment location aids in constraining the subsequent one. For joint margin delineation, first the boundary not affected by X-ray projection is extracted, followed by a backtrace process to obtain the actual joint margin. The accuracy of the quantified JSW is evaluated by comparison with the manually obtained ground truth. Results: Two of the 50 radiographs used for evaluation of the method did not yield a correct path through all three wrist joints. The delineated joint margins of the remaining 48 radiographs were used for JSW quantification. It was found that 90% of the joints had a JSW deviating less than 20% from the mean JSW of manual indications, with the mean JSW error less than 10%. Conclusion: The proposed method is able to automatically quantify the JSW of radiographic wrist joints reliably.

SIGNIFICANCE

The proposed method may aid clinical researchers to study the progression of wrist joint damage in RA studies.

Fully automated joint space width measurement and digital X-ray radiogrammetry in early RA

Objectives

To study fully automated digital joint space width (JSW) and bone mineral density (BMD) in relation to a conventional radiographic scoring method in early rheumatoid arthritis (eRA).

Methods

Radiographs scored by the modified Sharp van der Heijde score (SHS) in patients with eRA were acquired from the SWEdish FarmacOTherapy study. Fully automated JSW measurements of bilateral metacarpals 2, 3 and 4 were compared with the joint space narrowing (JSN) score in SHS. Multilevel mixed model statistics were applied to calculate the significance of the association between ΔJSW and ΔBMD over 1 year, and the JSW differences between damaged and undamaged joints as evaluated by the JSN.

Results

Based on 576 joints of 96 patients with eRA, a significant reduction from baseline to 1 year was observed in the JSW from 1.69 (±0.19) mm to 1.66 (±0.19) mm (p<0.01), and BMD from 0.583 (±0.068) g/cm² to 0.566 (±0.074) g/cm² (p<0.01). A significant positive association was observed between ΔJSW and ΔBMD over 1 year (p<0.0001). On an individual joint level, JSWs of undamaged (JSN=0) joints were wider than damaged (JSN>0) joints: 1.68 mm (95% CI 1.70 to 1.67) vs 1.54 mm (95% CI 1.63 to 1.46). Similarly the unadjusted multilevel model showed significant differences in JSW between undamaged (1.68 mm (95% CI 1.72 to 1.64)) and damaged joints (1.63 mm (95% CI 1.68 to 1.58)) (p=0.0048). This difference remained significant in the adjusted model: 1.66 mm (95% CI 1.70 to 1.61) vs 1.62 mm (95% CI 1.68 to 1.56) (p=0.042).

Conclusions

To measure the JSW with this fully automated digital tool may be useful as a quick and observer-independent application for evaluating cartilage damage in eRA.

Trial registration number

NCT00764725.

Performance of computer-based analysis using temporal subtraction to assess joint space narrowing progression in rheumatoid patients

Our computer-based method can detect the chronological change in joint space width between baseline and follow-up images as the joint space difference index (JSDI). The aim of this study was to verify the sensitivity and specificity of our computer-based method in assessment of joint space narrowing progression in rheumatoid patients. Twenty-seven patients (24 women and 3 men) with rheumatoid arthritis underwent radiography of the bilateral hand at baseline and at 1 year. The joint space narrowing (JSN) of a total of 252 metacarpophalangeal (MCP) joints and 229 carpal joints was assessed by our computer-based method, setting the Sharp/van der Heijde method as the gold standard. We constructed a receiver operating characteristic curve by using the Sharp/van der Heijde method as the gold standard and set the optimal cutoff on JSDI for MCP, carpal, and MCP/carpal joints. We then calculated the sensitivity and specificity for each cutoff in assessment of JSN progression. At the most discriminant cutoff, the sensitivity and specificity of the computer-based method for MCP joints was 78.6 versus 85.3 %, respectively (AUC = 0.837; P < 0.001). Carpal joints revealed a lower sensitivity and specificity with 64.7 and 86.8 % (AUC = 0.775; P < 0.001). Furthermore, the sensitivity and specificity for MCP/carpal joints was 71.0 versus 83.6 %, respectively (AUC = 0.778; P < 0.001). The computer-based method presented a reliable assessment of JSN progression with high sensitivity and specificity and may be useful in follow-up assessment of the joint damage in rheumatoid patients.

OARSI Clinical Trials Recommendations: Hand imaging in clinical trials in osteoarthritis

Tremendous advances have occurred in our understanding of the pathogenesis of hand osteoarthritis (OA) and these are beginning to be applied to trials targeted at modification of the disease course. The purpose of this expert opinion, consensus driven exercise is to provide detail on how one might use and apply hand imaging assessments in disease modifying clinical trials. It includes information on acquisition methods/techniques (including guidance on positioning for radiography, sequence/protocol recommendations/hardware for MRI); commonly encountered problems (including positioning, hardware and coil failures, sequences artifacts); quality assurance/control procedures; measurement methods; measurement performance (reliability, responsiveness, validity); recommendations for trials; and research recommendations.

Outcomes and costs of incorporating a multibiomarker disease activity test in the management of patients with rheumatoid arthritis

OBJECTIVE

The multibiomarker disease activity (MBDA) blood test has been clinically validated as a measure of disease activity in patients with RA. We aimed to estimate the effect of the MBDA test on physical function for patients with RA (based on HAQ), quality-adjusted life years and costs over 10 years.

METHODS

A decision analysis was conducted to quantify the effect of using the MBDA test on RA-related outcomes and costs to private payers and employers. Results of a clinical management study reporting changes to anti-rheumatic drug recommendations after use of the MBDA test informed clinical utility. The effect of treatment changes on HAQ was derived from 5 tight-control and 13 treatment-switch trials. Baseline HAQ scores and the HAQ score relationship with medical costs and quality of life were derived from published National Data Bank for Rheumatic Diseases data.

RESULTS

Use of the MBDA test is projected to improve HAQ scores by 0.09 units in year 1, declining to 0.02 units after 10 years. Over the 10 year time horizon, quality-adjusted life years increased by 0.08 years and costs decreased by US$457 (cost savings in disability-related medical costs, US$659; in productivity costs, US$2137). The most influential variable in the analysis was the effect of the MBDA test on clinician treatment recommendations and subsequent HAQ changes.

CONCLUSION

The MBDA test aids in the assessment of disease activity in patients with RA by changing treatment decisions, improving the functional status of patients and cost savings. Further validation is ongoing and future longitudinal studies are warranted.

Validity of a computer-assisted manual segmentation software to quantify wrist erosion volume using computed tomography scans in rheumatoid arthritis

BACKGROUND

To investigate the performance of conventional radiography (CR) for the detection of bone erosions of wrist in rheumatoid arthritis (RA) using multidetector computed tomography (CT) as the reference method and to evaluate the validity of a computer-assisted manual segmentation (outlining) technique to quantify erosion volume on CT scans.

METHODS

Twenty five RA patients and six controls underwent CT and radiographic evaluation of the dominant wrist on the same day. CT was performed by using a 64 GE light Speed VCT power. Wrists images were evaluated separately and scored for the presence of erosions according to the Outcome Measures in Rheumatology Rheumatoid Arthritis MRI Scoring System (RAMRIS) and the Sharp/van der Heijde scoring method. Measurements of bone erosion volumes were obtained using OsiriX medical imaging software. The mean value of the volumes of the CT bone erosions detected at two readings was used to calculate inter-rater agreement.

RESULTS

The overall sensitivity, specificity and accuracy of radiography for detecting erosions were 25.5%, 98.3% and 70.1%, respectively. Using computer-assisted manual segmentation (outlining) technique, erosion volume on CT measurements per subject was ranged from 0.001 cm³ to 2.01 cm³. Spearman's RAMRIS score of each wrist bones in all subjects (n = 25) were correlated with the total erosion volume on CT (p < 0.0001), with the ratio between erosion volume and the corresponding bone volume on a percentage basis (p < 0.0001). The total Sharp/van der Heijde erosion score of the all wrist bones was correlate with the RAMRIS score (p = 0.008). The intraclass correlation coefficients (ICC) for manual segmentation showed high agreement (ICC = 0.901).

CONCLUSIONS

Considering CT as the reference method, CR showed very low sensitivity. A close correlation with CT erosion volumes supports the OMERACT RAMRIS erosion score as a semiquantitative measure of joint damage in RA. Although the computer-assisted manual segmentation can be beneficial for diagnostic decision in cross-sectional CT examinations of the wrist in RA, this technique will require further evaluation in terms of responsiveness.

The usefulness of computer-aided joint space analysis in the assessment of rheumatoid arthritis

OBJECTIVE

Computer-aided joint space analysis (CAJSA) is a newly developed technique for the measurement of radiogeometrically detectable joint space widths of the metacarpal-phalangeal (JSD-MCP) and proximal-interphalangeal articulations (JSD-PIP). The aim of this study was to verify the sensitivity and specificity of these CAJSA measurements in the assessment of established RA.

METHODS

Four hundred and fifty-eight participants (248 healthy subjects, 210 RA patients) underwent computerized semi-automated measurements of the JSD-MCP and JSD-PIP articulations (CAJSA, Radiogrammetry Kit, Version 1.3.6) based on digitally performed radiographs. The Sharp joint space narrowing score was also performed to determine RA-related joint space narrowing.

RESULTS

The significant severity-dependent reduction for JSD-MCP was -44.0% and for JSD-PIP, -25.94% between Sharp scores 0 and 3. The sensitivity and specificity of JSD-MCP (total) was 88.1% versus 77.8%, respectively (AUC = 0.920; P < 0.001). Furthermore, JSD-PIP (total) revealed a lower sensitivity and specificity with 61.4% and 88.7% (AUC = 0.878; P < 0.001).

CONCLUSION

The CAJSA method presented a reliable assessment of disease-related joint space narrowing in patients suffering from RA with excellent sensitivity and specificity. By providing quantitative data, other scoring methods could be significantly improved, and thereby the accuracy of the diagnosis and a better therapeutic evaluation could be achieved.

Reliability and longitudinal validity of computer-assisted methods for measuring joint damage progression in subjects with rheumatoid arthritis

OBJECTIVE

To compare the metric properties of a computer-assisted erosion segmentation volume measurement with scoring using the Rheumatoid Arthritis Magnetic Resonance Imaging Score (RAMRIS) in a longitudinal cohort of patients with rheumatoid arthritis (RA).

METHODS

Thirty-two sets of baseline and 2-year followup magnetic resonance imaging (MRI) of metacarpal phalangeal 2-5 joints of patients with RA were scored using RAMRIS and segmented using OSIRIS software. The smallest detectable difference (SDD), standardized response mean (SRM), and paired t-test were used to evaluate the sensitivity to change. Eleven of the 32 patients' MRI were segmented by both readers to evaluate interreader agreement. The 28-joint Disease Activity Score (DAS28) and Sharp erosion scores further evaluated construct and longitudinal validity.

RESULTS

Reliability of erosion progression by computer-assisted volume measurement was superior to RAMRIS [intrareader interclass correlation coefficient (ICC) 0.97 (0.94-0.99) vs 0.52 (0.22-0.73)] and interreader ICC of volume measurement was 0.85 (0.53-0.96). Computer-assisted volume measurements identified 10 of 32 patients who progressed more than the SDD progression, whereas RAMRIS identified only 4 of 32 patients (p = 0.0013). By a paired t-test, however, all MRI measures progressed significantly over 2 years (irrespective of treatment arm) and there was little difference by SRM. Construct correlational validity of the MRI methods was 0.47-0.90 for status scores and 0.33-0.81 for progression. There was no relationship between the average DAS28 and erosion progression by any imaging method.

CONCLUSION

Computer-assisted measurement of erosion volume has good performance metrics. It had excellent intrareader and interreader reliability and was more sensitive to change than RAMRIS in this group of patients. www.ClinicalTrials.gov, NCT00451971.

Does digital X-ray radiogrammetry have a role in identifying patients at increased risk for joint destruction in early rheumatoid arthritis?

INTRODUCTION

The aim of this study was to investigate the role of hand bone mineral density (BMD) loss analyzed with digital X-ray radiogrammetry (DXR) in early rheumatoid arthritis (RA) as a predictor for progression of joint damage.

METHODS

In 379 patients with early RA, baseline and one-year hand BMD was measured with DXR and the hand bone loss (HBL) was analyzed using the smallest detectable change (HBLsdc) and tertiles (HBLtertiles). Joint damage in hands and feet were scored according to the Sharp van der Heijde (SHS) method at baseline and at one, two, five and eight years. At the same time-points Disease Activity Score (DAS28) was calculated and functional disability assessed. Rheumatoid factor (RF) and antibodies against cyclic citrullinated peptides (anti-CCP) were analyzed at baseline.

RESULTS

Sixty-six percent of the patients had hand BMD loss in the first year of RA determined by HBLsdc and 65% by HBLtertiles. Radiographic progression after two, five and eight years was associated with hand bone loss defined by HBLsdc. By HBLtertiles there were significant associations at all time-points except at eight years. The change in DXR at one year (ChDXR1yr) correlated significantly and inversely with the change in SHS (ChSHS) at two, five and eight years. Multivariate analysis showed that only change in SHS during the first year and the presence of anti-CCP were independent predictors of long-term progressive joint damage. If radiographic scores were not included, DXR-BMD loss was an independent predictor. Patients with great bone loss by HBLtertiles had significantly more often high disease activity after two years. However, neither bone loss by HBLsdc or HBLtertiles nor by ChDXR1yr was an independent predictor of remission after two, five and eight years.

CONCLUSIONS

This study confirms previous reports of an association of decrease in DXR-BMD during the first disease year with progression of radiographic joint damage over an extended period of time. This association was independent in a regression model only when radiological findings were excluded suggesting a possible predictive role of DXR-BMD in clinical practice when radiographic evaluation is not available. However, further studies are required before this can be established.

ON THE SEGMENTATION AND CLASSIFICATION OF HAND RADIOGRAPHS

This research is part of a wider project to build predictive models of bone age using hand radiograph images. We examine ways of finding the outline of a hand from an X-ray as the first stage in segmenting the image into constituent bones. We assess a variety of algorithms including contouring, which has not previously been used in this context. We introduce a novel ensemble algorithm for combining outlines using two voting schemes, a likelihood ratio test and dynamic time warping (DTW). Our goal is to minimize the human intervention required, hence we investigate alternative ways of training a classifier to determine whether an outline is in fact correct or not. We evaluate outlining and classification on a set of 1370 images. We conclude that ensembling with DTW improves performance of all outlining algorithms, that the contouring algorithm used with the DTW ensemble performs the best of those assessed, and that the most effective classifier of hand outlines assessed is a random forest applied to outlines transformed into principal components.

Automatic radiographic quantification of hand osteoarthritis; accuracy and sensitivity to change in joint space width in a phantom and cadaver study

OBJECTIVE

To validate a newly developed quantification method that automatically detects and quantifies the joint space width (JSW) in hand radiographs. Repeatability, accuracy and sensitivity to changes in JSW were determined. The influence of joint location and joint shape on the measurements was tested.

METHODS

A mechanical micrometer set-up was developed to define and adjust the true JSW in an acrylic phantom joint and in human cadaver-derived phalangeal joints. Radiographic measurements of the JSW were compared to the true JSW. Repeatability, systematic error (accuracy) and sensitivity (defined as the smallest detectable difference (SDD)) were determined. The influence of joint position on the JSW measurement was assessed by varying the location of the acrylic phantom on the X-ray detector with respect to the X-ray beam and the influence of joint shape was determined by using morphologically different human cadaver joints.

RESULTS

The mean systematic error was 0.052 mm in the phantom joint and 0.210 mm in the cadaver experiment. In the phantom experiments, the repeatability was high (SDD = 0.028 mm), but differed slightly between joint locations (p = 0.046), and a change in JSW of 0.037 mm could be detected. Dependent of the joint shape in the cadaver hand, a change in JSW between 0.018 and 0.047 mm could be detected.

CONCLUSIONS

The automatic quantification method is sensitive to small changes in JSW. Considering the published data of JSW decline in the normal and osteoarthritic population, the first signs of OA progression with this method can be detected within 1 or 2 years.

Reproducibility and influence of hand rotation on computer-aided joint space analysis

OBJECTIVES

Computer-aided joint space analysis (CAJSA) is a recently developed, semi-automated tool to aid in the measurement of joint space margins based on hand radiographs. The objective of this study is to verify the potential effect of hand rotation during X-ray imaging on the measurement of joint space width using CAJSA and to evaluate the reproducibility of the CAJSA technique in healthy subjects and in patients with rheumatoid arthritis.

METHODS

All joint space distance (JSD) measurements were performed using CAJSA-technology at the metacarpophalangeal articulation based on conventional and digital hand radiographs.

RESULTS

I. Hand rotation showed an effect on the reproducibility with CV ranging between 0.39% (angulation 1°) and 1.66% (angulation 19°). II. Regarding the overall reproducibility of the CAJSA method, the intra-radiograph reproducibility of JSD was calculated with CV=0.54% for conventional images and CV=0.38% for digital images. The inter-radiograph reproducibility error was observed with CV=0.66% (conventional images) and CV=0.63% (digital images). III. The study revealed a reproducibility for CAJSA measurements in RA ranging between 0.37% (JSD-MCP ring finger; van der Heijde-modification of the Sharp method score 1) to 1.37% (JSD-MCP index finger; van der Heijde-modification of the Sharp method score 3).

CONCLUSION

CAJSA measures JSD at the metacarpophalangeal articulation with high reproducibility in healthy subjects and in patients with differing stages of rheumatoid arthritis. Additional findings show that hand rotation during X-ray imaging has an impact on the CAJSA measurements; thus avoiding acquisition of hand radiographs with a rotation error of more than 15 degrees can be recommended.

A novel statistical morphometry imaging method for differentiating long bone geometry: methodological development and application with adolescent idiopathic scoliosis (AIS) patients

Non-invasive quantification of bone shape is crucial in orthopaedic research. The primary objective of this study was to develop an automated statistical morphometry method for comparing the cross-sectional images of normal and diseased bones. The secondary objective involved demonstrating the effectiveness of the proposed method in distinguishing AIS patients from normal controls. This framework is composed of bone segmentation followed by measurements of maximum and minimum bone diameters, inter-group and intra-group statistical morphometry, and statistical analysis of bone thickness. The proposed framework was applied to detect bone morphological abnormality in adolescent idiopathic scoliosis (AIS) patients. The forearm bones in cross-sectional peripheral quantitative computed tomography (pQCT) images from 23 AIS patients and 16 normal controls were analyzed. The radius outer contour was found to be rounder and the radius cortical bone was thinner in AIS patients compared to normal controls.

Patient repositioning reproducibility of joint space width measurements on hand radiographs

OBJECTIVE

Computer-based methods to measure radiographic joint space width (JSW) have the potential to improve the longitudinal assessment of rheumatoid arthritis (RA). The purpose of this report was to measure the long-term patient repositioning reproducibility of software-measured radiographic JSW.

METHODS

Patients underwent baseline and followup hand radiography examinations with a followup time of ≤3 years. To eliminate any JSW change due to real disease progression, the evaluation was performed on "unaffected" joints, defined as having JSW and erosion Sharp scores of 0 at both baseline and followup. The root mean square SD (RMSSD) and coefficient of variation (CV) were used as the reproducibility metrics.

RESULTS

The RMSSD was 0.14 mm (CV 10.5%) for all joints, 0.18 mm (CV 10.9%) for the metacarpophalangeal (MCP) joints, and 0.08 mm (CV 8.3%) for the proximal interphalangeal (PIP) joints. The distribution of JSW change was asymmetric, suggesting that narrowing due to RA progression occurred for several joints. A second analysis was performed, excluding joints where the loss of JSW was greater than 3 SDs. For this analysis, the RMSSD was 0.10 mm (CV 7.5%) for all joints, 0.12 mm (CV 7.3%) for the MCP joints, and 0.07 mm (CV 7.1%) for the PIP joints.

CONCLUSION

Repositioning reproducibility is very good but is likely to be a dominating factor compared to reader and software reproducibility. Additionally, further evidence is given that a software method is able to detect changes in some joints for which the Sharp score is insensitive.

Feasibility study of semi-automated measurements of finger joint space widths

The purpose of this study is to evaluate technical feasibility based on image capturing conditions (film-focus distance (FFD), film sensitivity, film brand, exposure level and tube voltage) that potentially alter radiographs and consequently may influence the semi-automated measurement of joint space distance (JSD) by computer-aided joint space analysis (CAJSA) in rheumatoid arthritis and osteoarthritis. The radiographs of a left hand (deceased man) were acquired under systematically changing image capturing conditions (exposure level: 4-8 mAs; FFD: 90-130 cm; film sensitivity: 200/400 and tube voltage: 40-52 kV with different image modalities: conventional radiographs, original digital radiographs, digital print-outs). All JSD-measurements were performed with the CAJSA-technology (Radiogrammetry Kit, Version 1.3.6; Sectra; Sweden) at the metacarpal-phalangeal articulation. JSD-analysis was not influenced by changes of FFD, exposure level, film sensitivity or film brand. JSD showed significant variation caused by tube voltage (conventional: CV = 1.913% for Agfa and CV = 2.448% for Kodak; digital: CV = 0.741% for Philips print-outs and CV = 0.620% with original digital images versus CV = 2.185% for Siemens print-outs and 0.951% with original digital images). Computer-aided joint space analysis for JSD-measurements is unaffected by the following image capturing parameters: film-focus distance, film sensitivity, film brand and exposure level. An influence of tube voltage was detected in a lesser extent for original digital images compared to the printed digital as well as conventional versions. Consequently, a standardized tube voltage is essential for accurate reproductions of CAJSA-measurements in rheumatoid arthritis and osteoarthritis.

Multiple computer-based methods of measuring joint space width can discriminate between treatment arms in the COBRA trial -- Update of an ongoing OMERACT project

Previously reported data on 5 computer-based programs for measurement of joint space width focusing on discriminating ability and reproducibility are updated, showing new data. Four of 5 different programs for measuring joint space width were more discriminating than observer scoring for change in narrowing in the 12 months interval. Three of 4 programs were more discriminating than observer scoring for the 0-18 month interval. The program that failed to discriminate in the 0-12 month interval was not the same program that failed in the 0-18 month interval. The committee agreed at an interim meeting in November 2007 that an important goal for computer-based measurement programs is a 90% success rate in making measurements of joint pairs in followup studies. This means that the same joint must be measured in images of both timepoints in order to assess change over time in serial radiographs. None of the programs met this 90% threshold, but 3 programs achieved 85%-90% success rate. Intraclass correlation coefficients for assessing change in joint space width in individual joints were 0.98 or 0.99 for 4 programs. The smallest detectable change was < 0.2 mm for 4 of the 5 programs, representing 29%-36% of the change within the 99th percentile of measurements.

Computational radiology in skeletal radiography

Recent years have brought rapid developments in computational image analysis in musculo-skeletal radiology. Meanwhile the algorithms have reached a maturity that makes initial clinical use feasible. Applications range from joint space measurement to erosion quantification, and from fracture detection to the assessment of alignment angles. Current results of computational image analysis in radiography are very promising, but some fundamental issues remain to be clarified, among which the definition of the optimal trade off between automatization and operator-dependency, the integration of these tools into clinical work flow and last not least the proof of incremental clinical benefit of these methods.

Assessement of rheumatic diseases with computational radiology: current status and future potential

In recent years, several computational image analysis methods to assess disease progression in rheumatic diseases were presented. This review article explains the basics of these methods as well as their potential application in rheumatic disease monitoring, it covers radiography, sonography as well as magnetic resonance imaging in quantitative analysis frameworks.

Implementation of Z-Scores as an Age- and Sex-independent Parameter for Estimating Joint Space Widths in Rheumatoid Arthritis

Objective.

To compare normative data of joint space distances (JSD) with the JSD of patients with rheumatoid arthritis (RA) as measured by computer-aided joint space analysis (CAJSA) at the metacarpophalangeal (MCP) articulations, and to differentiate age- and sex-related alterations from the disease-related joint space narrowing.

Methods.

In total, 256 healthy subjects and 248 patients with verified RA (following revised ACR criteria) underwent computerized semiautomated measurements of JSD (CAJSA, version 1.3.6) at the MCP articulation (JSD-MCP) based on digital radiographs. The Z-score, a comparative parameter that differentiates joint space alterations caused by RA-related cartilage destruction from age- and sex-related changes, was calculated.

Results.

Our data showed a relationship between measured joint space widths (MCP total and MCP thumb to little finger) and age for healthy subjects and also the RA group. The RA group revealed an age-related joint space narrowing that was surpassed by the RA-related narrowing of joint space widths classified by Sharp joint space narrowing score and resulting in smaller Z-scores for RA patients.

Conclusion.

The CAJSA technique seems to distinguish age-related JSD changes in healthy volunteers from RA-induced alterations. In addition the Z-score was also able to differentiate RA-dependent narrowing of JSD. Calculation of the Z-scores based on sex- and age-specific reference data may facilitate earlier identification of patients with RA, allowing initiation of a more optimal, individually adapted therapeutic strategy.

Automatic quantification of joint space narrowing and erosions in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic disease that affects and potentially destroys the joints of the appendicular skeleton. The precise and reproducible quantification of the progression of joint space narrowing and the erosive bone destructions caused by RA is crucial during treatment and in imaging biomarkers in clinical trials. Current manual scoring methods exhibit high interreader variability, even after intensive training, and thus, impede the efficient monitoring of the disease. We propose a fully automatic quantitative assessment of the radiographic changes that result from RA, to increase the accuracy, reproducibility, and speed of image interpretation. Initial joint location estimates are obtained by local linear mappings based on texture features. Bone contours are delineated by active shape models comprised of statistical models of bone shape and local texture. These models are refined by snakes which increase the accuracy and allow for a fitting of pathological deviations from the training population. The method then measures joint space widths and detects erosions on the bone contour. Joint space widths are measured with a coefficient of variation of 2%-7% for repeated measurements and erosion detection exhibits an area under the receiver operating characteristic (ROC) curve of 0.89. Model landmarks serve as a reference system along the contour. These landmarks enable the definition of joint regions and more specific follow-up monitoring. The automatic quantification allows for a remote analysis, relevant for multicenter clinical trials, and reduces the workload of clinical experts since parts of the process can be managed by nonexpert personnel.

Articular cartilage thickness of the humeral head: an anatomic study

This study determined the thickness of normal humeral head articular cartilage by anatomic cross section using computer-aided image analysis software. Sixteen adult cadaveric humeral heads were analyzed. Our findings reveal that the thickness of humeral articular cartilage is substantially thinner than articular cartilage found in the knee. The cartilage is thickest in the central portion of the head and becomes progressively thinner towards the periphery. Surgical techniques used to treat pathology in the glenohumeral joint, specifically thermal energy or mechanical debridement, may have deleterious effects on the relatively thin humeral articular cartilage.

Automatic quantification of osteoarthritis in hand radiographs: validation of a new method to measure joint space width

OBJECTIVE

To quantify osteoarthritis (OA) by automatic measurement of the joint space width (JSW) in the finger joints in hand radiographs, and validate this against semi-quantitative scoring.

METHODS

Forty subjects with primary OA at multiple sites were selected, 20 of which had no joint space narrowing (JSN) in the hand, and 20 had moderate to severe JSN. Conventional hand radiographs were acquired and the image set was split into a training and test set, for developing and evaluating the method, respectively. No methods to quantify the JSW automatically in OA exist; therefore a method was developed. First the semi-automatic image analysis program identifies the metacarpal and interphalangeal joints. Subsequently, the corresponding joint margins and measurement interval are determined. Finally the JSW is measured in millimetres as the average distance between the joint margins enclosed by the measurement interval. The success rates of different steps of the image analysis were evaluated and a comparison was made between the JSW and the semi-quantitative Osteoarthritis Research Society International score for JSN.

RESULTS

The success rates for the identification of the joint locations and margins were 97.5% and 64.2%, respectively. The JSW decreased with increasing semi-quantitative scores, but increased as the joint was severely damaged. The agreement between the JSW and the semi-quantitative score was confirmed by a significant correlation (R2=0.54 and P<0.01) while assessing hand pairs instead of individual joints.

CONCLUSION

The JSW measurement closely reflects semi-quantitative scoring of JSN. Therefore, this method is expected to offer a reproducible and accurate measurement of the JSW in OA.

Model-based erosion spotting and visualization in rheumatoid arthritis

RATIONALE AND OBJECTIVES

A method for the automatic detection and the visualization of erosions caused by rheumatoid arthritis is investigated. Erosion-enhanced viewing is a contribution to the computer-aided diagnosis of rheumatoid arthritis. It supports the clinician by providing the automatic marking of erosions and the visualization of any deviations from intact anatomy for a concise reviewing interface.

MATERIALS AND METHODS

A generative appearance model is used to capture the variability of intact bone and erosions. The algorithm marks erosions on hand radiographs using this model, and visualizes these erosions with the help of the residual appearance error after fitting the model built from intact bone texture. The algorithm was evaluated on 17 hand radiographs. The standard of reference was an annotation of the erosions by a musculoskeletal radiologist.

RESULTS

Detection results from the algorithm are reported for a set of 17 radiographs of moderately diseased hands. With a specificity of 84%, the detection of unequivocal erosions achieved a sensitivity of 85%. A receiver operating characteristic analysis yields an area under the curve of 0.92. The visualization provided a clear representation of the erosions as determined by two musculoskeletal radiologists.

CONCLUSION

The automatic spotting of erosions provides promising results, and the visualization of the deviation from healthy anatomy aids clinicians in the evaluation of the erosions and in the reviewing of automatic detection results.

Computerized digital imaging techniques provided by digital X-ray radiogrammetry as new diagnostic tool in rheumatoid arthritis

PURPOSE

Our study evaluates digital x-ray radiogrammetry (DXR) and Radiogrammetry Kit (RK) as a new diagnostic method for the measurement of disease-related osteoporosis including quantification of joint space narrowing dependent on the severity of rheumatoid arthritis (RA).

MATERIALS AND METHODS

A total of 172 unselected patients with RA underwent computerized measurements of bone mineral density (BMD) and metacarpal index (MCI) by DXR, as well as a semiautomated measurement of joint space distances at the metacarpal-phalangeal articulation (JSD-MCP 2-5), both were analyzed from plain radiographs of the nondominant hand.

RESULTS

Correlations between DXR-BMD and DXR-MCI vs. parameters of RK were all significant (0.34 < R < 0.61; p < 0.01). An expected negative association was observed between RK parameters and the different scoring methods (-0.27 < R < -0.59). The maximum relative decrease in BMD vs. MCI as measured by DXR between the highest and lowest RA severity group was -27.7% vs. -27.5% (p < 0.01) for the modified Larsen Score, whereas the minimal value of relative DXR-BMD and DXR-MCI reduction could be documented for the Sharp Erosion Score (-20.8% vs. -26.8%; p < 0.01). The relative reduction of mean JSD-MCP using RK significantly varied from -25.0% (Sharp Erosion Score) to -41.2% (modified Larsen Score). In addition, an excellent reproducibility of DXR and RK could be verified.

CONCLUSION

DXR in combination with RK could be a promising, widely available diagnostic tool to supplement the different scoring methods of RA with quantitative data, allowing an earlier and improved diagnosis and more precision in determining disease progression.

Computerized quantification of joint space narrowing and periarticular demineralization in patients with rheumatoid arthritis based on digital x-ray radiogrammetry

OBJECTIVES

The aim of our work was to evaluate digital x-ray radiogrammetry (DXR) for the quantification of disease-related periarticular demineralization and computerized analysis of joint space distances (JSDA) for the measurement of joint space narrowing as a new diagnostic method for the early detection of joint-associated alterations and for monitoring disease progression in patients with rheumatoid arthritis (RA).

MATERIALS AND METHODS

Digital radiographs in 313 patients with varying severity of RA were performed annually and assessed by 2 radiologists using modified Larsen and also the Sharp scores within an observation period of 3 years. The hand radiographs underwent measurements of bone mineral density (BMD) and metacarpal index (MCI) by DXR, as well as computerized JSDA at the metacarpal-phalangeal articulation (JSD-MCP) for a cross-sectional and longitudinal study design.

RESULTS

Both DXR-BMD (-29.6%; P < 0.01) and DXR-MCI (-31.0%; P < 0.01) revealed a notable reduction dependent on the severity of RA (from grade 1 to grade 5 of the modified Larsen score); the severity dependent decrease of mean JSD-MCP ranged from -31.9% (P < 0.01; Sharp erosion part) to -39.1% (P < 0.01) for the modified Larsen score. Over an observation period of 3 years, a significant decrease of DXR-BMD (-22.3%) and DXR-MCI (-23.3%) as well as JSD-MCP mean (-17.5%) was observed (P < 0.05), whereas an accentuated decline of DXR and JSDA parameters was verified for patients without disease-modifying antirheumatic drugs or methotrexate therapy.

CONCLUSION

Computerized analysis of hand radiographs by DXR and JSDA is a promising approach to assess the severity and to monitor the progression of RA because DXR and JSDA are timely able to measure periarticular demineralization and also narrowing of JSD-MCP dependent on the severity, the medical treatment and the course of RA.

Radiographic joint space width in the fingers of patients with rheumatoid arthritis of less than one year's duration

OBJECTIVE

To determine the radiographic joint space width (JSW) in undamaged metacarpophalangeal (MCP) and proximal interphalangeal (PIP) joints of patients with early rheumatoid arthritis (RA) and to identify important clinical determinants of JSW.

METHODS

Radiographs of patients with RA of <1 year's duration, from an early arthritis cohort at a tertiary care rheumatology clinic, were obtained. JSW was analyzed by joint, finger, age, sex, height, and a number of other clinically relevant variables. Multivariate analysis was also performed, to account for possible confounding between variables.

RESULTS

Thirty-eight patients were included in the study. We found that JSW was greater in the MCP joint than the PIP joint (P < 0.0001). JSW was significantly greater in men (P < 0.0001) and increased with increasing height (P < 0.003), but was not associated with age (P < 0.21). In multivariate analyses, sex was shown to be the most important predictor of JSW.

CONCLUSION

In patients with early RA, MCP and PIP JSW is significantly associated with sex and height. In studies of RA in which JSW measurements are included as an outcome, these differences may need to be accounted for in the analysis.

Application of explicit process of care measurement to rheumatoid arthritis: Moving from evidence to practice

OBJECTIVE

To construct quality measures with measurement validity and meaning for clinicians.

METHODS

We conducted a prospective cohort study of rates of change in disease-modifying antirheumatic drug (DMARD) and/or systemic corticosteroid drug or dose for 568 patients with rheumatoid arthritis (RA) across 6,159 clinical encounters within 12 months to examine how changes in clinical specifications change adherence.

RESULTS

Rates of DMARD change were sensitive to specifications regarding the intensity of disease activity (severe or moderate), duration of specified disease activity, and length of the observation period. Over 12 months, the proportions of 377 patients with severe disease activity observed for 1-month, 2-month, and 3-month time blocks who had a change in DMARD drug or dose were 36%, 57%, and 74%, respectively. Over 12 months, a change in DMARD drug or dose was observed for 44%, 50%, and 68% of 377 patients with severe disease within 3 months, 6 months, and 12 months, respectively, of the patient meeting criteria for severe disease activity. A change in DMARD drug or dose was observed for 21%, 23%, and 34% of 149 patients with moderate disease activity within 3, 6, and 12 months, respectively, of the patient meeting criteria for moderate disease activity.

CONCLUSION

Rates of pharmacologic interventions for patients with moderate and severe RA disease activity vary substantially by intensity and duration of disease activity and by duration of period for observing change. Lack of precision in explicit process criteria could substantially mislead comparisons of quality of care across comparison groups.

Performance of an automated computer-based scoring method to assess joint space narrowing in rheumatoid arthritis: A longitudinal study

OBJECTIVE

To compare the diagnostic performance of a computer-based method for measuring joint space width with the Sharp joint space narrowing (JSN) scoring method in patients with rheumatoid arthritis (RA).

METHODS

A random sample of patients with early RA, for whom sequential hand radiographs and Sharp scores were available, was selected from the National Data Bank for Rheumatic Diseases. Hand joint space width was measured using an automated, computer-based method in random order and with blinding for clinical information. We constructed a receiver operating characteristic curve and compared the diagnostic performance of the computer-based and Sharp methods based on the areas under the curve.

RESULTS

One hundred twenty-nine patients with early RA who underwent serial radiography were included. Changes in the computer-based and Sharp methods were highly correlated (r = 0.75, P < 0.001). The computer-based method was significantly more discriminant than the Sharp JSN subscale. The area under the curve of the computer-based method was 0.96 (95% confidence interval [95% CI] 0.94, 0.99) compared with 0.93 (95% CI 0.89, 0.96) for the Sharp subscale (P = 0.024). At the most discriminant cutoff, specificity of the computer-based method was 88.4% compared with 81.4% for the Sharp subscale (P = 0.11); sensitivity was 87.6% for the computer-based method compared with 82.2% for Sharp subscale (P = 0.19). The signal-to-noise ratio for the computer-based method was 83% compared with 70% for the Sharp subscale (P = 0.013).

CONCLUSION

The computer-based method for measuring joint space width is more discriminant than the semiquantitative Sharp JSN subscale.

Digital x-ray radiogrammetry combined with semiautomated analysis of joint space widths as a new diagnostic approach in rheumatoid arthritis: A cross-sectional and longitudinal study

OBJECTIVE

To evaluate digital x-ray radiogrammetry (DXR) and the Radiogrammetry Kit program as new diagnostic tools for quantifying disease-related periarticular osteoporosis and for measuring joint space narrowing according to the severity and duration of rheumatoid arthritis (RA).

METHODS

Using DXR, we performed computerized calculations of bone mineral density (BMD) and the metacarpal index (MCI) in 258 patients with active RA. Using the Radiogrammetry Kit program, we also performed semiautomated measurements of joint space width (JSW) at the second through the fifth metacarpophalangeal (MCP) joints in these patients.

RESULTS

All correlations between the different parameters of both techniques (BMD and the MCI as measured by DXR and MCP JSW as measured by the Radiogrammetry Kit) were significant (0.36 < or = R < or = 0.63; P < 0.01). As expected, a significant negative association was shown between the different MCP JSW results and the results of all scoring methods (-0.67 < or = R < or = -0.29). The BMD and the MCI measured by DXR both decreased significantly between Steinbrocker stage I and stage IV (by 32.7% and 36.6%, respectively; both P < 0.01). Reductions in the overall (mean) MCP JSW varied from 35.3% (Larsen score) to 52.9% (Steinbrocker stage). Over a period of 6 years, we observed relative decreases in BMD and the MCI as measured by DXR (32.1% and 33.3%, respectively), as well as in the overall (mean) MCP JSW (23.5%), and these were pronounced in early RA (duration <1 year). In addition, excellent reproducibility of DXR and Radiogrammetry Kit parameters was verified (coefficients of variation <1%).

CONCLUSION

DXR with the integrated Radiogrammetry Kit program could be a promising, widely available diagnostic tool for supplementing the different RA scoring methods with quantitative data, thus allowing an earlier and improved diagnosis of RA and more precision in determining disease progression.

Differential diagnosis of rheumatic diseases using conventional radiography

The approach to the differential diagnosis of rheumatic diseases using conventional radiography is systematic and module-oriented, which, with respect to future developments, forms the basis for computer-assisted diagnosis (CAD). The indications follow consensus-based referral criteria and attempts should be made to raise the evidence level of the recommendations. Investigation techniques have been improved in the last few years with the use of digital radiography. New imaging technologies may be available in the future that will achieve at least the same quality of images, while exposing patients to a significantly lower radiation dose. The interpretation of radiographical signs could be enhanced through a correlation with other imaging modalities. Computer-assisted techniques with image processing tools for automated measurements, lesion detection and in the form of expert systems are under development. With conventional radiography embedded in CAD systems, promising options will be available to enhance the differential diagnosis of rheumatic diseases.