DE102011088265A1 - Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten - Google Patents

Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten Download PDF

Info

Publication number
DE102011088265A1
DE102011088265A1 DE201110088265 DE102011088265A DE102011088265A1 DE 102011088265 A1 DE102011088265 A1 DE 102011088265A1 DE 201110088265 DE201110088265 DE 201110088265 DE 102011088265 A DE102011088265 A DE 102011088265A DE 102011088265 A1 DE102011088265 A1 DE 102011088265A1
Authority
DE
Germany
Prior art keywords
images
image
subtraction
ray
recording
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
Application number
DE201110088265
Other languages
English (en)
Other versions
DE102011088265B4 (de
Inventor
Christian Kruschel
Christoph Köhler
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Healthcare GmbH
Original Assignee
Siemens AG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Priority to DE102011088265.0A priority Critical patent/DE102011088265B4/de
Priority to US13/710,727 priority patent/US8768037B2/en
Publication of DE102011088265A1 publication Critical patent/DE102011088265A1/de
Application granted granted Critical
Publication of DE102011088265B4 publication Critical patent/DE102011088265B4/de
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5252Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data removing objects from field of view, e.g. removing patient table from a CT image
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/02Investigating 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/04Investigating 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20221Image fusion; Image merging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Veterinary Medicine (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

Verfahren zur Korrektur von aufgrund eines fest mit einem Röntgendetektor (17) verbundenen Streustrahlenrasters (19) auftretenden Bildartefakten bei der Aufnahme von dreidimensionalen Bilddatensätzen (11) mit einer Röntgeneinrichtung (14) mit einem C-Bogen (15), an dem sich gegenüberliegend ein Röntgenstrahler (16) und der Röntgendetektor (17) angeordnet sind, wobei aus unterschiedlichen Aufnahmegeometrien zweidimensionale Projektionsbilder (10) eines Objekts aufgenommen und zur Rekonstruktion des dreidimensionalen Bilddatensatzes (11) verwendet werden, umfassend folgenden Schritte: – Aufnahme jeweils eines Kalibrationsbildes (2) für jede Aufnahmegeometrie, – Ermittlung eines Durchschnittsbildes (4) aus allen aufgenommenen Kalibrationsbildern (2), – Ermittlung von Subtraktionsbildern (6) für alle Aufnahmegeometrien durch Subtraktion des Durchschnittsbildes (4) von den entsprechenden Kalibrationsbildern (2), – Glättung der Subtraktionsbilder (6) zur Entrauschung der Subtraktionsbilder (6), und – Speichern der Subtraktionsbilder (6) und des Durchschnittsbildes (4), wobei zur Korrektur von aufgenommenen Projektionsbildern (10) für jede Aufnahmegeometrie ein Korrekturbild durch Addition des der Aufnahmegeometrie zugeordneten Subtraktionsbildes (6) und des Durchschnittsbildes (4) ermittelt und das Korrekturbild zur Korrektur des Projektionsbildes (10) verwendet wird.

Description

  • Die Erfindung betrifft ein Verfahren zur Korrektur von aufgrund eines fest mit einem Röntgendetektor verbundenen Streustrahlenrasters auftretenden Bildartefakten bei der Aufnahme von dreidimensionalen Bilddatensätzen mit einer Röntgeneinrichtung mit einem C-Bogen, an dem sich gegenüberliegend ein Röntgenstrahler und der Röntgendetektor angeordnet sind, wobei aus unterschiedlichen Aufnahmegeometrien zweidimensionale Projektionsbilder eines Objekts aufgenommen und zur Rekonstruktion des dreidimensionalen Bilddatensatzes verwendet werden. Daneben betrifft die Erfindung eine Röntgeneinrichtung.
  • C-Bogen-Röntgeneinrichtungen, also Röntgeneinrichtungen, die einen C-Bogen aufweisen, an dem sich gegenüberliegend ein Röntgenstrahler und ein Röntgendetektor vorgesehen sind, sind im Stand der Technik bereits bekannt. Aufgrund wenigstens eines Bewegungsfreiheitsgrades des C-Bogens, insbesondere wenigstens einer Rotationsmöglichkeit, können der Röntgendetektor und der Röntgenstrahler so um ein aufzunehmendes Objekt, insbesondere einen Patienten, verfahren werden, dass verschiedene Zielgebiete unter unterschiedlichen Projektionswinkeln, also aus unterschiedlichen Aufnahmegeometrien, aufgenommen werden können.
  • Um Einflüsse durch Streustrahlung zu vermeiden, weisen derartige Röntgeneinrichtungen meist ein sogenanntes Streustrahlenraster auf, welches letztlich aus einem auf den Röntgendetektor aufgesetzten Gitter besteht. Bevorzugt werden dabei sogenannte fokussierte Streustrahlenraster verwendet, die fest mit dem Röntgendetektor verbunden sind und auf eine übliche Entfernung zwischen dem Röntgenstrahler und dem Röntgendetektor eingestellt sind, wie sie für gezielte, insbesondere hoch aufgelöste diagnostische zweidimensionale Röntgenbildaufnahmen, sogenannte Radiographie-Aufnahmen, benötigt wird.
  • Inzwischen wurde jedoch auch vorgeschlagen, mittels einer Röntgeneinrichtung mit einem C-Bogen dreidimensionale Bilddatensätze aufzunehmen, indem ähnlich der Computertomographie (CT) vorgegangen wird, daher auch die häufige Bezeichnung C-Bogen-CT. Dabei werden entlang einer Aufnahmetrajektorie, meist durch Rotation des C-Bogens realisiert, aus verschiedenen Aufnahmegeometrien, mithin unter verschiedenen Projektionsrichtungen, zweidimensionale Projektionsbilder aufgenommen, aus denen über übliche Verfahren dreidimensionale Bilddatensätze rekonstruiert werden können, beispielsweise mittels der iterativen Rekonstruktion oder der gefilterten Rückprojektion (FBP). Der dreidimensionale Bilddatensatz enthält dabei letztlich ein dreidimensionales Volumen von Schwächungskoeffizienten, das das Ziel beschreibt.
  • Das erwähnte Streustrahlenraster (englisch auch „anti scatter grid“) führt nicht nur zu einer Reduzierung der Streustrahlung, sondern auch der gewünschten Primärstrahlung, nachdem es eine auch die Primärstrahlung teilweise blockierende Struktur, insbesondere Gitterstruktur, aufweist. Dies führt zu Artefakten im Röntgen-Projektionsbild, die wiederum zu einer reduzierten Bildqualität im rekonstruierten Bilddatensatz führen. Um diese Artefakte zu reduzieren, ist es bekannt, eine pixelbasierte Gain-Korrektur mittels eines Korrekturbildes vorzunehmen, welches beispielsweise durch Kalibration erhalten werden kann. Dabei kann durch ein solches Korrekturbild nicht nur der Effekt der Gitterstrukturen des Streustrahlenrasters berücksichtigt werden, sondern auch ein gerade bei der C-Bogen-CT häufig auftretendes Problem, nachdem ein größerer Abstand zwischen dem Röntgenstrahler und dem Röntgendetektor verwendet wird, als der, für den das Streustrahlenraster fokussiert ist, so dass letztlich bei defokussiertem Streustrahlenraster aufgenommen wird, was zu einer ungleichen Intensitätsverteilung führt, insbesondere einer zu den Rändern des Detektors hin abnehmenden Intensität selbst bei homogener Ausleuchtung.
  • Insbesondere im Zusammenhang mit einer solchen Intensitätsschwankung aufgrund eines defokussierten Streustrahlenrasters ist auch ein weiterer nachteiliger Effekt bei der C-Bogen-CT zu betrachten, nämlich die Tatsache, dass aufgrund mechanischer Ungenauigkeiten der Röntgeneinrichtung zu verschiedenen Aufnahmezeitpunkten, also bei verschiedenen Aufnahmegeometrien, Abweichungen auftreten, insbesondere im Hinblick auf die Position des Fokus des Röntgenstrahlers relativ zum Röntgendetektor. Es kommt mithin zu Verschiebungen von Artefaktmustern und insbesondere auch der genannten Intensitätsverteilung, was zu deutlich sichtbaren Kanten im rekonstruierten Volumen führen kann, obwohl diese Effekte in den Projektionsbildern kaum feststellbar sind. Derartige mechanische Ungenauigkeiten der Röntgeneinrichtung basieren im Wesentlichen darauf, dass der C-Bogen nicht vollständig starr ist und daher beispielsweise in verschiedenen Stellungen des C-Bogens ein unterschiedlicher Einfluss der Schwerkraft auf die Röntgeneinrichtung vorliegt.
  • Die genannten mechanischen Ungenauigkeiten erlauben es letztlich nicht, eine globale Korrektur, mithin ein einziges Korrekturbild, für alle Aufnahmegeometrien zu verwenden, da dies wiederum zum Beibehalten von Artefakten oder zum Auftreten neuer Artefakte, wie beschrieben, führen kann.
  • Zur Lösung dieses Problems wurde vorgeschlagen, ein separates, pixelbasiertes Korrekturbild für jede Aufnahmegeometrie, mithin jedes aufzunehmende Projektionsbild, zu erstellen und vorzuhalten. Das Korrekturbild wird durch eine Kalibrierungsmessung bestimmt, wobei beispielsweise, um zu ähnlichen Intensitätswerten wie bei der Aufnahme eines Patienten zu gelangen, eine Kupferplatte in den Strahlengang eingebracht werden kann. Da das Rauschen ein kritischer Parameter in der Rekonstruktion ist, nachdem beispielsweise für Patienten das Rauschen direkt mit der Patientendosis korreliert ist, darf das Korrekturbild nur einen sehr geringen Rauschbeitrag enthalten. Um gute Korrekturbilder zu erhalten, werden daher im Stand der Technik eine Vielzahl von Kalibrationsbildern für jede Aufnahmegeometrie aufgenommen und gemittelt, was einen äußerst hohen Kalibrationsaufwand erfordert. Eine Entrauschung der Bilder durch eine einfache Glättung, beispielsweise einen Gauss-Filter, löst das Problem nicht, nachdem damit auch das Muster der Artefakte auf den Korrekturbildern zerstört wird, so dass die Korrektur die zerstörten Artefaktmuster nicht mehr beseitigen kann. Ein weiterer Nachteil des bekannten Verfahrens ist, dass eine sehr hohe Menge an Speicherplatz erforderlich ist, um Korrekturbilder für jede Aufnahmegeometrie vorzuhalten, denn im Rahmen einer Komprimierung würden auch die wesentlichen Strukturen verloren gehen.
  • Der Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren zur Korrektur von aufgenommenen Projektionsbildern im Hinblick auf durch ein Streustrahlenraster erfolgte Artefakte anzugeben, welches die mechanischen Unzulänglichkeiten hinreichend berücksichtigt, mit einem reduzierten Kalibrationsaufwand auskommt und insbesondere auch den Speicherplatzbedarf für Korrekturbilder reduziert.
  • Zur Lösung dieser Aufgabe sind bei einem Verfahren der eingangs genannten Art erfindungsgemäß folgende Schritte vorgesehen:
    • – Aufnahme jeweils eines Kalibrationsbildes für jede Aufnahmegeometrie,
    • – Ermittlung eines Durchschnittsbildes aus allen aufgenommenen Kalibrationsbildern,
    • – Ermittlung von Subtraktionsbildern für alle Aufnahmegeometrien durch Subtraktion des Durchschnittsbildes von den entsprechenden Kalibrationsbildern,
    • – Glättung der Subtraktionsbilder zur Entrauschung der Subtraktionsbilders, und
    • – Speichern der Subtraktionsbilder und des Durchschnittsbildes, wobei zur Korrektur von aufgenommenen Projektionsbildern für jede Aufnahmegeometrie ein Korrekturbild durch Addition des der Aufnahmegeometrie zugeordneten Subtraktionsbildes und des Durchschnittsbildes ermittelt und das Korrekturbild zur Korrektur des Projektionsbildes verwendet wird.
  • Erfindungsgemäß wird also das Problem adressiert, ein geglättetes, also entrauschtes Korrekturbild für jede Aufnahmegeometrie zu ermitteln, welches auch weiterhin das Artefaktmuster, welches durch das Streustrahlenraster ausgelöst wird, enthält. Hierzu wird zunächst ein einziges Kalibrationsbild für jede Aufnahmegeometrie aufgenommen. Diese einzelnen Kalibrationsbilder sind, wie bereits erwähnt wurde, verrauscht, so dass sie nicht unmittelbar als Korrekturbilder eingesetzt werden können.
  • Mithin wird erfindungsgemäß nun vorgeschlagen, zunächst ein Durchschnittsbild aller aufgenommenen Kalibrationsbilder zu errechnen. Nachdem das Streustrahlenraster fest mit dem Röntgendetektor verbunden ist, sind die entsprechenden unmittelbar durch das Streustrahlenraster ausgelösten Artefaktmuster mit hohen Frequenzen im Ortsraum immer gleich, so dass das Durchschnittsbild diese Information hoher Ortsfrequenzen des Artefaktmusters des Streustrahlenrasters enthält, während das Rauschen aufgrund der statistischen Betrachtung einer Vielzahl von Kalibrationsbildern reduziert ist. Die angesprochenen mechanischen Ungenauigkeiten der Röntgeneinrichtung tragen zu diesen Hochfrequenzinformationen aufgrund der nicht vorhandenen relativen Bewegung zwischen Streustrahlenraster und Röntgendetektor nicht bei.
  • Das Durchschnittsbild wird sodann von jedem aufgenommenen Kalibrationsbild abgezogen, um Subtraktionsbilder zu erhalten. Diese Subtraktionsbilder enthalten mithin das Artefaktmuster mit hohen Ortsfrequenzen nun nicht mehr und können daher geglättet werden. Zu korrigierende verbleibende Muster, beispielsweise eine Intensitätsverteilung aufgrund eines nicht für den Abstand zwischen Röntgendetektor und Röntgenstrahler bei der Aufnahme der Kalibrationsbilder fokussiertem Streustrahlenrasters ändern sich auf eher großen Abständen, weisen mithin niedrige Ortsfrequenzen auf, die auch bei einer Glättung, die erfindungsgemäß bevorzugt lokal sein kann, erhalten bleiben. Auf diese Weise ist es also möglich, die Subtraktionsbilder zu „entrauschen“ und sodann, insbesondere komprimiert, gemeinsam mit dem Durchschnittsbild abzuspeichern.
  • Um das tatsächlich zu verwendende Korrekturbild für eine gegebene Aufnahmegeometrie zu erhalten, muss nun lediglich – gegebenenfalls nach einer Dekompression des gespeicherten Subtraktionsbildes – das Durchschnittsbild wieder zu dem Subtraktionsbild addiert werden. Letztlich werden die von der konkreten Aufnahmegeometrie unabhängigen hochfrequenten Informationen der Artefaktmuster durch die Ermittlung des Durchschnittsbildes und dessen Subtraktion zum einen anhand des Durchschnittsbildes aufbewahrt, zum anderen der weiteren Betrachtung im Hinblick auf die Entrauschung entzogen, da erkannt wurde, dass die von der Aufnahmegeometrie abhängigen Artefakt-Informationen bei den niedrigen Ortsfrequenzen liegen. Mithin ist eine Entrauschung durch Glättung genau dann möglich, wenn die hochfrequenten Artefaktmuster nicht mehr verloren werden können.
  • Wie dargestellt wurde, bezieht sich das erfindungsgemäße Verfahren in seiner konkreten Problemlösung hauptsächlich auf Fälle, in denen bei der Aufnahme der Projektionsbilder ein defokussiertes Streustrahlenraster vorliegt und gleichzeitig mechanische Ungenauigkeiten, insbesondere ausgelöst aufgrund der nicht perfekten Starrheit des C-Bogens, insbesondere bezüglich der Schwerkraft. Mithin lässt sich das erfindungsgemäße Verfahren besonders vorteilhaft dann anwenden, wenn das Streustrahlenraster für eine andere Entfernung zwischen dem Röntgendetektor und dem Röntgenstrahler als die bei der Aufnahme der Projektionsbilder genutzte Entfernung fokussiert ist und/oder unter dem Einfluss der Schwerkraft eine Verschiebung des Röntgenstrahlers gegen den Röntgendetektor und mithin das Streustrahlenraster bei unterschiedlichen Positionen des C-Bogens auftritt.
  • Insgesamt ermöglicht es das erfindungsgemäße Verfahren bei einem deutlich reduzierten Kalibrierungsaufwand das Artefaktniveau in rekonstruierten Bilddatensätzen deutlich zu reduzieren, nachdem das Rauschen trotz Aufnahme nur eines einzigen Kalibrationsbildes deutlich reduziert werden kann.
  • Während es grundsätzlich denkbar ist, eine Glättung unter Verwendung eines Gauss-Filters in Betracht zu ziehen, ist diese jedoch nicht lokal und hat daher gegebenenfalls den Verlust von Artefaktmustern mittlerer oder auch niedriger Ortsfrequenzen zur Folge.
  • Daher kann im Rahmen der vorliegenden Erfindung in besonders bevorzugter Weise vorgesehen sein, dass ein auf Variation basierender, lokaler Glättungsalgorithmus zur Glättung der Subtraktionsbilder eingesetzt wird, wobei eine Minimierung der Summe einer Norm und eines insbesondere mit einem Gewichtungsfaktor eingehenden Regularisierungs-Funktionals vorgenommen wird. Auf diese Weise wird lokal gearbeitet, so dass auch Artefaktmuster mittlerer und niedrigerer Frequenzen auf verbesserte Weise erhalten bleiben, die durch die geänderten Aufnahmebedingungen in der jeweiligen Aufnahmegeometrie ausgelöst werden. Ein solcher auf Variation basierender, lokaler Glättungsalgorithmus kann mathematisch als ein Minimierungsproblem beschrieben werden, in Formeln als
    Figure 00070001
  • Darin bezeichnet U * / i das gewünschte geglättete Subtraktionsbild, Ui das Ausgangs-Subtraktionsbild und die Funktion ψ ein Regularisierungs-Funktional, das geeignet gewählt werden muss. λ ist ein Gewichtungsfaktor, der eine Feineinstellung ermöglicht und M der Minimierungsraum, also der Raum, aus dem die Lösungen gewählt werden.
  • Während es nun grundsätzlich möglich ist, eine übliche TV-Glättung (total variation smoothing) vorzunehmen, wird es erfindungsgemäß bevorzugt, dass eine Sobolev-Regularisierung verwendet wird, bei der der Sobolev-Raum W1,p als Lösungsraum für die Minimierung betrachtet wird. Die erwähnte Norm wird dabei als die dem Sobolev-Raum zugeordnete Norm gewählt. p wird mit besonderem Vorteil größer als 1 gewählt, aber möglichst nahe an 1, beispielsweise als p = 1.1, da sich hier gezeigt hat, dass auf diese Weise eine Glättung der Subtraktionsbilder bei hinreichender Aufrechterhaltung der Information über die Artefaktmuster auf hervorragende Weise möglich ist. Eine andere gute Wahl bezieht sich auf den Raum W1,1 mit Radon-Maßen.
  • Wie bereits erwähnt, ist ein weiteres bestehendes Problem bei heutigen Verfahren, insbesondere jenen, die aus mehreren Kalibrationsbildern gemittelte Korrekturbilder für jede Aufnahmegeometrie vorhalten, der hohe Speicherplatzbedarf. Mithin ist es wünschenswert, die Korrekturbilder (zumindest teilweise) komprimieren zu können, während die Artefaktmuster erhalten werden können. Nachdem erfindungsgemäß eine Aufteilung in die hochfrequenten Artefaktmuster in dem Durchschnittsbild und die Subtraktionsbilder gegeben ist, kann mit besonderem Vorteil vorgesehen sein, dass die Subtraktionsbilder vor dem Speichern mit einem Kompressionsalgorithmus komprimiert werden. Dabei muss darauf geachtet werden, dass ein eventuell auftretender Informationsverlust sich nicht auf die in den Subtraktionsbildern verbliebenen Artefaktmuster mittlerer oder höherer Ortsfrequenzen bezieht, so dass mit besonderem Vorteil vorgesehen sein kann, dass ein Wavelet-Kompressionsalgorithmus, insbesondere unter Verwendung eines CDF9/7- oder eines CDF5/3-Wavelets, verwendet wird. Insbesondere kann dabei ein Wavelet-Kompressionsverfahren mit weicher Schwellwertbetrachtung und variablem Schwellwert herangezogen werden. Eine weiche Schwellwertbetrachtung (soft thresholding) bedeutet letztlich, dass bei Erreichen des Schwellwerts nicht ein Sprung in der Kompressionsfunktion vorliegt, sondern der lineare Verlauf bei 0 beginnend unmittelbar nach dem Schwellwert beginnt, da hiermit gute Ergebnisse erreicht werden können. Ein variabler Schwellwert ist dahingehend zu verstehen, dass der Schwellwert selbst für Detailräume abhängig von den Ortsfrequenzen gewählt wird, mithin selbst variabel ist. Dieses Konzept wird im Englischen zusammenfassend auch als „Bayesian soft thresholding“ bezeichnet. Besonders geeignete Wavelets wurden bereits beschrieben, wobei das CDF5/3-Wavelet eine verlustfreie Kompression erlaubt, bei dem CDF9/7-Wavelet eine verlustbehaftete Kompression möglich ist.
  • Bei geeigneter Auswahl des Kompressionsverfahrens erlaubt die vorliegende Erfindung beispielsweise eine Reduzierung des benötigten Speicherplatzes auf etwa 7 %. Es ergibt sich also in jedem Fall eine deutliche Einsparung.
  • Neben dem Verfahren betrifft die vorliegende Erfindung auch eine Röntgeneinrichtung mit einem C-Bogen, an dem sich gegenüberliegend ein Röntgenstrahler und ein Röntgendetektor angeordnet sind, aufweisend ein fest mit dem Röntgendetektor verbundenes Streustrahlenraster und eine zur Durchführung des erfindungsgemäßen Verfahrens ausgebildete Steuereinrichtung. Sämtliche Ausführungen bezüglich des erfindungsgemäßen Verfahrens lassen sich analog auf die erfindungsgemäße Röntgeneinrichtung übertragen, mit welcher mithin ebenso die Vorteile der vorliegenden Erfindung erhalten werden können. Die Steuereinrichtung führt automatisch das erfindungsgemäße Verfahren aus, wofür beispielsweise geeignete Hardware- und/oder Softwarekomponenten vorgesehen sein können.
  • Weitere Vorteile und Einzelheiten der vorliegenden Erfindung ergeben sich aus den im Folgenden beschriebenen Ausführungsbeispielen sowie anhand der Zeichnung. Dabei zeigen:
  • 1 einen Ablaufplan des erfindungsgemäßen Verfahrens, und
  • 2 eine Prinzipskizze einer erfindungsgemäßen Röntgeneinrichtung
  • Im erfindungsgemäßen Verfahren, mithin auch im dargestellten Ausführungsbeispiel, geht es um eine Korrektur von aufgrund eines fest mit dem Röntgendetektor verbundenen Streustrahlenrasters auftretenden Artefakten, konkret also um eine Korrektur der Projektionsbilder, die für eine dreidimensionale Rekonstruktion genutzt werden, derart, dass Artefakte im dreidimensionalen Datensatz reduziert werden. Solche Artefakte, wie sie im hier dargestellten Ausführungsbeispiel betrachtet werden, haben letztlich verschiedene Ursachen. Neben den Artefaktmustern, die sich bereits in den Projektionsbildern durch Abschattungen durch ein Gitter des Streustrahlenrasters zwangsläufig ergeben, können zudem auch Identitätsschwankungen auftreten, wenn ein nicht auf die korrekte Entfernung zwischen dem Röntgenstrahler und dem Röntgendetektor fokussiertes Streustrahlenraster eingesetzt wird bzw. die Entfernung zwischen dem Röntgendetektor und dem Röntgenstrahler nicht so eingestellt wird, dass sich der Fokussierungsabstand ergibt. Diese auf die Abschattungen bezogenen Bildeffekte sind unabhängig von der konkreten Aufnahmegeometrie, wobei allerdings auch Effekte auftreten können, die von der Aufnahmegeometrie abhängig sind, beispielsweise von mechanischen Ungenauigkeiten der Röntgeneinrichtung, insbesondere des C-Bogens, ausgelöste Effekte. Diese beruhen üblicherweise auf einer relativen Verschiebung von dem Fokus des Röntgenstrahlers und dem Röntgendetektor mit dem Streustrahlenraster, beispielsweise aufgrund von Effekten der Schwerkraft. Durch eine Verschiebung des Fokus relativ zum Detektor tritt auch eine Verschiebung der oben genannten nicht homogenen Intensitätsverteilung auf, so dass bei deren Korrektur Kantenartefakte und dergleichen auftreten können.
  • Das hier dargestellte Ausführungsbeispiel des erfindungsgemäßen Verfahrens bezieht sich konkret auf eine Möglichkeit, mit geringem Messaufwand geeignete Korrekturbilder für die verschiedenen Aufnahmegeometrien zu erzeugen, die noch immer alle relevanten Artefaktmuster enthalten und effizient speicherbar sind.
  • Hierzu werden in einem Schritt 1, vgl. 1, zunächst Kalibrationsbilder 2 aufgenommen, wobei für jede hier betrachtete Aufnahmegeometrie nur ein einziges Kalibrationsbild 2 erforderlich ist. Beispielsweise kann dabei eine homogene Beleuchtung durch eine Kupferplatte sichergestellt werden.
  • In einem Schritt 3 werden nun alle diese zweidimensionalen Kalibrationsbilder 2 gemittelt, so dass ein einziges Durchschnittsbild 4 erhalten wird. Das Durchschnittsbild 4 enthält insbesondere die mit hohen Frequenzen im Ortsraum verbundenen Informationen zu den durch die Abschattung durch das Streustrahlenraster entstandenen Artefaktmustern, nachdem diese unabhängig von der Aufnahmegeometrie immer entstehen, da das Streustrahlenraster ortsfest mit dem Röntgendetektor verbunden ist.
  • In einem Schritt 5 werden nun Subtraktionsbilder 6 gebildet, indem von den Kalibrationsbildern 2 jeweils das Durchschnittsbild 4 abgezogen wird.
  • Die Subtraktionsbilder 6 werden in einem Schritt 7 nun geglättet, wobei ein lokales, auf Variation basierendes Glättungsverfahren verwendet wird. Konkret wird im hier dargestellten Ausführungsbeispiel dabei eine Sobolev-Regularisierung verwendet, wobei als Lösungsraum der Sobolev-Raum W1,p verwendet wird, hier beispielhaft mit p = 1.1, wobei selbstverständlich auch die entsprechende Norm verwendet wird. So ist es möglich, die Subtraktionsbilder 6 zu entrauschen, wobei gleichzeitig Artefaktmuster bzw. Effekte aufgrund des Streustrahlenrasters, beispielsweise eine verschobene Intensitätsverteilung, die bei mittleren und/oder niedrigen Ortsfrequenzen liegen, erhalten werden. Man erhält also als Ergebnis des Schrittes 7 entrauschte, also weniger verrauschte, Subtraktionsbilder 6.
  • Diese entrauschten Subtraktionsbilder 6 werden nun in einem Schritt 8 komprimiert, wobei eine wavelet-Kompression mit weicher Schwellwertbetrachtung (soft thresholding) und variablem Schwellwert eingesetzt wird. Dabei wird hier das CDF9/7-wavelet verwendet.
  • Die geglätteten, komprimierten Subtraktionsbilder 6 und das Durchschnittsbild 4 werden als ein Korrekturdatensatz 9 in einer Speichereinrichtung zur späteren Verwendung bei der Korrektur von Projektionsbildern aufbewahrt. Die Speichereinrichtung kann dabei Teil einer das Verfahren durchführenden Steuereinrichtung der Röntgeneinrichtung sein.
  • Zur Korrektur von bei unterschiedlichen Aufnahmegeometrien aufgenommenen Projektionsbildern 10, aus denen ein dreidimensionaler Bilddatensatz 11 rekonstruiert werden soll, werden nun in einem Schritt 12 zunächst Korrekturbilder ermittelt, indem die geglätteten, komprimierten Subtraktionsbilder 6 dekomprimiert werden und dann das Durchschnittsbild und das jeweilige geglättete, dekomprimierte Subtraktionsbild 6 für jede Aufnahmegeometrie, die benötigt wird, aufsummiert werden. Ebenso in Schritt 12 erfolgt die Korrektur der Projektionsbilder 10 unter Verwendung der so erhaltenen Korrekturbilder.
  • Unter Verwendung der korrigierten Projektionsbilder 10 erfolgt dann in einem Schritt 13 die Rekonstruktion des dreidimensionalen Bilddatensatzes 11, wobei beispielsweise ein Verfahren der gefilterten Rückprojektion oder der iterativen Rekonstruktion verwendet werden kann.
  • 2 zeigt eine Prinzipskizze einer erfindungsgemäßen Röntgeneinrichtung 14. Diese umfasst einen C-Bogen 15, an dem sich gegenüberliegend ein Röntgenstrahler 16 und ein Röntgendetektor 17 angeordnet sind. Der C-Bogen 15 kann dabei um mindestens eine Achse um eine Patientenliege 18, auf der ein hier nicht näher dargestellter Patient angeordnet werden kann, geschwenkt werden, um Röntgenbilder bzw. Projektionsbilder verschiedener Aufnahmegeometrien, mithin verschiedener Projektionsrichtungen, aufnehmen zu können. Um Streustrahlen effektiv zu mindern, ist auf dem Röntgendetektor 17 ein bezüglich diesem ortsfestes Streustrahlenraster 19 montiert.
  • Das Streustrahlenraster 19 ist auf eine bestimmte Entfernung zwischen dem Röntgendetektor 17 und dem Röntgenstrahler 16 hin fokussiert, wobei diese Entfernung ideal gewählt ist, um zweidimensionale Radiographie-Röntgenbilder, also zweidimensionale diagnostische Durchleuchtungsbilder, aufnehmen zu können. Sollen Projektionsbilder 10, aus denen ein dreidimensionaler Bilddatensatz 11 rekonstruiert werden soll, aufgenommen werden, wird jedoch meist die maximal mögliche Entfernung zwischen dem Röntgendetektor 17 und dem Röntgenstrahler 16 gewählt, so dass das Streustrahlenraster 19 für diese Entfernung nicht fokussiert ist.
  • Ferner kommt hinzu, dass bei der Verwendung des C-Bogens 15 mechanische Ungenauigkeiten auftreten, beispielsweise aufgrund dessen nicht komplett starrer Ausbildung und der unterschiedlich wirkenden Schwerkraft auf den Röntgenstrahler 16 und den Röntgendetektor 17 mit dem Streustrahlenraster 19. Dies führt zu einer leichten Verschiebung des Fokus des Röntgenstrahlers 16 im Bezug auf den Röntgendetektor 17 und das Streustrahlenraster 19, was sich auch in den Projektionsbildern als Effekt des Streustrahlenrasters 19 bemerkbar macht.
  • Die Röntgeneinrichtung 14 umfasst ferner zu ihrer Steuerung eine Steuereinrichtung 20, die zur Durchführung des erfindungsgemäßen Verfahrens ausgebildet ist. Dazu weist die Steuereinrichtung 20 insbesondere eine Speichereinrichtung 21 zum Speichern des Korrektur-Datensatzes 9 mit den geglätteten komprimierten Subtraktionsbildern 6 und dem Durchschnittsbild 4 auf.
  • Obwohl die Erfindung im Detail durch das bevorzugte Ausführungsbeispiel näher illustriert und beschrieben wurde, so ist die Erfindung nicht durch die offenbarten Beispiele eingeschränkt und andere Variationen können vom Fachmann hieraus abgeleitet werden, ohne den Schutzumfang der Erfindung zu verlassen.

Claims (7)

  1. Verfahren zur Korrektur von aufgrund eines fest mit einem Röntgendetektor (17) verbundenen Streustrahlenrasters (19) auftretenden Bildartefakten bei der Aufnahme von dreidimensionalen Bilddatensätzen (11) mit einer Röntgeneinrichtung (14) mit einem C-Bogen (15), an dem sich gegenüberliegend ein Röntgenstrahler (16) und der Röntgendetektor (17) angeordnet sind, wobei aus unterschiedlichen Aufnahmegeometrien zweidimensionale Projektionsbilder (10) eines Objekts aufgenommen und zur Rekonstruktion des dreidimensionalen Bilddatensatzes (11) verwendet werden, umfassend folgenden Schritte: – Aufnahme jeweils eines Kalibrationsbildes (2) für jede Aufnahmegeometrie, – Ermittlung eines Durchschnittsbildes (4) aus allen aufgenommenen Kalibrationsbildern (2), – Ermittlung von Subtraktionsbildern (6) für alle Aufnahmegeometrien durch Subtraktion des Durchschnittsbildes (4) von den entsprechenden Kalibrationsbildern (2), – Glättung der Subtraktionsbilder (6) zur Entrauschung der Subtraktionsbilder (6), und – Speichern der Subtraktionsbilder (6) und des Durchschnittsbildes (4), wobei zur Korrektur von aufgenommenen Projektionsbildern (10) für jede Aufnahmegeometrie ein Korrekturbild durch Addition des der Aufnahmegeometrie zugeordneten Subtraktionsbildes (6) und des Durchschnittsbildes (4) ermittelt und das Korrekturbild zur Korrektur des Projektionsbildes (10) verwendet wird.
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das Streustrahlenraster (19) für eine andere Entfernung zwischen dem Röntgenstrahler (16) und dem Röntgendetektor (17) als die bei der Aufnahme der Projektionsbilder (10) genutzte Entfernung fokussiert ist und/oder unter dem Einfluss der Schwerkraft eine Verschiebung des Röntgenstrahlers (16) gegen den Röntgendetektor (17) und mithin das Streustrahlenraster (19) bei unterschiedlichen Positionen des C-Bogens (15) auftritt.
  3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass ein auf Variation basierender, lokaler Glättungsalgorithmus zur Glättung der Subtraktionsbilder (6) eingesetzt wird, wobei eine Minimierung der Summe einer Norm und eines insbesondere mit einem Gewichtungsfaktor eingehenden Regularisierungs-Funktionals vorgenommen wird.
  4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass eine Sobolev-Regularisierung verwendet wird, bei der der Sobolev-Raum W1,p als Lösungsraum für die Minimierung betrachtet wird.
  5. Verfahren nach einem der vorangehenden Ansprüche, dadurch gekennzeichnet, dass die Subtraktionsbilder (6) vor dem Speichern mit einem Kompressionsalgorithmus komprimiert werden.
  6. Verfahren nach Anspruch 5, dadurch gekennzeichnet, dass ein Wavelet-Kompressionsalgorithmus, insbesondere unter Verwendung eines CDF9/7- oder eines CDF5/3-Wavelets, verwendet wird.
  7. Röntgeneinrichtung (14) mit einem C-Bogen (15), an dem sich gegenüberliegend ein Röntgenstrahler (16) und ein Röntgendetektor (17) angeordnet sind, aufweisend ein fest mit dem Röntgendetektor (17) verbundenes Streustrahlenraster (19) und eine zur Durchführung eines Verfahrens nach einem der vorangehenden Ansprüche ausgebildete Steuereinrichtung (20).
DE102011088265.0A 2011-12-12 2011-12-12 Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten Expired - Fee Related DE102011088265B4 (de)

Priority Applications (2)

Application Number Priority Date Filing Date Title
DE102011088265.0A DE102011088265B4 (de) 2011-12-12 2011-12-12 Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten
US13/710,727 US8768037B2 (en) 2011-12-12 2012-12-11 Method for correcting image artifacts occuring as a result of an anti scatter grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102011088265.0A DE102011088265B4 (de) 2011-12-12 2011-12-12 Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten

Publications (2)

Publication Number Publication Date
DE102011088265A1 true DE102011088265A1 (de) 2013-06-13
DE102011088265B4 DE102011088265B4 (de) 2017-10-19

Family

ID=48464418

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102011088265.0A Expired - Fee Related DE102011088265B4 (de) 2011-12-12 2011-12-12 Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten

Country Status (2)

Country Link
US (1) US8768037B2 (de)
DE (1) DE102011088265B4 (de)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102014206720A1 (de) * 2014-04-08 2015-10-08 Siemens Aktiengesellschaft Rauschreduktion in Tomogrammen
DE102016206559B3 (de) * 2016-04-19 2017-06-08 Siemens Healthcare Gmbh Verfahren zur Korrektur eines Röntgenbilds auf Effekte eines Streustrahlenrasters, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105335936B (zh) * 2014-06-19 2018-07-31 上海西门子医疗器械有限公司 X光机、x光图像中静态滤线栅伪影的去除方法及装置
US10755395B2 (en) * 2015-11-27 2020-08-25 Canon Medical Systems Corporation Dynamic image denoising using a sparse representation
CN108065944B (zh) * 2016-11-14 2021-10-22 上海西门子医疗器械有限公司 平板探测器与滤线栅的适配方法和设备及医疗成像系统
EP3427664A1 (de) 2017-07-13 2019-01-16 Koninklijke Philips N.V. Vorrichtung zur streuungskorrektur in einem röntgenbild und verfahren zur streuungskorrektur in einem röntgenbild
US10803555B2 (en) 2017-08-31 2020-10-13 Shanghai United Imaging Healthcare Co., Ltd. System and method for determining a trained neural network model for scattering correction
JP7233850B2 (ja) * 2018-04-25 2023-03-07 キヤノンメディカルシステムズ株式会社 医用画像処理装置、x線診断装置及び医用画像処理プログラム
CN111553960B (zh) * 2020-04-24 2023-04-07 重庆大学 一种基于投影均值图像的环状伪影快速校正方法
CN112288659B (zh) * 2020-12-18 2021-03-23 南京佗道医疗科技有限公司 一种透视图像矫正方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1081947A2 (de) * 1999-08-31 2001-03-07 General Electric Company Korrekturverfahren und Vorrichtung für digitale Röntgenbildaufnahmen
US6285781B1 (en) * 1997-03-27 2001-09-04 Canon Kabushiki Kaisha Radiographic apparatus and method

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010134295A1 (ja) * 2009-05-22 2010-11-25 株式会社島津製作所 同期型グリッドの箔影除去方法およびそれを用いた放射線撮影装置
US8363919B2 (en) * 2009-11-25 2013-01-29 Imaging Sciences International Llc Marker identification and processing in x-ray images
US8781243B2 (en) * 2011-01-07 2014-07-15 Wisconsin Alumni Research Foundation Method for constrained reconstruction of high signal-to-noise ratio images

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6285781B1 (en) * 1997-03-27 2001-09-04 Canon Kabushiki Kaisha Radiographic apparatus and method
EP1081947A2 (de) * 1999-08-31 2001-03-07 General Electric Company Korrekturverfahren und Vorrichtung für digitale Röntgenbildaufnahmen

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102014206720A1 (de) * 2014-04-08 2015-10-08 Siemens Aktiengesellschaft Rauschreduktion in Tomogrammen
US10282872B2 (en) 2014-04-08 2019-05-07 Siemens Healthcare Gmbh Noise reduction in tomograms
DE102016206559B3 (de) * 2016-04-19 2017-06-08 Siemens Healthcare Gmbh Verfahren zur Korrektur eines Röntgenbilds auf Effekte eines Streustrahlenrasters, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
US10159452B2 (en) 2016-04-19 2018-12-25 Siemens Healthcare Gmbh Correction of an x-ray image for effects of an anti-scatter grid

Also Published As

Publication number Publication date
US20130148786A1 (en) 2013-06-13
US8768037B2 (en) 2014-07-01
DE102011088265B4 (de) 2017-10-19

Similar Documents

Publication Publication Date Title
DE102011088265B4 (de) Verfahren zur Korrektur von aufgrund eines Streustrahlenrasters auftretenden Bildartefakten
DE102012204019B4 (de) Verfahren zur Reduzierung von Bewegungsartefakten
EP3236416B1 (de) Verfahren zur korrektur eines röntgenbilds auf effekte eines streustrahlenrasters, röntgeneinrichtung, computerprogramm und elektronisch lesbarer datenträger
WO2005124683A2 (de) Vorrichtung und verfahren für die streustrahlungskorrektur in der projektionsradiographie, insbesondere der mammographie
DE102012206714A1 (de) Verfahren, Recheneinheit, CT-System und C-Bogen-System zur Reduktion von Metallartefakten in CT-Bilddatensätzen
DE102009032059A1 (de) Sinogrammbearbeitung für die Metallartefaktreduktion in der Computertomographie
DE102006005803A1 (de) Verfahren zur Rauschreduktion in bildgebenden Verfahren
DE102010019632A1 (de) Verfahren zur Aufnahme und Rekonstruktion eines dreidimensionalen Bilddatensatzes und Röntgeneinrichtung
DE102005012654A1 (de) Verfahren und Computertomographie-System zur Erstellung tomographischer Aufnahmen eines Objektes
DE102007013570A1 (de) Verfahren zur Rauschverminderung in digitalen Bildern mit lokal unterschiedlichem und gerichtetem Rauschen
DE102011005715A1 (de) Verfahren zum Gewinnen eines von Spuren eines Metallobjektes befreiten 3D-Bilddatensatzes
DE102011006188A1 (de) Verfahren und Computertomographie-System zur Erstellung tomographischer Bilddarstellung mit mindestens zwei Strahler-Detektor-Systemen
DE102007061934A1 (de) Verfahren zur Qualitätssteigerung von computertomographischen Aufnahmeserien durch Projektionsdatenverarbeitung und CT-System mit Recheneinheit
DE102012212774A1 (de) Verfahren zur Korrektur von Metallartefakten und Röntgeneinrichtung
DE102006014629A1 (de) Verfahren zur Korrektur von Trunkierungsartefakten
DE102010026374A1 (de) Verfahren zur Rekonstruktion eines dreidimensionalen Bilddatensatzes und Röntgeneinrichtung
DE102013214689B4 (de) Verfahren zur Reduzierung von Artefakten in Bilddatensätzen und Recheneinrichtung
DE102007002417B4 (de) Verfahren zur Bestimmung einer Position für wenigstens eine halbtransparente Blende und zugehörige Einrichtung
DE102009019840A1 (de) Kontrastverstärkung von CT-Bildern mittels eines Multibandfilters
DE102013215807A1 (de) Verfahren zur Spiralaufnahme mit veränderlicher Tischgeschwindigkeit bei konstantem Pitch und Computertomographiegerät zur Durchführung eines solchen Verfahrens
DE102011075912A1 (de) Verfahren zum Bereitstellen eines 3D-Bilddatensatzes mit unterdrückten Bildartefakten und Computertomograph
DE102015224331B4 (de) Verfahren zur Helligkeitsanpassung eines Röntgenbildes, Röntgeneinrichtung, Computerprogramm und Datenträger
DE102019202878A1 (de) Verfahren zur Rauschreduzierung in Computertomographiedaten, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
DE102008048716B4 (de) Verfahren zur Aufhärtungskorrektur
DE102005003227A1 (de) Verfahren zur Reduzierung von Cupping-Artefakten in ConeBeam-CT-Bilddatensätzen

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHCARE GMBH, DE

Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE

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