DE102010026675A1 - Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt - Google Patents
Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt Download PDFInfo
- Publication number
- DE102010026675A1 DE102010026675A1 DE102010026675A DE102010026675A DE102010026675A1 DE 102010026675 A1 DE102010026675 A1 DE 102010026675A1 DE 102010026675 A DE102010026675 A DE 102010026675A DE 102010026675 A DE102010026675 A DE 102010026675A DE 102010026675 A1 DE102010026675 A1 DE 102010026675A1
- Authority
- DE
- Germany
- Prior art keywords
- images
- phase
- series
- registration
- determined
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000033001 locomotion Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000003384 imaging method Methods 0.000 title claims abstract description 10
- 238000004590 computer program Methods 0.000 title claims abstract description 6
- 238000006073 displacement reaction Methods 0.000 claims description 13
- 238000007408 cone-beam computed tomography Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 5
- 230000000241 respiratory effect Effects 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000004072 lung Anatomy 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- BUHVIAUBTBOHAG-FOYDDCNASA-N (2r,3r,4s,5r)-2-[6-[[2-(3,5-dimethoxyphenyl)-2-(2-methylphenyl)ethyl]amino]purin-9-yl]-5-(hydroxymethyl)oxolane-3,4-diol Chemical compound COC1=CC(OC)=CC(C(CNC=2C=3N=CN(C=3N=CN=2)[C@H]2[C@@H]([C@H](O)[C@@H](CO)O2)O)C=2C(=CC=CC=2)C)=C1 BUHVIAUBTBOHAG-FOYDDCNASA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000002594 fluoroscopy Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000036391 respiratory frequency Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
- A61B6/487—Diagnostic techniques involving generating temporal series of image data involving fluoroscopy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/262—Analysis of motion using transform domain methods, e.g. Fourier domain methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/486—Diagnostic techniques involving generating temporal series of image data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Multimedia (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Die Erfindung betrifft ein Verfahren zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie von einem periodisch bewegten Objekt, welche Bilderserie eine Vielzahl von sequentiellen angefertigten Abbildungen des Objekts umfasst, umfassend folgende Schritte: – Registrierung von verschiedenen Abbildungen der Bilderserie, sodass eine Deformation als Maß für die Veränderung, die zwischen den zueinander registrierten Abbildungen aufgetreten ist, ermittelt wird, – Bestimmung der Phase der Objektbewegung für wenigstens eine der Abbildungen der Bilderserie unter Verwendung der durch die Registrierung ermittelten Deformation. Weiterhin betrifft die Erfindung eine Vorrichtung mit einer Rechnereinheit, die zur Dürchführung eines entsprechenden Verfahrens ausgebildet ist, sowie ein entsprechendes Computerprogrammprodukt zur Implementierung eines derartigen Verfahrens und ein bildgebendes Gerät mit einer derartigen Rechnereinheit.
Description
- Die Erfindung betrifft ein Verfahren zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, in der ein sich bewegendes Objekt abgebildet ist. Weiterhin betrifft die Erfindung eine Vorrichtung zur Auswertung einer Bilderserie zur Bestimmung einer Phase einer Objektbewegung eines in der Bilderserie abgebildeten sich bewegenden Objekts. Die Erfindung betrifft ebenso ein entsprechendes bildgebendes Gerät und Computerprogrammprodukt. Derartige Verfahren bzw. Vorrichtungen finden Anwendung in der medizinischen Bildgebung, bei der es für die Rekonstruktion von dreidimensionalen Abbildungen darauf ankommt, in welcher Bewegungsphase sich ein sich bewegenden Objekt, wie beispielsweise eine Lunge, bei der Aufzeichnung der Rohbilddaten befunden hat.
- Ein seit Jahren bekanntes Verfahren der medizinischen Bildgebung ist die so genannte Cone-Beam-Computertomographie. Bei dieser Art von Bildgebung werden zweidimensionale Projektionen des Objekts mit Röntgenstrahlen angefertigt. Die Projektionen werden von verschiedenen Winkeln aus einem Winkelbereich von über 180° angefertigt. Aus den zweidimensionalen Projektionen kann dann eine dreidimensionale Abbildung des Objekts rekonstruiert werden.
- Es benötigt eine gewisse Zeit, bis alle zweidimensionalen Projektionen aus unterschiedlichen Richtungen angefertigt worden sind. Falls sich das abzubildende Objekt in dieser Zeit bewegt, wie beispielsweise eine Lunge, die eine quasiperiodische (Atem-)Bewegung ausführt, sind die einzelnen Projektionsbilddaten zu unterschiedlichen Phasen der Objektbewegung angefertigt. Wenn die Projektionen dann zur Rekonstruktion einer dreidimensionalen Abbildung verwendet werden, kann es aufgrund der inkonsistenten Projektionen zu Artefakten kommen, welche die Auswertung und/oder Befundung der dreidimensionalen Abbildung erschweren.
- Eine Möglichkeit, derartige Artefakte zu vermeiden, ist es, zur Rekonstruktion lediglich Projektionen zu verwenden, die einander bezüglich der Bewegungsphase des Objekts entsprechen. Beispielsweise können zur Rekonstruktion nur Projektionsbilder ein und derselben Atemlage verwendet werden.
- Eine Möglichkeit, die Atemlage festzustellen, ist die Verwendung eines Atemgurts. Diese Möglichkeit ist in der Schrift: Lars Dietrich et al, "Linac-integrated 4D cone beam CT: first experimental results", 2006, Phys. Med. Biol. 51, 2939 beschrieben.
- Die
US 7,349,564 offenbart ein Verfahren, das zur Feststellung der Atemlage anatomische Merkmale in den Bilddaten verwendet, beispielsweise die Position eines Zwerchfells. - Es ist die Aufgabe der Erfindung, ein Verfahren und eine Vorrichtung anzugeben, mit denen sich die Phase der Objektbewegung bei einer Bilderserie des sich bewegenden Objekts auf robuste Weise und/oder ohne externe Hilfsmittel festgestellt werden kann. Weiterhin ist es die Aufgabe der Erfindung ein entsprechendes bildgebendes Gerät und Computerprogrammprodukt anzugeben.
- Die Aufgabe der Erfindung wird durch die Merkmale der unabhängigen Ansprüche gelöst. Vorteilhafte Weiterbildungen der Erfindung finden sich in den Merkmalen der abhängigen Ansprüche.
- Das erfindungsgemäße Verfahren zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie von einem periodisch bewegten Objekt, welche Bilderserie eine Vielzahl von sequentiell angefertigten Abbildungen des Objekts umfasst, umfasst folgende Schritte:
- – Registrierung von verschiedenen Abbildungen der Bilderserie zueinander, sodass eine Deformation, insbesondere eine Verschiebung, als Maß für die Veränderung, die zwischen den zueinander registrierten Abbildungen aufgetreten ist, ermittelt wird,
- – Bestmmung der Phase der Objektbewegung für wenigstens eine der Abbildungen der Bilderserie unter Verwendung der durch die Registrierung ermittelten Deformation.
- Die Erfindung beruht auf der Erkenntnis, dass es vorteilhaft ist, wenn lediglich die Abbildungen selbst herangezogen werden, um die Phase der Objektbewegung für einzelne Abbildungen der Bilderserie zu bestimmen. Es wurde aber gleichzeitig erkannt, dass es mitunter problematisch sein kann, wenn man für diese Bestimmung ein anatomisches Merkmal in den Bilddaten isoliert betrachtet, wie beispielsweise die Position des Zwerchfells.
- Die Erfindung bietet hierfür eine einfache aber wirkungsvolle Lösung an. Grundlage für die Bestimmung der Phase der Objektbewegung bildet nun eine Registrierung, die zwischen verschiedenen Abbildungen der Bilderserie durchgeführt wird.
- Mit einer derartigen Registrierung lässt sich eine Deformation bzw. eine Verschiebung bestimmen, die ein Maß dafür darstellt, wie sich das Objekt zwischen den zwei zueinander registrierten Abbildungen verändert hat. Beispielsweise kann eine Verschiebung durch die Registrierung bestimmt werden. Die durch die Registrierung ermittelte Deformation kann dann dafür verwendet werden, um die Phase der Objektbewegung für eine und/oder für mehrere Abbildungen der Bilderserie daraus abzuleiten.
- Die Erfindung ist also vorteilhaft; eine Registrierung zwischen den Abbildungen lässt sich nämlich stets durchführen. Das Verfahren wird dadurch unempfindlich dagegen, dass in einzelnen Abbildungen nicht alle relevanten anatomischen Merkmale, auf denen eine merkmalsbasierte Phasenbestimmung beruht, sichtbar oder detektierbar sein müssen. Die Bestimmung der Phase der Objektbewegung ist damit auch möglich, wenn das anatomische Merkmal nicht in allen Projektionen sichtbar oder detektierbar ist.
- Die Registrierung kann beispielsweise eine lineare Registrierung sein, die auf besonders einfache Weise die Deformation zwischen zwei Abbildungen, insbesondere zwischen zwei aufeinander folgenden Abbildungen der Bilderserie bestimmt. Als Maß für die Stärke der Deformation kann die Verschiebung bestimmt werden, die bei einer Registrierung ermittelt wird, z. B. anhand des Verschiebevektor oder dessen Länge.
- Die Registrierung muss dabei nicht zwangsläufig auf den gesamten Bildinhalt der Abbildungen angewendet werden, beispielsweise kann es auch genügen, einen Teilbereich der Abbildungen zu definieren, auf den die Registrierung angewendet werden soll. Dieser Teilbereich kann beispielsweise durch eine Benutzereingabe ausgewählt werden.
- Die Registrierung kann bei je zwei aufeinander folgenden Abbildungen angewendet werden, d. h. zwei aufeinander folgende Abbildungen werden zueinander registriert, um jeweils eine Deformation als Maß für die Veränderung, die sukzessive von Abbildung zu Abbildung auftritt, zu ermitteln.
- Auf diese Weise kann für eine oder mehrere der Abbildungen der Bilderserie die Phasenlage ermittelt werden, beispielsweise die Atemlage oder der Zeitpunkt im Herzzyklus, an dem die Abbildung aufgezeichnet wurde.
- Die Abbildungen der Bilderserie kann eine Serie sequentieller Durchleuchtungsabbildungen bilden, die von einem zu untersuchenden Objekt angefertigt z. B. mittels Röntgenstrahlen angefertigt worden sind. Insbesondere können die Durchleuchtungsabbildungen Projektionen einer Cone-Beam-Computertomographie sein.
- Bei einer Cone-Beam-Computertomographie wirkt sich das Verfahren besonders vorteilhaft aus, da hierbei eine fehlerhafte Bestimmung der Phase der Objektbewegung bei der Rekonstruktion zu einer Artefakt-behafteten Rekonstruktion führt. Zudem bewirkt die Vielzahl von unterschiedlichen Richtungen, aus denen die Projektionen angefertigt werden, dass sich die Lage eines anatomischen Merkmals in den Projektionen üblicherweise deutlich ändert oder gar aus einzelnen Projektionen ”rutschen” kann. Dies kann eine merkmalsbasierte Phasenbestimmung erschweren oder fehlerbehaftet machen.
- Bei der Cone-Beam-Computertomographie kann insbesondere eine exzentrische Projektionsgeometrie verwendet werden. Insbesondere dann wirkt sich das Verfahren sehr vorteilhaft aus. Exzentrische Projektionsgeometrie bedeutet, dass das Rotationszentrum der CBCT nicht symmetrisch in dem zur Durchleuchtung ausgesendete Röntgen-Kegelstrahl liegt; der Röntgen-Kegelstrahl kann schräg seitlich am Zentrum vorbei gelenkt werden, um dann auf einen seitlich versetzten, exzentrischen Detektor zu treffen.
- Bedingt durch die exzentrische Projektionsgeometrie kann es besonders leicht auftreten, dass ein anatomisches Merkmal seine Position abhängig vom Projektionswinkel ändert oder gar ausgeblendet wird. Daher vermeidet das Registrierungsbasierte Verfahren bei einer exzentrischen Projektionsgeometrie besonders effizient die Nachteile einer merkmalsbasierten Rekonstruktion. Einfache Schwellwertverfahren könnten hingegen vermehrt zu falschen Ergebnissen führen.
- Ein weiterer Vorteil ist es, dass die Phase der Objektbewegung bereits aus wenigen aufeinander folgenden Abbildungen bestimmt werden kann. Es muss keine vollständige Periode der Objektbewegung abgewartet werden, um die Phasenlage in der Objektbewegung korrekt bestimmen zu können.
- Aus den Deformationen bzw. Verschiebungen, die zwischen Abbildungen der Bilderserie aufgetreten sind, kann ein Signal gebildet werden, das darauffolgend zur Bestimmung der Phase verwendet wird.
- Das Signal gibt also an, wie sehr sich die Abbildungen der Bilderserie im Lauf der Bilderserie verändern, also deformieren bzw. zueinander verschieben. Beispielsweise kann das Signal durch eine Addition der Verschiebungen erzeugt werden, die im Lauf der Bilderserie zwischen aufeinander folgenden Abbildungen aufgetreten sind.
- Die Bestimmung der Phase der Objektbewegung kann dann auf einfache Weise mithilfe des Signals erfolgen, beispielsweise indem eine Hilbert-Transformation auf das Signal angewendet wird. Die Hilbert-Transformation kann als eine Phasenverschiebung um 90° im Frequenzraum angesehen werden.
- Die Phase der Objektbewegung kann mithilfe der oben genannten Verfahren für einzelne, mehrere oder insbesondere für alle Abbildungen der Bilderserie jeweils bestimmt werden. Auf diese Weise lassen sich sämtliche Abbildungen der Bilderserie klassifizieren, beispielsweise können die Abbildungen nach der ihnen zugeordneten Phase der Objektbewegung sortiert werden.
- Für eine nachfolgende Rekonstruktion können dann lediglich Abbildungen, deren Phase der Objektbewegung in einem bestimmten Intervall liegt, verwendet werden. Das auf diese Weise rekonstruierte Bild weist weniger bewegungsinduzierte Artefakte auf.
- Die erfindungsgemäße Vorrichtung zur Auswertung einer Bilderserie von einem periodisch bewegten Objekt zur Bestimmung einer Phase einer Objektbewegung, welche Bilderserie eine Vielzahl von sequentiellen angefertigten Abbildungen des Objekts umfasst, umfasst eine Rechnereinheit, welche ausgebildet ist zur
- – Registrierung von verschiedenen Abbildungen der Bilderserie, sodass eine Deformation als Maß für die Veränderung, die zwischen den zueinander registrierten Abbildungen aufgetreten ist, ermittelbar ist, und zur
- – Bestimmung der Phase der Objektbewegung für wenigstens eine der Abbildungen der Bilderserie unter Verwendung der durch die Registrierung ermittelten Deformation.
- Die oben beschriebenen Verfahren können in der Rechnereinheit implementiert sein.
- Ein bildgebendes Gerät zur Abbildung eines Objekts weist eine Röntgenstrahlenquelle und einen Röntgendetektor zur Anfertigung der Bilderserie des Objekts auf. Es kann mit einer derartigen Vorrichtung ausgestattet sein.
- Das erfindungsgemäße Computerprogrammprodukt hat maschinenlesbare Anweisungen, die ein Verfahren nach einem der Ansprüche 1 bis 11 auf einer Rechnereinheit implementieren, wenn sie auf einer Rechnereinheit ausgeführt werden.
- Die vorangehende und die folgende Beschreibung der einzelnen Merkmale, deren Vorteile und deren Wirkungen bezieht sich sowohl auf die Vorrichtungskategorie als auch auf die Verfahrenskategorie, ohne dass dies im Einzelnen in jedem Fall explizit erwähnt ist; die dabei offenbarten Einzelmerkmale können auch in anderen als den gezeigten Kombinationen erfindungswesentlich sein.
- Ausführungsformen der Erfindung mit vorteilhaften Weiterbildungen gemäß den Merkmalen der abhängigen Ansprüche werden anhand der folgenden Zeichnung näher erläutert, ohne jedoch darauf beschränkt zu sein. Es zeigen:
-
1 eine schematische Darstellung eines Cone-Beam-Computertomographen mit exzentrischer Projektionsgeometrie, -
2 ein Ablaufdiagramm einer Ausführungsform des erfindungsgemäßen Verfahrens, -
3 ein Diagramm, dass die zwischen jeweils zwei Abbildungen aufgetretene Verschiebung über den Lauf der Bilderserie zeigt, -
4 ein Diagramm, dass das Signal zeigt, das aus Addition der Verschiebungen erzeugt wird, -
5 eine lokale Phase, die mittels Hilbert-Transformation aus dem Signal für die einzelnen Abbildungen der Bilderserie ermittelt wird, -
6 eine Darstellung einer graphischen Benutzerschnittstelle, mit der Einstellungen für das Verfahren festgelegt werden können. -
1 zeigt eine schematische Darstellung einer Cone-Beam-Computertomographie. Bei dieser Art von Abbildung werden Projektionsabbildungen11 eines zu untersuchenden Objekts13 angefertigt. - Eine Röntgenquelle
15 und ein Röntgendetektor17 rotieren um ein gemeinsames Rotationszentrum19 . Von der Röntgenquelle15 wird ein Röntgen-Kegelstrahl21 auf den Röntgendetektor17 gerichtet. Der Röntgendetektor17 befindet sich an exzentrischer Position. Dies bedeutet, dass das Rotationszentrum19 nicht zentral im Röntgen-Kegelstrahl21 liegt. - Durch Rotation der Röntgenquelle
15 und des Röntgendetektors17 werden sukzessive eine Vielzahl von Projektionsabbildungen11 erzeugt, die gemeinsam eine Bilderserie des zu untersuchenden Objekts13 bilden. - Aufgrund der Bewegung des Objekts
13 , angedeutet durch den Pfeil, kann es vorkommen, dass sich das Objekt13 bei einer Projektionsabbildung in einer anderen Bewegungsphase befindet als bei einer anderen Projektionsabbildung. Für die Rekonstruktion ist es vorteilhaft, nur Projektionsabbildungen zu verwenden, bei denen sich das Objekt13 in ähnlicher Bewegungsphase befunden hat. - Anhand von
2 bis5 wird erläutert, wie mithilfe einer Registrierung bestimmt werden kann, in welcher Bewegungsphase sich das Objekt bei Anfertigung einer Projektionsabbildung befunden hat. -
2 zeigt ein Ablaufdiagramm für eine Ausführungsform des erfindungsgemäßen Verfahrens. Zunächst werden die Projektionsabbildungen aufgezeichnet (Schritt31 ). - Anschließend werden je zwei aufeinander folgende Projektionsabbildungen zueinander registriert, beispielsweise durch eine Registrierung, um festzustellen, wie stark sich das zu untersuchende Objekt zwischen einer Projektionsabbildung und der darauf folgenden Projektionsabbildung verschoben hat (Schritt
33 ). Beispielsweise kann dies durch eine Optimierung der Summe der quadratischen Abweichungen (engl: ”Summed Squared Difference”) der Bild-Grauwerte bezüglich einer Verschiebung entlang einer vorgegebenen Achse, z. B. entlang einer Körperachse wie einer Körperlängsachse, erfolgen. Denkbar wäre auch eine lineare Registrierung. - Zwischen je zwei aufeinander folgenden Projektionsabbildungen kann die aufgetretene Verschiebung ermittelt werden (Schritt
35 ). - In einem nächsten Schritt werden die jeweils ermittelten Verschiebungen sukzessive addiert (Schritt
37 ). Hierdurch lässt sich ein Signal erzeugen, das angibt, wie sich das zu untersuchende Objekt im Laufe der Bilderserie insgesamt in Summe verschiebt. - Auf dieses Signal kann in einem Vorverarbeitungsschritt ein Frequenzfilter angewendet werden. Anhand von Glättungs- und/oder Kantenfilter oder einer Kombination der beiden Filter können Signalfrequenzen unterdrückt werden, die zu sehr von einer normalen Atemfrequenz abweichen.
- Dieses vorverarbeitete Signal wird anschließend einer diskreten Hilbert-Transformation unterworfen (Schritt
39 ). - Das Signal und die Hilbert-Transformierte des Signals können dann dazu verwendet werden, jedem Projektionsbild eine lokale Phase zuzuordnen (Schritt
41 ). - Wenn s(p) das Signal darstellt (p für Projektionsabbildungsnummer), kann ein analytisches Signal f(p) durch f(p) = s(p) + i Hs(p) erhalten werden, wobei Hs(p) die Hilbert-Transformierte von s(p) darstellt. Die komplexe Funktion f(p) lässt sich auch schreiben als f(p) = F(p)·exp(i ph(p))·ph(p) stellt dabei die lokale Phase dar, die dann den einzelnen Projektionen zugeordnet werden kann.
- Anhand der lokalen Phase können die Projektionsbilder dann klassifiziert werden (Schritt
43 ). - Eine nachfolgende Bildrekonstruktion kann mithilfe der klassifizierten Projektionsabbildungen durchgeführt werden (Schritt
45 ). -
3 zeigt ein Diagramm, das die zwischen zwei Projektionsabbildungen aufgetretene Verschiebung V, aufgetragen gegenüber der Projektionsabbildungsnummer p, zeigt. -
4 zeigt ein Diagramm, bei dem die addierten Verschiebungen V', aufgetragen gegenüber der Projektionsabbildungsnummer p, dargestellt werden. Einer Projektionsabbildungsnummer wird die Summe aller bis dahin aufgetretenen Verschiebungen zugeordnet. -
5 zeigt die lokale Phase ph, die für jede Projektionsabbildung mithilfe der Hilbert Transformation bestimmt werden kann, aufgetragen gegenüber der Projektionsabbildungsnummer p. Diese lokale Phase kann dann für eine robuste Zuordnung einer Projektionsabbildung zu einer Phase der Objektbewegung, beispielsweise zu einer Atemlage, verwendet werden. -
6 zeigt eine Darstellung einer graphischen Benutzerschnittstelle, mit der Einstellungen für das Verfahren festgelegt werden können. - Im linken Bilddrittel
61 werden alle Projektionsabbildungen dargestellt. Ein Anwender kann zudem den Bereich (ROI für ”region of interest”) festlegen, der dem Registrierungsalgorithmus zu Grunde gelegt werden soll. - Im rechten Bilddrittel
63 lassen sich die Anzahl der Phasenlage eines Bewegungszyklus festlegen, in die die Projektionsabbildungen sortiert bzw. klassifiziert werden sollen. Es wird angezeigt, wie viele Projektionsabbildungen für jede der Phasenlagen vorliegen. Ein Anwender kann zudem die Phasenlage auswählen. - Im mittleren Bilddrittel
65 lassen sich dann die Projektionsabbildungen der gewählten Phasenlage anzeigen. Eine Rekonstruktion kann mit den klassifizierten Projektionsabbildungen durchgeführt werden. - Bezugszeichenliste
-
- 11
- Projektionsabbildung
- 13
- Objekt
- 15
- Röntgenquelle
- 17
- Röntgendetektor
- 19
- Rotationszentrum
- 21
- Röntgen-Kegelstrahl
- 31
- Schritt
31 - 33
- Schritt
33 - 35
- Schritt
35 - 37
- Schritt
37 - 39
- Schritt
39 - 41
- Schritt
41 - 43
- Schritt
43 - 45
- Schritt
45 - 61
- linkes Bilddrittel
- 63
- rechtes Bilddrittel
- 65
- mittleres Bilddrittel
- ZITATE ENTHALTEN IN DER BESCHREIBUNG
- Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
- Zitierte Patentliteratur
-
- US 7349564 [0006]
- Zitierte Nicht-Patentliteratur
-
- Lars Dietrich et al, ”Linac-integrated 4D cone beam CT: first experimental results”, 2006, Phys. Med. Biol. 51, 2939 [0005]
Claims (15)
- Verfahren zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie von einem periodisch bewegten Objekt, welche Bilderserie eine Vielzahl von sequentiellen angefertigten Abbildungen des Objekts umfasst, umfassend folgende Schritte: – Registrierung von verschiedenen Abbildungen der Bilderserie, sodass eine Deformation als Maß für die Veränderung, die zwischen den zueinander registrierten Abbildungen aufgetreten ist, ermittelt wird, – Bestimmung der Phase der Objektbewegung für wenigstens eine der Abbildungen der Bilderserie unter Verwendung der durch die Registrierung ermittelten Deformation.
- Verfahren nach Anspruch 1, wobei bei der Registrierung eine Verschiebung als Maß für die Veränderung, die zwischen den zueinander registrierten Abbildungen aufgetreten ist, ermittelt wird.
- Verfahren nach Anspruch 1 oder 2, wobei die Bilderserie eine Serie sequentieller Durchleuchtungsabbildungen umfasst, die von einem zu untersuchenden Objekt angefertigt worden sind.
- Verfahren nach Anspruch 3, wobei die Durchleuchtungsabbildungen insbesondere Projektionen einer Cone-Beam-Computertomographie sind, welche insbesondere mit einer exzentrischen Projektionsgeometrie angefertigt worden sind.
- Verfahren nach einem der Ansprüche 1 bis 4, wobei bei der Registrierung je zwei aufeinander folgende Abbildungen zueinander registriert werden, um jeweils eine Deformation als Maß für die Veränderung, die sukzessive von Abbildung zu Abbildung auftritt, zu ermitteln.
- Verfahren nach einem der Ansprüche 1 bis 5, wobei bei der Bestimmung der Phase aus den Verschiebungen, die zwischen verschiedenen Abbildungen aufgetreten sind, ein Signal gebildet wird, das zur Bestimmung der Phase verwendet wird.
- Verfahren nach Anspruch 6, wobei das Signal durch Addition aufeinanderfolgender Verschiebungen erzeugt wird.
- Verfahren nach Anspruch 6 oder 7, wobei bei der Bestimmung der Phase unter Verwendung einer Hilbert-Transformation, die auf das Signal angewendet wird, bestimmt wird.
- Verfahren nach einem der Ansprüche 1 bis 8, wobei bei der Bestimmung der Phase für eine Vielzahl von Abbildungen der Bilderserie, insbesondere für alle Abbildungen der Bilderserie, jeweils die Phase der Objektbewegung bestimmt wird.
- Verfahren nach Anspruch 9, wobei die Abbildungen nach der ihnen zugeordneten Phase der Objektbewegung klassifiziert werden.
- Verfahren nach Anspruch 10, wobei eine Bildrekonstruktion lediglich Abbildungen, deren Phase der Objektbewegung in einem bestimmten Intervall liegt, verwendet werden.
- Vorrichtung zur Auswertung einer Bilderserie von einem periodisch bewegten Objekt zur Bestimmung einer Phase einer Objektbewegung, welche Bilderserie eine Vielzahl von sequentiellen angefertigten Abbildungen des Objekts umfasst, umfassend eine Rechnereinheit, welche ausgebildet ist zur Registrierung von verschiedenen Abbildungen der Bilderserie, sodass eine Deformation als Maß für die Veränderung, die zwischen den zueinander registrierten Abbildungen aufgetreten ist, ermittelbar ist, und Bestimmung der Phase der Objektbewegung für wenigstens eine der Abbildungen der Bilderserie unter Verwendung der durch die Registrierung ermittelten Deformation.
- Vorrichtung, wobei die Rechnereinheit ausgebildet ist zur Durchführung eines Verfahrens nach einem der Ansprüche 1 bis 11.
- Bildgebendes Gerät zur Abbildung eines Objekts, umfassend: eine Vorrichtung nach Anspruch 12 oder 13, und eine Röntgenstrahlenquelle und einen Röntgendetektor zur Anfertigung der Bilderserie des Objekts.
- Computerprogrammprodukt mit maschinenlesbaren Anweisungen, die ein Verfahren nach einem der Ansprüche 1 bis 11 auf einer Rechnereinheit implementieren, wenn sie auf einer Rechnereinheit ausgeführt werden.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102010026675.2A DE102010026675B4 (de) | 2010-07-09 | 2010-07-09 | Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt |
US13/179,002 US8515003B2 (en) | 2010-07-09 | 2011-07-08 | Determining a phase of an object movement in a series of images |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102010026675.2A DE102010026675B4 (de) | 2010-07-09 | 2010-07-09 | Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt |
Publications (2)
Publication Number | Publication Date |
---|---|
DE102010026675A1 true DE102010026675A1 (de) | 2012-01-12 |
DE102010026675B4 DE102010026675B4 (de) | 2019-07-25 |
Family
ID=45372462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE102010026675.2A Expired - Fee Related DE102010026675B4 (de) | 2010-07-09 | 2010-07-09 | Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt |
Country Status (2)
Country | Link |
---|---|
US (1) | US8515003B2 (de) |
DE (1) | DE102010026675B4 (de) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5226590B2 (ja) | 2009-04-02 | 2013-07-03 | キヤノン株式会社 | 画像解析装置、画像処理装置及び画像解析方法 |
JP5294956B2 (ja) * | 2009-04-08 | 2013-09-18 | キヤノン株式会社 | 画像処理装置及び画像処理装置の制御方法 |
US9526468B2 (en) * | 2014-09-09 | 2016-12-27 | General Electric Company | Multiple frame acquisition for exposure control in X-ray medical imagers |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004054443A1 (en) * | 2002-12-13 | 2004-07-01 | Koninklijke Philips Electronics N.V. | System and method for processing a series of image frames representing a cardiac cycle |
US20040136490A1 (en) * | 2002-07-23 | 2004-07-15 | Edic Peter Michael | Method and apparatus for correcting motion in image reconstruction |
US7349564B2 (en) | 2003-01-21 | 2008-03-25 | Elekta Ab (Publ) | Imaging internal structures |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6535570B2 (en) * | 1999-06-17 | 2003-03-18 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence Of Her Majesty's Canadian Government | Method for tracing organ motion and removing artifacts for computed tomography imaging systems |
US6421552B1 (en) * | 1999-12-27 | 2002-07-16 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for estimating cardiac motion using projection data |
US6901131B2 (en) * | 2001-12-28 | 2005-05-31 | General Electric Company | Methods and apparatus for computed tomography imaging |
US7062009B2 (en) * | 2002-09-12 | 2006-06-13 | Analogic Corporation | Helical interpolation for an asymmetric multi-slice scanner |
US7280630B2 (en) * | 2004-09-30 | 2007-10-09 | Wisconsin Alumni Research Foundation | Calculation of additional projection data from projection data acquired with a divergent beam |
WO2006058124A1 (en) * | 2004-11-24 | 2006-06-01 | Wisconsin Alumni Research Foundation | Fan-beam and cone-beam image reconstruction using filtered backprojection of differentiated projection data |
JP4777007B2 (ja) * | 2005-08-03 | 2011-09-21 | 東芝メディカルシステムズ株式会社 | X線コンピュータ断層撮影装置 |
US7940971B2 (en) * | 2006-07-24 | 2011-05-10 | Siemens Medical Solutions Usa, Inc. | System and method for coronary digital subtraction angiography |
JP2010069099A (ja) * | 2008-09-19 | 2010-04-02 | Toshiba Corp | 画像処理装置及びx線コンピュータ断層撮影装置 |
-
2010
- 2010-07-09 DE DE102010026675.2A patent/DE102010026675B4/de not_active Expired - Fee Related
-
2011
- 2011-07-08 US US13/179,002 patent/US8515003B2/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040136490A1 (en) * | 2002-07-23 | 2004-07-15 | Edic Peter Michael | Method and apparatus for correcting motion in image reconstruction |
WO2004054443A1 (en) * | 2002-12-13 | 2004-07-01 | Koninklijke Philips Electronics N.V. | System and method for processing a series of image frames representing a cardiac cycle |
US7349564B2 (en) | 2003-01-21 | 2008-03-25 | Elekta Ab (Publ) | Imaging internal structures |
Non-Patent Citations (2)
Title |
---|
Lars Dietrich et al, "Linac-integrated 4D cone beam CT: first experimental results", 2006, Phys. Med. Biol. 51, 2939 |
SUGIMOTO N. [u.a.]: Medical Processing in Collaboration with Medical Researchers, In: Informatics Research for Development of Knowledge Society Infrastructure, 2007, S. 77-86 * |
Also Published As
Publication number | Publication date |
---|---|
DE102010026675B4 (de) | 2019-07-25 |
US8515003B2 (en) | 2013-08-20 |
US20120008737A1 (en) | 2012-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102008028387B4 (de) | Tomographisches Bildrekonstruktionsverfahren zum Erzeugen eines Bildes von einem Untersuchungsobjekt und nach diesem Verfahren arbeitende bildgebende Einrichtung | |
DE102005049602B3 (de) | Verfahren und Vorrichtung zur Segmentierung zumindest einer Substanz in einem Röntgenbild | |
DE102009032059A1 (de) | Sinogrammbearbeitung für die Metallartefaktreduktion in der Computertomographie | |
DE102005049586A1 (de) | Verfahren zur Erzeugung von CT-Darstellungen in der Röntgen-Computertomographie | |
DE102009051384A1 (de) | Strahlaufhärtungskorrektur für CT-Perfusionsmessungen | |
DE102007061935A1 (de) | Verfahren zur Qualitätssteigerung von computertomographischen Aufnahmeserien durch Bildverarbeitung und CT-System mit Recheneinheit | |
DE102017200282B3 (de) | Verfahren zur Reduktion von Bildartefakten | |
DE102014200303B4 (de) | Verfahren zur Durchführung einer Positronenemissionstomographie in einer Hybridanlage und entsprechende Hybridanlage | |
DE102011083727B4 (de) | Verfahren zur Erzeugung eines rauschreduzierten CT-Bilddatensatzes, Rechensystem und CT-System | |
DE102010013360B4 (de) | Verfahren zur Rekonstruktion von Bilddaten eines zyklisch sich bewegenden Untersuchungsobjektes | |
DE102008048045A1 (de) | Verfahren zur Erzeugung von computertomographischen Bilddatensätzen eines Patienten in der Herz-CT bei einer Perfusionskontrolle unter Kontrastmittelapplikation | |
DE102012216652A1 (de) | Angiographisches Untersuchungsverfahren | |
DE102011076547A1 (de) | Verfahren zum Gewinnen eines 3D-Bilddatensatzes zu einem Bildobjekt | |
DE102010026675B4 (de) | Verfahren und Vorrichtung zur Bestimmung einer Phase einer Objektbewegung bei einer Bilderserie, bildgebendes Gerät und Computerprogrammprodukt | |
DE102008045449B4 (de) | Verfahren zur Erstellung von computertomographischen Aufnahmen eines Patienten mit metallischen Anteilen und Computersystem zur Durchführung des Verfahrens | |
DE102008038357B3 (de) | Verfahren zur Erzeugung von 2D-Schnittbildern aus 3D-Projektionsdaten, die mittels eines CT-Systems von einem metallische Anteile enthaltenden Untersuchungsobjekt erfasst wurden | |
DE102009007236A1 (de) | CT-Bildrekonstruktion eines sich bewegenden Untersuchungsobjektes | |
DE102017213479A1 (de) | Computertomographische Aufnahme mit verschiedenen Energieschwellensätzen | |
DE102008045633B4 (de) | Verfahren zur verbesserten Darstellung von Mehr-Energie-CT-Aufnahmen und zugehöriges Computersystem | |
DE102014217966A1 (de) | Adaptive Metallartefaktreduktion für die Computertomographie | |
DE102009010287B4 (de) | Verfahren und Bildrekonstruktionseinrichtung zur Rekonstruktion von Bilddaten | |
DE102012223745A1 (de) | CT-Bildrekonstruktion mit kantenerhaltender Filterung | |
EP4057229A1 (de) | Verarbeitungsvorrichtung zum erhalten einer lückenfüller-sinogramminformation, computertomograph, verfahren und computerprogramm | |
DE102011005554A1 (de) | Verfahren zur Reduzierung von Bewegungsartefakten bei der Dual-Energy-Computertomographie | |
DE102011075912A1 (de) | Verfahren zum Bereitstellen eines 3D-Bilddatensatzes mit unterdrückten Bildartefakten und Computertomograph |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R016 | Response to examination communication | ||
R079 | Amendment of ipc main class |
Free format text: PREVIOUS MAIN CLASS: A61B0019000000 Ipc: A61B0034000000 |
|
R081 | Change of applicant/patentee |
Owner name: SIEMENS HEALTHCARE GMBH, DE Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE |
|
R016 | Response to examination communication | ||
R018 | Grant decision by examination section/examining division | ||
R020 | Patent grant now final | ||
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |