Method for determining an examination point for the diaphanoscopic examination of a being and device for realizing the same

• US 6,311,083 B1
• Filed: 05/16/2000
• Issued: 10/30/2001
• Est. Priority Date: 07/21/1997
• Status: Expired due to Fees

• Current Assignee: Siemens AG
• Original Assignee: Siemens AG
• Inventors: Prapavat, Viravuth, Mueller, Gerhard-J., Abraham-Fuchs, Klaus, Beuthan, Juergen
• Primary Examiner(s): Lateef, Marvin M.
• Assistant Examiner(s): QADERI, RUNA S
• Application Number:
  US09/463,151
• Time in Patent Office: 532 Days
• Field of Search: 600/407, 600/310, 600/473, 600/476
• US Class Current: 600/407
• CPC Class Codes: A61B 5/0059 using light, e.g. diagnosis by transillumination, diascopy, fluorescence A61B5/0093, A61B5/1455, A61B5/02416 and A61B5/0261 take precedence
• CPC Class Codes: A61B 5/4528 Joints A61B5/4533, A61B5/4538 take precedence
• CPC Class Codes: G01N 21/49 within a body or fluid

Abstract

In a method and an apparatus for locating an examination size for diaphanoscopic examination of a living subject, a region of the living subject wherein the optimum examination location is suspected is sequentially transilluminated with radiation, preferably in a wavelength range of the optical tissue window, for registration of scattered light distributions in the form of location-related spread functions, particularly point spread functions, and at least one function-specific, location-related feature of each spread function is determined, based on which a position value defining the examination location is determined.

Claims

1. A method for locating an examination site for conducting a diaphanoscopic examination of a living subject, comprising the steps of:

successively transilluminating a region of an examination subject with radiation in a wavelength range of an optical tissue window of said region;

detecting a plurality of scattered light distributions, each represented as location-dependent spread functions, respectively for said transilluminations of said region;

for each of said spread functions, identifying at least one function-specific, location-dependent feature, and thereby obtaining a plurality of determined features; and

from said plurality of determined features, identifying a position value for an examination site within said region for subsequently conducting a diaphanoscopic examination.

2. A method as claimed in claim 1 comprising selecting said function-specific, location-dependent feature from the group consisting of a position value of a maximum and a position value of a center of gravity.

3. A method as claimed in claim 1 wherein the step of determining said examination site comprising averaging the function-specific, location-dependent features of the respective spread functions.

4. A method as claimed in claim 1 wherein the step of determining said examination site comprises weighting said function-specific, location-dependent features of the respective spread functions with respectively different weighting factors, to obtain a plurality of weighted features, and determining said position value of said examination site from said weighted features.

5. A method as claimed in claim 4 comprising determining said weighting factors respectively from at least one further function-specific feature of the spread function to which the weighting factor is applied.

6. A method as claimed in claim 5 wherein each of said spread functions has an overall radiation density and a standard deviation, and comprising determining each weighting factor by multiplying the overall radiation density and the standard deviation of the spread function to which the weighting factor is applied.

7. A method as claimed in claim 5 wherein said weighting factors consist of a highest weighting factor and a plurality of other weighting factors, and comprising the additional step of normalizing all of said other weighting factors relative to said highest weighting factor.

8. A method as claimed in claim 1 wherein said region is successively transilluminated in an irradiation plane, and comprising the additional step of back-transforming said position value into said irradiation plane.

9. A method as claimed in claim 8 comprising, for back transforming said position value into said irradiation plane, determining a further position value and back-transforming said position value into said irradiation plane by shifting said position value dependent on said further position value.

10. A method as claimed in claim 9 wherein said region has an irradiation side, at which said radiation is incoming, and a detection side, at which said radiation is detected after passing through said region, and comprising determining said further position value by the steps of:

identifying a spread function, among said plurality of spread functions, having a highest overall radiation intensity, as a highest-intensity spread function;

identifying a location of a maximum of said highest-intensity spread function at said irradiation side of said region, as a first location, and identifying a location of a maximum of said highest-intensity spread function at said detection side of said region, as a second location; and

determining said further position value as an offset between said first location and said second location.

11. A method as claimed in claim 9 wherein said region has an irradiation side, at which said radiation is incoming, and a detection side, at which said radiation is detected after passing through said region, and comprising determining said further position value by the steps of:

identifying a spread function, among said plurality of spread functions, having a highest overall radiation intensity, as a highest-intensity spread function;

identifying a location of a center of gravity of said highest-intensity spread function at said irradiation side of said region, as a first location, and identifying a location of a center of gravity of said highest-intensity spread function at said detection side of said region, as a second location; and

determining said further position value as an offset between said first location and said second location.

12. A method as claimed in claim 9 wherein said region has an irradiation side, at which said radiation is incoming, and a detection side, at which said radiation is detected after passing through said region, and wherein said further position value is determined by the steps of:

determining a plurality of weighting factors respectively for said function-specific, location-dependent features, said weighting factors including a highest weighting factor;

identifying a spread function having a highest weighting factor, as a highest weighting factor spread function;

identifying a location of a maximum of said highest weighting factor spread function at said irradiation side of said region, as a first location, and identifying a location of a maximum of said highest weighting factor spread function at said detection side of said region, as a second location; and

determining said further position value as an offset between said first location and said second location.

13. A method as claimed in claim 9 wherein said region has an irradiation side, at which said radiation is incoming, and a detection side, at which said radiation is detected after passing through said region, and wherein said further position value is determined by the steps of:

determining a plurality of weighting factors respectively for said function-specific, location-dependent features, said weighting factors including a highest weighting factor;

identifying a spread function having a highest weighting factor, as a highest weighting factor spread function;

identifying a location of a center of gravity of said highest weighting factor spread function at said irradiation side of said region, as a first location, and identifying a location of a center of gravity of said highest weighting factor spread function at said detection side of said region, as a second location; and

determining said further position value as an offset between said first location and said second location.

14. A method as claimed in claim 1 wherein each of said spread functions has a maximum, and comprising the step of smoothing each of said spread functions, at least in a region of the maximum thereof, before determining the respective function-specific location-dependent features of said spread functions.

15. A method as claimed in claim 1 comprising approximating said spread functions, with respective approximation functions, before determining said function-specific, location-dependent features of the respective spread functions.

16. A method as claimed in claim 15 comprising approximating said spread functions with Gaussian approximation functions.

17. A method as claimed in claim 1 comprising representing said scattered light distributions respectively as point spread functions.

18. An apparatus for locating an examination site for conducting a diaphanoscopic examination of a living subject, comprising the steps of:

a radiation unit for successively transilluminating a region of an examination subject with radiation in a wavelength range of an optical tissue window of said region;

a radiation detector disposed for detecting a plurality of scattered light distributions, each represented as location-dependent spread functions, respectively for said transilluminations of said region;

a computer supplied with said spread functions which, for each of said spread functions, identifies at least one function-specific, location-dependent feature, thereby obtaining a plurality of determined features, and which from said plurality of determined features, identifies a position value for an examination site within said region for subsequently conducting a diaphanoscopic examination.

19. An apparatus as claimed in claim 18 wherein said irradiation unit comprises a radiation source which emits said radiation at a single wavelength, and further comprising a positioning unit, operated by said control unit, for moving said radiation source relative to said region to successively transilluminate said region from a plurality of different positions.

20. An apparatus as claimed in claim 18 wherein said irradiation unit comprises a plurality of separate radiation sources, respectively operating at a same wavelength, disposed at a plurality of different positions relative to said region, and individually driven by said control unit.

21. An apparatus as claimed in claim 20 comprising a positioning unit, operated by said control unit, for moving at least one of said radiation sources dependent on said position value.

22. An apparatus as claimed in claim 20 wherein said radiation sources are respectively disposed in said radiation unit in a grid having a grid spacing, and wherein each of said radiation sources emits said radiation in a beam having beam diameter, and wherein a ratio of said grid spacing to said beam diameter is in a range between approximately 1:

1 through approximately 3;

1.

23. An apparatus as claimed in claim 22 wherein said ratio is 2:

1.

24. An apparatus as claimed in claim 22 wherein said grid spacing is in a range between approximately 250 μ

m through approximately 750 μ

m, and wherein said beam diameter is in a range between approximately 125 μ

m through approximately 5000 μ

m.

25. An apparatus as claimed in claim 24 wherein said grid spacing is approximately 500 μ

m.

26. An apparatus as claimed in claim 24 wherein said beam diameter is approximately 250 μ

m.

27. An apparatus as claimed in claim 22 wherein said grid spacing is in a range between approximately 125 μ

m through approximately 500 μ

m, and wherein said beam diameter is in a range between approximately 125 μ

m through approximately 500 μ

m.

28. An apparatus as claimed in claim 27 wherein said grid spacing is approximately 250 μ

m.

29. An apparatus as claimed in claim 27 wherein said beam diameter is approximately 250 μ

m.

30. An apparatus as claimed in claim 18 wherein said irradiation unit comprises a plurality of radiation sources respectively emitting radiation at different wavelengths in said wavelength range.