DE102012208764A1 - Method for reconstructing computed tomography image data sets of intensity measurements of flat panel detector, involves outputting and storing reconstructed tomographic representation of object under examination - Google Patents
Method for reconstructing computed tomography image data sets of intensity measurements of flat panel detector, involves outputting and storing reconstructed tomographic representation of object under examination Download PDFInfo
- Publication number
- DE102012208764A1 DE102012208764A1 DE201210208764 DE102012208764A DE102012208764A1 DE 102012208764 A1 DE102012208764 A1 DE 102012208764A1 DE 201210208764 DE201210208764 DE 201210208764 DE 102012208764 A DE102012208764 A DE 102012208764A DE 102012208764 A1 DE102012208764 A1 DE 102012208764A1
- Authority
- DE
- Germany
- Prior art keywords
- phantom
- values
- attenuation
- determined
- scanning
- 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.)
- Withdrawn
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/30—Accessories, mechanical or electrical features
- G01N2223/303—Accessories, mechanical or electrical features calibrating, standardising
-
- 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/401—Imaging image processing
Abstract
Description
Die Erfindung betrifft ein Verfahren und Röntgentomographiesystem zur Rekonstruktion von CT-Bilddatensätzen aus Intensitätsmessungen eines Flachdetektors, wobei ein zumindest bezüglich seiner räumlichen Schwächungskoeffizienten bekanntes Phantom durch ein Strahler-Flachdetektor-System aus einer Vielzahl von Projektionswinkeln abgetastet wird, und durch einen Vergleich von theoretisch berechneter Schwächungswerte durchlaufender Röntgenstrahlen und entsprechender gemessener Intensitätswerte für diese Röntgenstrahlen eine Zuordnung zwischen den Intensitätswerten und den theoretischen Schwächungswerten bestimmt wird, die bei einer späteren Abtastung eines Untersuchungsobjektes – unter Verwendung gleicher Geräteparameter – dazu verwendet wird, die gemessenen Intensitätswerte den richtigen Schwächungswerten zuzuordnen und eine Rekonstruktion einer tomographischen Darstellung mit räumlichen Schwächungskoeffizienten des Untersuchungsobjektes auszuführen. The invention relates to a method and X-ray tomography system for reconstructing CT image data sets from intensity measurements of a flat detector, wherein a phantom known at least with respect to its spatial attenuation coefficients is scanned by a radiator flat-detector system from a multiplicity of projection angles, and by a comparison of theoretically calculated attenuation values continuous X-rays and corresponding measured intensity values for these X-rays, an association between the intensity values and the theoretical attenuation values is determined, which is used in a later scan of an examination object - using the same device parameters - to assign the measured intensity values to the correct attenuation values and a reconstruction of a tomographic Representation with spatial attenuation coefficients of the examination object to execute.
Es ist allgemein bekannt, dass die Detektorsignale nicht unmittelbar zur Rekonstruktion artefaktfreier tomographischer Darstellungen verwendet werden können. Hierzu werden Schwächungswerte von Röntgenstrahlen, die ein Untersuchungsobjekt durchdringen, benötigt. Ein monochromer Strahl der Anfangsintensität wird beim Durchdringen eines Objektes gemäß der Formel I = I0exp(–∫µ(x)dx) geschwächt, wobei I der Strahlungsintensität nach der durchdrungenen Schicht, I0 der Anfangsintensität, µ(x) dem linearen Schwächungskoeffizienten an der Stelle x und d der Dicke der durchdrungenen Materie entspricht. Als logarithmierter Schwächungswert wird das Verhältnis ln(I/I0) bezeichnet, das dem Wert des Linienintegrals ∫µ(x)dx entspricht, wobei in der CT die Funktion µ(x) als tomographisches Graustufenbild dargestellt wird. It is well known that the detector signals can not be used directly to reconstruct artifact-free tomographic representations. For this purpose, attenuation values of X-rays penetrating an examination object are required. A monochromatic beam of the initial intensity is weakened when penetrating an object according to the formula I = I 0 exp (-∫μ (x) dx), where I is the radiation intensity after the penetrated layer, I 0 the initial intensity, μ (x) the linear attenuation coefficient at the point x and d corresponds to the thickness of the permeated matter. The logarithmic attenuation value is the ratio ln (I / I 0 ) which corresponds to the value of the line integral ∫μ (x) dx, the function μ (x) being represented in the CT as a tomographic greyscale image.
Wesentlich zur Rekonstruktion ist also die Kenntnis des Schwächungswertes, der durch eine ideale monochromatische Absorption im Untersuchungsobjekt entsteht. Essential to the reconstruction is thus the knowledge of the attenuation value, which results from an ideal monochromatic absorption in the examination object.
Es ist weiterhin bekannt, dass zur Übertragung der gemessenen Intensitäten einer multispektralen Röntgenstrahlung an einem Detektor auf einen entsprechenden monochromatischen Schwächungswert – basierend auf Modellbetrachtungen – Übertragungsfunktionen verwendet werden, die die gemessenen multispektralen Intensitäten auf einen monochromatischen Schwächungswert umrechnen. It is further known that for the transmission of the measured intensities of a multispectral X-radiation at a detector to a corresponding monochromatic attenuation value - based on model considerations - transfer functions are used which convert the measured multispectral intensities to a monochromatic attenuation value.
Es hat sich allerdings gezeigt, dass diese so ermittelten Übertragungsfunktionen in der Praxis nicht immer zu artefaktfreien CT-Darstellungen führen. However, it has been shown that these transfer functions determined in practice do not always lead to artifact-free CT images.
Es ist daher Aufgabe der Erfindung, ein verbessertes Verfahren zur Übertragung gemessener Intensitätswerte an einem Detektor zu Schwächungswerten zu finden, die für die Rekonstruktion artefaktfreier Darstellungen geeignet sind. It is therefore an object of the invention to find an improved method for transmitting measured intensity values at a detector to attenuation values which are suitable for the reconstruction of artefact-free representations.
Diese Aufgabe wird durch die Merkmale der unabhängigen Patentansprüche gelöst. Vorteilhafte Weiterbildungen der Erfindung sind Gegenstand untergeordneter Ansprüche. This object is solved by the features of the independent claims. Advantageous developments of the invention are the subject of the subordinate claims.
Der Erfinder hat erkannt, dass es möglich ist, auf der Basis von Intensitätsmessungen eines Phantoms durch Vergleich mit theoretisch berechneten Schwächungswerten jeweils sich entsprechender Strahlen unmittelbar Übertragungsfaktoren zwischen Intensitätswerten und Schwächungswerten – unter der Voraussetzung gleicher Geräteparameter bei der Phantomabtastung und der Abtastung des Untersuchungsobjektes – zu finden sind. The inventor has recognized that it is possible, on the basis of intensity measurements of a phantom, to find directly corresponding transmission factors between intensity values and attenuation values by comparison with theoretically calculated attenuation values of the corresponding beams - assuming the same device parameters in the phantom scan and the scan of the examination subject are.
Demgemäß schlägt der Erfinder ein Verfahren zur Rekonstruktion von CT-Bilddatensätzen aus Intensitätsmessungen eines Flachdetektors vor, welches die folgenden Verfahrensschritte aufweist:
- – Abtastung eines Phantoms mit bekannten Schwächungskoeffizienten durch ein Strahler-Flachdetektor-System aus einer Vielzahl von Projektionswinkeln und Aufzeichnung der gemessenen Intensitätswerte für alle Strahlen,
- – Bestimmung der vom Phantom besetzten Positionen,
- – Berechnung theoretischer Schwächungswerte für alle Strahlen aus der Abtastung auf der Basis der nun bestimmten Lage und Kontur, einschließlich der vorbekannten räumlichen Schwächungskoeffizienten des Phantoms,
- – Ermittlung einer Zuordnung, die jedem unter vorgegebenen Geräteparametern gemessenen Intensitätswert eines Strahls den theoretischen Schwächungswert dieses Strahls am Ort der Messung gegenüberstellt,
- – Abtastung eines Untersuchungsobjektes mit dem Strahler-Flachdetektor-System mit den gleichen Geräteparametern aus einer Vielzahl von Projektionswinkeln und Aufzeichnung der gemessenen Intensitätswerte für alle Strahlen,
- – Übertragung der gemessenen Intensitätswerte unter Verwendung der ermittelten Zuordnung auf Schwächungswerte,
- – Rekonstruktion einer Darstellung der räumlichen Schwächungskoeffizienten des Untersuchungsobjektes unter Verwendung der ermittelten Schwächungswerte,
- – Ausgabe oder Speicherung der rekonstruierten tomographischen Darstellung des Untersuchungsobjektes.
- Scanning a phantom with known attenuation coefficients by a radiator flat-detector system from a plurality of projection angles and recording the measured intensity values for all the beams,
- - determination of positions occupied by the phantom,
- Calculation of theoretical attenuation values for all beams from the scan on the basis of the now determined position and contour, including the known spatial attenuation coefficients of the phantom,
- Determination of an assignment which, for each intensity parameter of a beam measured under predetermined device parameters, compares the theoretical attenuation value of this beam at the location of the measurement,
- Scanning an object under examination with the emitter flat panel detector system with the same device parameters from a multiplicity of projection angles and recording the measured intensity values for all the beams,
- Transmission of the measured intensity values using the determined assignment to attenuation values,
- Reconstruction of a representation of the spatial attenuation coefficients of the examination object using the determined attenuation values,
- - Output or storage of the reconstructed tomographic representation of the examination object.
Zur Bestimmung der genauen Lage und der Kontur des Phantoms im Messfeld wird beispielsweise vorgeschlagen, dass die vom Phantom während der Abtastung besetzten Positionen durch die folgenden Verfahrensschritte bestimmt werden:
- – Rekonstruktion einer tomographischen Phantomdarstellung des abgetasteten Phantoms mit den bei der Abtastung gemessenen Intensitätswerten,
- – Segmentierung der tomographischen Phantomdarstellung bezüglich des Vorliegens oder Nicht-Vorliegens des Phantoms am jeweiligen Bildpixel,
- – Ersetzen der Werte der, das Phantom darstellenden Bildwerte durch die vorbekannten Schwächungswerte des Phantoms. Alternativ können die vom Phantom während der Abtastung besetzten Positionen im Messfeld auch bestimmt werden durch die folgenden Verfahrensschritte:
- – Verwendung eines Phantoms mit ausschließlich konvexer Kontur,
- – Bestimmung aller Randstrahlen in mindestens einer Abtastebene, welche das Phantom bei der Abtastung begrenzen,
- – Bestimmung der Kontur und Lage des Phantoms in der mindestens einen Abtastebene aus der von keinem Randstrahl durchdrungenen Fläche.
- Reconstruction of a tomographic phantom representation of the sampled phantom with the intensity values measured during the scan,
- Segmentation of the tomographic phantom representation with respect to the presence or absence of the phantom at the respective image pixel,
- Replacing the values of the image values representing the phantom with the previously known attenuation values of the phantom. Alternatively, the positions occupied by the phantom during the scan in the measurement field can also be determined by the following method steps:
- Use of a phantom with exclusively convex contour,
- Determination of all marginal rays in at least one scanning plane, which limit the phantom during scanning,
- - Determination of the contour and position of the phantom in the at least one scanning plane from the surface penetrated by no edge beam.
Weiterhin kann die Zuordnung zwischen gemessenen Intensitätswerten und theoretisch zu erwartenden Schwächungswerten durch Aufstellen einer LUT, also einer „look up“-Tabelle, erzeugt werden. Hierbei können in der LUT Mittelwerte von gemessenen Intensitätswerten gegenüber theoretischen Schwächungswerten abgelegt werden. Furthermore, the association between measured intensity values and theoretically expected attenuation values can be generated by setting up an LUT, ie a "look up" table. Here, mean values of measured intensity values can be stored in the LUT compared to theoretical attenuation values.
Andererseits besteht auch die Möglichkeit, die Zuordnung zwischen gemessenen Intensitätswerten und theoretisch zu erwartenden Schwächungswerten in Form einer mathematischen Funktion zu bestimmen. Bei dieser Variante kann als mathematische Funktion beispielsweise ein Polynom, eine Exponentialfunktion oder eine Spline-Funktion verwendet werden. On the other hand, it is also possible to determine the association between measured intensity values and theoretically expected attenuation values in the form of a mathematical function. In this variant, for example, a polynomial, an exponential function or a spline function can be used as a mathematical function.
Da es aufgrund der geometrischen Verhältnisse, einer nicht zu verhindernden Statistik der Messwerte und auch bedingt durch das Spektrum der verwendeten Röntgenstrahlung, bei der Abtastung des Phantoms zu Redundanzen und Unsicherheiten bei den Messwerten kommen kann, wird auch vorgeschlagen, zur Zuordnung ein „Least Square Fit“-Verfahren auszuführen, in dem die Parameter der mathematischen Funktion bestimmt werden. Since due to the geometric conditions, an unavoidable statistic of the measured values and also due to the spectrum of the X-ray radiation used, the sampling of the phantom can lead to redundancies and uncertainties in the measured values, it is also proposed to assign a Least Square Fit Method in which the parameters of the mathematical function are determined.
Weiterhin kann mit der zuvor bestimmten mathematischen Funktion wiederum eine eindeutige LUT erstellt werden, die zur Übertragung der beim Scan des Untersuchungsobjektes gemessenen Intensitätswerte auf die zur Rekonstruktion verwendeten Schwächungswerte verwendet wird. Furthermore, with the previously determined mathematical function, in turn, a unique LUT can be created, which is used to transmit the intensity values measured during the scan of the examination subject to the attenuation values used for the reconstruction.
Grundsätzlich basiert das vorliegende Verfahren darauf, dass das Röntgentomographiesystem zur Abtastung des Phantoms und zur Abtastung des Untersuchungsobjektes technisch gleich ist, also dass eine Übereinstimmung aller Geräteparameter vorliegt. Wird die LUT für unterschiedliche Randbedingungen, also unterschiedliche Geräte oder gleiche Geräte mit unterschiedlichen Einstellungen erzeugt, so sollten diese Parameter in der LUT vorliegen beziehungsweise je Parametersatz eine LUT vorliegen. Dabei können als Geräteparameter mindestens ein Parameter der folgenden Liste verwendet werden: Röhrenstrom, Röhrenspannung, Filtermaterial, Filterdicke, Detektortemperatur. In principle, the present method is based on the fact that the X-ray tomography system is technically the same for scanning the phantom and for scanning the examination object, ie that all device parameters match. If the LUT is generated for different boundary conditions, ie different devices or the same devices with different settings, these parameters should be present in the LUT or a LUT should be available for each parameter set. At least one parameter from the following list can be used as device parameters: tube current, tube voltage, filter material, filter thickness, detector temperature.
Neben dem erfindungsgemäßen Verfahren schlägt der Erfinder auch ein Röntgentomographiesystem mit mindestens einem Röntgenstrahler und einem gegenüberliegenden Flachdetektor, geeignet zur tomographischen Abtastung eines Untersuchungsobjektes aus einer Vielzahl von Projektionswinkeln, und einer Steuer- und Recheneinheit zur Steuerung des Röntgentomographiesystems und zur Erzeugung tomographischer Darstellungen des Untersuchungsobjektes, vor, wobei in einem Speicher der Steuer- und Recheneinheit ein Computerprogramm gespeichert ist, welches im Betrieb die Verfahrensschritte des voranstehenden beschriebenen Verfahrens ausführt. In addition to the method according to the invention, the inventor also proposes an X-ray tomography system having at least one X-ray emitter and an opposite flat detector suitable for tomographic scanning of an examination subject from a plurality of projection angles, and a control and processing unit for controlling the X-ray tomography system and for generating tomographic representations of the examination subject , wherein in a memory of the control and computing unit, a computer program is stored, which executes the method steps of the above described method during operation.
Im Folgenden wird die Erfindung mit Hilfe der Figuren näher beschrieben, wobei nur die zum Verständnis der Erfindung notwendigen Merkmale dargestellt sind. Es werden folgende Bezugszeichen und Variablen verwendet: A1: Aufnahme der akquirierten Intensitätsbilder eines Phantoms; A2: Rekonstruktion einer tomographische Abbildung des Phantoms mit einer genau bekannten Kontur; A3: Segmentierung des rekonstruierten Phantoms; A4: Berechnung der theoretisch zu erwartenden Schwächungswerte auf der Basis des segmentierten Phantoms mit bekannten Schwächungskoeffizienten; I/I0: Schwächungswerte; IM: Intensitätswert; P: Phantom, Prg1–Prgn: Computerprogramme; S: Strahlen; S1: Abtastung eines Phantoms mit bekannten Schwächungskoeffizienten; S2: Bestimmung der vom Phantom besetzten Positionen; S3: Berechnung theoretischer Schwächungswerte; S4: Ermittlung einer Zuordnung zwischen gemessenen Intensitätswert und dem theoretischen Schwächungswert; S5: Abtastung eines Untersuchungsobjektes; S6: Übertragung der gemessenen Intensitätswerte auf Schwächungswerte; S7: Rekonstruktion der Darstellung der räumlichen Schwächungskoeffizienten; S8: Ausgabe/Speicherung der rekonstruierten tomographischen Darstellung des Untersuchungsobjektes; µ(x): Schwächungskoeffizient;
Es zeigen im Einzelnen: They show in detail:
Die
Im Schritt S2 werden durch geeignete Maßnahmen die vom Phantom besetzten Positionen innerhalb des Messfeldes bestimmt. Beispielsweise kann dies durch eine physische Vermessung der Lage und Kontur des Phantoms geschehen. Vorteilhafter ist jedoch eine Bestimmung durch andere messtechnische Maßnahmen, wie eine Rekonstruktion mit Segmentierung oder die Bestimmung eines Satzes von Strahlen, die die sich an die Kontur des Phantoms anlegen. In step S2, the positions occupied by the phantom within the measuring field are determined by suitable measures. For example, this can be done by physically measuring the location and contour of the phantom. More advantageous, however, is a determination by other metrological measures, such as a reconstruction with segmentation or the determination of a set of rays that apply to the contour of the phantom.
Es folgt im Schritt S3 die Berechnung theoretischer Schwächungswerte I/I0 für alle Strahlen S, die auch bei der Abtastung des Phantoms verwendet wurden, wobei auf der Basis der zuvor bestimmten Lage und Kontur des Phantoms und dessen bekannter Schwächungskoeffizienten µ(x). Somit liegt für jeden Messstrahl der das Phantom durchdrungen hat und dessen Intensität nach der Schwächung durch das Phantom gemessen wurde ein korrespondierender Schwächungswert vor, der theoretisch beim Durchtritt des Strahls durch das Phantom entstehen sollte. It follows in step S3, the calculation of theoretical attenuation values I / I 0 for all the rays S, which were also used in the scanning of the phantom, based on the previously determined position and contour of the phantom and its known attenuation coefficient μ (x). Thus, for each measuring beam that has penetrated the phantom and whose intensity has been measured after attenuation by the phantom, there is a corresponding attenuation value that should theoretically arise as the beam passes through the phantom.
Im Schritt S4 wird nun auf der Basis der so ermittelten korrespondierenden Werte eine Zuordnung bestimmt, die jedem unter vorgegebenen Geräteparametern gemessenen Intensitätswert IM eines Strahls einen entsprechenden theoretischen Schwächungswert I/I0 dieses Strahls am Ort der Messung gegenüberstellt. Grundsätzlich können diese gegenüberstehenden gemessenen Intensitätswerte und theoretischen Schwächungswerte in einem Diagramm aufgetragen werden, so dass die – z.B. wegen der immer vorliegenden Messstatistik – vom idealen Zusammenhang abweichenden Unsicherheiten über eine Kurvenanpassung gemittelt werden können. Es liegt damit eine eindeutige Zuordnung zwischen den an einem Flachdetektor aufgenommenen Intensitätswerten und den Schwächungswerten vor. Hierbei ist zu bemerken, dass im Rahmen dieser Zuordnung zumindest teilweise auch automatisch Nebeneffekte, wie Streustrahlung und Strahlaufhärtung, korrigiert werden. In step S4, on the basis of the corresponding values determined in this way, an assignment is determined which compares a respective theoretical attenuation value I / I 0 of this beam at the location of the measurement to each intensity value I M of a beam measured under predetermined device parameters. In principle, these opposing measured intensity values and theoretical attenuation values can be plotted in a diagram, so that the uncertainties, which are deviating from the ideal relationship, for example because of the ever present measurement statistics, can be averaged over a curve fit. There is thus a clear association between the intensity values recorded on a flat detector and the attenuation values. It should be noted that, within the scope of this assignment, at least partially automatic side effects, such as scattered radiation and beam hardening, are corrected.
Zur eigentlichen Untersuchung eines Untersuchungsobjektes kann nun im Schritt S5 eine Abtastung dieses Untersuchungsobjektes mit dem Strahler-Flachdetektor-System mit den gleichen Geräteparametern aus einer Vielzahl von Projektionswinkeln und Aufzeichnung der gemessenen Intensitätswerte IM für alle Strahlen erfolgen. Zur Korrektur der Messwerte werden dann im Schritt S6 die gemessenen Intensitätswerte unter Verwendung der ermittelten Zuordnung auf die Schwächungswerte I/I0 übertragen. For the actual examination of an examination subject, scanning of this examination subject with the radiator flat-detector system with the same device parameters from a multiplicity of projection angles and recording of the measured intensity values I M for all the rays can now take place in step S5. In order to correct the measured values, the measured intensity values are then transferred to the attenuation values I / I 0 in step S 6 using the determined assignment.
Im Schritt S7 erfolgt die Rekonstruktion der Darstellung der räumlichen Schwächungskoeffizienten µ(x) des Untersuchungsobjektes unter Verwendung der ermittelten Schwächungswerte I/I0 in an sich bekannter Weise, wobei im Schritt S8 eine Ausgabe oder Speicherung der rekonstruierten tomographischen Darstellung des Untersuchungsobjektes erfolgt. In step S7, the representation of the spatial attenuation coefficients μ (x) of the examination object is reconstructed using the determined attenuation values I / I 0 in a manner known per se, with an output or storage of the reconstructed tomographic representation of the examination subject in step S8.
Die
Die
Die Steuerung des C-Bogen-Systems
Insgesamt wird also mit dieser Erfindung ein Verfahren und ein Röntgentomographiesystem zur Rekonstruktion von CT-Bilddatensätzen aus Intensitätsmessungen eines Flachdetektors vorgestellt, wobei durch Abtastung eines Phantoms mit bekannten Schwächungskoeffizienten durch ein Strahler-Flachdetektor-System eine Vielzahl von Intensitätswerten gemessen und diese mit theoretisch berechneten Schwächungswerten korreliert werden, so dass sich eine eindeutige Zuordnung zwischen Intensitätswerten und Schwächungswerten ergibt und nachfolgend bei der Abtastung eines Untersuchungsobjektes die gemessenen Intensitätswerte in Schwächungswerte transferiert werden, mit denen eine Rekonstruktion des Untersuchungsobjektes stattfindet. Overall, this invention thus provides a method and an X-ray tomography system for reconstructing CT image data sets from intensity measurements of a flat detector, wherein a plurality of intensity values are measured by sampling a phantom having known attenuation coefficients by a radiator flat-detector system and correlating these with theoretically calculated attenuation values so that there is an unambiguous association between intensity values and attenuation values, and subsequently, during the scanning of an examination object, the measured intensity values are transferred into attenuation values with which a reconstruction of the examination subject takes place.
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. Although the invention has been further illustrated and described in detail by the preferred embodiment, the invention is not limited by the disclosed examples, and other variations can be derived therefrom by those skilled in the art without departing from the scope of the invention.
BezugszeichenlisteLIST OF REFERENCE NUMBERS
- A1 A1
- Aufnahme der akquirierten Intensitätsbilder eines PhantomsRecording the acquired intensity images of a phantom
- A2 A2
- Rekonstruktion einer tomographische Abbildung des Phantoms mit einer genau bekannten KonturReconstruction of a tomographic image of the phantom with a well-known contour
- A3 A3
- Segmentierung des rekonstruierten Phantoms Segmentation of the reconstructed phantom
- A4 A4
- Berechnung der theoretisch zu erwartenden Schwächungswerte auf der Basis des segmentierten Phantoms mit bekannten SchwächungskoeffizientenCalculation of the theoretically expected attenuation values on the basis of the segmented phantom with known attenuation coefficients
- I/I0 I / I 0
- Schwächungswerte attenuation values
- IM I M
- Intensitätswert intensity value
- P P
- Phantom phantom
- Prg1–Prgn Prg 1 prr n
- Computerprogramme computer programs
- S S
- Strahlen radiate
- S1 S1
- Abtastung eines Phantoms mit bekannten SchwächungskoeffizientenScanning a phantom with known attenuation coefficients
- S2 S2
- Bestimmung der vom Phantom besetzten PositionenDetermination of positions occupied by the phantom
- S3 S3
- Berechnung theoretischer Schwächungswerte Calculation of theoretical attenuation values
- S4 S4
- Ermittlung einer Zuordnung zwischen gemessenen Intensitätswert und dem theoretischen SchwächungswertDetermining an association between the measured intensity value and the theoretical attenuation value
- S5 S5
- Abtastung eines Untersuchungsobjektes Scanning of an examination object
- S6 S6
- Übertragung der gemessenen Intensitätswerte auf SchwächungswerteTransmission of the measured intensity values to attenuation values
- S7 S7
- Rekonstruktion der Darstellung der räumlichen SchwächungskoeffizientenReconstruction of the representation of the spatial attenuation coefficients
- S8 S8
- Ausgabe/Speicherung der rekonstruierten tomographischen Darstellung des UntersuchungsobjektesOutput / storage of the reconstructed tomographic representation of the examination object
- µ(x)μ (x)
- Schwächungskoeffizient attenuation coefficient
- 11
- C-Bogen-System C-arm system
- 22
- Röntgenstrahler X-ray
- 33
- Flachdetektor FPD
- 44
- Patient patient
- 5 5
- Computer computer
- 6 6
- Gehäuse casing
- 7 7
- C-Bogen C-arm
- 8 8th
- Patientenliege patient support
Claims (11)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE201210208764 DE102012208764A1 (en) | 2012-05-24 | 2012-05-24 | Method for reconstructing computed tomography image data sets of intensity measurements of flat panel detector, involves outputting and storing reconstructed tomographic representation of object under examination |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE201210208764 DE102012208764A1 (en) | 2012-05-24 | 2012-05-24 | Method for reconstructing computed tomography image data sets of intensity measurements of flat panel detector, involves outputting and storing reconstructed tomographic representation of object under examination |
Publications (1)
Publication Number | Publication Date |
---|---|
DE102012208764A1 true DE102012208764A1 (en) | 2013-11-28 |
Family
ID=49547044
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE201210208764 Withdrawn DE102012208764A1 (en) | 2012-05-24 | 2012-05-24 | Method for reconstructing computed tomography image data sets of intensity measurements of flat panel detector, involves outputting and storing reconstructed tomographic representation of object under examination |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE102012208764A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111789624A (en) * | 2020-06-29 | 2020-10-20 | 上海联影医疗科技有限公司 | Defocused radiation measuring method, defocused radiation measuring device, computer equipment and readable storage medium |
CN117770863A (en) * | 2024-02-27 | 2024-03-29 | 赛诺威盛科技(北京)股份有限公司 | Z-direction alignment method and device for CT detector |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5774519A (en) * | 1997-01-30 | 1998-06-30 | Analogic Corporation | Method of and apparatus for calibration of CT scanners |
-
2012
- 2012-05-24 DE DE201210208764 patent/DE102012208764A1/en not_active Withdrawn
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5774519A (en) * | 1997-01-30 | 1998-06-30 | Analogic Corporation | Method of and apparatus for calibration of CT scanners |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111789624A (en) * | 2020-06-29 | 2020-10-20 | 上海联影医疗科技有限公司 | Defocused radiation measuring method, defocused radiation measuring device, computer equipment and readable storage medium |
CN111789624B (en) * | 2020-06-29 | 2023-08-08 | 上海联影医疗科技股份有限公司 | Defocused radiation measurement method, apparatus, computer device and readable storage medium |
CN117770863A (en) * | 2024-02-27 | 2024-03-29 | 赛诺威盛科技(北京)股份有限公司 | Z-direction alignment method and device for CT detector |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102016203257B4 (en) | Generating contrast-enhanced image data based on multi-energy X-ray imaging | |
DE102011004120B4 (en) | Method, image data record processing device, X-ray system and computer program for the correction of image data of an examination object | |
DE102012204980B4 (en) | Method for the reconstruction of CT images with scatter correction, in particular for dual source CT devices | |
DE102013200337B4 (en) | Method, computer tomograph and computer program product for determining intensity values of an X-ray radiation for dose modulation | |
DE102016207437B4 (en) | Spectral-independent determination of calcium deposits in blood vessels | |
DE102011055616A1 (en) | System and method for breast imaging using X-ray computed tomography | |
DE102012217555A1 (en) | Method and Computed Tomography System for Determining Bone Mineral Density Values | |
DE102012105560A1 (en) | Method and system for correcting scattering in X-ray imaging | |
EP1764040A2 (en) | Method for radiologic 3D imaging with reduced artefacts, medical imaging device and method of establishing a treatment plan | |
EP3236416B1 (en) | Method for adjusting an x-ray image for effects of an anti-scatter grid, x-ray device, computer program and electronically readable data carrier | |
DE102015207107A1 (en) | Method for generating a virtual X-ray projection on the basis of an image data set obtained by means of an X-ray image recording device, computer program, data carrier and X-ray image recording device | |
DE102012217569A1 (en) | Automatic determination of a spectral distribution of X-rays of a number of X-ray sources | |
DE10141346A1 (en) | Procedure for the acquisition of measurement data with a computer tomograph | |
DE102016223831B3 (en) | Method for providing quantitative CT image data acquired from a measurement field by means of a dual-source CT device, as well as arithmetic unit, dual-source CT device and computer program product | |
DE102011083854B4 (en) | Time-resolved tomosynthesis imaging | |
DE102011003857B4 (en) | Method for adapting a gray value window, arithmetic unit, X-ray device and data carrier | |
DE102007014829B3 (en) | Method for scattered radiation correction in imaging X-ray devices and X-ray imaging system | |
DE102015206363B3 (en) | Filter arrangement for CT system with multiple X-ray sources | |
DE102012212774A1 (en) | Method for correcting metal artifacts and x-ray equipment | |
DE102018221691A1 (en) | Customized generation of virtual image data based on multi-energy X-ray imaging | |
DE102011005161B4 (en) | Method, image data set reconstruction device, X-ray system and computer program for artifact correction | |
DE60214022T2 (en) | METHOD FOR REDUCING ARTEFACT IN OBJECT IMAGES | |
DE102011079496A1 (en) | Method for generating combined tomographic emission- and transmission representation of breathing and cyclically moving patients with multiple rest phases, involves receiving emission-detection data from emission incidents in patients | |
DE102017213479A1 (en) | Computed tomographic image with different energy threshold sets | |
DE102012208764A1 (en) | Method for reconstructing computed tomography image data sets of intensity measurements of flat panel detector, involves outputting and storing reconstructed tomographic representation of object under examination |
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 |
|
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |