Method and system for reducing grid line artifacts in X-ray image

Patent metadata

Field	Value
Publication number	US-10297010-B2
Application number	US-201715459701-A
Country	US
Kind code	B2
Filing date	Mar 15, 2017
Priority date	Mar 17, 2016
Publication date	May 21, 2019
Grant date	May 21, 2019

Abstract

A method for reducing grid line artifacts in an X-ray image is disclosed, which includes acquiring an X-ray image by scanning an object, wherein the X-ray image comprises grid line artifacts; decomposing the X-ray image into a high frequency image and a low frequency image, wherein the high frequency image comprises the grid line artifacts; filtering the high frequency image to reduce the grid line artifacts in the high frequency image so as to obtain a filtered high frequency image; and combining the filtered high frequency image with the low frequency image to reconstruct an output image. A system adopting the above method is also disclosed.

First claim

What is claimed is:

1. A method for reducing grid line artifacts in an X-ray image, comprising:
acquiring an X-ray image by scanning an object, wherein the X-ray image comprises grid line artifacts;
decomposing the X-ray image into a high frequency image and a low frequency image, wherein the high frequency image comprises the grid line artifacts;
filtering the high frequency image to reduce the grid line artifacts in the high frequency image so as to obtain a filtered high frequency image, said filtering comprising:
identifying a strong edge in the high frequency image;
separating information of the strong edge from the high frequency image to create a separated high frequency image, said separating comprising suppressing the strong edge to a lower amplitude to obtain a suppressed high frequency image, and obtaining and storing an image difference between the high frequency image with the identified strong edge and the suppressed high frequency image; and
filtering the separated high frequency image to reduce the grid line artifacts so as to obtain the filtered high frequency image, said filtering comprising filtering the suppressed high frequency image; and

combining the filtered high frequency image with the low frequency image to reconstruct an output image, said reconstructing comprising recovering the information of the strong edge in the filtered high frequency image.

2. The method of claim 1 , wherein decomposing the X-ray image into the high frequency image and the low frequency image comprises:
applying a high-pass filter to the X-ray image so as to obtain the high frequency image; and
subtracting the high frequency image from the X-ray image to obtain the low frequency image.

3. The method of claim 1 , wherein decomposing the X-ray image into the high frequency image and the low frequency image comprises:
applying a low-pass filter to the X-ray image so as to obtain the low frequency image; and
subtracting the low frequency image from the X-ray image to obtain the high frequency image.

4. The method of claim 1 , wherein recovering the information of the strong edge in the filtered high frequency image comprises: adding the stored image difference to the filtered high frequency image.

5. The method of claim 1 , wherein identifying the strong edge comprises:
comparing amplitudes of individual pixels in the high frequency image with a predetermined amplitude threshold; and
regarding regions of pixels whose amplitudes are higher than the predetermined amplitude threshold as the strong edge.

6. The method of claim 5 , wherein identifying the strong edge further comprises:
determining whether the regions of pixels higher than the predetermined amplitude threshold are continuous; and
regarding continuous regions of pixels as the strong edge.

7. The method of claim 1 , wherein filtering the high frequency image to reduce the grid line artifacts comprises:
performing a Fourier transform to convert the high frequency image in a spatial domain to a frequency domain;
applying a band-rejection filter to the high frequency image in the frequency domain so as to reduce the grid line artifacts; and
performing an inverse Fourier transform to convert a filtered high frequency image in the frequency domain back to the spatial domain.

8. A system for reducing grid line artifacts in an X-ray image, comprising:
an image acquisition device configured to scan an object to acquire an X-ray image which comprises grid line artifacts; and
a processor comprising:
an image decomposition module configured to decompose the X-ray image into a high frequency image and a low frequency image, wherein the high frequency image comprises the grid line artifacts;
a filter module configured to filter the high frequency image to reduce the grid line artifacts in the high frequency image so as to obtain a filtered high frequency image;
an image reconstruction module configured to combine the filtered high frequency image with the low frequency image to reconstruct an output image;
an identification module configured to identify a strong edge in the high frequency image; and
a separation module configured to separate information of the strong edge from the high frequency image to create a separated high frequency image, said separation module comprising a suppression module configured to suppress the strong edge to a lower amplitude to obtain a suppressed high frequency image, and a difference image acquiring module configured to obtain and store an image difference between the high frequency image with the identified strong edge and the suppressed high frequency image,
wherein the filter module filters the separated high frequency image to reduce the grid line artifacts so as to obtain the filtered high frequency image, the image reconstruction module is further configured to recover the information of the strong edge in the filtered high frequency image, and filter module filters the suppressed high frequency image to obtain the filtered high frequency image.

9. The system of claim 8 , wherein the image decomposition module comprises:
a high-pass filter configured to filter the X-ray image to obtain the high frequency image; and
a subtractor configured to subtract the high frequency image from the X-ray image to obtain the low frequency image.

10. The system of claim 8 , wherein the image decomposition module comprises:
a low-pass filter configured to filter the X-ray image to obtain the low frequency image; and
a subtractor configured to subtract the low frequency image from the X-ray image to obtain the high frequency image.

11. The system of claim 8 , wherein the image reconstruction module is further configured to add the stored image difference to the filtered high frequency image.

12. The system of claim 8 , wherein the identification module is configured to:
compare amplitudes of individual pixels in the high frequency image with a predetermined amplitude threshold; and
regard regions of pixels whose amplitudes are higher than the predetermined amplitude threshold as the strong edge.

13. The system of claim 12 , wherein the identification module is further configured to:
determining whether the regions of pixels higher than the predetermined amplitude threshold are continuous; and
regarding continuous regions of pixels as the strong edge.

14. The system of claim 8 , wherein the filter module comprises a Fourier transform module, a band-rejection filter and an inverse Fourier transform module, wherein
the Fourier transform module is configured to perform a Fourier transform to convert the high frequency image in a spatial domain to a frequency domain;
the filter module is configured to filter the high frequency image in the frequency domain to reduce the grid line artifacts; and
the inverse Fourier transform module is configured to perform an inverse Fourier transform to convert a filtered high frequency image in the frequency domain back to the spatial domain.

15. The method of claim 1 , further comprising displaying the output image on a display.

16. The system of claim 8 , further comprising a display to display the output image.

Assignees

• Gen Electric

Inventors

• Zhang Zhaoxia
• Tao Kun
• Lai Hao
• Yan Ming
• Kang Han

Classifications

• G06T5/10

  using non-spatial domain filtering · CPC title

• G06T2207/20224

  Image subtraction · CPC title

• G06T2207/10116

  X-ray image · CPC title

• A61B6/5258

  involving detection or reduction of artifacts or noise · CPC title

• G06T5/50 Primary

  using two or more images, e.g. averaging or subtraction · CPC title