X-ray image processing

• US 7,315,640 B1
• Filed: 02/21/2000
• Issued: 01/01/2008
• Est. Priority Date: 03/01/1999
• Status: Expired due to Fees

• Current Assignee: Oxford University Innovation Limited (University of Oxford)
• Original Assignee: Mirada Solutions Ltd. (Siemens AG)
• Inventors: Highnam, Ralph, Brady, John Michael
• Primary Examiner(s): Dang; Duy M.
• Application Number:
  US09/914,460
• Time in Patent Office: 2,871 Days
• Field of Search: 382/128, 382/130, 382/131, 382/132, 382/169, 378/98.4, 378/207
• US Class Current: 382/132
• CPC Class Codes: G06T 2207/10116 X-ray image
• CPC Class Codes: G06T 2207/30068 Mammography; Breast
• CPC Class Codes: G06T 5/70 Denoising; Smoothing
• CPC Class Codes: G06T 5/73 Deblurring; Sharpening

Abstract

A method of enhancing and normalizing x-ray images, particularly mammograms, by correcting the image for digitizer blur, glare from the intensifying screen and the anode-heel effect. The method also allows the calculation of the compressed thickness of the imaged breast and calculation of the contribution to the mammograms of the extra focal radiation. The correction of the image for glare from the intensifying screen allows the detection of noise, such as film shot noise, in the image, and in particular the differentiation between such noise and microcalcifications.

Claims

1. In a method of x-ray imaging using an intensifying screen to receive x-rays and emit light to be recorded on an x-ray film, a method of calculating from a film image density an energy which was imparted to the intensifying screen, comprising steps of:

calibrating the film and intensifying screen by measuring a response of the film and intensifying screen to a plurality of different intensities of received x-rays;

fitting a theoretical model of an expected response to the measured response, and using the fitted theoretical model to calculate the imparted energy from image density,

wherein the theoretical model is a serpentine curve of the form;

x²y+a²y−

b²x=0,

where a and b are constants, x is a logarithm of the energy imparted to the intensifying screen and y is the image density.

2. In a method of x-ray imaging using an intensifying screen to receive x-rays and emit light to be recorded on an x-ray film, a method of calculating from a film image density an energy which was imparted to the intensifying screen, comprising steps of:

calibrating the film and intensifying screen by measuring a response of the film and intensifying screen to a plurality of different intensities of received x-rays;

fitting a theoretical model of an expected response to the measured response, and using the fitted theoretical model to calculate the imparted energy from image density,

wherein the response of the film and intensifying screen to a plurality of different intensities of received x-rays is measured by exposing the film-screen combination to x-rays through a lucite step wedge and measuring image density produced by the exposing through different steps of the wedge.

3. Apparatus for calculating from image densities of an x-ray film image the energy which was imparted to an intensifying screen used to receive x-rays and emit light to be recorded on the x-ray film, the apparatus comprising:

fitting means for fitting a theoretical model of the expected response to the x-ray film and intensifying screen to a plurality of different intensities of received x-rays to a measured response, and

calculation means for using the fitted theoretical model to calculate the imparted energy from image density,

wherein the theoretical model is a serpentine curve of the form;

x²y+a²y−

b²x=0,

where a and b are constants, x is a logarithm of the energy imparted to the intensifying screen and y is the image density.

4. A computer program storage medium readable by a computer system and encoded with a computer program for controlling a computer to calculate from image densities of an x-ray film the energy which was imparted to an intensifying screen used to receive x-rays and emit light recorded on the x-ray film, by a method comprising steps of:

fitting a theoretical model of an expected response of the x-ray film and intensifying screen to a plurality of different intensities of received x-rays to a measured response, and using the fitted theoretical model to calculate the imparted energy from image density,

wherein the theoretical model is a serpentine curve of the form;

x²y+a²y−

b²x=0,

where a and b are constants, x is a logarithm of the energy imparted to the intensifying screen and y is the image density.

5. A method of calculating from a mammogram a compressed thickness of an imaged breast, comprising a step of delimiting in the mammogram a region corresponding to a part of the breast which is compressed from a region corresponding to an uncompressed breast edge by detecting smoothness of curves of equal intensity in the mammogram.

6. A method according to claim 5, further comprising steps of detecting a smooth equal intensity curve which is positioned furthest from the breast edge, calculating from intensities along that curve a thickness of fat in the breast which gives those intensities, and calculating the thickness of the compressed breast from that thickness of fat.

7. A method according to claim 5, comprising as an initial step a conversion of the mammogram into an h_int representation representing thicknesses of fat and interesting tissue in regions of the breast contributing to each point in the mammogram.

8. A method according to claim 7, wherein the conversion is performed on the basis of an underestimate of the compressed breast thickness, the conversion being iterated with successively changing estimates of the compressed breast thickness until a smooth delimitation between the region corresponding to the part of the breast which is compressed and the region corresponding to the uncompressed breast edge is found.

9. A method according to claim 5, wherein a predefined smoothness threshold is set to detect smooth curves.

10. A method of calculating a contribution to a mammogram of extra-focal radiation by calculating the compressed thickness of the imaged breast according to the method of claim 5, calculating from an x-ray exposure and an estimate of scattering of radiation an expected intensity in the mammogram along the delimitation between the region corresponding to the part of the breast which is compressed and the region corresponding to the uncompressed breast edge, and comparing the expected intensity to an actual intensity in the mammogram.

11. A computer program storage medium readable by a computer system and encoded with a computer program for controlling a computer to calculate a contribution to a mammogram of extra-focal radiation by calculating the compressed thickness of the imaged breast according to the method of claim 5, calculating from an x-ray exposure and an estimate of scattering of radiation an expected intensity in the mammogram along the delimitation between the region corresponding to the part of the breast which is compressed and the region corresponding to the uncompressed breast edge, and comparing the expected intensity to an actual intensity in the mammogram.

12. Apparatus for calculating from a mammogram a compressed thickness of an imaged breast, comprising a display for displaying a mammogram and means for delimiting in the mammogram a region corresponding to a part of the breast which is compressed from a region corresponding to an uncompressed breast edge by detecting smoothness of curves of equal intensity in the mammogram.

13. Apparatus according to claim 12, further comprising detection means for detecting a smooth equal intensity curve which is positioned furthest from the breast edge, and calculating means for calculating from intensities along that curve a thickness of fat in the breast which gives those intensities, and for calculating the thickness of the compressed breast from that thickness of fat.

14. Apparatus according to claim 12, further comprising conversion means for converting the mammogram into an h_int representation representing thicknesses of fat and interesting tissue in regions of the breast contributing to each point in the mammogram.

15. Apparatus according to claim 14, wherein the conversion means first underestimates the compressed breast thickness, and iterates the conversion with successively changing estimates of the compressed breast thickness until a smooth delimitation between the region corresponding to the part of the breast which is compressed and the region corresponding to the uncompressed breast edge is found.

16. Apparatus according to claim 12, wherein a predefined smoothness threshold is set to detect smooth curves.

17. Apparatus for calculating a contribution to a mammogram of extra-focal radiation comprising apparatus according to claim 12 for calculating the compressed thickness of the imaged breast, and further comprising expected intensity calculating means for calculating from an x-ray exposure and an estimate of scattering of radiation an expected intensity in the mammogram along the delimitation between the region corresponding to the part of the breast which is compressed and the region corresponding to the uncompressed breast edge, and comparing means for comparing the expected intensity to an actual intensity in the mammogram.

18. A computer program storage medium readable by a computer system and encoded with a computer program for controlling a computer to calculate from a mammogram a compressed thickness of an imaged breast by a method comprising a step of delimiting in the mammogram a region corresponding to a part of the breast which is compressed from a region corresponding to an uncompressed breast edge by detecting smoothness of curves of equal intensity in the mammogram.

19. A computer program storage medium according to claim 18, wherein the method further comprises steps of detecting a smooth equal intensity curve which is positioned furthest from the breast edge, calculating from intensities along that curve a thickness of fat in the breast which gives those intensities, and calculating the thickness of the compressed breast from that thickness of fat.

20. A computer program storage medium according to claim 18, wherein the method further comprises as an initial step a conversion of the mammogram into an h_int representation representing thicknesses of fat and interesting tissue in regions of the breast contributing to each point in the mammogram.

21. A computer program storage medium according to claim 20, wherein the conversion is performed on the basis of an underestimate of the compressed breast thickness, the conversion being iterated with successively changing estimates of the compressed breast thickness until a smooth delimitation between the region corresponding to the part of the breast which is compressed and the region corresponding to the uncompressed breast edge is found.

22. A computer program storage medium according to claim 18, wherein a predefined smoothness threshold is set to detect smooth curves.

23. Apparatus for calculating from image densities of an x-ray film image the energy which was imparted to an intensifying screen used to receive x-rays and emit light to be recorded on the x-ray film, the apparatus comprising a processor controlled in accordance with a control program to fit a theoretical model of the expected response of the x-ray film and intensifying screen to a plurality of different intensities of received x-rays to a measured response, and to use the fitted theoretical model to calculate the imparted energy from image density,

wherein the theoretical model is a serpentine curve of the form:

x²y+a²y−

b²x=0,

where a and b are constants, x is a logarithm of the energy imparted to the intensifying screen and y is the image density.

24. Apparatus for calculating from a mammogram a compressed thickness of an imaged breast, the apparatus comprising a processor controlled in accordance with a control program to delimit in the mammogram a region corresponding to a part of the breast which is compressed from a region corresponding to an uncompressed breast edge by detecting smoothness of curves of equal intensity in the mammogram.