Method and system for propagation of myocardial infarction from delayed enhanced cardiac imaging to cine magnetic resonance imaging using hybrid image registration

• US 8,682,054 B2
• Filed: 11/15/2011
• Issued: 03/25/2014
• Est. Priority Date: 11/15/2010
• Status: Active Grant

• Current Assignee: Siemens Healthineers AG (Siemens AG)
• Original Assignee: Siemens AG
• Inventors: Xue, Hui, Liu, Yixun, Guetter, Christoph, Jolly, Marie-Pierre, Ghring, Jens
• Primary Examiner(s): Patel, Nirav G
• Application Number:
  US13/296,487
• Publication Number:
  US 20120121154A1
• Time in Patent Office: 861 Days
• Field of Search: None
• US Class Current: 382/131
• CPC Class Codes: G06T 2207/10088 Magnetic resonance imaging [MRI]
• CPC Class Codes: G06T 2207/30048 Heart; Cardiac
• CPC Class Codes: G06T 7/33 using feature-based methods

Abstract

A method and system for propagation of myocardial infarction from delayed enhanced magnetic resonance imaging (DE-MRI) to cine MRI is disclosed. A reference frame is selected in a cine MRI sequence. Deformation fields are calculated within the cine MRI sequence to register the frames of the cine MRI sequence to the reference frame. A DE-MRI image having an infarction region is registered to the reference frame of the cine MRI sequence. The DE-MRI image may be registered to the infarction region using a hybrid registration algorithm that unifies both intensity and feature points into a single cost function. Infarction information in the DE-MRI image is then propagated cardiac phases of the frames in the cine MRI sequence based on the registration of the DE-MRI image to the reference frame and the plurality of deformation fields calculated within the cine MRI sequence.

Claims

1. A method for propagating at least a portion of a static image over a cardiac cycle, comprising:

selecting a reference frame of a cine image sequence comprising a plurality of frames;

calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame;

registering the static image to the reference frame of the cine image sequence; and

propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence based on a deformation field resulting from registering the static image to the reference frame and the plurality of deformation fields calculated within the cine image sequence.

2. The method of claim 1, wherein the step of selecting a reference frame of a cine image sequence comprises:

selecting one of the plurality of frames in the cine image sequence that is most similar to the static image.

3. The method of claim 2, wherein the step of selecting one of the plurality of frames in the cine image sequence that is most similar to the static image comprises:

selecting one of the plurality of frames in the cine image sequence corresponding to a cardiac phase closest to a cardiac phase of the static image based on trigger times associated with the frames of the cine image sequence.

4. The method of claim 2, wherein the step of selecting one of the plurality of frames in the cine image sequence that is most similar to the static image comprises:

calculating a cross-correlation value between each frame of the cine image sequence and the static image; and

selecting the frame with the largest cross-correlation value as the reference frame.

5. The method of claim 1, wherein the step of calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame comprises:

registering each frame of the cine image sequence other than the reference frame to the reference frame using variational non-rigid registration with an inverse consistent constraint.

6. The method of claim 1, wherein the step of calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame comprises:

calculating a plurality of inverse consistent deformation fields to register the frames of the cine image sequence other than the reference frame to the reference frame.

7. The method of claim 1, wherein the step of registering the static image to the reference frame of the cine image sequence comprises:

registering the static image to the reference frame using a hybrid registration algorithm that unifies intensity-based and point-based similarity into a single cost function.

8. The method of claim 7, wherein the cost function includes a feature point matching term and an intensity matching term.

9. The method of claim 1, wherein the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

propagating the portion of the static image to a phase of the reference frame using the deformation field resulting from registering the static image to the reference frame; and

propagating the portion of the static image from the phase of the reference frame to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

10. The method of claim 1, wherein the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

deforming the static image using the deformation field resulting from registering the static image to the reference frame; and

propagating the deformed static image to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

11. The method of claim 1, wherein the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

deforming a delineated contour corresponding to a region of interest in the static image using the deformation field resulting from registering the static image to the reference frame; and

propagating the deformed contour to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

12. The method of claim 1, wherein the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

generating a region-of-interest in the reference frame by superimposing the portion of the static image onto the reference frame using the deformation field resulting from registering the static image to the reference frame; and

propagating the region-of-interest from the reference frame to each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

13. An apparatus for propagating at least a portion of static image over a cardiac cycle, comprising:

means for selecting a reference frame of a cine image sequence comprising a plurality of frames;

means for calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame;

means for registering a static image to the reference frame of the cine image sequence; and

means for propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence based on a deformation field resulting from registering the static image to the reference frame and the plurality of deformation fields calculated within the cine image sequence.

14. The apparatus of claim 13, wherein the means for selecting a reference frame of a cine image sequence comprises:

means for selecting one of the plurality of frames in the cine image sequence that is most similar to the static image.

15. The apparatus of claim 13, wherein the means for calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame comprises:

means for registering each frame of the cine image sequence other than the reference frame to the reference frame using variational non-rigid registration with an inverse consistent constraint.

16. The apparatus of claim 13, wherein the means for calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame comprises:

means for calculating a plurality of inverse consistent deformation fields to register the frames of the cine image sequence other than the reference frame to the reference frame.

17. The apparatus of claim 13, wherein the means for registering a static image to the reference frame of the cine image sequence comprises:

means for registering the static image to the reference frame using a hybrid registration algorithm that unifies intensity-based and point-based similarity into a single cost function.

18. The apparatus of claim 17, wherein the cost function includes a feature point matching term and an intensity matching term.

19. The apparatus of claim 13, wherein the means for propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

means for propagating the portion of the static image to a phase of the reference frame using the deformation field resulting from registering the static image to the reference frame; and

means for propagating the portion of the static image from the phase of the reference frame to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

20. The apparatus of claim 13, wherein the means for propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

means for deforming the static image using the deformation field resulting from registering the static image to the reference frame; and

means for propagating the deformed static image to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

21. The apparatus of claim 13, wherein the means for propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

means for deforming a delineated contour corresponding to a region of interest in the static image using the deformation field resulting from registering the static image to the reference frame; and

means for propagating the deformed contour to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

22. The apparatus of claim 13, wherein the means for propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprises:

means for generating a region-of-interest in the reference frame by superimposing the portion of the static image onto the reference frame using the deformation field resulting from registering the static image to the reference frame; and

means for propagating the region-of-interest from the reference frame to each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

23. A non-transitory computer readable medium encoded with computer executable instructions for propagating at least a portion of a static image over a cardiac cycle, the computer executable instructions defining steps comprising:

selecting a reference frame of a cine image sequence comprising a plurality of frames;

calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame;

registering a static image to the reference frame of the cine image sequence; and

propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence based on a deformation field resulting from registering the static image to the reference frame and the plurality of deformation fields calculated within the cine image sequence.

24. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of selecting a reference frame of a cine image sequence comprise computer executable instructions defining the step of:

selecting one of the plurality of frames in the cine image sequence that is most similar to the static image.

25. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame comprise computer executable instructions defining the step of:

registering each frame of the cine image sequence other than the reference frame to the reference frame using variational non-rigid registration with an inverse consistent constraint.

26. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of calculating a plurality of deformation fields within the cine image sequence to register the frames of the cine image sequence to the reference frame comprise computer executable instructions defining the step of:

calculating a plurality of inverse consistent deformation fields to register the frames of the cine image sequence other than the reference frame to the reference frame.

27. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of registering a static image to the reference frame of the cine MRI sequence comprise computer executable instructions defining the step of:

registering the static image to the reference frame using a hybrid registration algorithm that unifies intensity-based and point-based similarity into a single cost function.

28. The non-transitory computer readable medium of claim 27, wherein the cost function includes a feature point matching term and an intensity matching term.

29. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprise computer executable instructions defining the steps of:

propagating the portion of the static image to a phase of the reference frame using the deformation field resulting from registering the static image to the reference frame; and

propagating the portion of the static image from the phase of the reference frame to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

30. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprise computer executable instructions defining the steps of:

deforming the static image using the deformation field resulting from registering the static image to the reference frame; and

propagating the deformed static image to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

31. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprise computer executable instructions defining the steps of:

deforming a delineated contour corresponding to a region of interest in the static image using the deformation field resulting from registering the static image to the reference frame; and

propagating the deformed contour to a phase associated with each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.

32. The non-transitory computer readable medium of claim 23, wherein the computer executable instructions defining the step of propagating at least a portion of the static image to a plurality of cardiac phases corresponding to the plurality of frames in the cine image sequence comprise computer executable instructions defining the steps of:

generating a region-of-interest in the reference frame by superimposing the portion of the static image onto the reference frame using the deformation field resulting from registering the static image to the reference frame; and

propagating the region-of-interest from the reference frame to each frame of the cine image sequence other than the reference frame using an inverse of a respective one of the plurality of deformation fields calculated within the cine image sequence.