Method and system for the computerized radiographic analysis of bone

US 6,205,348 B1
Filed: 04/26/1999
Issued: 03/20/2001
Est. Priority Date: 11/29/1993
Status: Expired due to Fees

First Claim

Patent Images

1. A system for analyzing bone, comprising:

a bone image acquisition device;

an image processing circuit connected to said image acquisition device;

a bone texture measure circuit connected to said image processing circuit, wherein said texture measure circuit comprises first means for determining a root mean square variation and a first moment of a power spectrum and for determining at least one of an angular dependence of said variation and an angular dependence of said first moment.

View all claims

0 Assignments

Timeline View

Assignment View

0 Petitions

Accused Products

Abstract

A computerized method and system for the radiographic analysis of bone structure and risk of future fracture with or without the measurement of bone mass. Techniques including texture analysis for use in quantitating the bone structure and risk of future fracture. The texture analysis of the bone structure incorporates directionality information, for example in terms of the angular dependence of the RMS variation and first moment of the power spectrum of a ROI in the bony region of interest. The system also includes using dual energy imaging in order to obtain measures of both bone mass and bone structure with one exam. Specific applications are given for the analysis of regions within the vertebral bodies on conventional spine radiographs. Techniques include novel features that characterize the power spectrum of the bone structure and allow extraction of directionality features with which to characterize the spatial distribution and thickness of the bone trabeculae. These features are then merged using artificial neural networks in order to yield a likelihood of risk of future fracture. In addition, a method and system is presented in which dual-energy imaging techniques are used to yield measures of both bone mass and bone structure with one low-dose radiographic examination; thus, making the system desirable for screening (for osteoporosis and risk of future fracture).

121 Citations

34 Claims

1. A system for analyzing bone, comprising:
- a bone image acquisition device;
  
  an image processing circuit connected to said image acquisition device;
  
  a bone texture measure circuit connected to said image processing circuit, wherein said texture measure circuit comprises first means for determining a root mean square variation and a first moment of a power spectrum and for determining at least one of an angular dependence of said variation and an angular dependence of said first moment.
- View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24)
- - 2. A system as recited in claim 1, comprising:
3. A system as recited in claim 2, wherein said analysis circuit comprises:
- a neural network adapted to yield a likelihood of fracture.
4. A system as recited in claim 2, comprising:
- a superimposing circuit connected to said analysis circuit; and
  
  a display.
5. A system as recited in claim 1, wherein said processing circuit comprises:
- a region of interest selection circuit connected to said acquisition device; and
  
  a trend correction circuit connected to said selection circuit.
6. A system as recited in claim 1, wherein said first means comprises:
- means for determining said root mean square variation as R, where;
  
  $R = \sqrt{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} V^{2} (u, v) | F (u, v) |^{2} \partial u \partial v}$ means for determining said first moment as M, where;
  
  $M = \frac{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} \sqrt{u^{2} + v^{2}} V^{2} (u, v) {\langle F (u, v) \rangle}^{2} \partial u \partial v}{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} V^{2} (u, v) {\langle F (u, v) \rangle}^{2} \partial u \partial v}$ where;
  
  V(u,v) is a filter function; and
  
  F(u,v) is a power spectrum.
7. A system as recited in claim 1, wherein said first means comprises:
- means for determining an angular dependence R_θ of said variation as;
  
  $R_{E} (θ_{1} ≦ θ < θ_{2}) = \sqrt{\underset{m}{Σ} \underset{n}{Σ} {\langle F_{m, n} \rangle}^{2}}$ $for θ_{1} ≦ \tan^{- 1} (\frac{n}{m}) < θ_{2}$ determining an angular dependence M_θ of said first moment as;
  
  $M_{E} (θ_{1} ≦ θ < θ_{2}) = \frac{\underset{m}{Σ} \underset{n}{Σ} \sqrt{m^{2} + n^{2}} {\langle F_{m, n} \rangle}^{2}}{\underset{m}{Σ} \underset{n}{Σ} {\langle F_{m, n} \rangle}^{2}}$ $for θ_{1} ≦ \tan^{- 1} (\frac{n}{m}) < θ_{2} .$
8. A system as recited in claim 1, wherein said texture measure circuit comprises:
- means for dividing said power spectrum into a number of sectors; and
  
  means for determining said angular dependence of said variation and said angular dependence of said first moment for each of said sectors.
9. A system as recited in claim 7, wherein said texture measure circuit comprises:
- means for determining a minimum of said angular dependence of said first moment.
10. A system as recited in claim 1, wherein said first means comprises:
- means for determining a root mean square variation of said power spectrum;
  
  means for determining an angular dependence of said variation;
  
  means for determining a maximum of said angular dependence;
  
  means for determining a minimum of said angular dependence;
  
  means for determining a standard deviation of said variation;
  
  means for determining a relative standard deviation of said variation;
  
  means for determining a first moment of said power spectrum;
  
  means for determining an angular dependence of said first moment;
  
  means for determining a maximum of said angular dependence;
  
  means for determining a minimum of said angular dependence;
  
  means for determining a standard deviation of said first moment; and
  
  means for determining a relative standard deviation of said first moment.
11. A system as recited in claim 10, comprising:
- an analysis circuit connected to said texture measure circuit and adapted to receive selected of said variation, said angular dependence of said variation, said maximum of said angular dependence, said minimum of said angular dependence, said standard deviation of said variation, said relative standard deviation of said variation said first moment, said angular dependence of said first moment, said maximum of said angular dependence, said minimum of said angular dependence, said standard deviation of said first moment, said relative standard deviation of said first moment, and yield a likelihood of fracture.
12. A system as recited in claim 1, comprising:
- an analysis circuit connected to said texture measure circuit comprising;
  
  means for inputting at least one texture measure determined by said texture measure circuit; and
  
  means for determining a likelihood of risk of fracture using said at least one texture measure.
13. A system as recited in claim 1, wherein said texture measure circuit comprises means for determining texture measures selected from at least one of:
- a root means square variation of said power spectrum;
  
  a first angular dependence of said variation;
  
  a first maximum of said angular dependence;
  
  a first minimum of said angular dependence;
  
  a first standard deviation of said variation;
  
  a first relative standard deviation of said variation;
  
  a first moment of said power spectrum;
  
  a second angular dependence of said first moment;
  
  a second maximum of said angular dependence;
  
  a second minimum of said angular dependence;
  
  a second standard deviation of said first moment; and
  
  a second relative standard deviation of said first moment.
14. A system as recited in claim 13, comprising:
- an analysis circuit connected to said texture measure circuit comprising;
  
  means for merging said texture measures; and
  
  means for determining a likelihood of risk of fracture using said texture measures.
15. A system as recited in claim 1, further comprising:
- a bone mass circuit connected to said region of interest selection circuit.
16. A system as recited in claim 15, wherein said bone mass circuit comprises:
- an integration circuit; and
  
  a weighted sum circuit connected to said integration circuit.
17. A system as recited in claim 15, comprising:
- an analysis circuit connected to said bone mass circuit and said texture measure circuit.
18. A system as recited in claim 15, wherein said bone mass circuit comprises:
- means for integrating a region of interest in a first image obtained by said image acquisition device and for integrating a second region of interest in a second image obtained by said image acquisition device; and
  
  means for determining bone mass using integrated first and second regions of interest.
19. A system as recited in claim 18, wherein said means for determining bone mass comprises a weighted sum circuit.
20. A system as recited in claim 15, wherein said bone mass circuit comprises:
- a integration circuit integrating a first region of interest in a first image to obtain a first integration value, said first image being obtained from said image acquisition device at a first energy level, and integrating a second region of interest in a second image to obtain a second integration value, said second image being obtained from said image acquisition device at a second energy level higher than said first energy level; and
  
  a weighted sum circuit determining first and second weighted sums of said first and second integration values.
21. A system as recited in claim 1, wherein said texture measure circuit comprises:
- means for performing multi-fractal analysis.
22. A system as recited in claim 21, wherein said means comprises:
- means for determining a surface area of a region of interest in an image, obtained from said image acquisition device, at each of a plurality of levels of resolution;
  
  means for determining a slope using said surface areas and said levels of resolution; and
  
  means for determining a fractal dimension using said slope.
23. A system as recited in claim 22, comprising:
- means for determining said fractal dimension as D=2−
  
  H, where H is said slope.
24. A system as recited in claim 22, comprising:
- means for determining a first surface area of said region of interest based upon said pixels;
  
  means for selectively combining said pixels to obtain combined pixels;
  
  means for determining a second surface area based upon said combined pixels; and
  
  means for determining a fractal dimension based upon said first and second surface areas.

25. A system for analyzing bone, comprising:
- a bone image acquisition device;
  
  an image processing circuit connected to said image acquisition device;
  
  a bone texture measure circuit connected to said image processing circuit, wherein said texture measure circuit comprises;
  
  a Fourier transform circuit;
  
  a power spectrum circuit connected to said Fourier transform circuit; and
  
  a first circuit adapted to determine at least one texture measure from an angular dependence of at least one of a root mean square variation and a first moment of a power spectrum determined using said power spectrum circuit.
- View Dependent Claims (26, 27, 28, 29, 30, 31, 32, 33, 34)
- - 26. A system as recited in claim 25, comprising:
27. A system as recited in claim 25, wherein said analysis circuit comprises:
- a neural network adapted to yield a likelihood of fracture.
28. A system as recited in claim 25, comprising:
- a superimposing circuit connected to said analysis circuit; and
  
  a display.
29. A system as recited in claim 25, wherein first circuit is adapted to determine said root mean square variation as R, where:
- $R = \sqrt{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} V^{2} (u, v) | F (u, v) |^{2} \partial u \partial v}$ and to determine said first moment as M, where;
  
  $M = \frac{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} \sqrt{u^{2} + v^{2}} V^{2} (u, v) {\langle F (u, v) \rangle}^{2} \partial u \partial v}{\int_{- \infty}^{\infty} \int_{- \infty}^{\infty} V^{2} (u, v) {\langle F (u, v) \rangle}^{2} \partial u \partial v}$ where;
  
  V(u,v) is a filter function; and
  
  F(u,v) is a power spectrum.
30. A system as recited in claim 25, wherein said first circuit is further adapted to determine an angular dependence R_θ
- of said variation as;
  
  $R_{Θ} (Θ_{1} ≦ Θ < Θ_{2}) = \sqrt{\underset{m}{Σ} \underset{n}{Σ} {\langle F_{m, n} \rangle}^{2}}$ $for Θ_{1} ≦ \tan^{- 1} (\frac{n}{m}) < Θ_{2}$ and to determine an angular dependence M_θof said first moment as;
  
  $M_{Θ} (Θ_{1} ≦ Θ < Θ_{2}) = \frac{\underset{m}{Σ} \underset{n}{Σ} \sqrt{m^{2} + n^{2}} {\langle F_{m, n} \rangle}^{2}}{\underset{m}{Σ} \underset{n}{Σ} {\langle F_{m, n} \rangle}^{2}}$ $for Θ_{1} ≦ \tan^{- 1} (\frac{n}{m}) < Θ_{2} .$
31. A system as recited in claim 25, wherein said first circuit is further adapted to divide said power spectrum into a number of sectors, and to determine said angular dependence of said variation and said angular dependence of said first moment for each of said sectors.
32. A system as recited in claim 25, wherein first circuit is further adapted to determine a minimum of said angular dependence of said first moment.
33. A system as recited in claim 25, wherein said first circuit is adapted to determine:
- a root mean square variation of said power spectrum;
  
  an angular dependence of said variation;
  
  a maximum of said angular dependence;
  
  a minimum of said angular dependence;
  
  a standard deviation of said variation;
  
  a relative standard deviation of said variation;
  
  a first moment of said power spectrum;
  
  an angular dependence of said first moment;
  
  a maximum of said angular dependence;
  
  a minimum of said angular dependence;
  
  a standard deviation of said first moment; and
  
  a relative standard deviation of said first moment.
34. A system as recited in claim 33, comprising:
- an analysis circuit connected to said texture measure circuit and adapted to receive selected of said variation, said angular dependence of said variation, said maximum of said angular dependence, said minimum of said angular dependence, said standard deviation of said variation, said relative standard deviation of said variation said first moment, said angular dependence of said first moment, said maximum of said angular dependence, said minimum of said angular dependence, said standard deviation of said first moment, said relative standard deviation of said first moment, and yield a likelihood of fracture.

Specification

Resources

Litigation Campaign Assessment

Current Assignee
ARCH Development Corporation (University of Chicago)
Original Assignee
ARCH Development Corporation (University of Chicago)
Inventors
Giger, Maryellen L., Doi, Kunio
Primary Examiner(s)
Casler, Brian L.

Application Number

US09/298,852
Time in Patent Office

694 Days
Field of Search

600/407, 600/410, 600/437, 128/920, 128/916, 382/128, 382/131, 382/132, 382/190, 382/195, 382/108
US Class Current

600/407
CPC Class Codes

A61B 6/4042   K-edge filters

A61B 6/482   involving multiple energy i...

A61B 6/505   for diagnosis of bone

A61B 6/583   using calibration phantoms

G06T 2207/10116   X-ray image

G06T 2207/30008   Bone

G06T 2207/30012   Spine; Backbone

G06T 7/0012   Biomedical image inspection

Method and system for the computerized radiographic analysis of bone

First Claim

0 Assignments

0 Petitions

Accused Products

Abstract

121 Citations

34 Claims

Specification

Use Cases

Quick Links

Others

Method and system for the computerized radiographic analysis of bone

First Claim

0 Assignments

Subscription Required

Subscription Required

0 Petitions

Subscription Required

Accused Products

Subscription Required

Abstract

121 Citations

34 Claims

Specification

Subscription Required

Use Cases

Quick Links

Others