Guided Noise Reduction with Streak Removal for High Speed C-Arm CT

• US 20150279084A1
• Filed: 04/25/2014
• Published: 10/01/2015
• Est. Priority Date: 04/01/2014
• Status: Abandoned Application

• Current Assignee: Siemens AG
• Original Assignee: Yu Deuerling-Zheng
• Inventors: Deuerling-Zheng, Yu, Manhart, Michael
• Application Number:
  US14/262,277
• Publication Number:
  US 20150279084A1
• Time in Patent Office: Days
• US Class Current: 1/1
• CPC Class Codes: G06T 2207/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/30016 Brain
• CPC Class Codes: G06T 2207/30101 Blood vessel; Artery; Vein; Vascular
• CPC Class Codes: G06T 5/50 using two or more images, e.g. averaging or subtraction
• CPC Class Codes: G06T 5/70 Denoising; Smoothing

Abstract

In order to minimize streak artifacts within reconstructions of perfusion maps generated from moving C-arm acquisitions, a streak reduction method includes threshold analysis, time-contrast curve analysis, and total variation analysis. One or more mask volumes and a plurality of contrast agent enhanced volumes are generated based on a plurality of projections generated using a moving C-arm X-ray device. A maximum contrast attenuation volume is generated based on the plurality of contrast agent enhanced volumes. Voxels are identified as streaks based on analyses applied to the one or more mask volumes, the plurality of contrast agent enhanced volumes, and the maximum contrast attenuation volume, respectively.

Claims

1. A method for artifact removal within image data, the method comprising:

identifying, by a processor, a first three dimensional (3D) dataset, the first 3D dataset comprising image data representing an object without a contrast agent;

identifying, by the processor, a plurality of second 3D datasets, each second 3D dataset of the plurality of second 3D datasets comprising a plurality of voxels representing the object with the contrast agent;

generating, by the processor, a third 3D dataset based on the plurality of second 3D datasets;

segmenting, by the processor, a subset of data from the first 3D dataset, the segmented subset of data corresponding to a subset of voxels of the plurality of voxels;

generating, by the processor, time attenuation curves (TACs) based on the plurality of second 3D datasets;

identifying, by the processor, one or more artifacts within the third 3D dataset based on one or more thresholds and the generated TACs for voxels corresponding to locations of voxels of the subset; and

removing, by the processor, the one or more identified artifacts from the third 3D dataset.

2. The method of claim 1, wherein the identifying of the one or more artifacts comprises classifying a first portion of voxels and a second portion of voxels within the subset of voxels of the third 3D dataset based on at least two radiodensity thresholds.

3. The method of claim 2, wherein the classifying is based on a difference between a peak value and a value from which an increase to the peak value begins for each of the TACs.

4. The method of claim 3, wherein the classifying is based on a total variation threshold.

5. The method of claim 3, wherein the first portion of voxels is classified as data representing the one or more artifacts, and the second portion of voxels is classified as data representing vessels.

6. The method of claim 5, further comprising combining the first portion of voxels and the second portion of voxels, the combining comprising performing a dilation operation on the second portion of voxels, and performing an erosion operating and a dilation operation on the first portion of voxels.

7. The method of claim 1, wherein the removing comprises smoothing with a truncated Gaussian kernel averaging.

8. The method of claim 1, wherein the segmenting comprises categorizing data of the first 3D dataset based on two radiodensity thresholds, the categorized data comprising a first portion of data, a second portion of data, and a third portion of data, and

wherein the first portion of data represents air, the second portion of data represents bone, and the third portion of data represents brain tissue.

9. The method of claim 8, wherein the segmented subset of data from the first 3D dataset corresponds to the categorized third portion of data

10. The method of claim 1, further comprising generating a plurality of fourth 3D datasets, the generating of the plurality of fourth 3D datasets comprising subtracting the first 3D dataset from each second 3D dataset of the plurality of second 3D datasets,

wherein generating the third 3D dataset comprises generating the third 3D dataset based on a maximum contrast attenuation across the plurality of fourth 3D datasets for each voxel of the plurality of voxels.

11. In a non-transitory computer-readable storage medium that stores instructions executable by one or more processors to identify and remove artifacts within image data, the instructions comprising:

generating a first three dimensional (3D) dataset, the first 3D dataset comprising image data representing an object without a contrast agent;

generating a plurality of second 3D datasets, each second 3D dataset of the plurality of second 3D datasets comprising a plurality of voxels representing the object with the contrast agent;

generating a third 3D dataset based on a maximum over the plurality of second 3D datasets for each voxel of the plurality of voxels;

generating time attenuation curves (TACs) based on the plurality of second 3D datasets;

identifying one or more artifacts within a subset of voxels of the third 3D dataset based on one or more thresholds and the generated TACs; and

at least partially removing the one or more identified artifacts from the third 3D dataset.

12. The non-transitory computer-readable storage medium of claim 11, wherein the instructions further comprise registering the plurality of second 3D datasets with the first 3D dataset.

13. The non-transitory computer-readable storage medium of claim 12, wherein the instructions further comprise generating subtracted 3D datasets, the generating of the subtracted 3D datasets comprising subtracting the first 3D dataset from each second 3D dataset of the plurality of second 3D datasets,

wherein generating the third 3D dataset comprises generating the third 3D dataset based on the maximum contrast attenuation over the subtracted 3D datasets for each voxel of the plurality of voxels.

14. The non-transitory computer-readable storage medium of claim 13, wherein the instructions further comprise:

generating an initial third 3D dataset based on maximum contrast attenuation over the plurality of subtracted 3D datasets for each voxel of the plurality of voxels;

denoising the initial third 3D dataset;

denoising the subtracted 3D datasets; and

generating an updated third 3D dataset based on maximum contrast attenuation over the plurality of denoised subtracted 3D datasets for each voxel of the plurality of voxels,

wherein the updated third 3D dataset corresponds to the third 3D dataset.

15. The non-transitory computer-readable storage medium of claim 14, wherein the denoising of the initial third 3D dataset comprises bilateral filtering, and

wherein the denoising of the subtracted 3D datasets comprises bilateral filtering, the bilateral filtering of the subtracted 3D datasets being based on the initial third 3D dataset.

16. A system for artifact removal within computed tomography (CT) image data, the system comprising:

a processor configured to;

identify a first three dimensional (3D) dataset, the first 3D dataset comprising image data representing an object without a contrast agent;

identify a plurality of second 3D datasets, each second 3D dataset of the plurality of second 3D datasets comprising a plurality of voxels representing the object with the contrast agent;

generate a third 3D dataset based on a maximum over the plurality of second 3D datasets for each voxel of the plurality of voxels;

segment a subset of data from the first 3D dataset, the segmented subset of data corresponding to a subset of voxels of the plurality of voxels;

generate time attenuation curves (TACs) based on the plurality of second 3D datasets;

identify one or more artifacts within the third 3D dataset based on one or more thresholds and the generated TACs for voxels corresponding to locations of voxels of the subset; and

remove the one or more identified artifacts from the third 3D dataset;

a memory operatively connected to the processor and configured to store the third 3D dataset from which the one or more identified artifacts have been removed; and

a display operatively connected to the memory and configured to display the third 3D dataset from which the one or more identified artifacts have been removed.

17. The system of claim 16, further comprising a C-arm X-ray device comprising a C-arm and an X-ray detector attached to the C-arm,

wherein the X-ray detector is configured to:

generate a plurality of first projections over an angular range of the C-arm; and

generate a plurality of second projections over a plurality of angular ranges of the C-arm, and

wherein the processor is configured to;

reconstruct the first 3D dataset based on the plurality of generated first projections; and

reconstruct the plurality of second 3D datasets based on the plurality of generated second projections.

18. The system of claim 17, wherein the angular range and each angular range of the plurality of angular ranges is 200°

.

19. The system of claim 18, wherein the plurality of angular ranges comprises multiple consecutive rotations of the C-arm, and

wherein the multiple consecutive rotations alternate forward and backward, respectively.

20. The system of claim 17, wherein the C-arm is operable to rotate up to 120°

/s.

21. A method for artifact removal within computed tomography (CT) image data, the method comprising:

generating, by a processor, a 3D data set representing maximum contrast attenuation over time of a contrast agent in an object;

generating, by the processor, time attenuation curves (TACs) of the contrast agent in the object for voxels corresponding to the 3D data set;

identifying, by the processor, one or more streak artifacts in the 3D data set based on one or more thresholds and the generated TACs; and

removing, by the processor, the one or more identified artifacts from the 3D dataset.