Cartilage curvature
View Patent ↗A method for the analysis of three dimensional scan data representing an articular cartilage is provided to extract a quantitative parameter indicative of joint pathology. A measure of local curvature of the cartilage is determined within a region of interest. The value of the quantitative parameter of this joint derived from this measure is compared with the value of a similar quantitative parameter previously established in respect of healthy joints and/or joints characterized by a pathology.
1. A method for the analysis of three dimensional scan data representing an articular cartilage to extract a quantitative parameter indicative of joint pathology, which method comprises determining a measure of the local curvature of the cartilage within a region of interest, and further comprising comparing the measured value of said quantitative parameter for the joint with values of the quantitative parameter previously established in respect of joints characterized by a pathology.
2. The method as claimed in claim 1 , wherein the method comprises determining a measure of the 3D local curvature at a scale arranged to result in a measure of curvature above 40 m −1 .
3. The method as claimed in claim 2 , wherein said scale is arranged to result in a measure of curvature between 40 m −1 and 200 m −1 .
4. The A method as claimed in claim 1 , wherein said region of interest lies in a load bearing area of said cartilage.
5. The method as claimed in claim 1 , wherein the region of interest is selected to lie in an interior region of a cartilage sheet separated from the edges of said sheet.
6. The method as claimed in claim 1 , wherein said cartilage is a knee cartilage.
7. The method as claimed in claim 4 , wherein said knee cartilage is a tibial cartilage.
8. The method as claimed in claim 5 , wherein said tibial cartilage is a medial tibial cartilage.
9. The method as claimed in claim 1 , wherein a comparison of said quantitative parameter is made both with values of the quantitative parameter previously established in respect of healthy joints and with values of said quantitative parameter previously established in respect of joints characterized by a pathology.
10. The method as claimed in claim 6 , wherein said pathology is osteoarthritis, rheumatoid arthritis, villonodular synovitis, Lipoid dematoarthritis (Multicentric reticulohistiocytosis), Enteropathic arthritis, hemophilia (intraarticular bleeding), Gout, Familial Mediterranean fever, Pseudogout, Ochronotic arthropathy, Secondary osteoarthritis, Syphilis (tabes dorsalis), Pyrogenic arthritis, Tuberculous arthritis or Fungal arthritis.
11. The method as claimed in claim 1 , wherein said three dimensional scan data is produced by magnetic resonance imaging (MRI).
12. The method as claimed in claim 8 , wherein said quantitative parameter is derived from changes between adjacent locations within the region of interest in directions defined with respect to normals to a selected cartilage surface or interior contour.
13. The method as claimed in claim 12 , wherein the changes in the directions of normals are weighted by the inverse distance to the local centre point.
14. The method as claimed in claim 1 , wherein said quantitative parameter represents an average of said changes in direction over said region of interest.
15. The method as claimed in claim 12 , wherein said quantitative parameter is quantitatively related to the average value over the region of interest of ‘curvature (p)’ given for each point p in the region of interest by the expression:
curvature
(
p
)
=
1
N
∑
n
∈
N
a
cos
(
normal
(
p
)
·
normal
(
n
)
)
reflecting the mean angle between the surface normal at surface point p with a neighborhood set N and the normals at the neighboring surface points n.
16. The method as claimed in claim 12 , wherein points p are sampled at a spacing of from 0.2 to 2 mm in evaluating said average value of curvature(p).
17. The method as claimed in claim 15 , wherein said neighborhood set N of neighboring points for any point P occupies an area of from 1 mm 2 to 4 mm 2 .
18. The method as claimed in claim 9 , wherein said quantitative parameter is derived from changes in mean curvature flow over the region of interest taken from the maximal and minimal curvature associated with corresponding directions.
19. The method as claimed in claim 15 , wherein said mean curvature flow is calculated by iteratively calculating the average of the absolute value of mean curvature at an articular surface.
20. The method as claimed in claim 1 , wherein said region of interest is from 1 cm 2 to 8 cm 2 in area.
21. The method as claimed in claim 1 , further comprising combining said measure of local curvature with one or more of (a) a measure of cartilage volume, (b) a measure of cartilage thickness, and (c) a measure of cartilage smoothness, to provide an improved discrimination between said combined measure when representing joint pathology and equivalent measures derived from healthy joints.
22. The method as claimed in claim 1 , further comprising analyzing said scan data to perform an automatic segmentation of image data representing cartilage from bone and other background prior to extracting said quantitative parameter.
23. The method as claimed in claim 1 , wherein a value for said quantitative parameter which is higher than that established for healthy joints is taken as indicative of probable joint pathology.
24. The method as claimed in claim 1 , further comprising normalizing the position of the knee prior to segmentation.