Method for processing tomosynthesis acquisitions in order to obtain a representation of the contents of an organ

• US 9,842,415 B2
• Filed: 12/22/2010
• Issued: 12/12/2017
• Est. Priority Date: 12/22/2009
• Status: Active Grant

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Bernard, Sylvain, Iordache, Razvan, Souchay, Henri
• Primary Examiner(s): Kao, Glen
• Application Number:
  US12/975,583
• Publication Number:
  US 20110150178A1
• Time in Patent Office: 2,547 Days
• Field of Search: 378 21- 26, 382131
• CPC Class Codes: G06T 11/008 Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
• CPC Class Codes: G06T 2207/10112 Digital tomosynthesis [DTS]
• CPC Class Codes: G06T 2211/436 Limited angle

Abstract

A method for processing tomosynthesis images of an object of interest, using an imaging system, the imaging system comprising an X-ray source positioned facing a detector on which the object of interest is positioned. With the method of an embodiment of the invention, it is possible to display a three-dimensional (3D) reconstruction slice, as well as a two-dimensional (2D) image of the object of interest.

Claims

1. A method of processing tomosynthesis images of an object of interest, using an imaging system, the imaging system comprising an X-ray source positioned facing a detector on which the object of interest is positioned, the processing comprising the following steps:

acquiring a plurality of 2D projection images of the object of interest according to a plurality of orientations located with respect to a perpendicular to the detector, a zero orientation being the closest to the perpendicular;

applying a filter on the acquired 2D projection images so as to obtain filtered projection images of the object of interest;

reconstructing slices of the object of interest from at least two of the filtered projection images, the reconstructed slices being a filtered reconstructed volume of the object of interest;

reprojecting the reconstructed slices according to a determined orientation direction to obtain an intermediate 2D image of the object of interest; and

using a pixel-to-pixel linear combination of the intermediate 2D image obtained by the reprojection and at least one of the acquired plurality of 2D projection images corresponding to the determined orientation to obtain a final 2D image of the object of interest.

2. The method according to claim 1, wherein the reprojection of the reconstructed slices is a Sorted Intensity Pixel (SIP) reprojection in the determined orientation direction, the SIP reprojection comprising, for each pixel of the intermediate 2D image:

determining, in the filtered reconstructed volume, an intensity voxel calculated by using sorting of the values of pixels along the ray from the source to the pixel of the intermediate 2D image; and

storing, in a memory unit of the imaging system, an identifier of the reconstructed slice in which is found the voxel from the sorting, so as to have in the memory unit, depth information connecting each pixel of the intermediate 2D image to the associated reconstructed slice from which this pixel stems.

3. The method according to claim 1, wherein the reprojection of the reconstructed slices is a Maximum Intensity Pixel (MIP) reprojection in the determined orientation direction, the MIP reprojection comprising, for each pixel of the intermediate 2D image:

determining, in the filtered reconstructed volume, the maximum intensity voxel along the ray from the source to the pixel of the intermediate 2D image; and

storing, in a memory unit of the imaging system, an identifier of the reconstructed slice in which the maximum intensity voxel is found, so as to have in the memory unit, depth information connecting each pixel of the intermediate 2D image to the associated reconstructed slice from which this pixel stems.

4. The method according to claim 1, wherein the reprojecting of the reconstructed slices is performed by sorting of values of pixels according to their intensity.

5. The method according to claim 4, wherein a voxel having the highest probability of belonging to a lesion along the ray from the source to the pixel is selected.

6. The method according to claim 5, comprising a step for determining a volume, obtained using a system for automatic detection of lesions, wherein each voxel represents a probability of belonging to a lesion.

7. A display method for displaying a set of radiography images of the object of interest by a unit for displaying radiography images, radiography images being the reconstructed slices of the object of interest and a final 2D image of the object of interest obtained by a method of processing tomosynthesis images according to claim 1, the method of processing tomosynthesis images further including the step of:

storing depth information (argMax′

) in a memory unit, the display method comprising the steps of;

displaying the final 2D image on the display unit,

selecting a pixel from the final 2D image;

searching, among the depth information (argMax′

) stored in the memory unit, for an identifier of the reconstructed slice associated with the selected pixel;

adjusting the size of the reconstructed slice associated with the selected pixel to the size of the final 2D image; and

displaying by means of the display unit, the adjusted reconstructed slice associated with the selected pixel.

8. The method of claim 1 further comprising:

displaying the final 2D image obtained from the combination of the tomosynthesis projection image and the reprojection image;

selecting a pixel from the final 2D image;

searching among a depth information for an identifier of a reconstructed slice associated with the selected pixel; and

displaying the reconstructed slice.

9. The method of claim 8, wherein the size of the reconstructed slice is adjusted to the size of the final 2D image.

10. The method of claim 9, wherein an adjustment factor for the adjusting depends on the height of the reconstructed slice.

11. The method of claim 8, wherein the size of the final 2D image is adjusted to the size of the reconstructed slice.

12. The method of claim 11, wherein an adjustment factor for the adjusting depends on the height of the reconstructed slice.

13. The method of claim 8, wherein the final 2D image and the reconstructed slice are displayed concurrently.

14. The method of claim 8 further comprising toggling between the reconstructed slice and the final 2D image.

15. The method of claim 8, wherein the depth information is stored in a memory unit.

16. A method of displaying radiography images of an object of interest using a display unit and a processing unit comprising:

displaying a final 2D image;

selecting a pixel from the final 2D image, wherein the final 2D image is obtained from the combination of a 2D projection acquired in a determined orientation and a reprojection of reconstructed slices in the corresponding determined orientation;

searching among a depth information for an identifier of a reconstructed slice associated with the selected pixel;

adjusting the size of the identified reconstructed slice to be the same size as the displayed final 2D image; and

displaying the adjusted reconstructed slice.

17. The method of claim 16, wherein the final 2D image and the reconstructed slice of the object of interest were processed using an imaging system, the imaging system comprising an X-ray source positioned facing a detector on which the object of interest is positioned, the processing comprising the following steps:

acquiring a plurality of 2D projection images of the object of interest according to a plurality of orientations located with respect to a perpendicular to the detector where at least one of the plurality of orientations includes the determined orientation, a zero orientation being the closest to the perpendicular;

applying a filter on the acquired 2D projection images so as to obtain filtered projection images of the object of interest;

reconstructing the slices of the object of interest from at least two of the filtered projection images, the reconstructed slices being a filtered reconstructed volume of the object of interest;

obtaining an intermediate 2D image of the object of interest by reprojecting the reconstructed slices according to the determined orientation; and

using a pixel-to-pixel linear combination of the intermediate 2D image obtained by the reprojection and the at least one of the acquired plurality of 2D projection images corresponding to the determined orientation to obtain the final 2D image of the object of interest.