DE102013203410A1 - Method for producing flat or curved image plane associated reconstruction image of object i.e. organ of patient, by image pickup device, involves correcting back projected image data using filter kernel according to plane and trajectory - Google Patents
Method for producing flat or curved image plane associated reconstruction image of object i.e. organ of patient, by image pickup device, involves correcting back projected image data using filter kernel according to plane and trajectory Download PDFInfo
- Publication number
- DE102013203410A1 DE102013203410A1 DE201310203410 DE102013203410A DE102013203410A1 DE 102013203410 A1 DE102013203410 A1 DE 102013203410A1 DE 201310203410 DE201310203410 DE 201310203410 DE 102013203410 A DE102013203410 A DE 102013203410A DE 102013203410 A1 DE102013203410 A1 DE 102013203410A1
- Authority
- DE
- Germany
- Prior art keywords
- image
- recording
- trajectory
- reconstruction
- image plane
- 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
-
- 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]
-
- 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
-
- 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
-
- 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
-
- 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
-
- 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/416—Exact reconstruction
-
- 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]
Abstract
Description
Die Erfindung betrifft ein Verfahren zur Erzeugung eines einer flachen oder gekrümmten Bildebene zugeordneten Rekonstruktionsbildes eines Objekts aus einer Mehrzahl von zweidimensionalen, unter unterschiedlichen Projektionsrichtungen mit einer entlang einer Aufnahmetrajektorie relativ zu dem Objekt bewegten Aufnahmeanordnung einer Bildaufnahmeeinrichtung aufgenommenen zweidimensionalen Projektionsbildern. Daneben betrifft die Erfindung eine Bildaufnahmeeinrichtung. The invention relates to a method for generating a reconstructed image of an object assigned to a flat or curved image plane from a plurality of two-dimensional projection images recorded under different projection directions with a recording arrangement of an image recording device moved along a recording trajectory relative to the object. In addition, the invention relates to an image pickup device.
Die Ermittlung dreidimensionaler Informationen eines Objekts aus zweidimensionalen Projektionsbildern ist im Stand der Technik bereits weithin bekannt. Dabei werden häufig dedizierte medizinische Bildgebungssysteme eingesetzt, um die Darstellung anatomischer Strukturen im Raum zu ermöglichen. Häufig werden aus den zweidimensionalen Projektionsbildern dreidimensionale Bilddatensätze rekonstruiert, wobei es jedoch auch bekannt ist, Bildinformationen, beispielsweise Schwächungswerte, für bestimmte Bildebenen im dreidimensionalen Raum aus zweidimensionalen Projektionsbildern zu rekonstruieren. Ein einer Bildebene zugeordnetes Rekonstruktionsbild gibt also insbesondere Schwächungswerte bzw. Objektdichtewerte in der Bildebene wider. The determination of three-dimensional information of an object from two-dimensional projection images is already widely known in the prior art. Dedicated medical imaging systems are often used to facilitate the visualization of anatomical structures in space. Frequently, three-dimensional image data sets are reconstructed from the two-dimensional projection images, although it is also known to reconstruct image information, for example attenuation values, for certain image planes in three-dimensional space from two-dimensional projection images. A reconstruction image assigned to an image plane therefore reflects in particular attenuation values or object density values in the image plane.
Eine in Bezug auf hinreichende Projektionsdaten artefaktfreie Rekonstruktion dreidimensionaler Bilddatensätze, oft auch als „vollständige Rekonstruktion“ bezeichnet, ist nur dann möglich, wenn die sogenannte Tuy-Bedingung erfüllt ist. Beispielsweise in der Computertomographie (CT) erfordert dies eine Datenaufnahme aus unterschiedlichen Projektionsrichtungen bei einer Rotation um 180° plus den Kegelstrahlwinkel, in der Praxis also beispielsweise 220°. An artifact-free reconstruction of three-dimensional image data sets, often referred to as "complete reconstruction", is possible only if the so-called Tuy condition is met. For example, in computed tomography (CT), this requires a data acquisition from different projection directions with a rotation of 180 ° plus the cone angle, in practice, for example, 220 °.
Insbesondere für die dreidimensionale Bildrekonstruktion in der interventionellen Radiologie, beispielsweise unter Verwendung von Röntgeneinrichtungen mit einem C-Arm, stellt die Datenakquisition über ein solches Winkelintervall eine starke Einschränkung für den Anwender dar. Die Aufnahmeanordnung, in diesem Fall umfassend eine Strahlungsquelle und einen Röntgendetektor, muss eine dedizierte Drehbewegung über etwa 220° oder dergleichen ausführen. Insbesondere im Hinblick auf die Bedienbarkeit sind dabei starke Beeinträchtigungen des „medical workflow“ gegeben. Particularly for three-dimensional image reconstruction in interventional radiology, for example using X-ray equipment with a C-arm, data acquisition over such an angular interval represents a severe limitation for the user. The acquisition arrangement, in this case comprising a radiation source and an X-ray detector, must perform a dedicated rotational movement over about 220 ° or the like. In particular with regard to the operability, there are severe impairments of the "medical workflow".
Bildaufnahmeeinrichtungen, die aus zweidimensionalen Projektionsbildern Informationen im dreidimensionalen Raum wiedergebende Rekonstruktionsbilder erzeugen, werden häufig auch als tomographische Bildaufnahmesysteme bezeichnet. Imaging devices that generate reconstructive images representing two-dimensional projection images in three-dimensional space are also often referred to as tomographic image acquisition systems.
Um den Workflow in Teilbereichen zu verbessern, wurden Sonderformen der tomographischen Bildgebung vorgeschlagen, die die Rekonstruktion dann artefaktbehafteter dreidimensionaler Rekonstruktionsbilderstapel oder Rekonstruktionsbilder erlauben. Ein Beispiel hierfür ist die Tomosynthese, bei der die Akquisition von Projektionsbildern über einen kleineren Winkelbereich durchgeführt wird. Häufig bleibt dabei der Detektor ortsfest und nur die Strahlungsquelle führt die notwendige Drehbewegung aus. Hieraus werden dann einer Bildebene zugeordnete, insbesondere also Schwächungswerte oder Objektdichten in der Bildebene wiedergebende, Rekonstruktionsbilder erzeugt und wiedergegeben. In order to improve the workflow in subareas, special forms of tomographic imaging have been proposed, which allow the reconstruction of then artifact-laden three-dimensional reconstruction image stacks or reconstruction images. An example of this is tomosynthesis, in which the acquisition of projection images is carried out over a smaller angular range. Frequently, the detector remains stationary and only the radiation source performs the necessary rotational movement. This then generates and reproduces reconstruction images assigned to an image plane, in particular thus attenuation values or object densities in the image plane.
Der Erfindung liegt daher die Aufgabe zugrunde, eine Möglichkeit zur Ermittlung möglichst artefaktfreier Rekonstruktionsbilder bei einer Nutzung von Aufnahmetrajektorien, die nicht die Tuy-Bedingung erfüllen, anzugeben. The invention is therefore based on the object of specifying a possibility for determining artifact-free reconstruction images when using acquisition trajectories that do not fulfill the Tuy condition.
Zur Lösung dieser Aufgabe ist bei einem Verfahren der eingangs genannten Art erfindungsgemäß vorgesehen, dass eine nicht die Tuy-Bedingung zur vollständigen Rekonstruktion erfüllende Aufnahmetrajektorie verwendet wird, wobei das Rekonstruktionsbild durch Rückprojektion für innerhalb der Bildebene liegende Punkte und Korrektur der rückprojizierten Bilddaten unter Verwendung eines von der Bildebene und der Aufnahmetrajektorie abhängigen Tiefpassfilterkerns ermittelt wird. To solve this problem, it is provided according to the invention in a method of the type mentioned above that a reception trajectory not satisfying the Tuy condition for complete reconstruction is used, the reconstruction image being backprojected for points lying within the image plane and correction of the backprojected image data using a the image plane and the receiving trajectory dependent low-pass filter kernel is determined.
Es ist noch immer eine Relativbewegung zwischen Objekt und Aufnahmeanordnung vorgesehen, mithin eine bestimmte Aufnahmetrajektorie, die jedoch nicht mehr die Tuy-Bedingung erfüllen muss. Dabei sei bereits an dieser Stelle deutlich hervorgehoben, dass nicht zwangsläufig die Aufnahmeanordnung bewegt sein muss, sondern auch das Objekt bewegt werden kann, um die Aufnahmetrajektorie zu realisieren. Beispielsweise kann, falls das Objekt ein Patient oder in einem Patienten vorhanden ist, eine Patientenliege genutzt werden, die bewegbar ausgestaltet ist. Dies ist insbesondere vorteilhaft bei der Anwendung der vorliegenden Erfindung auf nuklearmedizinische Bildaufnahmeverfahren, beispielsweise PET oder SPECT, wobei beispielsweise bei der PET der Patient durch die ringförmige Detektoranordnung gefahren werden kann, so dass sich insgesamt eine Aufnahmetrajektorie ergibt. Hauptsächlich findet die Erfindung jedoch Anwendung auf Röntgeneinrichtungen. There is still a relative movement between object and recording arrangement provided, thus a certain recording trajectory, but no longer has to fulfill the Tuy condition. It should be emphasized at this point that not necessarily the recording arrangement must be moved, but also the object can be moved to realize the recording trajectory. For example, if the object is a patient or in a patient, a patient bed can be used, which is designed to be movable. This is particularly advantageous in the application of the present invention to nuclear medicine image recording methods, for example PET or SPECT, wherein, for example, in PET, the patient can be driven through the annular detector arrangement, so that overall results in a recording trajectory. Mainly, however, the invention finds application to X-ray equipment.
Mithin kann zweckmäßigerweise vorgesehen sein, dass als Bildaufnahmeeinrichtung eine Röntgeneinrichtung verwendet wird, die eine eine Strahlungsquelle und einen, insbesondere mit der Strahlungsquelle bewegten, Detektor umfassende Aufnahmeanordnung aufweist, wobei die Aufnahmetrajektorie für die Strahlungsquelle definiert ist. Dabei sei darauf hingewiesen, dass selbst bei einer Bewegung des Patienten der Patient bzw. das Objekt als ruhend angenommen werden kann, mithin eine Aufnahmetrajektorie für die Strahlungsquelle definiert werden kann. Die Röntgeneinrichtung kann dabei beispielsweise eine Röntgeneinrichtung mit einem C-Bogen sein, an dem sich gegenüberliegend die Strahlungsquelle und der Detektor angeordnet sind. Selbstverständlich lässt sich das Verfahren jedoch auch auf Computertomographie-Bildaufnahmeeinrichtungen anwenden. Consequently, it may be expediently provided that an X-ray device is used as the image recording device, which has a recording device comprising a radiation source and a detector, in particular moving with the radiation source, the recording trajectory being defined for the radiation source. It should be noted that even with a movement of the patient, the patient or the object can be assumed to be at rest, thus a recording trajectory for the radiation source can be defined. The X-ray device can be, for example, an X-ray device with a C-arm, on which the radiation source and the detector are arranged opposite each other. Of course, however, the method can also be applied to computed tomography imagers.
Die der Erfindung zugrunde liegende Bewegung sei für das Beispiel einer geradlinigen Bewegung der Aufnahmeanordnung senkrecht zu einem Zentralstrahl bei einer flachen Bildebene, die insbesondere parallel zu einer Detektorfläche verläuft, genauer erläutert. Die Aufnahmeanordnung führt also eine geradlinige Bewegung in eine bestimmte Richtung, beispielsweise eine z-Richtung, aus, wobei während der Bewegung an verschiedenen Positionen entlang der Richtung Projektionsbilder des darzustellenden Objekts aufgenommen werden. Diese Projektionsdaten sind durch die entsprechenden Detektorelemente (k, l) und die z-Position der Aufnahmeanordnung kennzeichnet, lassen sich also als p(z, k, l) schreiben. The movement on which the invention is based is explained in more detail for the example of a rectilinear movement of the receiving arrangement perpendicular to a central beam at a flat image plane, which runs in particular parallel to a detector surface. The recording arrangement thus executes a rectilinear movement in a specific direction, for example a z-direction, during which movement projection images of the object to be displayed are recorded at different positions along the direction. These projection data are characterized by the corresponding detector elements (k, l) and the z-position of the recording arrangement, can thus be written as p (z, k, l).
Werden nun die Projektionsdaten p(z, k, l) aus den unterschiedlichen z-Positionen in eine beliebige Ebene rückprojiziert, wobei selbstverständlich die durch die Bildebene und die Aufnahmetrajektorie beschriebene Aufnahmegeometrie berücksichtigt wird, erfolgt eine adäquate Rückprojektion der Strukturen im Objekt nur für exakt diese Bildebene genau, also „scharf“. Objektstrukturen aus anderen, zur Bildebene parallelen Ebenen werden im Rekonstruktionsbild geometrisch verschmiert, was sich mathematisch durch einen Tiefpassfilter beschreiben lässt. Diese Tiefpassfilterung der Projektionsdaten für Objektstrukturen außerhalb der zu rekonstruierenden Bildebene führt zu einer Dämpfung des entsprechenden Bildkontrastes für diese Objektstrukturen. Lediglich für die Bildebene selber erfolgt eine korrekte Darstellung der Objektstrukturen bei maximalem Bildkontrast. If now the projection data p (z, k, l) are backprojected from the different z-positions into an arbitrary plane, whereby, of course, the imaging geometry described by the image plane and the recording trajectory is taken into account, an adequate backprojection of the structures in the object takes place only for exactly these Image plane exactly, so "sharp". Object structures from other planes parallel to the image plane are blurred geometrically in the reconstruction image, which can be mathematically described by a low-pass filter. This low-pass filtering of the projection data for object structures outside the image plane to be reconstructed leads to a damping of the corresponding image contrast for these object structures. Only for the image plane itself is a correct representation of the object structures at maximum image contrast.
Aus dieser Erkenntnis, die der vorliegenden Erfindung zugrunde liegt, folgt jedoch, dass auch eine Korrektur möglich ist, wenn die Tiefpassfilterung, die sich also aus der Aufnahmegeometrie bestimmt, letztlich wieder umgekehrt wird. Die Idee ist es, die tiefpassgefilterten Bildbeiträge aus anderen Ebenen als der Bildebene zu eliminieren, so dass ein scharfes, möglichst artefaktfreies Rekonstruktionsbild der ausgewählten Bildebene ermittelt wird. Dabei gibt es verschiedene Möglichkeiten, den Tiefpassfilterkern, der die beschriebene Tiefpassfilterung beschreibt, einzusetzen, worauf im Folgenden noch näher eingegangen werden wird. From this finding, which is the basis of the present invention, it follows, however, that a correction is possible if the low-pass filtering, which thus determined from the recording geometry, is ultimately reversed again. The idea is to eliminate the low-pass filtered image contributions from planes other than the image plane, so that a sharp, as artifact-free reconstruction image of the selected image plane is determined. There are various possibilities for using the low-pass filter kernel, which describes the described low-pass filtering, which will be discussed in more detail below.
Dabei bringt das erfindungsgemäße Verfahren eine Vielzahl von Vorteilen mit sich. So ist die Rekonstruktion von dreidimensionalen Objektstrukturen auch bei nicht die Tuy-Bedingung erfüllenden Aufnahmetrajektorien möglich, insbesondere also auch bei geradlinigen Bewegungstrajektorien. Die Aufnahmetrajektorie kann so gewählt werden, dass Rotationsbewegungen, die eine stark störende Wirkung auf den Workflow haben können, entfallen. In Ausführungsformen der Erfindung ist eine kurze geradlinige Bewegung der Aufnahmeanordnung relativ zu dem Objekt bereits ausreichend, um die dreidimensionale Struktur darzustellen. Hieraus ergibt sich auch, dass insgesamt, nachdem die mechanische Bewegung verringert werden kann, die Datenakquisition in einem sehr kurzen Zeitfenster erfolgen kann, so dass sich eine kürzere Aufnahmezeit ergibt und die zeitliche Auflösung der Rekonstruktionsbilder verbessert wird. The method according to the invention brings a large number of advantages. Thus, the reconstruction of three-dimensional object structures is also possible in the case of non-Tuy-conditional recording trajectories, especially in the case of linear motion trajectories. The recording trajectory can be chosen so that rotational movements, which can have a strong disruptive effect on the workflow, are eliminated. In embodiments of the invention, a short rectilinear motion of the receiver assembly relative to the object is already sufficient to represent the three-dimensional structure. It also follows that overall, after the mechanical movement can be reduced, the data acquisition can be done in a very short time window, resulting in a shorter recording time and the temporal resolution of the reconstruction images is improved.
Wie bereits dargelegt wurde, sieht eine besonders vorteilhafte Ausführungsform der vorliegenden Erfindung vor, dass die Aufnahmetrajektorie wenigstens teilweise, bevorzugt vollständig, eine Gerade ist. Gerade Bewegungen lassen sich deutlich einfacher realisieren als Drehbewegungen oder sonstige Kurven entlang der Aufnahmetrajektorie, so dass eine deutliche Vereinfachung der benötigten mechanischen Bewegung gegeben ist, die auch weniger hinderlich auf den Zugang zum Objekt und den sonstigen medizinischen Workflow wirkt. As already stated, a particularly advantageous embodiment of the present invention provides that the recording trajectory is at least partially, preferably completely, a straight line. Straight movements can be realized much easier than turning movements or other curves along the recording trajectory, so that a clear simplification of the required mechanical movement is given, which also acts less hindering on the access to the object and the other medical workflow.
Besonders zweckmäßig ist die Verwendung geradliniger Aufnahmetrajektorien dann, wenn auch die Bildebene flach, also krümmungsfrei, ist. The use of rectilinear photography trajectories is particularly expedient if the image plane is also flat, that is, free of curvature.
Wie bereits dargelegt wurde, sind im Rahmen der vorliegenden Erfindung jedoch auch gekrümmte Bildebenen denkbar. Dann kann vorgesehen sein, dass die Aufnahmetrajektorie wenigstens teilweise so gewählt wird, dass ihre Tangente parallel zu einer Tangente der insbesondere gekrümmten Bildebene an einem über den Zentralstrahl verbundenen Ort liegt. Wird mithin beispielsweise die Aufnahmetrajektorie über die Strahlungsquelle definiert, existiert an jeder Aufnahmeposition entlang der Aufnahmetrajektorie ein Zentralstrahl, der die Zuordnung von Punkten auf der Aufnahmetrajektorie zu Punkten auf der Bildebene erlaubt. Wird nun die Aufnahmetrajektorie so gewählt, dass sie auf die beschriebene Art parallel zur Bildebene läuft, ergibt sich eine besonders einfache Wahl des Tiefpassfilterkerns, worauf im Folgenden noch näher eingegangen werden wird. Es sei darauf hingewiesen, dass es auch äußerst zweckmäßig ist, wenn bei einer flachen Bildebene die Aufnahmetrajektorie parallel zur Bildebene gewählt wird. As has already been explained, however, curved image planes are also conceivable within the scope of the present invention. Then it can be provided that the recording trajectory is at least partially selected so that their tangent parallel to a tangent of the particular curved image plane at one over the central beam connected location is located. If, for example, the recording trajectory is defined via the radiation source, a central beam exists at each recording position along the recording trajectory, which allows the assignment of points on the recording trajectory to points on the image plane. Now, if the recording trajectory is chosen so that it runs parallel to the image plane in the described manner, results in a particularly simple choice of the low-pass filter core, which will be discussed in more detail below. It should be noted that it is also extremely useful if the recording trajectory is selected parallel to the image plane in the case of a flat image plane.
Es ist jedoch auch denkbar, dass eine wenigstens teilweise senkrecht zu der ebenen Bildebene verlaufende Aufnahmetrajektorie gewählt wird, zumindest dann, wenn die Kegelstrahlgeometrie verwendet wird. In dieser Ausführungsform bewegt sich die Aufnahmeanordnung relativ bezüglich des Objektes entlang einer geradlinigen Aufnahmetrajektorie senkrecht zur Bildebene und senkrecht zur Detektorebene. Bei einer Röntgeneinrichtung wird also der Abstand von der Strahlungsquelle zum Objekt geändert. Die notwendige geometrische Verschmierung von Objektstrukturen außerhalb der zu rekonstruierenden Bildebene entsteht in diesem Fall durch Veränderung der Kegelstrahlgeometrie. Dabei sei an dieser Stelle nochmals angemerkt, dass eine beliebige Kombination der beschriebenen Bewegungstrajektorien selbstverständlich auch möglich ist, der Einfachheit halber hier jedoch nicht im Detail beschrieben sein soll. However, it is also conceivable that a recording trajectory extending at least partially perpendicular to the flat image plane is selected, at least when the cone beam geometry is used. In this embodiment, the pickup assembly moves relative to the object along a rectilinear photography trajectory perpendicular to the image plane and perpendicular to the detector plane. In the case of an X-ray device, therefore, the distance from the radiation source to the object is changed. The necessary geometric smearing of object structures outside the image plane to be reconstructed arises in this case by changing the cone beam geometry. It should be noted again at this point that any combination of the described movement trajectories is of course also possible, for the sake of simplicity, however, should not be described in detail here.
Äußerst zweckmäßig ist es ferner, wenn Rekonstruktionsbilder für mehrere parallele Bildebenen ermittelt werden. So kann beispielsweise eine feste Anzahl von N Rekonstruktionsbildern in N parallelen Bildebenen rekonstruiert werden, wobei jede Bildrekonstruktion die dedizierte Aufnahmegeometrie, wie beschrieben, berücksichtigt. Ein auf diese Art erzeugter Rekonstruktionsbildstapel lässt sich in jedem gängigen medizinischen Bildbetrachtungsprogramm anzeigen, wo auch eine Nachverarbeitung stattfinden kann. It is also extremely expedient if reconstruction images for a plurality of parallel image planes are determined. For example, a fixed number of N reconstruction images may be reconstructed in N parallel image planes, with each image reconstruction taking into account the dedicated acquisition geometry as described. A reconstruction image stack generated in this way can be displayed in any conventional medical image viewer, where post-processing can also take place.
In Weiterbildung der vorliegenden Erfindung kann vorgesehen sein, dass, insbesondere bei dem Krümmungsprofil der Bildebene folgenden Aufnahmetrajektorien, als Tiefpassfilterkern eine Rechteckfunktion verwendet wird. Bei Untersuchungen im Rahmen der vorliegenden Erfindung hat sich gezeigt, dass dann, wenn die Aufnahmetrajektorie parallel zur Bildebene verläuft, die auftretende Tiefpassfilterung von Objektstrukturen anderer Bildebenen als der rekonstruierten Bildebene durch eine Rechteckfunktion, in der Mathematik als rect(z) bezeichnet, verwendet werden kann. Findet eine Faltung mit dem Tiefpassfilterkern statt, ist der wesentliche, von der Aufnahmegeometrie, also der Aufnahmetrajektorie und der Bildebene, abhängige Parameter die Breite der Rechteckfunktion, welche beispielsweise, worauf noch näher eingegangen werden wird, in einem Optimierungsverfahren oder durch theoretische Berechnungen ermittelt werden kann. In a development of the present invention, provision can be made for a rectangular function to be used as the low-pass filter kernel, in particular in the case of the absorption profile of the image plane following the recording trajectories. Investigations in the context of the present invention have shown that, when the recording trajectory is parallel to the image plane, the low-pass filtering of object structures of image planes other than the reconstructed image plane can be used by a rectangular function in which mathematics is called rect (z) , If a convolution occurs with the low-pass filter kernel, the essential parameter dependent on the recording geometry, that is to say the recording trajectory and the image plane, is the width of the rectangular function, which, for example, as will be discussed in more detail, can be determined in an optimization method or by theoretical calculations ,
Dabei sei an dieser Stelle nochmals angemerkt, dass die Aufnahmetrajektorie und die Bildebene zwar die Aufnahmegeometrie definiert, welche letztlich anzeigt, welche durch den Detektor gemessenen Strahlen, insbesondere Röntgenstrahlen, wo die Bildebene durchtreten haben. Dabei sind selbstverständlich auch weitere Größen, die die Aufnahmeanordnung betreffen, relevant, beispielsweise im Fall der Kegelstrahlgeometrie der Öffnungswinkel des Kegels, der Abstand von der Strahlungsquelle zum Detektor und dergleichen. Derartige Größen und ihr Einfluss auf die Aufnahmegeometrie sind im Stand der Technik jedoch weithin bekannt, so dass hier und in der weiteren Beschreibung davon ausgegangen wird, dass diese auch entsprechend Berücksichtigung finden. It should be noted again at this point that the recording trajectory and the image plane, although defining the recording geometry, which ultimately indicates which measured by the detector beams, in particular X-rays where the image plane have passed. In this case, of course, other sizes that affect the receiving arrangement, relevant, for example, in the case of the cone beam geometry, the opening angle of the cone, the distance from the radiation source to the detector and the like. However, such quantities and their influence on the geometry of the recording are well known in the prior art, so that it is assumed here and in the further description that these are also taken into account accordingly.
Eine besonders vorteilhafte Ausführungsform der vorliegenden Erfindung sieht vor, dass zur Korrektur von den rückprojizierten Bilddaten die mit dem Tiefpassfilter gefalteten rückprojizierten Bilddaten derselben Bildebene oder wenigstens einer zu der Bildebene parallelen Ebene, insbesondere bei Ermittlung von Rekonstruktionsbildern in mehreren parallelen Bildebenen jeder anderen Bildebene, subtrahiert werden. Damit sind konkrete Korrekturvorschriften gegeben, die im vorliegenden Fall letztlich zwei grundsätzliche, alternative Wege wiedergeben. A particularly advantageous embodiment of the present invention provides that the backprojected image data folded with the low-pass filter are subtracted from the backprojected image data for the same image plane or at least one plane parallel to the image plane, in particular when determining reconstruction images in several parallel image planes of every other image plane , This gives concrete correction rules, which in the present case ultimately reflect two fundamental, alternative ways.
In einem ersten Fall werden mithin die rückprojizierten Bilddaten in der zu rekonstruierenden Bildebene mit dem Tiefpassfilterkern gefaltet und von den rückprojizierten Bilddaten der gerade betrachteten Bildebene, also derselben Bildebene, abgezogen. Im Endeffekt wird dabei ein Hochpassfilter auf ein rekonstruiertes Bild einer bestimmten Bildebene angewendet. Beispielsweise kann dabei die Faltung mit einer Rechteckfunktion vorgesehen sein, welche die Einflüsse aus anderen Objektebenen gemäß der Aufnahmegeometrie berücksichtigt, so dass die tiefpassgefilterten Bildbeiträge aus anderen Ebenen des Objekts eliminiert werden. In a first case, therefore, the backprojected image data in the image plane to be reconstructed is folded with the low-pass filter kernel and subtracted from the backprojected image data of the currently viewed image plane, ie the same image plane. In effect, a high-pass filter is applied to a reconstructed image of a given image plane. For example, the convolution may be provided with a rectangular function, which takes into account the influences from other object planes according to the acquisition geometry, so that the low-pass filtered image contributions from other planes of the object are eliminated.
Ein anderer Ansatz sieht vor, dass eine feste Anzahl von Rekonstruktionsbildern in parallelen Bildebenen rekonstruiert werden. Dabei soll ein initial rekonstruiertes Bild unter Berücksichtigung der tiefpassfilternden Wirkung des Rekonstruktionsverfahrens durch initial rekonstruierte Bilder anderer Bildebenen korrigiert werden. Dabei ist vorgesehen, dass zunächst tiefpassgefilterte rückprojizierte Bilddaten für alle Bildebenen ermittelt werden, beispielsweise durch Faltung mit einer Rechteckfunktion als Tiefpassfilterkern. Um nun das korrigierte Rekonstruktionsbild einer bestimmten Bildebene n zu erhalten, werden von den rückprojizierten (unkorrigierten) Bilddaten dieser Bildebene n alle gefilterten rückprojizierten Bilddaten aller anderen Bildebenen m ungleich n abgezogen. Another approach is to reconstruct a fixed number of reconstruction images in parallel image planes. In this case, an initially reconstructed image is corrected, taking into account the low-pass filtering effect of the reconstruction process, by means of initially reconstructed images of other image planes become. It is provided that initially low-pass filtered back-projected image data for all image planes are determined, for example, by convolution with a rectangular function as a low-pass filter core. In order to obtain the corrected reconstruction image of a specific image plane n, all filtered backprojected image data of all other image planes m not equal to n are subtracted from the backprojected (uncorrected) image data of this image plane n.
Dies alles sei in einer konkreten Ausführungsform im Folgenden näher erläutert. Dabei wird, was im Rahmen der vorliegenden Erfindung auch grundsätzlich möglich ist, ein Rückprojektionsverfahren eingesetzt, welches lediglich die Aufnahmegeometrie bezüglich der zu rekonstruierenden Bildebene n berücksichtigt und weitere auf die Rückprojektion bezogene Rechenoperationen auf den Projektionsdaten vermeidet, so dass beispielsweise keine Fourier-Transformation durchgeführt wird. Werden die Projektionsdaten der Projektionsbilder wiederum als p(z, k, l) geschrieben und ist die Aufnahmetrajektorie eine Gerade entlang der zur Detektorebene und zur Bildebene parallelen Richtung z, lässt sich schreiben: All this will be explained in more detail below in a concrete embodiment. In this case, what is fundamentally possible in the context of the present invention, a backprojection method is used, which takes into account only the recording geometry with respect to the image plane to be reconstructed n and avoids further on the back projection related arithmetic operations on the projection data, so that, for example, no Fourier transform is performed , If the projection data of the projection images are again written as p (z, k, l) and if the recording trajectory is a straight line along the direction z parallel to the detector plane and to the image plane, it is possible to write:
Hierin bezeichnet x die weitere, neben z vorhandene Koordinate in der Bildebene n, In,nativ(x, z) die rückprojizierten Bilddaten (noch unkorrigiert), rSDD ist der Abstand zwischen der Strahlungsquelle und dem Detektor, NDet beschreibt die Gesamtanzahl der Detektorelemente, ΔkDet bezeichnet die Detektorauflösung in k-Richtung, ΔlDet ist die Detektorauflösung in l-Richtung und rRecon,n bezeichnet den Abstand der Bildebene von der Strahlungsquelle. Dabei sei an dieser Stelle noch angemerkt, dass sich aus der Rekonstruktionsvorschrift die offensichtliche Notwendigkeit ergibt, dass die Aufnahme der Projektionsdaten mit einem Sampling-Inkrement in z-Richtung Δz << (ΔkDet·Nk Det) erfolgen muss, das deutlich kleiner als die räumliche Detektorausdehnung in k-Richtung ist. Allgemein bedeutet dies, dass die Schrittweite zwischen der Aufnahme zweier Projektionsbilder entlang der Aufnahmetrajektorie kleiner als ein Zehntel der Detektorausdehnung in der Richtung der Aufnahmetrajektorie sein sollte. Here x denotes the further coordinate present in the image plane n, I n, native (x, z) the backprojected image data (still uncorrected), r SDD is the distance between the radiation source and the detector, N det describes the total number of Detector elements, Δk Det denotes the detector resolution in the k-direction, Δl Det is the detector resolution in the l-direction and r Recon, n denotes the distance of the image plane from the radiation source. It should be noted at this point that arises from the reconstruction rule the obvious need that the recording of the projection data with a sampling increment in the z-direction Δz << (Δk Det · N k Det ) must be made, which is significantly smaller than is the spatial detector extent in the k direction. Generally, this means that the step size between the acquisition of two projection images along the acquisition trajectory should be less than one tenth of the detector extent in the direction of the acquisition trajectory.
Wird ein derartiges Rekonstruktionsverfahren eingesetzt, äußert sich die Tiefpassfilterung als Tiefpassfilterkern einer Rechteckfunktion. Entsprechend ist dann, wenn die Korrektur in derselben Bildebene erfolgen soll, das korrigierte Rekonstruktionsbild In,corr zu berechnen zu
Dabei deutet das Zeichen „x“ die Faltung an. Die Breite der Rechteckfunktion rect(z) ergibt sich in Abhängigkeit der konkreten Aufnahmegeometrie und kann aus Berechnungen und/oder in einem Optimierungsverfahren ermittelt werden. The sign "x" indicates the convolution. The width of the rectangular function rect (z) depends on the specific recording geometry and can be determined from calculations and / or in an optimization method.
Soll die Korrektur, wie oben beschrieben, durch Beachtung der rückprojizierten Bilddaten in anderen, parallelen Bildebenen erfolgen, wird mithin eine feste Anzahl von n > 1 Rekonstruktionsbildern ermittelt. Bezeichnet man die noch unkorrigierten, rückprojizierten Bilddaten gemäß Formel (1), also das initial rekonstruierte Bild, wieder mit In,nativ, werden in dieser Ausführungsform der vorliegenden Erfindung zunächst die tiefpassgefilterten Bildbeiträge aus Bildebenen m ungleich n in der Bildebene n abgeschätzt. Diese ergeben sich zu
Diese entsprechenden Bildbeiträge In,m der Bildebenen m sollen entsprechend eliminiert werden, so dass das korrigierte Rekonstruktionsbild In,corr bereitgestellt wird. These corresponding image contributions I n, m of the image planes m should be eliminated accordingly, so that the corrected reconstruction image I n, corr is provided.
Dabei sei an dieser Stelle noch angemerkt, dass sich die Formel (1) auch leicht für eine Aufnahmegeometrie umschreiben lässt, in der eine senkrecht zu der ebenen Bildebene, welche wiederum parallel zur Detektorebene verläuft, verlaufende Aufnahmetrajektorie gewählt wird. Dann ergibt sich zur Ermittlung der rückprojizierten Bilddaten It should also be noted at this point that the formula (1) can also be easily described for a recording geometry in which a perpendicular to the flat image plane, which in turn runs parallel to the detector plane, running recording trajectory is selected. Then results in the determination of the backprojected image data
In einer zweckmäßigen Weiterbildung kann vorgesehen sein, dass die Korrektur iterativ angewandt wird, indem neue Rekonstruktionsbilder durch Anwendung des Tiefpassfilters auf bisherige Rekonstruktionsbilder und entsprechende Subtraktion vom dem ursprünglichen rückprojizierten Bilddaten ermittelt werden. Das bedeutet, es ergibt sich, insbesondere im Fall der Formeln (3) und (4), ein iteratives Multi-Ebenen-Korrekturverfahren. Es wird eine feste Anzahl von N Bildern in N parallelen Bildebenen rekonstruiert und gemäß (3) und (4) korrigiert. Die korrigierten Rekonstruktionsbilder werden wieder an den Korrekturalgorithmus gegeben und erneut verwendet, um die tiefpassgefilterten Bildbeiträge in den anderen Bildebenen gemäß (3) abzuschätzen und gemäß (4) vom initial rekonstruierten Bild abzuziehen. Ein derartiger iterativer Ansatz kann mit einer festen Zahl an Iterationsschritten erfolgen, denkbar ist es jedoch auch, ein anderes Abbruchkriterium zu verwenden, beispielsweise nach Art einer Optimierungsschwelle oder dergleichen. In an expedient development it can be provided that the correction is applied iteratively by determining new reconstruction images by applying the low-pass filter to previous reconstruction images and corresponding subtraction from the original backprojected image data. That is, in the case of formulas (3) and (4) in particular, it produces an iterative multi-level correction method. A fixed number of N pictures in N parallel picture planes are reconstructed and corrected according to (3) and (4). The corrected reconstruction images are again applied to the correction algorithm and reused to estimate the low-pass filtered image contributions in the other image planes according to (3) and subtract from the initial reconstructed image according to (4). Such an iterative approach can be done with a fixed number of iteration steps, but it is also conceivable to use a different termination criterion, for example in the manner of an optimization threshold or the like.
Wie bereits angedeutet wurde, kann wenigstens ein die konkrete Form des Tiefpassfilterkerns beschreibender Filterparameter, beispielsweise die Breite einer Rechteckfunktion, in einem Optimierungsverfahren angepasst werden. Dabei können beispielsweise Kriterien der Bildqualität berücksichtigt werden, insbesondere im Hinblick auf den Kontrast bestimmter Strukturen. Auch die Verwendung von Vergleichsbildern und dergleichen ist denkbar. Es sei an dieser Stelle nochmals hervorgehoben, dass die konkrete Wahl des Tiefpassfilterkerns stark von der Aufnahmegeometrie, insbesondere also der Lage der Bildebene und der Aufnahmetrajektorie gemeinsam mit den Eigenschaften der Aufnahmeanordnung an sich, abhängig ist. Für bestimmte Sachverhalte, beispielsweise dem Krümmungsprofil der Bildebene folgende Aufnahmetrajektorien und dergleichen, bei denen die Detektorebene auch parallel zur Tange bzw. Bildebene liegt, ergibt sich, wie Untersuchungen des Erfinders belegt haben, als Tiefpassfilterkern eine Rechteckfunktion, die zu einer hervorragenden Korrektur verwendet werden kann. Für andere, insbesondere komplexere Aufnahmegeometrien und Kombinationen von Aufnahmegeometrien, beispielsweise teilweise senkrecht und teilweise parallel zur Bildebene verlaufende Aufnahmetrajektorien, ergeben sich gegebenenfalls andere, sinnvoll zu verwendende Tiefpassfilterkerne. Gezeigt hat sich jedoch, dass der Tiefpassfilterkern der Rechteckfunktion auch dann anwendbar ist, wenn die Aufnahmetrajektorie senkrecht zu der ebenen, zur Detektorebene parallelen Bildebene verläuft. As already indicated, at least one filter parameter describing the specific form of the low-pass filter kernel, for example the width of a rectangular function, can be adapted in an optimization method. In this case, for example, criteria of image quality can be taken into account, in particular with regard to the contrast of certain structures. The use of comparative images and the like is conceivable. It should again be emphasized here that the specific choice of the low-pass filter kernel is strongly dependent on the recording geometry, in particular therefore the position of the image plane and the recording trajectory together with the properties of the recording arrangement per se. For certain situations, such as the curvature profile of the image plane following recording trajectories and the like, in which the detector plane is also parallel to the tangent or image plane results, as investigations of the inventor have proven, as a low-pass filter kernel a rectangular function that can be used to an excellent correction , For other, in particular more complex recording geometries and combinations of recording geometries, for example, partially perpendicular and partially parallel to the image plane extending recording trajectories, there may be other, useful to use low-pass filter cores. However, it has been shown that the low-pass filter kernel of the rectangular function can also be used if the recording trajectory is perpendicular to the flat image plane parallel to the detector plane.
In einer weiteren, besonders vorteilhaften Ausführungsform der vorliegenden Erfindung kann vorgesehen sein, dass die Aufnahmetrajektorie und/oder die Bildebene in Abhängigkeit von Vorwissen über das aufzunehmende Objekt definiert werden. Es ist also denkbar, die Bildebene, aber auch die Aufnahmetrajektorie, mithin die gesamte Aufnahmegeometrie, gezielt so zu wählen, dass sich für das aufzunehmende Objekt eine optimale Bildqualität ergibt. Dabei kann vorgesehen sein, dass das Vorwissen aus einem Planungsbilddatensatz gewonnen wird, jedoch sind auch grundsätzlich andere Möglichkeiten denkbar. Planungsbilddatensätze, aus denen Eigenschaften des aufzunehmenden Objekts abgeleitet werden können, welche dann weitere Aufnahmeparameter beeinflussen, sind grundsätzlich bekannt. Insbesondere kann die Geometrie des aufzunehmenden Objekts betrachtet werden. Betrachtet man die der Erfindung zugrunde liegende Erkenntnis, dass ein Rückprojektionsverfahren, beispielsweise gemäß (1) oder (5), in der Wirkung einer Tiefpassfilterung in der Richtung der Aufnahmetrajektorie, insbesondere also in z-Richtung, entspricht, ergibt sich, dass eine optimale Darstellung eines Objekts in einer beliebigen Bildebene erfolgen kann, wenn sich die Vorzugsrichtung des Objekts in der Bildebene, aber senkrecht zur parallel zur Bildebene verlaufenden Aufnahmetrajektorie erstreckt. Es kann also mit besonderem Vorteil vorgesehen sein, dass die Aufnahmetrajektorie und die Bildebene so gewählt werden, dass eine Vorzugsrichtung des Objekts, insbesondere eine Richtung größter Ausdehnung des Objekts, in der Bildebene liegt und die Aufnahmetrajektorie senkrecht zu der Vorzugsrichtung und parallel zu der Bildebene verläuft. Dies bedeutet nichts anderes, als dass bei einer zweidimensionalen Modulationstransferfunktion des Objekts sich die Niederfrequenzsignalanteile senkrecht zur Aufnahmetrajektorie ausbilden. Für solche Trajektorien kann dann ein maximaler Kontrast zu Objekten aus anderen Bildebenen gewährleistet werden. In a further, particularly advantageous embodiment of the present invention, it can be provided that the recording trajectory and / or the image plane are defined as a function of prior knowledge of the object to be recorded. It is therefore conceivable to selectively select the image plane, but also the recording trajectory, and consequently the entire recording geometry, in such a way that optimum picture quality results for the object to be recorded. It can be provided that the prior knowledge is obtained from a planning image data record, but in principle, other possibilities are conceivable. Planning image data records from which properties of the object to be recorded can be derived, which then influence further acquisition parameters, are known in principle. In particular, the geometry of the object to be recorded can be considered. If one considers the knowledge underlying the invention that a backprojection method, for example according to (1) or (5), corresponds to the effect of low-pass filtering in the direction of the recording trajectory, ie in particular in the z-direction, this results in an optimum representation An object can take place in any image plane if the preferred direction of the object extends in the image plane but perpendicular to the recording trajectory parallel to the image plane. It can therefore be provided with particular advantage that the recording trajectory and the image plane are selected so that a preferred direction of the object, in particular a direction of maximum extension of the object, in the image plane and the receiving trajectory is perpendicular to the preferred direction and parallel to the image plane , This means nothing else than that in a two-dimensional modulation transfer function of the object, the low-frequency signal components perpendicular to the recording trajectory form. For such trajectories, a maximum contrast to objects from other image planes can then be guaranteed.
Wird als Bildaufnahmeeinrichtung, wie oben bereits dargelegt, eine Röntgeneinrichtung verwendet, die eine Strahlungsquelle und einen Detektor als Aufnahmeanordnung umfasst, ergeben sich weitere vorteilhafte Möglichkeiten, wenn mehr als eine Aufnahmeanordnung an der Röntgeneinrichtung vorgesehen sind. If an X-ray device which comprises a radiation source and a detector as a recording device is used as the image recording device, as already explained above, further advantageous possibilities result if more than one recording device is provided on the X-ray device.
So kann vorgesehen sein, dass eine Röntgeneinrichtung mit wenigstens zwei Aufnahmeanordnungen verwendet wird, die jeweils entlang einer ihnen zugeordneten Aufnahmetrajektorie Projektionsbilder aufnehmen. Beispiele für solche Röntgeneinrichtungen sind sogenannte Biplan-Systeme, insbesondere Biplan-C-Bogen-Systeme. Dann ist es möglich, die vorgestellten Berechnungsvorschriften und Korrekturverfahren entsprechend zu erweitern, so dass eine spezifische Berücksichtigung der unterschiedlichen Aufnahmeanordnungen und deren Lage zum Objekt erfolgt, wobei mehrere vorteilhafte Ausprägungen denkbar sind. Thus, it can be provided that an X-ray device with at least two recording arrangements is used, each of which receives projection images along an associated recording trajectory. Examples of such X-ray devices are so-called biplane systems, in particular biplane C-arm systems. It is then possible to correspondingly expand the presented calculation rules and correction methods, so that a specific consideration of the different recording arrangements and their position relative to the object takes place, whereby several advantageous features are conceivable.
So kann vorgesehen sein, dass das Rekonstruktionsbild unter Berücksichtigung von Projektionsbildern beider Aufnahmeanordnungen ermittelt wird. Sind beispielsweise die Aufnahmeanordnungen um einen bestimmten Winkel, beispielsweise 30°, gegeneinander verkippt, kann dies in die Rückprojektions- und Korrekturmöglichkeiten so einfließen, dass für eine bestimmte Bildebene Projektionsbilder beider Aufnahmeanordnungen Berücksichtigung finden können. Thus it can be provided that the reconstruction image is determined taking into account projection images of both recording arrangements. If, for example, the recording arrangements are tilted relative to one another by a certain angle, for example 30 °, this can be incorporated into the backprojection and correction possibilities in such a way that projection images of both recording arrangements can be taken into consideration for a specific image plane.
Besonders zweckmäßig ist es jedoch, wenn das Vorliegen von Projektionsbildern mehrerer Aufnahmeanordnungen dahingehend genutzt wird, auch dreidimensionale Volumenbilder des Objekts erzeugen zu können. So kann vorgesehen sein, dass aus Projektionsbildern beider Aufnahmeanordnungen jeweils denselben bestimmten Raumbereich abdeckende Rekonstruktionsbilder ermittelt werden, wonach aus den Rekonstruktionsbildern, insbesondere durch Mittelung von in einem gemeinsamen Voxel liegenden Bilddaten der Rekonstruktionsbilder, ein dreidimensionaler Rekonstruktionsbilddatensatz des Raumbereichs ermittelt wird. Dies sei anhand eines Beispiels näher erläutert. Wird ein Biplan-System mit zwei zueinander senkrecht ausgerichteten Aufnahmeanordnungen verwendet, ist es denkbar, einen quaderförmigen Raumbereich zu definieren, so dass für einen Projektionsbildersatz Rekonstruktionsbilder paralleler Bildebenen in der einen Hauptrichtung des Quaders ermittelt werden können, aus dem anderen Projektionsbildersatz Rekonstruktionsbilder in parallelen, zu den ersten Bildebenen senkrechten zweiten Bildebenen in einer zur ersten Hauptrichtung senkrechten Hauptrichtung ermittelt werden. Werden dann, beispielsweise an Schnittpunkten der ersten und zweiten Bildebenen, Voxel definiert, können die Bilddaten der Rekonstruktionsbilder innerhalb eines Voxels, gegebenenfalls auch gewichtet, zu einem einzigen Bilddatum des dreidimensionalen Rekonstruktionsbilddatensatzes zusammengeführt werden, beispielsweise durch Mittelung. So kann vorteilhaft ein dreidimensionaler, den Raumbereich wiedergebender Rekonstruktionsbilddatensatz rekonstruiert werden. However, it is particularly expedient if the presence of projection images of a plurality of recording arrangements is used to be able to also produce three-dimensional volume images of the object. It can thus be provided that reconstruction images covering the same specific spatial region are determined from projection images of both recording arrangements, according to which a three-dimensional reconstruction image data record of the spatial region is determined from the reconstruction images, in particular by averaging image data of the reconstruction images lying in a common voxel. This will be explained in more detail with reference to an example. If a biplanar system with two mutually vertically aligned recording arrangements is used, it is conceivable to define a cuboid spatial area, so that reconstruction images of parallel image planes in one main direction of the cuboid can be determined for a projection image set, and reconstruction images in parallel, from the other projection image set the first image planes perpendicular second image planes are determined in a main direction perpendicular to the first main direction main direction. If voxels are then defined, for example at intersections of the first and second image planes, the image data of the reconstruction images within a voxel, optionally also weighted, can be merged into a single image datum of the three-dimensional reconstruction image dataset, for example by averaging. Thus, advantageously, a three-dimensional reconstruction image data record representing the spatial region can be reconstructed.
Wie im konkreten Beispiel durch die Formeln (1) und (5) schon dargestellt, ist es allgemeiner gesprochen im erfindungsgemäßen Verfahren zweckmäßig, wenn zur Ermittlung der rückprojizierten Bilddaten für einen Bildpunkt des Rekonstruktionsbildes alle Beiträge aus den Projektionsbildern für diesen Bildpunkt aufsummiert werden. Letztlich werden also, wie aus den Formeln (1) und (5) bereits ersichtlich, die Strahlen von der Strahlungsquelle zum Detektor berücksichtigt, die durch den zu rekonstruierenden Bildpunkt des Rekonstruktionsbildes laufen, indem sie aufsummiert werden. Weitere Manipulationen im Rahmen der Rückprojektion sieht das erfindungsgemäße Verfahren dann nicht vor, insbesondere keine weitere Fourier-Transformation oder dergleichen. As already shown in the concrete example by the formulas (1) and (5), it is generally more appropriate in the method according to the invention if all contributions from the projection images for this pixel are summed to determine the backprojected image data for a pixel of the reconstruction image. Ultimately, therefore, as already apparent from formulas (1) and (5), the rays from the radiation source to the detector are taken into account, which run through the pixel of the reconstruction image to be reconstructed by summing them up. Further manipulations in the context of rear projection then does not provide the method according to the invention, in particular no further Fourier transformation or the like.
Neben dem Verfahren betrifft die vorliegende Erfindung auch eine Bildaufnahmeeinrichtung, die eine Aufnahmeanordnung und eine zur Durchführung des erfindungsgemäßen Verfahrens ausgebildete Steuereinrichtung aufweist. Sämtliche Ausführungen bezüglich des erfindungsgemäßen Verfahrens lassen sich analog auf die erfindungsgemäße Bildaufnahmeeinrichtung übertragen, so dass auch mit dieser die Vorteile der Erfindung erhalten werden können. In addition to the method, the present invention also relates to an image recording device which has a recording device and a control device designed to carry out the method according to the invention. All statements relating to the method according to the invention can be analogously transferred to the image recording device according to the invention, so that the advantages of the invention can be obtained with this.
Dabei handelt es sich mithin insbesondere um eine Röntgeneinrichtung, deren Aufnahmeanordnung, bevorzugt an einem C-Bogen befestigt, durch eine Strahlungsquelle und einen zugehörigen Detektor gebildet wird, der gemeinsam mit der Strahlungsquelle bewegbar sein kann. Wird eine Biplan-Röntgeneinrichtung verwendet, lassen sich auch die entsprechenden Ausgestaltungen des erfindungsgemäßen Verfahrens realisieren. In particular, it is therefore an X-ray device whose receiving arrangement, preferably attached to a C-arm, is formed by a radiation source and an associated detector, which can be movable together with the radiation source. If a biplane X-ray device is used, the corresponding embodiments of the method according to the invention can also be realized.
Weitere Vorteile und Einzelheiten der vorliegenden Erfindung ergeben sich aus den im Folgenden beschriebenen Ausführungsbeispielen sowie anhand der Zeichnungen. Dabei zeigen: Further advantages and details of the present invention will become apparent from the embodiments described below and with reference to the drawings. Showing:
Die hier dargestellten Ausführungsbeispiele beziehen sich hauptsächlich auf Röntgeneinrichtungen als Bildaufnahmeeinrichtungen und betrachten flache (also ebene) Bildebenen, die parallel zum Detektor der Röntgeneinrichtung verlaufen. Selbstverständlich sind auch Ausführungsbeispiele mit gekrümmten Bildebenen und für andere Bildaufnahmeeinrichtungen tomographischer Art, beispielsweise PET-Einrichtungen, denkbar, wie dies in der allgemeinen Beschreibung dargelegt wurde. The embodiments shown here relate mainly to X-ray devices as image recording devices and consider flat (ie flat) image planes that are parallel to the detector of the X-ray device. Of course, embodiments with curved image planes and other tomographic imaging devices, such as PET devices, are also conceivable, as set out in the general description.
Betrachtet man nun also wiederum
Entsprechend dieser Erkenntnis soll in einem Schritt
An dieser Stelle sei im Übrigen angemerkt, dass es grundsätzlich auch denkbar wäre, eine Korrektur nicht durchzuführen, sondern einem Betrachter unmittelbar ein unkorrigiertes Rekonstruktionsbild, beispielsweise das Rekonstruktionsbild
Um die Bildqualität jedoch auch weiterhin zu verbessern, wird die Korrektur im zweiten Ausführungsbeispiel iterativ durchgeführt, das bedeutet, in einem Schritt
Ist die Abbruchbedingung in Schritt
An dieser Stelle sei angemerkt, dass, insbesondere zur Bestimmung der Vorzugsrichtung
Eine zweckmäßige weitere Ausgestaltung des bezüglich
Wie aus
Die hier dargestellten Ausführungsbeispiele betreffen also hauptsächlich ein Verfahren zur Erzeugung eines einer flachen Bildebene zugeordneten Rekonstruktionsbildes eines Objekts aus einer Mehrzahl von zweidimensionalen, unter unterschiedlichen Projektionsrichtungen mit einer eine relativ zu dem Objekt entlang einer Aufnahmetrajektorie bewegten Strahlungsquelle und einen insbesondere mit der Strahlungsquelle bewegungsgekoppelten Detektor aufweisenden Aufnahmeanordnung einer Röntgeneinrichtung aufgenommenen Projektionsbildern, welches sich dadurch auszeichnet, dass eine nicht die Tuy-Bedingung erfüllende, parallel zu der Bildebene und der Empfangsfläche des Detektors (Detektorfläche/Detektorebene) oder senkrecht zu der Bildebene und der Empfangsfläche des Detektors verlaufende Aufnahmetrajektorie verwendet wird, welche gerade ist, wobei das Rekonstruktionsbild durch Rückprojektion für jeden zu rekonstruierenden, innerhalb der Bildebene liegenden Punkt, wobei rückprojizierte Bilddaten ermittelt werden, und Korrektur der rückprojizierten Bilddaten unter Verwendung eines von der Bildebene und der Aufnahmetrajektorie abhängigen Tiefpassfilterkerns, insbesondere einer Rechteckfunktion, ermittelt wird. Insbesondere werden dabei zur Ermittlung der rückprojizierten Bilddaten für einen Bildpunkt des Rekonstruktionsbildes alle Beiträge aus den Projektionsbildern für diesen Bildpunkt aufsummiert, und die Schrittweite zwischen der Aufnahme zweier Projektionsbilder entlang der Aufnahmetrajektorie sollte kleiner als ein Zehntel, idealerweise kleiner als ein Hundertstel, der Detektorausdehnung in der Richtung der Aufnahmetrajektorie sein. Bevorzugt wird die Aufnahmetrajektorie so gewählt, dass sie senkrecht zu der Vorzugsrichtung des Objekts liegt. Thus, the exemplary embodiments illustrated here relate mainly to a method for generating a reconstruction image of an object from a plurality of two-dimensional image sources associated with a flat image plane, with a radiation source moving relative to the object along a recording trajectory and a recording arrangement that is in particular motion-coupled to the radiation source An X-ray device recorded projection images, which is characterized in that a non-Tuy condition, parallel to the image plane and the receiving surface of the detector (detector surface / detector plane) or perpendicular to the image plane and the receiving surface of the detector extending recording trajectory is used, which straight where the reconstruction image is by backprojection for each point to be reconstructed within the image plane, with backprojected B ilddaten be determined, and correction of the backprojected image data using a dependent of the image plane and the recording trajectory low-pass filter core, in particular a rectangular function, is determined. In particular, all contributions from the projection images for this pixel are summed up to determine the backprojected image data for a pixel of the reconstruction image, and the increment between the acquisition of two projection images along the acquisition trajectory should be less than one tenth, ideally smaller than one hundredth, of the detector extent in the Be direction of the recording trajectory. Preferably, the recording trajectory is chosen so that it is perpendicular to the preferred direction of the object.
Wie jedoch dargestellt wurde, sind auch andere Ausgestaltungen der vorliegenden Erfindung, insbesondere mit gekrümmten Bildebenen und bei anderen Bildaufnahmeeinrichtungen denkbar. In jedem Fall sei nochmals hervorgehoben, dass die Aufnahmetrajektorie durch relative Bewegung zwischen dem Objekt und der Strahlungsquelle bzw. Aufnahmeanordnung definiert wird, mithin auch das Objekt bewegt werden kann. As has been shown, however, other embodiments of the present invention, in particular with curved image planes and in other image recording devices are conceivable. In any case, it should be emphasized once again that the recording trajectory is defined by relative movement between the object and the radiation source or recording arrangement, and consequently the object can also be moved.
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.
Claims (19)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE201310203410 DE102013203410A1 (en) | 2013-02-28 | 2013-02-28 | Method for producing flat or curved image plane associated reconstruction image of object i.e. organ of patient, by image pickup device, involves correcting back projected image data using filter kernel according to plane and trajectory |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE201310203410 DE102013203410A1 (en) | 2013-02-28 | 2013-02-28 | Method for producing flat or curved image plane associated reconstruction image of object i.e. organ of patient, by image pickup device, involves correcting back projected image data using filter kernel according to plane and trajectory |
Publications (1)
Publication Number | Publication Date |
---|---|
DE102013203410A1 true DE102013203410A1 (en) | 2014-08-28 |
Family
ID=51349513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE201310203410 Withdrawn DE102013203410A1 (en) | 2013-02-28 | 2013-02-28 | Method for producing flat or curved image plane associated reconstruction image of object i.e. organ of patient, by image pickup device, involves correcting back projected image data using filter kernel according to plane and trajectory |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE102013203410A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102015224176A1 (en) * | 2015-12-03 | 2017-06-08 | Siemens Healthcare Gmbh | Tomography system and method for generating a sequence of volume images of a vascular system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5473655A (en) * | 1994-07-21 | 1995-12-05 | General Electric Company | Artifact reduction by z-dependent filtration of three-dimensional cone beam data |
US8189735B2 (en) * | 2010-07-22 | 2012-05-29 | General Electric Company | System and method for reconstruction of X-ray images |
DE102011077389A1 (en) * | 2011-06-10 | 2012-12-13 | Siemens Aktiengesellschaft | X-ray image receiving device, has X-ray tube and X-ray detector arranged at X-ray c-arms, and evaluating device for computing three-dimensional X-ray image data set from two-dimensional data sets obtained by X-ray c-arms |
-
2013
- 2013-02-28 DE DE201310203410 patent/DE102013203410A1/en not_active Withdrawn
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5473655A (en) * | 1994-07-21 | 1995-12-05 | General Electric Company | Artifact reduction by z-dependent filtration of three-dimensional cone beam data |
US8189735B2 (en) * | 2010-07-22 | 2012-05-29 | General Electric Company | System and method for reconstruction of X-ray images |
DE102011077389A1 (en) * | 2011-06-10 | 2012-12-13 | Siemens Aktiengesellschaft | X-ray image receiving device, has X-ray tube and X-ray detector arranged at X-ray c-arms, and evaluating device for computing three-dimensional X-ray image data set from two-dimensional data sets obtained by X-ray c-arms |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102015224176A1 (en) * | 2015-12-03 | 2017-06-08 | Siemens Healthcare Gmbh | Tomography system and method for generating a sequence of volume images of a vascular system |
US10702220B2 (en) | 2015-12-03 | 2020-07-07 | Siemens Healthcare Gmbh | Tomography system and method for producing a series of volume images of a vascular system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102012204019B4 (en) | Method for reducing motion artifacts | |
DE102012207629B4 (en) | CT image reconstruction in the extended measuring field | |
DE102013218819B3 (en) | Method of reducing artifacts in an image data set and X-ray device | |
DE102012204977B4 (en) | Method for iterative image reconstruction for bi-modal CT data and corresponding computation unit, computed tomography system, computer program and data carrier | |
DE102009032059A1 (en) | Sinogram processing for metal artifact reduction in computed tomography | |
DE602004004877T2 (en) | SYSTEM AND METHOD FOR EXACT IMAGE RECONSTRUCTION FOR SPIRAL BEAM COMPUTER TERTOMOGRAPHY WITH REDUNDANT DATA | |
DE102011087337B4 (en) | A method of reconstructing a reconstruction data set containing two-dimensional virtual X-ray images | |
DE102011083643A1 (en) | Method, computer system and CT system for determining a motion field and for motion-compensating reconstruction with this motion field | |
DE102010034099B4 (en) | Iterative image filtering with anisotropic noise model for a CT image | |
DE102009039987A1 (en) | Iterative CT image filter for noise reduction | |
DE102010006585A1 (en) | CT image reconstruction in the extended measuring field | |
DE102010019632A1 (en) | Method for recording and reconstructing a three-dimensional image data set and x-ray device | |
DE102011083647A1 (en) | Motion compensated computer tomography-image data set creating method, involves reconstructing final image data set using movement-compensating reconstruction method based on reconstruction algorithm and movement field | |
DE102012220028A1 (en) | Angiographic examination procedure | |
DE112015002802T5 (en) | IMPROVED IMAGE RECONSTRUCTION FOR A VOLUME BASED ON PROJECTION RECORDS | |
DE102011005715A1 (en) | Method for obtaining a 3D image data set freed from traces of a metal object | |
DE102005044653A1 (en) | Method and device for reconstructing a three-dimensional image volume from two-dimensional projection images | |
DE102007061934A1 (en) | Method for improving the quality of computer tomographic acquisition series by projection data processing and CT system with arithmetic unit | |
DE102017200930A1 (en) | Method for determining a three-dimensional image data record | |
DE102008031530A1 (en) | Method for correcting truncation artifacts | |
DE102013200329B4 (en) | Method and device for misalignment correction for imaging methods | |
DE102009007236A1 (en) | CT image reconstruction of a moving examination object | |
DE102012205245A1 (en) | Angiographic investigation method of organ of patient, involves changing resolution such that large volume data acquisition takes place in low resolution and small volume data acquisition in high resolution | |
DE102011081167A1 (en) | Rear projection of a projection image data set with depth-dependent filtering | |
DE102012217613A1 (en) | Method for artifact reduction in specific-dimensional image data set, involves determining corrected projection data, in which for each pixel of projection image, correction data is determined |
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 |