Method, computer system and CT system for determining a motion field and for motion-compensated reconstruction using said motion field

Patent metadata

Field	Value
Publication number	US-9060733-B2
Application number	US-201213626184-A
Country	US
Kind code	B2
Filing date	Sep 25, 2012
Priority date	Sep 28, 2011
Publication date	Jun 23, 2015
Grant date	Jun 23, 2015

Abstract

A method is disclosed for determining a motion field in an image space of a CT image data set for a partially and cyclically moving object, wherein a projection data set from a CT system covering a motion phase and a projection angle range is recorded, and an iterative determination of the motion field by minimizing a cost function from CT image data sets which are motion-compensated and reconstructed making use of motion fields and which have a first image resolution, with a first analytical reconstruction algorithm, the cost function being based on raw data or transformed into the image space. A method for generating a motion-compensated CT image data set for a partially and cyclically moving object is also disclosed, along with a computer system for image reconstruction and a CT system.

First claim

What is claimed is:

1. A method for determining a motion field, including a plurality of location-specific motion vectors in an image space of a CT image data set for a partially and cyclically moving object, the method comprising:
detecting or transferring a projection data set from a computed tomography system covering a motion phase and a projection angle range, to enable direct reconstruction of a CT image data set;
iteratively determining the motion field by minimizing a cost function from CT image data sets which are motion-compensated and reconstructed making use of motion fields and which have a first image resolution, with a first analytical reconstruction algorithm, the cost function being based on raw data or transformed into the image space; and
at least one of storing and outputting the motion field;
wherein, as the raw data-based cost function, the equation min arg L=∥p in −A( , )f( , )∥ 2 is used.

2. The method of claim 1 , wherein the projection angle range is 180° plus the fan-beam angle of a ray beam used.

3. The method of claim 2 , wherein the analytical reconstruction method is a method from the following list:
FDK-reconstruction methods (FDK=Feldmann-Davis-Kress),
Clack-Defrise-reconstruction methods,
reconstruction methods based on Hilbert transformations,
reconstruction methods based on Fourier transforms, or
Radon methods.

4. The method of claim 1 , wherein the analytical reconstruction method is a method from the following list:
FDK-reconstruction methods (FDK=Feldmann-Davis-Kress),
Clack-Defrise-reconstruction methods,
reconstruction methods based on Hilbert transformations,
reconstruction methods based on Fourier transforms, or
Radon methods.

5. The method of claim 4 , wherein, as the image property to be optimized, at least one of the following image properties is used:
entropy,
L 2 -norm of the difference between the reconstructed raw data and the “synthetic” data determined with motion-sensitive forward projection,
total variation/total fluctuation, and
compressibility.

6. The method of claim 1 , wherein, as the image property to be optimized, at least one of the following image properties is used:
entropy,
L 2 -norm of the difference between the reconstructed raw data and the “synthetic” data determined with motion-sensitive forward projection,
total variation/total fluctuation, and
compressibility.

7. The method of claim 1 , wherein the detector data for generating the used projection data sets are collected from a plurality of motion cycles.

8. The method of claim 1 , wherein the motion field is calculated over only a partial region of the object.

9. The method of claim 1 , wherein the detector data for creating the projection data sets used is collectable from one or a plurality of motion cycles.

10. The method of claim 1 , wherein the detector data for creating the projection data sets used is collectable from one or a plurality of motion cycles.

11. A non-transitory computer readable medium including program segments for, when executed on a computer device, causing the computer device to implement the method of claim 1 .

12. A method for generating a motion-compensated CT image data set for a partially and cyclically moving object, the method comprising:
detecting or transferring a projection data set from a computed tomography system covering a motion phase and a projection angle range, to enable the direct reconstruction of a CT image data set;
determining a motion field according to the method of claim 1 ;
reconstructing a final CT image data set which has a second image resolution using a motion-compensated reconstruction method based on a second reconstruction algorithm and the motion field; and
at least one of storing the final CT image data set and outputting the final CT image data set to an image reproduction system.

13. The method of claim 12 , wherein the first image resolution is relatively lower than the second image resolution.

14. The method of claim 13 , wherein the second reconstruction algorithm differs from the first reconstruction algorithm.

15. The method of claim 12 , wherein the second reconstruction algorithm differs from the first reconstruction algorithm.

16. The method of claim 15 , wherein the second reconstruction algorithm is an analytical reconstruction algorithm.

17. The method of claim 15 , wherein the second reconstruction algorithm is an iterative reconstruction algorithm.

18. The method of claim 15 , wherein the second reconstruction algorithm is a non-analytical reconstruction algorithm.

19. A non-transitory computer readable medium including program segments for, when executed on a computer device, causing the computer device to implement the method of claim 12 .

20. A computer system, comprising:
a memory unit, configured to store computer programs; and
a processor configured to carry out the stored computer programs, at least one of the computer programs being stored in the memory unit, the at least one computer program carrying out, operation of the computer system, at least the following:
detecting or transferring a projection data set from a computed tomography system covering a motion phase and a projection angle range, to enable direct reconstruction of a CT image data set;
iteratively determining a motion field by minimizing a cost function from CT image data sets which are motion-compensated and reconstructed making use of motion fields and which have a first image resolution, with a first analytical reconstruction algorithm, the cost function being based on raw data or transformed into the image space; and
at least one of storing and outputting the motion field,
wherein, as the raw data-based cost function, the equation min arg L=∥p in −A( , )f( , )∥ 2 is used.

21. A CT system comprising the computer system of claim 20 .

22. A computer system for image reconstruction, comprising:
a memory unit, configured to store computer programs; and
a processor configured to carry out the stored computer programs, at least one of the computer programs being stored in the memory unit, the at least one computer program carrying out, operation of the computer system, at least the following:
detecting or transferring a projection data set from a computed tomography system covering a motion phase and a projection angle range, to enable the direct reconstruction of a CT image data set;
determining a motion field by at least,
detecting or transferring a projection data set from a computed tomography system covering a motion phase and a projection angle range, to enable direct reconstruction of a CT image data set,
iteratively determining a motion field by minimizing a cost function from CT image data sets which are motion-compensated and reconstructed making use of motion fields and which have a first image resolution, with a first analytical reconstruction algorithm, the cost function being based on raw data or transformed into the image space, and
at least one of storing and outputting the motion field;

reconstructing a final CT image data set which has a second image resolution using a motion-compensated reconstruction method based on a second reconstruction algorithm and the motion field; and
at least one of storing the final CT image data set and outputting the final CT image data set to an image reproduction system,
wherein, as the raw data-based cost function, the equation min arg L=∥p in −A( , )f( , )∥ 2 is used.

23. A

Assignees

• Siemens Ag

Inventors

• Bruder Herbert
• Rohkohl Christopher

Classifications

• A61B6/5205

  involving processing of raw data to produce diagnostic data · CPC title

• A61B6/503

  for diagnosis of the heart · CPC title

• A61B6/5264

  due to motion · CPC title

• A61B6/032

  Transmission computed tomography [CT] · CPC title

• A61B6/507

  for determination of haemodynamic parameters, e.g. perfusion CT · CPC title