DE102007026115A1 - Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers - Google Patents
Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers Download PDFInfo
- Publication number
- DE102007026115A1 DE102007026115A1 DE102007026115A DE102007026115A DE102007026115A1 DE 102007026115 A1 DE102007026115 A1 DE 102007026115A1 DE 102007026115 A DE102007026115 A DE 102007026115A DE 102007026115 A DE102007026115 A DE 102007026115A DE 102007026115 A1 DE102007026115 A1 DE 102007026115A1
- Authority
- DE
- Germany
- Prior art keywords
- projection
- detector
- data
- ray
- projections
- 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
- 238000000034 method Methods 0.000 title claims description 17
- 238000001914 filtration Methods 0.000 claims abstract description 19
- 238000013507 mapping Methods 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 description 2
- 238000000926 separation method Methods 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
- 210000000746 body region Anatomy 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000002244 precipitate Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/40—Imaging
- G01N2223/419—Imaging computed tomograph
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/432—Truncation
Landscapes
- Theoretical Computer Science (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Pulmonology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Immunology (AREA)
- Radiology & Medical Imaging (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
In der DE 102006041033.5 ist beschrieben, dass zur Erzeugung einer 3-D-Rekonstruktion eines besonders großen Körpers, der durch eine einzelne Projektion nicht einnehmbar ist, bei jeder Stellung eines Röntgen-C-Bogen zumindest zwei Projektionen aufgenommen werden, die gemeinsam den Körper erfassen. Aus den beiden Projektionen werden Daten auf einen virtuellen Detektor projiziert, und diese werden zur gefilterten Rückprojektion verwendet. Ausgegangen wird hierbei davon, dass die reale Quelle ortsunveränderlich bleibt und sich nur der Detektor bewegt. Für den Fall, dass reale Quellen Q1 und Q2 bei den beiden Projektionen nicht zusammenfallen, wird ein virtueller Detektor D1' bzw. D2' lediglich zum Durchführen einer großräumigen Filterung verwendet. Anschließend wird wieder zu zwei unabhängigen Projektionen zurückgekehrt, die einzeln bei der gefilterten Rückprojektion zum Erzeugen der 3-D-Rekonstruktion verwendet werden.
Description
- Die Erfindung betrifft ein Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers.
- Eine 3D-Rekonstruktion ist ein Datensatz, bei dem einzelnen Volumenelementen (Voxeln) eines Körpers Zahlenwerte (Grauwerte) zugeordnet sind. In der Gesamtschau erhält man eine Aussage über die Struktur des Körpers. Üblicherweise werden 3D-Rekonstruktionen dadurch erzeugt, dass eine Folge von zweidimensionalen Abbildungen (Projektionen) aufgenommen wird, und zwar aus unterschiedlichen Blickwinkeln. Die unterschiedlichen Blickwinkel werden dadurch erzeugt, dass ein Röntgen-C-Bogen, an dem eine Röntgenstrahlenquelle und ein Röntgendetektor angeordnet sind, gedreht wird, d. h. in unterschiedliche Stellungen verbracht wird. Jede Projektion (also jede Aufnahme) wird nun ins Dreidimensionale rückprojiziert. Hierbei wird üblicherweise davon ausgegangen, dass entlang eines Projektionsstrahls von der Röntgenstrahlungsquelle zum Röntgendetektor gleichmäßige Beiträge zur Abschwächung des Röntgenstrahls gegeben sind. Dadurch, dass dies anhand einer Vielzahl von Projektionen durchgeführt wird, erhält man Strukturen in der 3D-Rekonstruktion.
- Ein Problem bei der 3D-Rekonstruktion bei größeren Körpern ist es, dass der Röntgendetektor aufgrund seiner eingeschränkten Größe lediglich einen Strahlenfächer erfasst, der den interessierenden Körper nicht vollständig überstreicht. Mit anderen Worten gibt es Körperregionen, die sich in der Abbildung auf dem Röntgendetektor nicht niederschlagen.
- In der
DE 10 2006 041 033.5 , die nach dem Anmeldetag der vorliegenden Anmeldung veröffentlicht wird, ist beschrieben, dass zur Lösung des Problems zu breiter Körper bei jeder Stellung des Röntgen-C-Bogens zwei (oder auch mehr) Projektionen aufgenommen werden, wobei von Projektion zu Projektion eine Rotation der Fokus-Detektoranordnung stattfindet. Es wird hierbei davon ausgegangen, dass die Röntgenstrahlungsquelle auf demselben Platz verbleibt. - Die beiden Projektionen werden mittels einer so genannten Homografie auf einen virtuellen Detektor abgebildet, so dass man eine kombinierte, größere Aufnahme erhält, die das gesamte Objekt abbildet. Diese kombinierte Aufnahme wird anschließend mittels eines herkömmlichen Rekonstruktionsverfahrens (z. B. unter Verwendung des so genannten Feldkampalgorithmus) rekonstruiert. Üblicherweise wird ein Filter verwendet, der sich über horizontale Linien im Detektor erstreckt. Dies kann ein Rampenfilter oder ein Hilbertfilter sein.
- Implizite Voraussetzung bei der Erzeugung des virtuellen Detektors in dem vorangemeldeten Verfahren ist es, dass sich die Röntgenstrahlungsquelle zwischen den beiden Projektionen nicht translatorisch bewegt. Steuerungstechnische Gründe, die mit dem Aufbau des Röntgen-C-Bogens und seiner Beweglichkeit zusammenhängen, können jedoch bewirken, dass die Quellposition für die beiden Projektionen nicht mehr identisch ist, sondern dass es unterschiedliche Quellpositionen für die beiden Projektionen gibt. Damit können die Projektionen nicht mehr artefaktfrei überlagert werden. Die Artefakte treten in der schließlich erzeugten 3D-Rekonstruktion des Körpers auf.
- Es ist Aufgabe der Erfindung, einen Weg aufzuzeigen, wie bei dem Verfahren aus der
DE 10 2006 041 033.5 Artefakte der genannten Art verhindert werden können. - Die Aufgabe wird durch ein Verfahren gemäß Patentanspruch 1 gelöst.
- Erfindungsgemäß umfasst das Verfahren somit die Schritte:
- a) Aufnehmen von Röntgen-Projektionspaaren jeweils für eine Vielzahl von Stellungen eines Röntgen-C-Bogens, an dem eine Röntgenstrahlungsquelle und ein Röntgendetektor angeordnet sind, wobei bei möglichst unveränderter Lage der Röntgenstrahlungsquelle die erste Projektion des Projektionspaars bei einer ersten Relativstellung von Röntgenstrahlungsquelle und Röntgendetektor und die zweite Projektion des Projektionspaars bei einer zweiten Relativstellung von Röntgenstrahlungsquelle und Röntgendetektor aufgenommen wird,
- b) Durchführen folgender Schritte zu jedem Projektionspaar:
- b1) Abbilden der Projektionsdaten der ersten Projektion auf einen virtuellen Detektor und der Projektionsdaten der zweiten Projektion auf denselben virtuellen Detektor,
- b2) Durchführen einer Filterung derart, dass auf die Projektionsdaten der ersten Projektion und die Projektionsdaten der zweiten Projektion gleichzeitig zurückgehende Abbildungsdaten erzeugt werden,
- b3) Umkehr der Abbildung der Projektionsdaten zum Erzeugen einer ersten, durch die Filterung abgeänderten Projektion und einer zweiten, durch die Filterung abgeänderten Projektion,
- b4) Verwenden beider abgeänderter Projektionen unabhängig voneinander bei einer gefilterten Rückprojektion zum Erzeugen der 3D-Rekonstruktion.
- Während in der
DE 10 2006 041 033.5 die auf den virtuellen Detektor projizierten Daten rückprojiziert werden, also beide Projektionen untrennbar gemeinsam in die 3D-Rekonstruktion eingehen, werden die Projektionen vorliegend getrennt voneinander bei der gefilterten Rückprojektion eingesetzt. Es ist jedoch nicht so, dass diese Projektionen getrennt voneinander verwendet werden, so wie sie ursprünglich aufgenommen sind. Da die Filterung langreichweitig sein soll und sich über die gesamte Breite des abzubildenden Körpers erstrecken soll, werden die Projektionen auf dem virtuellen Detektor zusammengefasst. Es kann dann eine langreichweitige Filterung wie beispielsweise eine Rampenfilterung oder eine Hilbertfilterung erfolgen derart, dass zumindest ein Teil der Daten (Grauwerte von einzelnen Pixeln der digitalen Röntgen-Projektion) so gestaltet ist, dass er nicht aus einer einzi gen der beiden Projektionen hervorgehen hätte können, sondern dass Daten aus beiden Projektionen eingegangen sind. Um wie gewünscht die beiden Projektionen dennoch getrennt voneinander bei der gefilterten Rückprojektion einzusetzen, wird nach der Filterung wieder eine Trennung vollzogen. - Es ist möglich, einen einzigen virtuellen Detektor zu definierten, auf den beide Projektionen mittels Homografie abgebildet werden. Nach Durchführen der Filterung kann dann die Homografie umgekehrt werden, und man erhält die geänderten Projektionen.
- Es hat sich als zweckmäßig erwiesen, zur Durchführung der Filterung zu jeder einzelnen der beiden Projektionen des Projektionspaars gesondert einen virtuellen Detektor zu verwenden. So kann der virtuelle Detektor dem realen Detektor einer der beiden Projektionen entsprechen, z. B. der ersten Projektion. Die genannte Abbildung der Projektionsdaten (Grauwerte der Pixel der Projektion) der ersten Projektion ist dann nichts anderes als die Einheitsabbildung. Die Abbildung der Projektionsdaten der zweiten Projektion erfolgt dann durch Berechnen eines Verbindungsstrahls zwischen der (gegebenenfalls durch Aufnehmen von Kalibrierungsbildern) ermittelten realen Quelle, wie sie bei der ersten Projektion vorgegeben ist, und dem realen Detektor, wie er bei der zweiten Projektion vorgegeben ist. Es werden dann dem realen Detektor (also seinen Bildpunkten) bei der zweiten Projektion zugeordnete Projektionsdaten (also Grauwerte) auf den Schnittpunkt des Verbindungsstrahls mit dem virtuellen Detektor gesetzt. Unter realem Detektor und realer Quelle wird ein geometrisches Konstrukt verstanden, das die geometrischen Relativbeziehungen zwischen Quelle und Detektor, wie sie bei den Projektionen vorlagen, korrekt wiedergibt, so dass der Strahlengang reproduzierbar ist.
- Der virtuelle Detektor ist somit nichts anderes als ein erweiterter Detektor. Zu den Projektionsdaten auf dem ursprünglichen Detektor treten lediglich am Rand Daten von der zwei ten Projektion hinzu. Nach dem Filtern muss dann lediglich der hinzugekommene Anteil weggelassen werden, d. h. diejenigen Datenorte auf dem virtuellen Detektor (also Bildpunkte), die Schnittpunkte des Verbindungsstrahls der Verbindung aus realer Quelle der ersten Projektion mit dem realen Detektor bei der zweiten Projektion waren, werden anschließend wieder weggelassen.
- Naturgemäß erhält man durch diesen Schritt eine erste Projektion, die unter Verwendung von Daten aus der zweiten Projektion gefiltert wurde. Um umgekehrt eine zweite Projektion zu erhalten, die unter Verwendung von Daten aus der ersten Projektion gefiltert wurde, können die Schritte b1) bis b3) bei Definition der Abbildung wie soeben erwähnt zusätzlich unter Vertauschung der Rollen von erster und zweiter Projektion wiederholt werden. Mit anderen Worten entspricht der virtuelle Detektor dann dem realen Detektor bei der zweiten Projektion, für die die Einheitsabbildung verwendet wird. Die reale Quelle bei der zweiten Projektion wird mit dem Detektor bei der ersten Projektion verbunden, und auf den Schnittpunkten des Verbindungsstrahls mit dem virtuellen Detektor, der der vergrößerte reale Detektor bei der zweiten Projektion ist, werden die Daten von der ersten Projektion eingetragen. Dann wird gefiltert, und man erhält die gefilterte zweite Projektion durch Weglassen der Datenorte (Bildpunkte), die als Schnittpunkt definiert waren.
- Bei geeigneter Wahl der Filterung sind die Daten aus der ersten Projektion und aus der zweiten Projektion in gewisser Weise miteinander korreliert. Es ist somit nicht unbedingt notwendig, die beiden Projektionen bei der gefilterten Rückprojektion so zu behandeln, als handle es sich um Projektionen von unterschiedlichen Stellungen des Röntgen-C-Bogens. Es kann berücksichtigt werden, dass die Projektion derselben Stellung des Röntgen-C-Bogens zugeordnet sind, indem gewisse Gewichtungen vorgenommen werden. Bei der gefilterten Rückprojektion können die Projektionsdaten beider Projektionen nach einer Gewichtungsfunktion gewichtet werden, die in Abhängig keit von den Bildpunktdaten eines virtuellen Detektors definiert ist. Damit wird wieder ein wenig an den virtuellen Detektor aus der
DE 10 2006 041 033.5 angeknüpft. - Da die Gewichtungsfunktion nicht unbedingt gewährleistet, dass auf einzelne Volumenelemente zurückgehende Daten gleich gewichtet sind, kann bei der 3D-Rekonstruktion eine voxelweise Normierung im Anschluss an Schritt b4) durchgeführt werden.
- Nachfolgend wird eine bevorzugte Ausführungsform der Erfindung unter Bezug auf die Zeichnung beschrieben, wobei
-
1 schematisch die Anordnung aus realer Röntgenstrahlungsquelle und Röntgendetektor bei dem der Erfindung zugrunde liegenden Verfahren veranschaulicht und -
2 für die Anordnung aus1 einen virtuellen Detektor und eine Gewichtungsfunktion veranschaulicht. - Zur Erzeugung einer 3D-Rekonstruktion eines Körpers wird eine Folge von digitalen Röntgen-Aufnahmen aufgenommen. Bei einer vorbestimmten Stellung eines Röntgen-C-Bogens werden zwei derartige digitale Röntgen-Projektionen erzeugt. Es erfolgt eine Relativdrehung der Einheit aus Röntgenstrahlungsquelle und Röntgendetektor des C-Bogens bei vorbestimmter Stellung des Röntgen-C-Bogens zwischen der Aufnahme der ersten Projektion des Projektionspaares und der der zweiten Projektion des Projektionspaars. Vorgesehen ist, dass die Röntgenstrahlungsquelle an demselben Ort verbleibt. In der Realität ergibt es sich, dass die Röntgenstrahlungsquelle ihren Ort zumindest leicht ändert. In
1 ist der Ort der Röntgenstrahlungsquelle für die erste Aufnahme mit Q1 bezeichnet und der Ort für die Röntgenstrahlungsquelle für die zweite Aufnahme mit Q2 bezeichnet. Der Verlauf des Röntgendetektors bei der ersten Projektion ist mit D1 bezeichnet, und der Verlauf des Röntgendetektors bei der zweiten Projektion mit D2. Den geometrischen Verlauf des Röntgendetektors kann man als „realen Detektor zur Projektion" bezeichnen. In derDE 10 2006 041 033.5 ist beschrieben, dass aus dem realen Detektor ein virtueller Detektor erzeugbar ist, wobei der virtuelle Detektor die Daten beider Projektionen aufnimmt. In dieser Druckschrift ist jedoch davon ausgegangen, dass Q1 gleich Q2 ist, was in der Realität nicht immer gegeben ist. Für den Fall, dass Q1 ungleich Q2, wie in1 gezeigt, muss das Verfahren aus derDE 10 2006 041 033.5 etwas abgewandelt werden:
Es werden zwei virtuelle Detektoren D1' und D2' definiert. Vorliegend sind diese beiden virtuellen Detektoren nichts anderes als Erweiterungen der realen Detektoren D1 und D2, wobei der Einfachheit halber in1 eine parallele Versetzung des Detektors gezeigt ist. Dem virtuellen Detektor D1' werden zunächst sämtliche Grauwerte der Projektion zugeordnet, die man mit dem realen Detektor D1 erhalten hat. Der virtuelle Detektor D1' erstreckt sich jedoch etwas über einen etwas größeren Bereich, damit ein größerer Körper, der mit einer einzigen Projektion allein nicht abzubilden ist, datenmäßig vollständig erfasst wird. Auf den virtuellen Detektor D1' sollen nun Daten aus der zweiten Projektion abgebildet werden. Hierzu wird eine Vielzahl von Strahlen definiert, die Bildpunkte auf dem realen Detektor D2 mit der Quelle Q1 verbindet. Beispielhaft ist der Strahl S gezeigt. Der Strahl S verbindet einen Bildpunkt B des realen Detektors D2 mit der realen Quelle Q1. Dem Bildpunkt B ist ein bestimmter Grauwert zugeordnet. Dieser Grauwert wird nun von dem realen Detektor D2 auf den virtuellen Detektor D1' übertragen: Der Datenpunkt DP enthält denselben Grauwert wie der Punkt B. In analoger Weise werden dem virtuellen Detektor D2' Daten zugeordnet, nämlich im größten Teil des virtuellen Detektors D2' die Daten aus dem realen Detektor D2, und in dem in1 linken Bereich von D2' Daten aus D1 durch Bilden von Projektionsstrahlen von der realen Quelle Q2 zum realen Detektor D1. Abermals werden die Daten bestimmter Bildpunkte dem Schnitt punkt der Verbindungslinie zugeordnet, so wie dem Schnittpunkt DP die Daten des Bildpunkts B zugeordnet werden. - Nachdem nun für den virtuellen Detektor D1' den einzelnen Datenpunkten (Pixeln) Grauwerte zugeordnet sind, kann eine großräumige Filterung der Daten erfolgen. Es kann ein so genannter Rampenfilter verwendet werden, der einzelne horizontale Detektorzeilen als Ganzes verarbeitet. Dadurch werden gleichzeitig Daten, die auf den realen Detektor D1 (also die erste Projektion) zurückgehen, und Daten, die auf den realen Detektor D2 (also auf die zweite Projektion) zurückgehen, gemeinsam in einem Filterschritt verarbeitet. Dies bewirkt, dass im in
1 linken Bereich des virtuellen Detektors D1' Datenpunkten Grauwerte zugeordnet sind, die durch die Grauwerte vom rechten Bereich beeinflusst sind und umgekehrt. Zur weiteren Vorgehensweise werden nun die anhand des Strahls S und ähnlicher Strahlen konstruierten Grauwerte vom virtuellen Detektor D1' weggelassen, und analog wird dies bei D2' durchgeführt. Der virtuelle Detektor D1' wird so gekürzt, dass sich wieder nichts anderes als der reale Detektor D1 ergibt, und der virtuelle Detektor D2' wird so gekürzt, dass sich nichts anderes als der reale Detektor D2 ergibt. Kein Detektor D1 und D2 für sich genügt, um einen bestimmten abzubildenden Körper in seiner ganzen Breite abzubilden. Das Verfahren hat dazu geführt, dass man ausgehend von einer ersten und einer zweiten Projektion eine Abbildung auf einen virtuellen Detektor D1' und D2' erzeugt hat, dann eine Filterung durchgeführt hat und anschließend wieder eine Trennung der Daten vorgenommen hat, so dass man schließlich die erste Projektion in abgewandelter Form, nämlich gefiltert, erhält und die zweite Projektion ebenfalls in abgewandelter Form, nämlich ebenfalls gefiltert, erhält. Die beiden Projektionen können nun im Rahmen einer gefilterten Rückprojektion zur Erzeugung einer 3D-Rekonstruktion verwendet werden. - Bei der Rückprojektion geht es darum, einen breiten Körper, der durch die einzelnen realen Detektoren D1 und D2 nicht vollständig abbildbar ist, insgesamt vollständig abzubilden.
- In den Randbereichen werden lediglich Daten aus dem Detektor D1 alleine bzw. dem Detektor D2 alleine verwendet. Die Detektoren D1 und D2 haben jedoch einen Überlappbereich. Für diesen können beide Projektionen gleichzeitig verwendet werden. Würden die realen Quellen Q1 und Q2 zusammenfallen, kann man einfach die beiden Projektionen vor der Rückprojektion jeweils mit einer konstanten Gewichtsfunktion multiplizieren. Die Gewichtsfunktion kann als Überblendungsfunktion auf einen virtuellen Detektor definiert werden. Vorliegend fallen die Quellen Q1 und Q2 nicht zusammen. Dies hat folgenden Effekt: Ein von der Quelle Q2 ausgehender Röntgenstrahl, der auf dem realen Detektor D2 auftrifft, durchläuft die Punkte P1 und P2 gleichzeitig. Eine für den Bildpunkt B' auf dem Detektor D2 definiertes Gewicht betrifft somit gleichzeitig Volumenelemente um den Punkt P1 und Volumenelemente um den Punkt P2. Hingegen liegen aus Sicht der Röntgenquelle Q1 die Punkte P1 und P2 nicht auf demselben Strahl: Der Strahl S1, der durch den Punkt P1 verläuft, trifft auf einem anderen Bildpunkt B1 des realen Detektors D1 auf, als der Strahl S2, der durch den Punkt P2 läuft, der nämlich auf dem Bildpunkt B2 auftrifft. Bei einer Gewichtung würden somit die Volumenelemente um den Punkt P1 und um den Punkt P2 herum unterschiedliche Gewichte erhalten.
- Vorliegend wird zur Erzeugung einer passenden Gewichtung für den Überlappbereich eine Abbildung (Homografie) auf einen virtuellen Detektor D' verwendet. Für den virtuellen Detektor D' wird eine Gewichtungsfunktion G1, die den Bildpunkten des realen Detektors D1 zugeordnet ist, und eine Gewichtungsfunktion G2, die den Bildpunkten des realen Detektors D2 zugeordnet ist, verwendet. Durch das Auseinanderfallen der Quellen Q1 und Q2 kann es sein, dass die Summe der Gewichte bei der Rückprojektion nicht eins ist. Es wird dann eine voxelweise Normierung des Bildes durchgeführt, bei der ortsabhängige Gewichte verwendet werden.
- Bei der vorliegenden Erfindung werden Artefakte, die auf das Auseinanderfallen der realen Quellen Q1 und Q2 zurückzuführen ist, dadurch vermieden, dass die beiden Projektionen unabhängig voneinander bei der gefilterten Rückprojektion verwendet werden. Die Grauwerte der Projektionen werden jedoch einer Filterung unterzogen, bei der die jeweils andere Projektion berücksichtigt wird, so dass langreichweitige Effekte in einem abzubildenden Körper berücksichtigt werden, der durch eine Projektion alleine nicht vollständig abgebildet ist.
- 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
-
- - DE 102006041033 [0004, 0007, 0010, 0015, 0020, 0020]
Claims (5)
- Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers, mit den Schritten: a) Aufnehmen von digitalen Röntgen-Projektionspaaren jeweils für eine Vielzahl von Stellungen eines Röntgen-C-Bogens, an dem eine Röntgenstrahlungsquelle (Q1, Q2) und ein Röntgendetektor (D1, D2) angeordnet sind, wobei bei möglichst unveränderter Lage der Röntgenstrahlungsquelle (Q1, Q2) die erste Projektion des Projektionspaars bei einer ersten Relativstellung von Röntgenstrahlungsquelle (Q1) und Röntgendetektor (D1) und die zweite Projektion des Projektionspaars bei einer zweiten Relativstellung von Röntgenstrahlungsquelle (Q2) und Röntgendetektor (D2) aufgenommen wird, b) Durchführen folgender Schritte zu jedem Projektionspaar: b1) Abbilden der Projektionsdaten der ersten Projektion auf einen virtuellen Detektor (D1', D2') und der Projektionsdaten der zweiten Projektion auf denselben virtuellen Detektor (D1', D2'), b2) Durchführen einer Filterung derart, dass auf die Projektionsdaten der ersten Projektion und die Projektionsdaten der zweiten Projektion zugleich zurückgehende Abbildungsdaten erzeugt werden, b3) Umkehr der Abbildung der Projektionsdaten zum Erzeugen einer ersten, durch die Filterung abgeänderten Projektion und einer zweiten, durch die Filterung abgeänderten Projektion, b4) Verwenden beider abgeänderter Projektionen unabhängig voneinander bei einer gefilterten Rückprojektion zum Erzeugen der 3D-Rekonstruktion.
- Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass in den Schritten b1) bis b3) der virtuelle Detektor (D1') dem realen Detektor (D1) bei der ersten Projektion entspricht und die Abbildung bei den Projektionsdaten der ersten Projektion die Einheitsabbildung ist und die Abbildung bei den Projektionsdaten der zweiten Pro jektion durch Berechnen eines Verbindungsstrahls (S) zwischen der ermittelten realen Quelle (Q1) bei der ersten Projektion und dem realen Detektor (D2) bei der zweiten Projektion und ferner durch Setzen von entsprechenden, dem realen Detektor (D2) bei der zweiten Projektion zugeordneten Projektionsdaten auf den Schnittpunkt (DP) des Verbindungsstrahls (S) mit dem virtuellen Detektor (D1') erfolgt, und dass die Umkehr der Abbildung lediglich das Weglassen der auf die zweite Projektion zurückgehenden Datenorte (DP) umfasst.
- Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass die Schritte b1) bis b3) zusätzlich unter Vertauschung der Rollen von erster und zweiter Projektion wiederholt werden.
- Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass bei der gefilterten Rückprojektion die Projektionsdaten beider Projektionen nach Gewichtungsfunktionen (G1, G2) gewichtet werden.
- Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass der gefilterten Rückprojektion eine voxelweise Normierung der 3D-Rekonstruktion folgt.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102007026115.4A DE102007026115B4 (de) | 2007-06-05 | 2007-06-05 | Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers |
US12/156,932 US7742557B2 (en) | 2007-06-05 | 2008-06-05 | Method for generating a 3D reconstruction of a body |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102007026115.4A DE102007026115B4 (de) | 2007-06-05 | 2007-06-05 | Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers |
Publications (2)
Publication Number | Publication Date |
---|---|
DE102007026115A1 true DE102007026115A1 (de) | 2008-12-18 |
DE102007026115B4 DE102007026115B4 (de) | 2017-10-12 |
Family
ID=39985990
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE102007026115.4A Expired - Fee Related DE102007026115B4 (de) | 2007-06-05 | 2007-06-05 | Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers |
Country Status (2)
Country | Link |
---|---|
US (1) | US7742557B2 (de) |
DE (1) | DE102007026115B4 (de) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007042333A1 (de) | 2007-09-06 | 2009-03-12 | Siemens Ag | Verfahren zum Ermitteln einer Abbildungsvorschrift und Verfahren zum Erzeugen einer 3D-Rekonstruktion |
DE102009017437A1 (de) | 2009-04-15 | 2010-10-28 | Siemens Aktiengesellschaft | Verfahren zum Ermitteln einer Mehrzahl von Abbildungsvorschriften, insbesondere zu unterschiedlichen Stellungen eines Röntgenstrahlendetektors bezüglich einer Röntgenstrahlenquelle |
DE102011006991A1 (de) | 2011-04-07 | 2012-10-11 | Siemens Aktiengesellschaft | Röntgenverfahren und Röntgenvorrichtung zum Zusammenfügen von Röntgenaufnahmen und Ermitteln dreidimensionaler Volumendaten |
DE102020209703A1 (de) | 2020-07-31 | 2022-02-03 | Siemens Healthcare Gmbh | Bereitstellen eines 3D-Bilddatensatzes eines Untersuchungsobjekts |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10670740B2 (en) | 2012-02-14 | 2020-06-02 | American Science And Engineering, Inc. | Spectral discrimination using wavelength-shifting fiber-coupled scintillation detectors |
US10070828B2 (en) | 2013-03-05 | 2018-09-11 | Nview Medical Inc. | Imaging systems and related apparatus and methods |
US10846860B2 (en) | 2013-03-05 | 2020-11-24 | Nview Medical Inc. | Systems and methods for x-ray tomosynthesis image reconstruction |
US9171365B2 (en) * | 2013-11-29 | 2015-10-27 | Kabushiki Kaisha Toshiba | Distance driven computation balancing |
US9986983B2 (en) | 2014-10-31 | 2018-06-05 | Covidien Lp | Computed tomography enhanced fluoroscopic system, device, and method of utilizing the same |
JP6746603B2 (ja) | 2015-03-20 | 2020-08-26 | ラピスカン システムズ、インコーポレイテッド | 手持ち式携帯型後方散乱検査システム |
US10716525B2 (en) | 2015-08-06 | 2020-07-21 | Covidien Lp | System and method for navigating to target and performing procedure on target utilizing fluoroscopic-based local three dimensional volume reconstruction |
US10702226B2 (en) | 2015-08-06 | 2020-07-07 | Covidien Lp | System and method for local three dimensional volume reconstruction using a standard fluoroscope |
US10674982B2 (en) | 2015-08-06 | 2020-06-09 | Covidien Lp | System and method for local three dimensional volume reconstruction using a standard fluoroscope |
US11172895B2 (en) | 2015-12-07 | 2021-11-16 | Covidien Lp | Visualization, navigation, and planning with electromagnetic navigation bronchoscopy and cone beam computed tomography integrated |
US11051886B2 (en) | 2016-09-27 | 2021-07-06 | Covidien Lp | Systems and methods for performing a surgical navigation procedure |
US10699448B2 (en) | 2017-06-29 | 2020-06-30 | Covidien Lp | System and method for identifying, marking and navigating to a target using real time two dimensional fluoroscopic data |
CN111373448B (zh) | 2017-09-22 | 2023-12-08 | 尼维医疗公司 | 使用机器学习正则化器的图像重建 |
US10893843B2 (en) | 2017-10-10 | 2021-01-19 | Covidien Lp | System and method for identifying and marking a target in a fluoroscopic three-dimensional reconstruction |
US10905498B2 (en) | 2018-02-08 | 2021-02-02 | Covidien Lp | System and method for catheter detection in fluoroscopic images and updating displayed position of catheter |
US10893842B2 (en) | 2018-02-08 | 2021-01-19 | Covidien Lp | System and method for pose estimation of an imaging device and for determining the location of a medical device with respect to a target |
US10830911B2 (en) | 2018-06-20 | 2020-11-10 | American Science And Engineering, Inc. | Wavelength-shifting sheet-coupled scintillation detectors |
US11340361B1 (en) | 2020-11-23 | 2022-05-24 | American Science And Engineering, Inc. | Wireless transmission detector panel for an X-ray scanner |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006041033A1 (de) | 2006-09-01 | 2008-03-13 | Siemens Ag | Verfahren zur Rekonstruktion eines dreidimensionalen Bildvolumens und Röntgengeräte |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5319693A (en) * | 1992-12-30 | 1994-06-07 | General Electric Company | Three dimensional computerized tomography scanning configuration for imaging large objects with smaller area detectors |
US5740224A (en) * | 1994-09-27 | 1998-04-14 | University Of Delaware | Cone beam synthetic arrays in three-dimensional computerized tomography |
US6236704B1 (en) * | 1999-06-30 | 2001-05-22 | Siemens Corporate Research, Inc. | Method and apparatus using a virtual detector for three-dimensional reconstruction from x-ray images |
JP2002263094A (ja) * | 2001-03-06 | 2002-09-17 | Shimadzu Corp | X線透視撮影装置 |
US7106825B2 (en) * | 2002-08-21 | 2006-09-12 | Breakaway Imaging, Llc | Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system |
US6944260B2 (en) * | 2003-11-11 | 2005-09-13 | Ge Medical Systems Global Technology Company, Llc | Methods and apparatus for artifact reduction in computed tomography imaging systems |
US8218720B2 (en) * | 2007-03-12 | 2012-07-10 | Varian Medical Systems, Inc. | Method and apparatus to facilitate reconstructing an image using fan-beam data |
-
2007
- 2007-06-05 DE DE102007026115.4A patent/DE102007026115B4/de not_active Expired - Fee Related
-
2008
- 2008-06-05 US US12/156,932 patent/US7742557B2/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102006041033A1 (de) | 2006-09-01 | 2008-03-13 | Siemens Ag | Verfahren zur Rekonstruktion eines dreidimensionalen Bildvolumens und Röntgengeräte |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007042333A1 (de) | 2007-09-06 | 2009-03-12 | Siemens Ag | Verfahren zum Ermitteln einer Abbildungsvorschrift und Verfahren zum Erzeugen einer 3D-Rekonstruktion |
DE102009017437A1 (de) | 2009-04-15 | 2010-10-28 | Siemens Aktiengesellschaft | Verfahren zum Ermitteln einer Mehrzahl von Abbildungsvorschriften, insbesondere zu unterschiedlichen Stellungen eines Röntgenstrahlendetektors bezüglich einer Röntgenstrahlenquelle |
DE102011006991A1 (de) | 2011-04-07 | 2012-10-11 | Siemens Aktiengesellschaft | Röntgenverfahren und Röntgenvorrichtung zum Zusammenfügen von Röntgenaufnahmen und Ermitteln dreidimensionaler Volumendaten |
US8693622B2 (en) | 2011-04-07 | 2014-04-08 | Siemens Aktiengesellschaft | X-ray method and X-ray system for merging X-ray images and determining three-dimensional volume data |
DE102011006991B4 (de) | 2011-04-07 | 2018-04-05 | Siemens Healthcare Gmbh | Röntgenverfahren und Röntgenvorrichtung zum Zusammenfügen von Röntgenaufnahmen und Ermitteln dreidimensionaler Volumendaten |
DE102020209703A1 (de) | 2020-07-31 | 2022-02-03 | Siemens Healthcare Gmbh | Bereitstellen eines 3D-Bilddatensatzes eines Untersuchungsobjekts |
US11717239B2 (en) | 2020-07-31 | 2023-08-08 | Siemens Healthcare Gmbh | Providing a 3D image dataset of an examination object |
Also Published As
Publication number | Publication date |
---|---|
US20080304617A1 (en) | 2008-12-11 |
US7742557B2 (en) | 2010-06-22 |
DE102007026115B4 (de) | 2017-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102007026115A1 (de) | Verfahren zum Erzeugen einer 3D-Rekonstruktion eines Körpers | |
DE102006045423B4 (de) | 07.09.07Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur | |
DE102005012654B4 (de) | Verfahren und Computertomographie-System zur Erstellung tomographischer Aufnahmen eines Objektes | |
DE10100572A1 (de) | Verfahren zur Darstellung des Blutflusses in einem Gefäßbaum | |
DE102010053458B4 (de) | Verfahren zur Aufbereitung von IR-Bildern und korrespondierende Wärmebildkamera | |
DE102019215460A1 (de) | Verfahren und Vorrichtung zur Rauschreduktion von Bildaufnahmen | |
EP3407299B1 (de) | Verfahren und vorrichtung zur findung einer positionierung und zur datenfusion | |
DE3726595A1 (de) | Vorrichtung zur korrektur von zwischen verschiedenen ansichten stattfindender bewegung in einem abbildungssystem sowie verfahren dafuer | |
DE19945636A1 (de) | Verfahren und Vorrichtung zum dreidimensionalen CT-Abbilden eines interessierenden Bereichs in einem Objekt | |
DE102012204019A1 (de) | Verfahren zur Reduzierung von Bewegungsartefakten | |
DE102005051620A1 (de) | Verfahren zur Rekonstruktion einer tomographischen Darstellung eines Objektes | |
DE102011086456A1 (de) | Rekonstruktion von Bilddaten | |
DE102016219887A1 (de) | Verfahren und System zur Nutzung von Messdaten | |
DE102006005804A1 (de) | Verfahren zur Rauschreduktion in tomographischen Bilddatensätzen | |
DE102008052690B4 (de) | Vorrichtung und Verfahren zur Rekonstruktion und Visualisierung von Projektionsdaten | |
DE102011075917A1 (de) | Verfahren zum Bereitstellen eines 3D-Bilddatensatzes mit unterdrückten Messfeldüberschreitungsartefakten und Computertomograph | |
EP1498851A1 (de) | Bestimmung einer dreidimensionalen Form eines Körpers, insbesondere einer anatomischen Struktur, aus zweidimensionalen Projektionsbildern | |
DE102008038357B3 (de) | Verfahren zur Erzeugung von 2D-Schnittbildern aus 3D-Projektionsdaten, die mittels eines CT-Systems von einem metallische Anteile enthaltenden Untersuchungsobjekt erfasst wurden | |
EP3579193B1 (de) | Verfahren zur ermittlung von bildwerten in markierten pixeln wenigstens eines projektionsbildes, röntgeneinrichtung, computerprogramm und elektronisch lesbarer datenträger | |
DE102010029187A1 (de) | Verfahren zur Rekonstruktion eines einer Schnittebene durch ein aufgenommenes Objekt entsprechenden zweidimensionalen Schichtbildes und Röntgeneinrichtung | |
DE102012217613A1 (de) | Verfahren zur Reduzierung von Artefakten bei der Rekonstruktion von Bilddatensätzen aus Projektionsbildern | |
DE102015007934A1 (de) | Verfahren und Computerprogrammprodukt zum Erzeugen eines artefaktreduzierten Voxeldatensatzes | |
DE10307331B4 (de) | Bildgebendes Verfahren zur rechnergestützten Auswertung computer-tomographischer Messungen durch direkte iterative Rekonstruktion | |
DE3420575A1 (de) | Reprojektionsanordnung | |
DE102008051157B4 (de) | Verfahren zum Erhalten einer 3D-Rekonstruktion eines Objekts |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
OP8 | Request for examination as to paragraph 44 patent law | ||
8120 | Willingness to grant licences paragraph 23 | ||
R085 | Willingness to licence withdrawn | ||
R085 | Willingness to licence withdrawn |
Effective date: 20140911 |
|
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 | ||
R079 | Amendment of ipc main class |
Free format text: PREVIOUS MAIN CLASS: G01N0023060000 Ipc: G01N0023046000 |
|
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |