PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Skeletal Radiol. Author manuscript; available in PMC 2011 January 10.
Published in final edited form as:
PMCID: PMC3018074
NIHMSID: NIHMS246384

Magnetic resonance image segmentation using semi-automated software for quantification of knee articular cartilage—initial evaluation of a technique for paired scans

M. H. Brem
Department of Radiology, Brigham and Women's Hospital, 75 Francis St., Boston, MA 02115, USA
Division of Orthopaedic and Trauma Surgery, Department of Surgery, Friedrich-Alexander-University Erlangen Nuremberg, Erlangen, Germany
P. K. Lang, G. Neumann, and P. M. Schlechtweg
Department of Radiology, Brigham and Women's Hospital, 75 Francis St., Boston, MA 02115, USA
E. Schneider
SciTrials, LLC, Rocky River, OH, USA
The Cleveland Clinic, Imaging Institute, Cleveland, OH, USA
R. Jackson and J. Yu
Diabetes and Metabolism and Radiology, Department of Endocrinology, The Ohio State University, Columbus, OH, USA
C. B. Eaton
Memorial Hospital of Rhode Island, Center for Primary Care and Prevention and the Warren Alpert Medical School of Brown University, Providence, RI, USA
F. F. Hennig
Division of Orthopaedic and Trauma Surgery, Department of Surgery, Friedrich-Alexander-University Erlangen Nuremberg, Erlangen, Germany

Abstract

Purpose

Software-based image analysis is important for studies of cartilage changes in knee osteoarthritis (OA). This study describes an evaluation of a semi-automated cartilage segmentation software tool capable of quantifying paired images for potential use in longitudinal studies of knee OA. We describe the methodology behind the analysis and demonstrate its use by determination of test–retest analysis precision of duplicate knee magnetic resonance imaging (MRI) data sets.

Methods

Test–retest knee MR images of 12 subjects with a range of knee health were evaluated from the Osteoarthritis Initiative (OAI) pilot MR study. Each subject was removed from the magnet between the two scans. The 3D DESS (sagittal, 0.456 mm×0.365 mm, 0.7 mm slice thickness, TR 16.5 ms, TE 4.7 ms) images were obtained on a 3-T Siemens Trio MR system with a USA Instruments quadrature transmit–receive extremity coil. Segmentation of one 3D-image series was first performed and then the corresponding retest series was segmented by viewing both image series concurrently in two adjacent windows. After manual registration of the series, the first segmentation cartilage outline served as an initial estimate for the second segmentation. We evaluated morphometric measures of the bone and cartilage surface area (tAB and AC), cartilage volume (VC), and mean thickness (ThC.me) for medial/lateral tibia (MT/LT), total femur (F) and patella (P). Test–retest reproducibility was assessed using the root-mean square coefficient of variation (RMS CV%).

Results

For the paired analyses, RMS CV % ranged from 0.9% to 1.2% for VC, from 0.3% to 0.7% for AC, from 0.6% to 2.7% for tAB and 0.8% to 1.5% for ThC.me.

Conclusion

Paired image analysis improved the measurement precision of cartilage segmentation. Our results are in agreement with other publications supporting the use of paired analysis for longitudinal studies of knee OA.

Keywords: Osteoarthritis, Knee, Cartilage, MRI, Segmentation

Osteoarthritis (OA) of the knee is a prevalent disease with a large social and economic cost [1, 2]. Due to the increased lifespan of the population and shifting demographics, the number of individuals with OA is increasing. A 2003 study estimated the cost of arthritis to the United States economy to be over $116 billion in 1997 dollars [1] and the economic problem is of similar magnitude in European countries [2]. New drugs and treatments for OA require evaluation by methods that best reflect disease progression. Radiographic measures of joint space provide only an indirect measurement of cartilage loss and do not incorporate other soft tissue changes such as those in the menisci, and show significant loss only when joint damage is already present [3, 4]. MRI offers a more precise tool for assessing cartilage degradation [5, 6], since cartilage can be visualized directly.

MRI data sets can be assessed visually using qualitative or semi-quantitative methods whereby a reader evaluates the image features according to a consistent set of rules or by comparison to representative cases from an atlas [3, 7-9]. MRI data sets are inherently digital, therefore, it is also possible to use software-based methods to evaluate the images and produce outcome measures for clinical studies. A common class of these software applications is the segmentation algorithm. The goal of segmentation software is to outline the relevant structure in space and identify which pixels or voxels belong to the object. For some clinical applications it is often critical to characterize the shape of the segmented structure to aid in diagnosis. In mammography, for example, analysis software is also employed to determine the level of spiculation for a lesion [10]. For knee OA assessment, the goal of the segmentation software is to outline the articular cartilage in three dimensions and calculate the volume, thickness, and surface area or other biomarkers. There are different approaches and software tools published in the literature to achieve 3D segmentation of the articular cartilage and this is an active area of research [5, 11-16].

A “fully automated” application would accurately delineate the cartilage on each slice in the data set without the need for any reader or observer to guide the process. Such an automated software tool could be considered the ultimate goal of these endeavors; however such a tool is not currently available. Although there has been some effort at creating fully automated segmentation software [17], normal anatomical variation and different degrees of disease status produce a wide range of image appearance for which the application of a single software approach is unlikely to be 100% successful. Therefore, a “semi-automated” method, which uses a skilled reader in combination with image-analysis software, may offer the most practical approach to cartilage segmentation.

We previously published a validation study of a semi-automated software tool which uses a hybrid approach [15] similar to other methods [11]. This tool can measure the cartilage volume (VC), cartilage and bone surface area (AC and tAB), as well as mean and minimum thickness (ThC. me and ThC.min). Additional imaging-based biomarkers can also be computed from the segmented images [18]. With semi-automated analysis, an expert reader with anatomic knowledge is required to guide the software. The software has numerous automated image-processing tools to increase the speed and precision of the cartilage segmentation. We present here the results of a second component of software validation. The initial validation study investigated the re-measurement reproducibility of cartilage volume and thickness when analyzed in a completely blinded manner [15]. While the software tool as reported [15] could be used to evaluate longitudinal studies of knee MRI, it has been previously shown that paired image analysis will reduce the measurement precision errors and this approach has become the standard approach used for longitudinal clinical trials [19]. We hypothesized that a semi-automated software method, which permitted the reader to simultaneously evaluate images series from multiple time points of a study, would increase the precision of the segmentation. To this end, we have modified the existing software tool to permit a reading of paired data sets. In this study, we present a validation of the method using cross-sectional test–retest knee MRI data.

Our study used images from the pilot MR study for the Osteoarthritis Initiative (OAI). The OAI is jointly sponsored by the National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS), several other Institutes at the National Institutes of Health (NIH), and the pharmaceutical industry. The OAI enrolled 4,794 subjects and will acquire MRI exams of both knees on every subject at baseline as well as at a minimum of four annual follow-up exams, providing nearly 48,000 individual knee exams. In addition to the OAI, there are other large-cohort studies and clinical trials using knee MRI that will require fast and low-cost tools to evaluate the substantial number of subjects [20-22]. The cost of a study can be reduced by using a faster method that is equally accurate, reproducible, and sensitive to longitudinal change as manual assessments since a major expense of these studies is quantitative evaluation.

Materials and methods

Data and images were obtained as part of the OAI pilot MR study. The pilot study was performed at the Ohio State University and Memorial Hospital of Rhode Island under review and approval of the local institutional review boards. MRI exams were performed on 12 knees from ten adult subjects (three men and seven women). The participants were selected according to the OAI study design [23] and were participants in other OAI pilot MR studies [15, 19, 24, 25]. All ten subjects met the requirements for the incident cohort of the OAI. The six knees (five subjects) classified as having symptomatic knee OA had pain, aching, or stiffness on the majority of days within a month in the past 12 months, and had OA diagnosed by a physician. Another three knees (three subjects) had only infrequent knee pain, aching, or stiffness over the prior year. Seven of the ten subjects in this pilot study also participated in the OAI (incidence cohort) and underwent fixed-flexion posterior-anterior (P/A) radiography [26] of both knees. The last three knees (two subjects) were asymptomatic and had no prior diagnosis of OA and did not participate in the OAI. Based on the OAI knee radiographs, five subjects had Kellgren–Lawrence grade (KLG) 1, one subject had KLG 2, and one subject had KLG 3. The mean age of the participants was 52.2 years (range, 45–73 years) and the mean body mass index (BMI) 28.2 kg/m2 (range: 21.8–34.6 kg/m2). Each subject was imaged twice to permit measurement of repositioning reproducibility. The subjects were removed from the scanner and table, allowed to walk for about 10 min, and were then repositioned and reimaged.

The images were obtained with a Siemens (Erlangen, Germany) Trio 3-Tesla (3T) MR system and a USA Instruments (Aurora, OH) quadrature transmit/receive extremity coil using a 3D sagittal dual echo in the steady state (DESS) acquisition with water excitation. The slice thickness was 0.7 mm, the in-plane pixel size was 0.365 mm×0.456 mm, TR was 16.5 ms, TE was 4.7 ms, and each acquisition produced 160 images.

The semi-automated software initially reported performed cartilage segmentation on a single knee 3D MRI acquisition [15] using an “active contour model” algorithm [27] to refine the segmented margin on each slice of the 3D image series. The active contour model step serves to provide amoreprecise and objectivesegmentationbymaking automatic modifications to the computer-determined cartilage margins until the full contour best matched the true anatomical edge. The steps are demonstrated in Fig. 1 and below. Segmentation of each successive slice is performed by copying the final segmentation from the previous slice and applying the active contour algorithm. Manual correction of mis-identified cartilage boundaries can be performed at each step. The manual correction was done by two readers with more than 3 years experience of cartilage segmentation.

Fig. 1
Demonstration of articular cartilage segmentation. In a, one mid-condyle slice (N)is shown. In b, the semi-automated edge tracking algorithm produces the final segmentation. c Shows the subsequent slice (N+1) on which the segmented contours from the previous ...

With the new image-analysis method, a matched pair of image sets is analyzed using the previously described technique with some modifications made to improve speed, convenience, and objectivity. The reader first performs a complete segmentation of one 3D image series using the method described above. Once completed, the software then displays a single slice from each acquisition of the pair. A tool is provided that allows the user to scroll through the images from the second 3D acquisition until the displayed images for both scans best match (one dimensional image alignment, without slice interpolation). From this point forward, the image series are synchronized so the reader can scroll through both scans simultaneously. Assuming identical slice thickness, scan prescription, and knee positioning, the two scans will continue to match as the reader scrolls though the images. In practice, minor adjustments are made by the reader during this process to provide the best agreement between the pair of images. Segmentation of the second scan is then accomplished by first copying the correct segmentation from the matched slice of the initial image set and then applying the active contour model. As above, the software tool then applies the initial segmentation to each successive slice, the active contour algorithm performs the correction, and the reader adjusts or corrects any anomalous boundaries.

After both knees have been segmented, the software counts the number of pixels contained within each contour, multiplies by the pixel size, and sums the results from all slices. This value represents the volume of articular cartilage (VC). The bone surface area (tAB) is the cartilage boundary adjacent to the bone; the cartilage surface area (AC) is the opposite cartilage boundary. The cartilage thickness (ThC) is computed by the distance perpendicular from the bone surface to the cartilage surface. After these computations were performed, the root-mean square coefficient of variation (RMS CV%) was used as a metric to quantify the test–retest reproducibility [28].

Results

For all segmentations, the paired analyses RMS CV% ranged from 0.9% to 1.2% for VC, from 0.3% to 0.7% for AC, from 0.6% to 2.7% for tAB and 0.8% to 1.5% for ThC. me (Table 1).

Table 1
Paired reproducibility (RMS CV%) for all subjects. MT, LT, P, and F refer to the medial tibial plateau, lateral tibial plateau, patella, and femur

We also examined the reproducibility for the knees with clinical diagnosis of OA and those from the other volunteers independently. For knees with clinical OA, the RMS CV% ranged from 0.7% to 1.2% for VC, from 0.2% to 1.0% for AC, from 0.7% to 2.8% for tAB and 0.7% to 1.8% for ThC.me (Table 2). For the other knees, the RMS CV% ranged from 0.8% to 1.3% for VC, from 0.3% to 0.6% for AC, from 0.6% to 2.6% for tAB and 0.8% to 1.2% for ThC.me (Table 3). There was no significant difference in reproducibility between the two groups which indicates that the presence of osteophytes and thinning cartilage should not be problematic for analysis.

Table 2
Paired reproducibility (RMS CV%) for knees with clinically defined OA. MT, LT, P, and F refer to the medial tibial plateau, lateral tibial plateau, patella, and femur
Table 3
Paired reproducibility (RMS CV%) for knees without clinically defined OA

Discussion

The goal of our study was to evaluate a software tool to segment test–retest 3D MRI series in a pair-wise manner for use in longitudinal studies of knee OA. We found excellent reproducibility for all cartilage plates. In addition, no substantial difference between the knees with clinical OA and the other knees was found which suggests that the method is not influenced by obscured or ambiguous cartilage margins, or osteophytes that are often present in OA. Our re-measurement reproducibility results compare favorably with other published studies using a paired-reading method [24, 29-32]. Eckstein et al. performed a similar paired analysis of VC, AC, and tAB of the MT and LT using the identical data set [24]. The RMS CV% values are similar and ranged from 1.2% to 2.2% [24] compared to 0.5% to 2.7% for this study.

As with unpaired image segmentation [15], the paired analysis software uses a semi-automated active contour algorithm to trace the cartilage borders. For segmentation of a single 3D image series, the active contour is applied to successive slices in a scan, while for the paired readings, the active contour is applied to matched-slices from different exams of the same subject. For both unpaired and paired image analysis, the semi-automated software algorithm is fast and robust since each segmentation begins with an accurate initial estimate. The active contour model also performs the segmentation according to a mathematically rigorous algorithm, and minimizes the amount of reader correction and hence the impact of their judgment, which improves the reproducibility.

In a previous reproducibility study of similar data using unpaired analysis and our segmentation software, the intrareader reproducibility was considerably higher than the results reported here [15]. This is expected since prior work has shown unpaired precision to be much larger than that of paired segmentation analysis [29, 33]. We have also found reader judgement to be a dominant factor for cartilage segmentation reproducibility and that substantial improvement is possible once this factor is mitigated. Improved reader reproducibility could be accomplished by additional training as well as by increasing the level of software automation. To the extent that the new software is more automated, the results will be less influenced by the reader, and the reproducibility should improve. Although the sensitivity of our segmentation algorithm has not yet been tested on longitudinal MRI exams, the work presented in this report is the first step to validate the methodology. Once the rate of longitudinal cartilage change has been quantified, the results presented here can be used for power calculations for designing clinical studies of OA using MRI.

In longitudinal studies of knee OA, the active contour refinement step should produce a segmentation that reflects any local cartilage changes between scans since it is designed with sufficient flexibility to vary the results and conform to the true bone and cartilage contours. Since the majority of the cartilage morphometry is likely to remain stable, using a paired segmentation as a starting point is an efficient method to segment each slice. Measures derived from a total cartilage plate or even a subregion may not reflect progression precisely since small changes are likely to be insignificant compared to the error associated with the measurement [34]. Therefore, quantifying changes for each defect individually may prove to the most favourable technique to observe progression. In future longitudinal studies, we anticipate using 3D image registration (spatial matching of baseline and follow-up exams) to examine changes in the region adjacent to a cartilage defect.

There was evidence of systematic bias in several of the morphometry measures when comparing the baseline to follow-up values. As an example, we show Fig. 2, a Bland Altman plot for the cartilage surface area (AC). This systematic change in the measurement demonstrates a potential shortcoming of the method, which suggests further improvement to the re-measurement reproducibility may be possible once the source of the bias is understood and corrected.

Fig. 2
Bland–Altman plot of the difference between the two segmentation results of the cartilage surface area (AC) of the medial tibial (MT) plateau. The plot demonstrates the systematic overestimate for the second segmentation. Such systematic biases ...

In agreement with prior studies, we found paired analysis improves the segmentation reproducibility compared to separate independent (unpaired) analysis. In most instances, reader judgement error decreases the reproducibility, especially in cases with cartilage defects and other abnormalities found over the course of OA [15]. Using paired analysis the reader can directly visualize the previous findings and the error can be reduced resulting from judgement of cartilage borders, by making a more consistent evaluation.

In contrast with most prior studies, we found no difference in reproducibility between the knees with and without clinically defined OA [15, 29, 32]. This is anticipated because the paired analysis technique allows the reader can directly visualize the previous findings and systematically select appropriate contours to accommodate osteophytes growth.

One limitation of this study is the relatively small number of subjects; it is difficult to draw definite conclusions since even a single case can skew the results. A further limitation is the use of only one scanner type and pulse sequence, but is not anticipated to be problematic or to require changes to the approach or to the software. However, the applicability of the analysis software to different MR system vendors and pulse sequences should also be evaluated. We have tested the cartilage segmentation technique using high quality data, from a well standardized study, with thin slices and high in-plane spatial resolution. All cartilage segmentation is potentially less precise on lower quality data as has been shown in other studies [30, 35]. Another limitation is that the results from this study were not validated using histological or arthroscopy as this was a pilot study with volunteers and no surgical intervention was included. Since each reader evaluated a separate portion of the data set, we were not able to measure the inter-reader reproducibility. In this study, we did not examine cartilage subregions, which have been shown to be more sensitive to change than the total cartilage plates [19]

The next step in validation of this paired analysis software approach is to perform measurements OA progression rates on longitudinal image data with the reader blinded to time point. Additional development of the core image-processing algorithms should enable further increase of the level of automation, thereby producing a quicker, more objective and robust tool which is free of systematic bias. Furthermore, we plan to add features to the graphical user interface to permit segmentation of more than two visits of a longitudinal study and to implement 3D image registration to track individual cartilage lesions or defects.

Conclusion

Software algorithms and user interfaces allow for fast and reproducible analysis of paired scans for articular cartilage segmentation. Such methods are necessary tools for use in large longitudinal studies of knee OA and are generally applicable to analyze other anatomic features. In addition to VC, AC, tAB, and ThC, it is also possible to measure other image-based metrics using segmented MRI data sets such as those described in Eckstein et al. [18]

Acknowledgments

The Osteoarthritis Initiative (OAI) and this pilot study are conducted and supported by the National Institute of Arthritis and Musculoskeletal and Skin Diseases in collaboration with the OAI Investigators and Consultants. This manuscript has been reviewed by the OAI Publications committee for scientific content and data interpretation. The research reported in this article was supported in part by contracts N01-AR-2-2261, N01-AR-2-2262 and N01-AR-2-2258 from NIAMS. Support for this project was also provided by a contract with the NIAMS intramural program.

This work was also supported by a contract with the NIH/NIAMS intramural program. We would like to thank Raphaela Goldbach-Mansky of the NIH/NIAMS Intramural Research Program for her support in developing early versions of the software.

NIAMS funded this work in part (contracts N01-AR-2-2261, N01-AR-2-2262 and N01-AR-2-2258).

Footnotes

Conflict of Interest Statement ES, RJ, JY, and CBE received direct salary support or had fee for service contracts associated with the OAI. In particular:

ES is the principal of SciTrials, LLC, is the NIAMS OAI Technical Advisor and is under contract to NIAMS for this purpose; RJ and JU are at The Ohio State University that is under contract (N01-AR-2-2261) as a clinical center for the OAI; CBE is at the Memorial Hospital of Rhode Island that is under contract (N01-AR-2-2262) as a clinical center for the OAI.

References

1. Cisternas M, Yelin E, Trupin L, Murphy L, Helmick CG. Direct and indirect costs of arthritis and other rheumatic conditions - United States, 1997. Report. Centers for Disease Control and Prevention; Atlanta: Nov 21, 2003. 2003 Contract No.: Document Number.
2. Pientka L. Arthrose als Volkskrankheit. Klin Forsch. 2000;6(Suppl 2):2–3.
3. Link TM, Steinbach LS, Ghosh S, Ries M, Lu Y, Lane N, et al. Osteoarthritis: MR imaging findings in different stages of disease and correlation with clinical findings. Radiology. 2003;226(2):373–381. [PubMed]
4. Kornaat PR, Bloem JL, Ceulemans RY, Riyazi N, Rosendaal FR, Nelissen RG, et al. Osteoarthritis of the knee: association between clinical features and MR imaging findings. Radiology. 2006;239(3):811–817. [PubMed]
5. Peterfy CG, van Dijke CF, Janzen DL, Gluer CC, Namba R, Majumdar S, et al. Quantification of articular cartilage in the knee with pulsed saturation transfer subtraction and fat-suppressed MR imaging: optimization and validation. Radiology. 1994;192(2):485–491. [PubMed]
6. Biswal S, Hastie T, Andriacchi TP, Bergman GA, Dillingham MF, Lang P. Risk factors for progressive cartilage loss in the knee: a longitudinal magnetic resonance imaging study in forty-three patients. Arthritis Rheum. 2002;46(11):2884–2892. [PubMed]
7. Yoshioka H, Stevens K, Genovese M, Dillingham MF, Lang P. Articular cartilage of knee: normal patterns at MR imaging that mimic disease in healthy subjects and patients with osteoarthritis. Radiology. 2004;231(1):31–38. [PubMed]
8. Peterfy CG, Guermazi A, Zaim S, Tirman PF, Miaux Y, White D, et al. Whole-Organ Magnetic Resonance Imaging Score (WORMS) of the knee in osteoarthritis. Osteoarthr Cartil. 2004;12(3):177–190. [PubMed]
9. Hunter DJ, Lo GH, Gale D, Grainger AJ, Guermazi A, Conaghan PG. The reliability of a new scoring system for knee osteoarthritis MRI and the validity of bone marrow lesion assessment: BLOKS (Boston Leeds Osteoarthritis Knee Score) Ann Rheum Dis. 2008;67(2):206–211. [PubMed]
10. Baker JA, Kornguth PJ, Lo JY, Williford ME, Floyd CE., Jr Breast cancer: prediction with artificial neural network based on BIRADS standardized lexicon. Radiology. 1995;196(3):817–822. [PubMed]
11. Stammberger T, Eckstein F, Michaelis M, Englmeier KH, Reiser M. Interobserver reproducibility of quantitative cartilage measurements: comparison of B-spline snakes and manual segmentation. Magn Reson Imaging. 1999;17(7):1033–1042. [PubMed]
12. Solloway S, Hutchinson CE, Waterton JC, Taylor CJ. The use of active shape models for making thickness measurements of articular cartilage from MR images. Magn Reson Med. 1997;37(6):943–952. [PubMed]
13. Kshirsagar AA, Watson PJ, Tyler JA, Hall LD. Measurement of localized cartilage volume and thickness of human knee joints by computer analysis of three-dimensional magnetic resonance images. Invest Radiol. 1998;33(5):289–299. [PubMed]
14. McWalter EJ, Wirth W, Siebert M, von Eisenhart-Rothe RM, Hudelmaier M, Wilson DR, et al. Use of novel interactive input devices for segmentation of articular cartilage from magnetic resonance images. Osteoarthr Cartil. 2005;13(1):48–53. [PubMed]
15. Duryea J, Neumann G, Brem MH, Koh W, Noorbakhsh F, Jackson RD, et al. Novel fast semi-automated software to segment cartilage for knee MR acquisitions. Osteoarthritis Cartilage. 2007;15(5):487–492. [PubMed]
16. Raynauld JP, Kauffmann C, Beaudoin G, Berthiaume MJ, de Guise JA, Bloch DA, et al. Reliability of a quantification imaging system using magnetic resonance images to measure cartilage thickness and volume in human normal and osteoarthritic knees. Osteoarthr Cartil. 2003;11(5):351–360. [PubMed]
17. Fripp J, Crozier S, Warfield SK, Ourselin S. Automatic segmentation of the bone and extraction of the bone-cartilage interface from magnetic resonance images of the knee. Phys Med Biol. 2007;52(6):1617–1631. [PubMed]
18. Eckstein F, Ateshian G, Burgkart R, Burstein D, Cicuttini F, Dardzinski B, et al. Proposal for a nomenclature for magnetic resonance imaging based measures of articular cartilage in osteoarthritis. Osteoarthr Cartil. 2006;14(10):974–983. [PubMed]
19. Wirth W, Hellio Le Graverand MP, Wyman BT, Maschek S, Hudelmaier M, Hitzl W, et al. Regional analysis of femorotibial cartilage loss in a subsample from the Osteoarthritis Initiative progression subcohort. Osteoarthr Cartil. 2008 Sep 11; (in press) [PMC free article] [PubMed]
20. Hunter DJ, Zhang YQ, Niu JB, Felson DT, Kwoh K, Newman A, et al. Patella malalignment, pain and patellofemoral progression: the Health ABC Study. Osteoarthritis Cartilage. 2007;15(10):1120–1127. [PMC free article] [PubMed]
21. Neogi T, Nevitt M, Niu J, LaValley MP, Hunter DJ, Terkeltaub R, et al. Lack of association between chondrocalcinosis and increased risk of cartilage loss in knees with osteoarthritis: results of two prospective longitudinal magnetic resonance imaging studies. Arthritis Rheum. 2006;54(6):1822–1828. [PubMed]
22. Neogi T, Booth SL, Zhang YQ, Jacques PF, Terkeltaub R, Aliabadi P, et al. Low vitamin K status is associated with osteoarthritis in the hand and knee. Arthritis Rheum. 2006;54(4):1255–1261. [PubMed]
24. Eckstein F, Kunz M, Hudelmaier M, Jackson R, Yu J, Eaton CB, et al. Impact of coil design on the contrast-to-noise ratio, precision, and consistency of quantitative cartilage morphometry at 3 Tesla: a pilot study for the osteoarthritis initiative. Magn Reson Med. 2007;57(2):448–454. [PubMed]
25. Eckstein F, Maschek S, Wirth W, Hudelmaier M, Hitzl W, Wyman B, et al. One year change of knee cartilage morphology in the first release of participants from the Osteoarthritis Initiative Progression Subcohort—association with sex, Body Mass Index, symptoms, and radiographic OA status. Ann Rheum Dis. 2008 Jul 7; (in press) [PMC free article] [PubMed]
26. Kothari M, Guermazi A, von Ingersleben G, Miaux Y, Sieffert M, Block JE, et al. Fixed-flexion radiography of the knee provides reproducible joint space width measurements in osteoarthritis. Eur Radiol. 2004;14(9):1568–1573. [PubMed]
27. Kass M, Witkin A, Terzopoulos D. Snakes: active contour models. Int J Comput Vis. 1987;1(4):321–331.
28. Gluer CC, Blake G, Lu Y, Blunt BA, Jergas M, Genant HK. Accurate assessment of precision errors: how to measure the reproducibility of bone densitometry techniques. Osteoporos Int. 1995;5(4):262–270. [PubMed]
29. Eckstein F, Heudorfer L, Faber SC, Burgkart R, Englmeier KH, Reiser M. Long-term and resegmentation precision of quantitative cartilage MR imaging (qMRI) Osteoarthr Cartil. 2002;10(12):922–928. [PubMed]
30. Eckstein F, Charles HC, Buck RJ, Kraus VB, Remmers AE, Hudelmaier M, et al. Accuracy and precision of quantitative assessment of cartilage morphology by magnetic resonance imaging at 3.0T. Arthritis Rheum. 2005;52(10):3132–3136. [PubMed]
31. Glaser C, Burgkart R, Kutschera A, Englmeier KH, Reiser M, Eckstein F. Femoro-tibial cartilage metrics from coronal MR image data: technique, test-retest reproducibility, and findings in osteoarthritis. Magn Reson Med. 2003;50(6):1229–1236. [PubMed]
32. Brem MH, Pauser J, Yoshioka H, Brenning A, Stratmann J, Hennig FF, et al. Longitudinal in vivo reproducibility of cartilage volume and surface in osteoarthritis of the knee. Skeletal Radiol. 2007;36(4):315–320. [PubMed]
33. Eckstein FCF, Raynauld JP, Waterton JC, Peterfy CG, editors. Magnetic resonance imaging (MRI) of articular cartilage in knee osteoarthritis (OA): Morphological assessment. Osteoarthr Cartil. 2006 [PubMed]
34. Koo S, Gold GE, Andriacchi TP. Consideration in measuring cartilage thickness using MRI: factors including reproducibility and accuracy. Osteoarthr Cartil. 2005;13:782–789. [PubMed]
35. Eckstein F, Buck RJ, Burstein D, Charles HC, Crim J, Hudelmaier M, et al. Precision of 3.0 Tesla quantitative magnetic resonance imaging of cartilage morphology in a multicentre clinical trial. Ann Rheum Dis. 2008;67(12):1683–1688. [PubMed]