EP2770911A1 - Verfahren zur erzeugung optimierter tomografie-aufnahmen - Google Patents

Verfahren zur erzeugung optimierter tomografie-aufnahmen

Info

Publication number
EP2770911A1
EP2770911A1 EP12784529.5A EP12784529A EP2770911A1 EP 2770911 A1 EP2770911 A1 EP 2770911A1 EP 12784529 A EP12784529 A EP 12784529A EP 2770911 A1 EP2770911 A1 EP 2770911A1
Authority
EP
European Patent Office
Prior art keywords
data set
values
optimized
time
measurement
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
Application number
EP12784529.5A
Other languages
English (en)
French (fr)
Inventor
Frank-Detlef Scholle
Joachim Hütter
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Life Molecular Imaging SA
Original Assignee
Piramal Imaging SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Piramal Imaging SA filed Critical Piramal Imaging SA
Publication of EP2770911A1 publication Critical patent/EP2770911A1/de
Withdrawn legal-status Critical Current

Links

Classifications

    • 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/508—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Definitions

  • tomographic methods allow the generation of sectional images and three-dimensional representations (3D images).
  • a sectional image reflects the internal structures of the examined body as they were after cutting out a thin layer.
  • a 3D representation shows how the examined structures are spatially present.
  • CT computed tomography
  • X-ray absorption profiles of the body to be examined are generated from many directions. From these absorption profiles, the degree of absorption can then be calculated for each volume element of the body and sectional images and 3D representations can be constructed.
  • positron emission tomography is based on the detection of a variety of annihilation events.
  • the more events that are registered the higher the number of data used for reconstruction and the higher the signal-to-noise ratio.
  • the number of events can in principle be influenced both by the amount of tracer administered and the duration of the scan.
  • the burden of the body with radioactive substances should, however, be kept as low as possible in order to avoid side effects. To minimize side effects, the amount of tracer administered should therefore be minimized.
  • the extension of the scan duration is also limited.
  • the examined body area should not move during the recording, since movements in the recordings lead to a false representation of the tracer distribution.
  • Motionless persistence is a burden on the patient. Some movements, such as respiratory movements and cardiac muscle movements, can not be avoided when taking measurements on living organisms.
  • restrict factors such as the half-life of the radioactive isotopes of the tracer and / or the degradation of the tracer in the body its temporal detectability and / or informative value
  • the aim of the development is to provide a tracer that delivers specific biochemical information to the body under investigation with a high signal-to-noise ratio and low body stress. Any increase in the signal-to-noise ratio resulting from improvements in the measurement and containment technique would be a valuable contribution that could help minimize the burden on the body of a tracer.
  • a first subject of the present invention is a method for generating optimized tomography images, comprising at least the steps:
  • the method according to the invention comprises at least the following steps: a) provision of a data record which represents an area in a body during a measuring time, the representation of the body area in the data set being subdivided into a plurality of discrete partial areas, the measuring time in the data set being in a plurality is subdivided by discrete measurement intervals, wherein each subregion is assigned a discrete structure value for each measurement interval; b) establishing boundary conditions about the expected time course of a structure size in the region of the body during the measurement time; c) calculating optimized structure values for each individual subarea on the basis of structural values of the individual subarea at temporally successive measurement intervals taking into account the boundary conditions; d) Output of an optimized data set which represents the body or an area in the body at arbitrary times within the measuring time and which is based on the optimized structure values.
  • the method according to the invention generates from a first data set, which represents an area in a body during a measuring time, a second, optimized data set, which represents an area in the body during freely selectable times within the measuring time
  • representations of the body area can be generated at freely selectable times within the measuring time
  • Step b) can take place before or after step a), i. the designation of the steps with a) and b) does not necessarily mean that first step a) and then step b) takes place
  • the boundary conditions determine the laws that follow the temporal course of the structure size in the area of the body.
  • the temporal course of the structure size is not arbitrary but it inevitably follows the laws, for example, by the anatomy, morphology and / or physiology of the body area and use a tracer or contrast agent are determined by the physical and chemical properties of the tracer or contrast agent. For example, it is extremely unlikely that the absorbance in the computer tomography of a patient increases and decreases oscillatory as a structure size after a single application of a contrast agent
  • step c) of the inventive method optimized structural values are calculated for each individual subsection.
  • Step c) requires the presence of a first data set and of boundary conditions, so that step c) can take place only after the steps a) and b).
  • the calculation is based on the measured structural values and taking into account the boundary conditions.
  • measured structural values are correlated to temporally successive measurement intervals
  • the sections must contain at least one measuring interval. At z. As the computed tomography or magnetic resonance imaging this is to be considered in the measurement of the data set.
  • c2) Means of the structure values within each section, if there is more than one measurement time range in the selected temporal section. Alternatively, instead of the averaging in a section, a corresponding data record with the time length of the considered section can be reconstructed, as for example possible in the case of PET. c3) fitting a compensation curve into the averaged structure values, the compensation curve providing optimized structure values
  • the size of the sections is adapted to the existing measured structure values In the areas of the measuring time in which large changes in the structure values are to be found, the sections are shorter than in the areas of the measuring time in which the structure values are less from one measuring interval to the next measuring interval change strongly.
  • the decisive factor is therefore the first derivation of the structure values according to time. The larger this is, the shorter the sections are.
  • each section is inversely proportional to the amount of first derivative of the structure values by time.
  • the sections can be chosen such that two sections each adjoin one another; it is also conceivable to design the sections in such a way that two or more sections each overlap.
  • the sections are designed such that in each case two temporally successive sections overlap in their edge regions. In a particularly preferred embodiment, two temporally successive sections overlap each in one edge point
  • Averaging is understood to mean the formation of known mathematical mean values such as, for example, the arithmetic or geometric or harmonic or quadratic mean or weighted mean.
  • the choice of the respective mean value depends above all on the considered structure size and the existing boundary conditions.
  • the arithmetic mean is formed.
  • the average values are preferably assigned to the middle of the respective time segment, so that a mean value curve results which represents the average structural values as a function of time.
  • a mean value curve results which represents the average structural values as a function of time.
  • a compensation curve is fitted in the mean value curve.
  • the compensation curve is selected on the basis of the boundary conditions set up in step b) of the method according to the invention.
  • the compensation curve is adjusted so that the deviations between the mean value curve and the compensation curve are as small as possible. It is also a weighted adjustment conceivable. Weighting means that the compensation curve in the area of the higher-weighted structure values may have a smaller deviation from the mean value curve than in the area of the lower-weighted structure values.
  • spline functions are suitable as compensation curves.
  • apart from recirculation peaks for example, a global maximum for the application of a tracer or contrast agent and, if appropriate, in each case a local maximum at zJB. present extravasation, leakage in tumors, specific or unspecific enrichment in the mathematical function.
  • the beginning of the curve can be extrapolated with the aid of the slope of the first two mean values.
  • the compensation curve provides optimized structure values at arbitrary times within the measurement interval, since the compensation curve represents a continuous time curve and does not consist of discrete values
  • the result is a data set with optimized structure values at freely selectable times within the measurement interval. Due to the boundary conditions taken into account, in the optimized data record that has been taken there is information which makes it possible to specifically highlight or suppress morphological and / or physiological structures within the data record.
  • a mathematical model is used to calculate the optimized feature width in step c).
  • This embodiment of the method according to the invention comprises the following steps: c1) providing a mathematical model which describes the temporal behavior of the structure value in the regions of the body; c2) for each subarea: adapting at least one parameter of the model to the measured structure values and determining a model function which optimally reproduces the temporal course of the measured structure values as a result of a mathematical optimization method, the model function being optimized
  • the mathematical model represents the boundary conditions that have been set up in step b) of the method according to the invention.
  • auxiliary means such as e.g. a tracer or contrast agent - preferably a single or multi-compartment model.
  • the considered Kötper Scheme is considered as a built up of one or more compartments body.
  • a compartment in the model is used for each temporal change of the structure value. For example, after a bolus application in the blood of a patient, a tracer distributes itself in a manner and speed characteristic of the patient and the tracer and is gradually eliminated and, if necessary, metabolized.
  • a compartment in the model function is to be provided.
  • various mathematical methods can be used. For example, a model function can be obtained by solving the differential equations that can be established for the model, as is the case in pharmacokinetic modeling.
  • model function can also be obtained by simulating the time evolution of the considered structure values over the measurement time.
  • mode functional parameters By varying the mode functional parameters, a mathematical adaptation of the model function to the temporal behavior of the structure values is possible.
  • the determination of a model function by adaptation to a mathematical model is preferably carried out in the inventive method using the simulation approach.
  • the result is a model function that optimally reproduces the temporal behavior of the structure values in the mathematical sense.
  • the model function provides optimized structure values at arbitrary times within the measurement interval the model function represents a continuous time curve and does not consist of discrete values. For each subarea of the scanned body, a set of optimized parameters, which indicates the influence of each compartment on the temporal course of the structure value, results from the named method variant
  • the contrasting of the vascular system in the output data set can be suppressed or highlighted as required.
  • the result of the model adaptation is thus a data set with optimized structure values and a data set with associated mode parameters, with which the optimized data record can be output in various variants useful for understanding the examination data.
  • step d) of the method according to the invention the output of an optimized data set takes place.
  • the optimized data set represents an area in the examined body.
  • the range in step d) coincides with the range in step a).
  • the area in step d) represents only a partial area of the area from step a). It is conceivable that partial areas have been discarded in the course of or following the calculation of the optimized structural values in step c) or by a motion correction , This applies in particular to edge regions of the data set that may not spatially coincide due to motion in all measurement time intervals.
  • step d) can take place only after step c).
  • the optimized data set may be output in the form of one or more two- or three-dimensional representations of the area of the body on a screen or as an expression. Likewise, it is conceivable that the output takes place on a data carrier in the form of machine-readable data.
  • the optimized data set which has been produced by means of the method according to the invention is likewise an object of the present invention.
  • a further subject of the present invention is a computer program product with program code stored on a machine-readable carrier for carrying out the method according to the invention on a computer.
  • the method according to the invention is suitable for optimizing all known 3D images or tomography images, for example for optimizing SPECT, PET, CT or MRT images, or measurement data from a 3D or 4D ultrasound method or optical tomography ( See relevant literature such as: Ashok Kliurana, Nirvikar Dahiya: 3D & 4D Vitrasound - A Text and Atlas, Jaypee Brothers Medical Pubiishers (P) Ltd., 2004, R.
  • Movements which occur during the measuring time in the scanned body or in subregions of the scanned body are reduced by the method according to the invention in many cases, which is particularly advantageous in highly noisy data sets image blurring, as they are unavoidable in static recordings with only one record per total measurement time are reduced by the method according to the invention and the spatial resolution is closer to the physically possible resolution of the scanning device.
  • the structure values per section located in the different temporal sections are averaged and corrected according to the selected boundary conditions for a main maximum and at most a secondary maximum, if necessary, in the amount of the value.
  • the slightly higher average of the penultimate Section (Minute 44-52) down to the mean of the third last section (minute 36-44) down, since there may be no further maximum in the curve due to the boundary conditions except the much larger secondary maximum at less than 20 minutes.
  • FIGS. 2 to 4 show, by way of example, a section of a measurement data set in the anatomically customary planes.
  • FIG. 2 shows the data set without processing using the method according to the invention.
  • the noise reduction with the method according to the invention can be seen in FIG. 3 on the basis of structures which can be easily recognized and substantially fewer individual spots.
  • FIG. 4 the structure recognizable in FIG. 3 is confirmed.
  • the dataset illustrated in FIG. 4 does not allow any conclusions to be drawn about the kinetics of the tracer distribution in the scanbody, in contrast to the dataset from FIG. figure description
  • Figure 1 Representation of an exemplary time course of the tracer concentration during an in vivo PET scan in a discrete portion of a PET data set
  • the section bars in FIG. 1b are each entered at the level of the value obtained from the section average.
  • the start of the PET scan was done immediately after application of the tracer
  • Figure 2 representation of the anatomical views
  • the scan was performed on a C nomolgus monkey after application of a thrombus tracer from PET tracer research with a small animal PET scanner. Shown is the measurement data record number 28 of 60 successive scans without noise reduction by the inventive method. The measurement duration of each measurement data set is 1 minute. The measurement of all data records took place one after the other without a break.
  • the planes for the views shown are identical to those of Figure 3a-c and Figure 4a-c.
  • the crosses recognizable in the figures represent the cursor position in the computer program product according to the invention with which the figures were created.
  • Figure 3 representation of the anatomical views
  • the scan was performed on a cynomolgus monkey after application of a thrombus tracer from PET tracer research with a small animal PET scanner. Shown is the measurement data set number 28 of 60 successive scans after application of the inventive method. The measurement duration of each measurement data set is 1 minute. The measurement of all data records took place one after the other without a break.
  • the planes for the views shown are identical to those of Figure 2a-c and Figure 4a-c.
  • the crosses recognizable in the figures represent the cursor position in the computer program product according to the invention with which the figures were created.
  • the scan was performed on a cynomolgus monkey after application of a thrombus tracer from PET tracer research with a small animal PET scanner. Shown is the averaging of all 60 individual data sets that were scanned during the total measurement time. The measurement duration of each measurement data set is 1 minute. The measurement of all data records took place one after the other without a break. The individual data sets have not been processed by the method according to the invention.
  • the planes for the views shown are identical to those of Figure 2a-c and Figure 3a-c.
  • the crosses recognizable in the figures represent the cursor position in the Computeiprograinm Cool invention, with which the figures were created.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Nuclear Medicine (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Image Processing (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Image Analysis (AREA)

Abstract

Die vorliegende Erfindung betrifft das technische Gebiet der bildgebenden Verfahren insbesondere zu diagnostischen Zwecken. Gegenstand der vorliegenden Erfindung ist ein Verfahren zur Erzeugung von optimierten Tomografie-Autnahmen, ein Computerprogrammprodukt zur Ausführung des erfindungsgemäßen Verfahrens auf einem Computer sowie die mittels des erfindungsgemäßen Verfahrens erzeugten optimierten Aufnahmen.

Description

Verfahren zur Erzeugung optimierter Tomografie-Aufnahmen
Die vorliegende Erfindung betrifft das technische Gebiet der bildgebenden Verfahren insbesondere zu diagnostischen Zwecken. Gegenstand der vorliegenden Erfindung ist ein Verfahren zur Erzeugung von optimierten Tomografie-Aufnahmen, ein Computerprogrammprodukt zur Ausführung des erfmdungs gemäßen Verfahrens auf einem Computer sowie die mittels des erfindungs gemäßen Verfahrens erzeugten optimierten Aufnahmen.
In der heutigen Medizin kommen verschiedene bildgebende Verfahren zum Einsatz, um anatomische und funktionelle Strukturen bei lebenden Menschen oder Tieren sichtbar zu machen und so den Gesundheitszustand zu beurteilen. Im Gegensatz zu den Projektionsverfahren wie beispielsweise der gewöhnlichen Röntgenaufnahme, bei der sich Strukturen, die im Strahlengang der Röntgenstrahlen hintereinander liegen, in der Aufnahme überlagern, erlauben tomo grafische Verfahren die Erzeugung von Schnittbildern und dreidimensionalen Darstellungen (3D-Bilder). Ein Schnittbild gibt die inneren Strukturen des untersuchten Körpers so wieder, wie sie nach dem Herausschneiden einer dünnen Schicht vorlägen. Eine 3D-Darstellung zeigt, wie die untersuchten Strukturen räumlich vorliegen.
Bei der Computertomografie (CT) werden beispielsweise Röntgenabsorptionsprofiie des zu untersuchenden Körpers aus vielen Richtungen erzeugt. Aus diesen Absorptionsprofilen können dann für jedes Volumenelement des Körpers der Absorptionsgrad errechnet und Schnittbilder sowie 3D-Darstellungen konstruiert werden.
Während mittels Computertomografie, gegebenenfalls unter Einsatz von Kontrastmitteln, die morphologische/anatomische Struktur eines Körpers dargestellt werden kann, erlaubt beispielsweise die Positronen-Emissions-Tomografie (PET) die Darstellung von biochemischen Funktionalitäten eines Organismus. Bei der PET wird hierzu ein radioaktiv markierter Tracer in den Körper eines Patienten appliziert. Der Tracer bindet selektiv an bestimmte Biomoleküle und durch Aufnahme der vom Tracer abgegebenen Strahlimg kann die Aktivität der Biomoleküle im Körper sichtbar gemacht werden. Nach der Verabreichung eines Tracers dauert es einige Zeit, bis der Tracer eine gewünschte Verteilung im Körper erreicht hat. Üblicherweise wird der Tracer intravenös verabreicht und gelangt somit über die Blutbahn an das gewünschte Ziel. Ein Teil der verabreichten Tracer- Moleküle bindet spezifisch an die gewünschten Zielbereiche, ein anderer Teil verteilt sich unspezifisch. Um Tomografie-Aufhahmen mit einem hohen S ignal-Rausch- Verhältnis zu erhalten, ist es oftmals sinnvoll nach der Verabreichung mit den Aufnahmen zu warten, bis ein Großteil der unspezifisch bindenden oder verteilten Tracer-Moleküle den zu untersuchenden Körperbereich wieder verlassen hat, da die nicht spezifisch bindenden Tracer- Moleküle zum Hintergrundsignal in den PET-Aumahmen beitragen.
In Abhängigkeit vom verwendeten Tracer und den physiologischen Parametern des untersuchten Patienten gibt es ein Zeitfenster nach der Verabreichung des Tracers und vor seinem Abtransport aus dem betrachteten Körperbereich oder seinem metabolischen Abbau, in dem ein optimales Signal-Rausch-Verhältnis erzielt werden kann.
Die Aufnahme von PET-Scans nimmt eine gewisse Zeit in Anspruch, denn die Positronen- Emissions-Tomografie basiert auf der Detektion einer Vielzahl von Annihilaitonsereignissen. Je mehr Ereignisse registriert werden, desto höher ist die Zahl der zur Rekonstruktion verwendeten Daten und desto höher ist das S ignal-Raus ch- Verhältnis . Die Zahl der Ereignisse kann prinzipiell sowohl durch die Menge des verabreichten Tracers als auch die Scan-Dauer beeinflusst werden.
Die Belastung des Körpers mit radioaktiven Substanzen sollte jedoch möglichst gering gehalten werden, um Nebenwirkungen zu vermeiden. Zur Minimierung von Nebenwirkungen sollte die Menge des verabreichten Tracers daher möglichst gering gehalten werden.
Auch der Ausdehnung der Scan-Dauer sind Grenzen gesetzt. Zum einen sollte sich der untersuchte Körperbereich wahrend der Aufnahme nicht bewegen, da Bewegungen in den Aufnahmen zu einer falschen Darstellung der Tracer-Verteilung führen. Das regungslose Verharren stellt für den Patienten jedoch eine Belastung dar. Manche Bewegungen wie beispielsweise Atembewegungen und Bewegungen des Herzmuskels können bei Messungen an lebenden Organismen nicht vermieden werden. Zum anderen beschränken Faktoren wie die Halbwertszeit der radioaktiven Isotope der Tracer und/oder der Abbau des Tracers im Körper seine zeitliche Detektierbariceit und/oder Aussagekraft
Bei der Entwicklung eines neuen Tracers spielen viele unterschiedliche Faktoren eine Rolle, Ziel der Entwicklung ist die Bereitstellung eines Tracers, der eine spezifische biochemische Information Ober den untersuchten Körper mit hohem Signal-Rausch- Verhältnis und bei geringer Belastung des Körpers liefert. Dabei wäre jede Erhöhung des Signal-Rausch- Verhältnisses, die durch Verbesserungen in der Mess- und Aufhahmetechnik zustande kommt, ein wertvoller Beitrag, der dazu fuhren kann, die Belastung des Körpers durch einen Tracer zu minimieren.
Die obigen Überlegungen gelten in ahnlicher Weise auch für weitere Tomo grafic- Verfahren, insbesondere für solche Verfahren, bei denen Hilfsmittel zur Signalerzeugung oder Signalverstärkung wie beispielsweise Tracer, Kontrastmittel oder Fluoreszenzfarbstoffe dem zu untersuchenden Körper verabreicht werden.
Es wäre wünschenswert, Tomografie-Aufoahmen mit einem hohen Signal-Rausch-Verhältnis erzeugen zu können, wobei die Belastung des untersuchten Patienten sowohl in Hinblick auf die dem Körper ausgesetzte Strahlendosis und/oder die Menge an appliziertem Hilfsmittel als auch in Hinblick auf die Untersuchungsdauer zu minimieren ist
Die bisherigen Betrachtungen bezogen sich vornehmlich auf die Erzeugung von statischen
Momentaufnahmen anatomischer und/oder funktioneller Strukturen.
Sie gelten jedoch in besonderem Maße auch für die zeitliche Verfolgung von Vorgängen in einem Körper, wobei Körper hier sowohl den Körper eines Menschen oder Tieres als auch ein lebloses Objekt wie beispielsweise ein Messphantom oder eine Materialprobe umfasst Bei der Erzeugung von Aufnahmen, die das dynamische Verhalten eines applizierten Hilfsmittels in einem Körperbereich darstellen, werden Messungen am betrachteten Körperbereich über einen längeren Zeitraum vorgenommen. Hieraus lassen sich wertvolle Informationen zum zeitlichen Ablauf von physiologischen Prozessen gewinnen.
Anschließend werden die Messdaten in mehrere Zeitbereiche unterteilt, die Signalintensitäten in jedem Volumenelement für die einzelnen Zeitbereiche ermittelt und eine Signalintensitäts- Zeit-Kurve erstellt Hierbei tritt das Problem auf, dass die Unterteilung der Gesamtmesszeit in zunehmend kürzere Abschnitte zwar zu einer zunehmend höheren zeitlichen Auflösung führt, die Verkürzung der Zeitbereiche jedoch ein stärker verrauschtes Signal zur Folge hat.
Man erhält also entweder eine hohe räumliche Auflösung mit geringem Rauschen bei geringer oder fehlender zeitlicher Information oder eine hohe zeitliche Auflösung bei geringer räumlicher Auflösung.
Es wäre daher wünschenswert, den durch die Erhöhung der zeitlichen Auflösung bedingten Verlust an räumlichem Auflösevermögen zumindest teilweise kompensieren zu können.
Erfindungsgemäß werden die genannten Aufgaben durch die Verknüpfung der räumlichen Messdaten mit dazugehörigen zeitlichen Informationen unter Berücksichtigung physiologischer Randbedingungen gelöst
Ein erster Gegenstand der vorliegenden Erfindung ist ein Verfahren zur Erzeugung optimierter Tomografie- Aufnahmen, mindestens umfassend die Schritte:
a) Bereitstellen eines Datensatzes, der einen Bereich im Körper eines Patienten während einer Messzeit repräsentiert, wobei die Repräsentation des Körperbereichs im Datensatz in eine Mehrzahl von diskreten Teilbereichen unterteilt ist, wobei die Messzeit im Datensatz in eine Mehrzahl von diskreten Messintervallen unterteilt ist, wobei jedem Teilbereich zu jedem Messintervall ein diskreter Stnikturwert zugeordnet ist;
b) Aufstellen von Randbedingungen über den zu erwartenden zeitlichen Verlauf einer Strukturgröße in dem Bereich des Körpers während der Messzeit; c) Berechnen von optimierten Strukturwerten für jeden einzelnen Teilbereich auf Basis von Strukturwerten des einzelnen Teilbereichs zu zeitlich aufeinander folgenden Messintervallen unter Berücksichtigung der Randbedingungen; d) Ausgabe eines optimierten Datensatzes, der einen Bereich im Körper zu beliebig wählbaren Zeitpunkten während der Messzeit repräsentiert und der auf den optimierten Strukturwerten basiert Unter einer Tomo grafie- Aufnahme wird ein Satz an Daten verstanden, der einen Bereich in einem Körper während einer Zeitspanne repräsentiert Der Begriff Tomografie- Aufnahme soll nicht auf Schnittbildcr beschrankt sein sondern soll auch Datensätze umfassen, die einen Körperbereich dreidimensional darstellen. Die Repräsentation des Körperbereichs erfolgt auf Basis einer Strukturgröße und entsprechenden Strukturwerten, die weiter unten naher beschrieben sind.
Das erfindungsgemäße Verfahren umfasst mindestens die folgenden Schritte: a) Bereitstellen eines Datensatzes, der einen Bereich in einem Körper während einer Messzeit repräsentiert, wobei die Repräsentation des Körperbereichs im Datensatz in eine Mehrzahl von diskreten Teilbereichen unterteilt ist, wobei die Messzeit im Datensatz in eine Mehrzahl von diskreten Messintervallen unterteilt ist, wobei jedem Teilbereich zu jedem Messintervall ein diskreter Strukturwert zugeordnet ist; b) Aufstellen von Randbedingungen über den zu erwartenden zeitlichen Verlauf einer Strukturgröße in dem Bereich des Körpers während der Messzeit; c) Berechnen von optimierten Strukturwerten für jeden einzelnen Teilbereich auf Basis von Strukturwerten des einzelnen Teilbereichs zu zeitlich aufeinander folgenden Messintervallen unter Berücksichtigung der Randbedingungen; d) Ausgabe eines optimierten Datensatzes, der den Körper oder einen Bereich im Körper zu beliebig wählbaren Zeitpunkten innerhalb der Messzeit repräsentiert und der auf den optimierten Strukturwerten basiert.
Das erfindungsgemäße Verfahren erzeugt aus einem ersten Datensatz, der einen Bereich in einem Körper während einer Messzeit repräsentiert, einen zweiten, optimierten Datensatz, der einen Bereich in dem Körper während frei wählbarer Zeitpunkte innerhalb der Messzeit repräsentiert
Der zweite, optimierte Datensatz ist durch folgende Punkte gekennzeichnet: der Rauschanteil ist gegenüber dem ersten Datensatz reduziert,
Bildunscharfen, wie sie in zeitlich längeren Scans unumgänglich sind, sind reduziert, und die räumliche Auflösung liegt näher an der physikalischen Auflösung des
Scangerätes,
Verschiebungen, Stauchungen, Dehnungen, Rotationen, usw., die im ersten Datensatz während der Messzeit enthalten sein können, sind in der Regel vermindert,
es lassen sich Darstellungen des Körperbereichs zu frei wählbaren Zeitpunkten innerhalb der Messzeit erzeugen,
es lassen sich gezielt morphologische und/oder physiologische Funktionen hervorheben oder unterdrücken.
Der erste Datensatz resultiert aus Messungen, die an einem menschlichen oder tierischen oder einem sonstigen Köper vorgenommen worden sind. Vorzugsweise sind die Messungen an einem lebenden Organismus vorgenommen worden.
Der erste Datensatz ist beispielweise eine Sequenz von PET-Rekonstruktionen, von CT- Aufhahmen, von Magnetresonanztomogiafie-Aufhahmen (MRT-Aufhahmen) oder vergleichbaren Aufnahmen. Jede einzelne Aufnahme ist innerhalb eines Messintervalls entstanden. Die Sequenz zeigt die Aumahmen zu aufeinander folgenden Zeitabschnitten bzw. Messintervallen. Die Begriffe„Sequenz" und„zeitliche Abfolge" werden hier synonym verwendet.
Alle Messintervalle zusammengefasst ergeben die Messzeit
Beim ersten und zweiten Datensatz kann es sich um eine dreidimensionale Repräsentation handeln. Es kann sich jedoch auch um eine zweidimensionale Repräsentation, also um ein Schnittbild, handeln. Unabhängig davon, ob es sich um eine zwei- oder dreidimensionale Repräsentation handelt, wird im Folgenden auch von der Repräsentation eines räumlichen Bereichs gesprochen. Die Repräsentation des räumlichen Bereichs im Datensatz ist quantisiert, das heißt, der räumliche Bereich ist in eine diskrete Anzahl von Teilbereichen (Flächenelemente oder Volumenelemente) aufgeteilt, wobei jeder einzelne Teilbereich durch seine Koordinaten im Raum charakterisiert ist Die Koordinaten im Raum sollten sich idealerweise während der Messzeit nicht ändern. Sie ändern sich dann nicht, wenn der Bereich des Körpers bei der Aufnahme der Messwerte zur Erzeugung des ersten Datensatzes während der Messzeit in Bezug zur Messvorrichtung nicht bewegt worden ist. Zunächst wird der Einfachheit halber angenommen, dass während der Messzeit weder eine Bewegung des Bereichs noch Bewegungen innerhalb des Bereichs des Körpers stattgefunden haben, so dass die Koordinaten der einzelnen Teilbereiche wahrend der Messzeit konstant sind.
Den einzelnen Teilbereichen ist zu jedem Messintervall ein Strukturwert zugeordnet Die Strukturwerte charakterisieren den Zustand des Teilbereichs innerhalb des betrachteten Messintervalls. Der Zustand eines jeden Teilbereichs ist durch eine Reihe von Größen bestimmt Mindestens eine Größe, die hier als Strukturgröße bezeichnet wird, wird in dem erfindungsgemäßen Verfahren betrachtet Es ist auch denkbar, mehrere Größen zu betrachten. Strukturgrößen können beispielsweise Größen wie Röntgenabsorption (CT), Anzahl der Zerfallsereignisse pro Zeit (PET), MR-Relaxationszeiten usw. sein. Um die obigen Definitionen näher zu verdeutlichen, seien beispielhaft die Computertomografie und die Positronen-Emis sions-Tomo grafie angeführt. Computertomografische Aufnahmen sind räumliche Datensätze, die aus einer diskreten Anzahl von Volumenelementen aufgebaut sind, wobei jedes einzelne Volumenelement durch seine Koordinaten im Raum und einen Absorptionswert charakterisiert ist Üblicherweise stellt der Absorptionsweit eine Graustufe dar, wobei beispielsweise „schwarz" den niedrigsten Absorptionsgrad (Graustufe 0) und„weiß" den höchsten Absorptionsgrad (z.B. bei 100 Graustufen die Graustufe 99) darstellt Dadurch lassen sich die räumlichen Datensätze bildlich darstellen. Die Strukturgröße, die im Fall der CT betrachtet wird, ist der Absorptionsgrad des Gewebes für Röntgenstrahlung.
Im Fall von PET werden die Zerfalle der verwendeten Radionuklide über der Messzeit detektiert. Die räumlichen Datensätze können dann für beliebige, die gesamte Messzeit unterteilende Zeitintervalle, rekonstruiert werden. Jedes einzelne Volumenelement ist hierbei durch seine Koordinaten im Raum und eine Zerfallsrate charakterisiert. Das erfindungsgemäße Verfahren erfordert mehrere räumliche Datensätze, die jeweils den Zustand des untersuchten Körperbereichs in einem zeitlichen Abstand voneinander darstellen. Der zeitliche Abstand voneinander kann gleichbleibend oder variabel sein; wichtig ist, dass der zeitliche Abstand voneinander und die Zeitdauer bei den einzelnen Datensätzen bekannt sind. Des Weiteren sind die zeitlichen Abstände und die Zeitdauern entweder bei der Messung oder, wie im Fall von PET, bei der Rekonstruktion so zu wählen, dass die interessierenden zeitlichen Änderungen des betrachteten Strukturwertes zeitlich aufgelöst werden. Die zeitlichen Abstände und die Zeitdauern sollten also kleiner sein als die betrachteten zeitlichen Änderungen des Strukturwertes.
Schritt a) des erfindungsgemäßen Verfahrens stellt die Bereitstellung eines ersten Datensatzes dar. Da dieser Datensatz aus Messungen resultiert, d.h. empirisch generiert wurde, weist er einen Rauschanteil auf. Insbesondere PET-Aufhahmen weisen einen signifikanten Rauschanteil aufgrund der Statistik der Zerfallsereignisse auf, der umso höher ist, je kürzer der Zeitabschnitt ist, währenddessen Annihilationsereignisse registriert werden, um eine PET-Aufnahme zu erzeugen.
Die Reduktion des Rauschanteils gelingt erfindungsgemäß durch Verknüpfen der räumlichen Messdaten mit den zugehörigen zeitlichen Informationen unter Berücksichtigung physiologischer Randbedingungen.
Diese Randbedingungen werden in Schritt b) des erfindungs gemäßen Verfahrens aufgestellt Schritt b) kann zeitlich vor oder nach Schritt a) erfolgen, d.h. die Bezeichnung der Schritte mit a) und b) bedeutet nicht zwangsläufig, dass zuerst Schritt a) und dann Schritt b) erfolgt
Die Randbedingungen legen die Gesetzmäßigkeiten fest, denen der zeitliche Verlauf der Strukturgröße in dem Bereich des Körpers folgt Der zeitliche Verlauf der Strukturgröße ist nicht beliebig sondern er folgt zwangsläufig den Gesetzmäßigkeiten, die beispielsweise durch die Anatomie, Morphologie und/oder Physiologie des Körperbereichs und bei Einsatz eines Tracers oder Kontrastmittels durch die physikalischen und chemischen Eigenschaften des Tracers oder Kontrastmittels festgelegt sind. So ist es beispielsweise äußerst unwahrscheinlich, dass der Absorptionsgrad bei der Computertomografie eines Patienten als Strukturgröße nach einmaliger Applikation eines Kontrastmittels oszillatorisch zu- und abnimmt
Wird ein Tracer oder Kontrastmittel verabreicht, so wird dieses in den betrachteten Körperbereich eintreten und diesen nach einer Verweilzeit wieder verlassen. Sieht man von Rezirkulationspeaks ab, sollte die messtechnische Verfolgung des Tracers oder Kontrastmittels daher einen Signalanstieg mit anschließendem Signalabfall zeigen (Hauptmaximum). Zusätzlich kann es jeweils höchstens zu einem weiteren Signalanstieg mit anschließendem Signalabfall aufgrund von z.B. Extravasation, Leakage in Tumoren, spezifischer oder unspezifischer Anreicherung kommen (Nebenmaxium), wobei das Nebenmaximum dem Hauptmaximum zeitlich nachgelagert ist
Die Randbedingungen legen demnach fest, in welchen Grenzen sich ein Strukturweit bewegen kann und welche zeitlichen Änderungen des Strukturwerts mit Naturgesetzen vereinbar sind. Randbedingungen können beispielsweise sein:
Zeitkonstante des Tracers oder Kontrastmittels in der betrachteten Spezies für die Verdünnung im Blutvolumen nach Applikation
Zeitkonstante des Tracers oder Kontrastmittels in der betrachteten Spezies für die Elimination aus dem Blut
- Typische Zeitverläufe für die Konzentration eines Tracers oder Kontrastmittels.
Beispielsweise kann es nach Applikation des Tracers oder Kontrastmittels nur einen
Signalanstieg mit anschließendem Abfall im Gefäßanteil in vivo geben sowie zusätzlich jeweils höchstens einem Anstieg und Abfall aufgrund z.B.
Extravasation (wenn Tracer oder Kontrastmittel klein genug sind, um Gefäßwände zu durchdringen), Leakage in Tumoren, spezifische oder unspezifische
Anreicherung, usw.
Diese Zeitverläufe können auch durch eine pharmakokinetische Modellfuiiktion beschrieben werden. In Schritt c) des erfmdungsgemäßen Verfahrens werden für jeden einzelnen Teilbereich optimierte Strukturwerte berechnet Schritt c) erfordert das Vorhandensein eines ersten Datensatzes und von Randbedingungen, so dass Schritt c) zeitlich nur nach den Schritten a) und b) erfolgen kann. Die Berechnung erfolgt auf Basis der gemessenen Strukturwerte und unter Berücksichtigung der Randbedingungen. Für die Berechnung der optimierten Strukturwerte werden gemessene Strukturwerte zu zeitlich aufeinanderfolgenden Messintervallen miteinander in Beziehung gesetzt
Die Berechnung kann auf verschiedene Weisen durchgeführt werden. Nachfolgend werden zwei bevorzugte Ausführungsformen näher beschrieben. /. Abschnittsweise Glättung
In einer ersten bevorzugten AusfÜhrungsform des erfindungsgemäßen Verfahrens werden ftir jeden einzelnen Teilbereich die folgenden mathematischen Operationen durchgeführt: cl) Unterteilen der Messzeit in eine Mehrzahl von Abschnitten, wobei die einzelnen Abschnitte umso kürzer sind, je größer die Änderung der Strukturwerte in einem
Bereich der Messzeit ist Die Abschnitte müssen dabei mindestens ein Messintervall enthalten. Bei z. B. der Computertomografie oder der Magnetresonanztomografie ist dieses bei der Messung des Datensatzes zu berücksichtigen. c2) Mitteln der Strukturwerte innerhalb eines jeden Abschnitts, sofern mehr als ein Messzeitbereich im gewählten zeitlichen Abschnitt liegt Alternativ kann statt der Mittelung in einem Abschnitt auch ein entsprechender Datensatz mit der zeitlichen Länge des betrachteten Abschnitts rekonstruiert werden, wie beispielsweise im Fall von PET möglich. c3) Einpassen einer Ausgleichskurve in die gemittelten Strukturwerte, wobei die Ausgleichskurve optimierte Strukturwerte liefert
Die Schritte cl) bis c3) erfolgen nacheinander in der angegebenen Reihenfolge. In Figur 1 wird die Berechnung bildlich verdeutlicht und in dem unten beschriebenen Beispiel näher erläutert.
Die Größe der Abschnitte wird an die vorliegenden gemessenen Strukturwerte angepasst In den Bereichen der Messzeit, in denen große Änderungen der Strukturwerte zu verzeichnen sind, sind die Abschnitte kürzer als in den Bereichen der Messzeit, in denen sich die Strukturwerte von einem Messintervall zum nächsten Messintervall weniger stark ändern. Ausschlaggebend ist demnach die erste Ableitung der Strukturwerte nach der Zeit Je größer diese ist, desto kürzer sind die Abschnitte.
Vorzugsweise ist die Größe eines jeden Abschnitts umgekehrt proportional zum Betrag der ersten Ableitung der Strukturwerte nach der Zeit. Die Abschnitte können so gewählt werden, dass jeweils zwei Abschnitte aneinander angrenzen; ebenso ist es denkbar, die Abschnitte so zu gestalten, dass jeweils zwei oder mehr Abschnitte sich überlappen. Vorzugsweise werden die Abschnitte so gestaltet, dass sich jeweils zwei zeitlich aufeinander folgende Abschnitte in ihren Randbereichen überlappen. In einer besonders bevorzugten AusfÜhrungsform überlappen sich zwei zeitlich aufeinander folgende Abschnitte in jeweils einem Randpunkt
Sobald die Abschnitte festgelegt sind, erfolgt eine Mittelung der innerhalb eines jeden Abschnittes liegenden Strukturwerte. Als Mittelung wird die Bildung bekannter mathematischer Mittelwerte wie beispielsweise des arithmetischen oder geometrischen oder harmonischen oder quadratischen Mittelwerts oder gewichteten Mittels verstanden. Die Wahl des jeweiligen Mittelwertes richtet sich vor allem nach der betrachteten Strukturgroße und den vorliegenden Randbedingungen. Üblicherweise wird der arithmetische Mittelwert gebildet.
Die Mittelwerte werden vorzugsweise der Mitte des jeweiligen Zeitabschnitts zugeordnet, so dass sich eine Mittelwertkurve ergibt, welche die mittleren Strukturwerte als Funktion der Zeit darstellt Es ist aber auch denkbar, die Mittelwerte jeweils dem ersten oder letzten oder einem anderen Zeitpunkt des entsprechenden Zeitabschnitts zuzuordnen.
In die Mittelwertkurve wird eine Ausgleichskurve eingepasst Die Ausgleichskurve wird auf Basis der Randbedingungen, die in Schritt b) des erfindungsgemäßen Verfahrens aufgestellt worden sind, gewählt Die Ausgleichskurve wird so eingepasst, dass die Abweichungen zwischen der Mittelwertkurve und der Ausgleichskurve möglichst gering sind. Es ist auch eine gewichtete Anpassung denkbar. Unter Gewichtung wird verstanden, dass die Ausgleichskurve im Bereich der höher gewichteten Strukturwerte eine geringere Abweichung von der Mittelwertkurve haben darf als im Bereich der niedriger gewichteten Strukturwerte. Als Ausgleichskurven eignen sich beispielsweise Spline-Funktionen. Je nach Randbedingungen wird abgesehen von Rezirkulationspeaks beispielsweise ein globales Maximum für die Applikation eines Tracers oder Kontrastmittels und gegebenenfalls jeweils ein lokales Maximum bei zJB. vorliegender Extravasation, Leakage in Tumoren, spezifischer oder unspezifischer Anreicherung in der mathematischen Funktion zugelassen.
Hierbei ist besonderes Augenmerk auf den Beginn der Ausgleichskurve zu legen. Da direkt nach Applikation eines Tracers oder Kontrastmittels schnelle Änderungen von hohen Signalwerten auftreten können, ist bei der Wahl der Ausgleichskurvenberechnung darauf zu achten, dass die Ausgleichskurve für die Zeitpunkte vor dem mittleren ersten Zeitabschnitt die Strukturwertentwicklung sinnvoll widerspiegelt
Beispielsweise kann in einer einfachen Variante der Kurvenbeginn mit Hilfe der Steigung der ersten beiden Mittelwerte extrapoliert werden.
Für die Einpassung der Ausgleichskurve lassen sich die dem Mathematiker bekannten mathematischen Optimierungsverfahren verwenden (siehe z.B.: J. A. Snyman: Practical Mathematical Optimization; Springer- Verlag 2005 / C. Daniel et al. : Fitting equations to data; 2nd ed., Wiley 1980 / P. Dierckx: Curve and Surface Fitting with Splines, Oxford Science Publications 1996).
Die Ausgleichskurve stellt optimierte Strukturwerte zu beliebigen Zeitpunkten innerhalb des Messintervalls bereit, da die Ausgleichskurve eine kontinuierliche Zeitkurve darstellt und nicht aus diskreten Werten besteht
Das Ergebnis ist also ein Datensatz mit optimierten Strukturwerten zu frei wählbaren Zeitpunkten innerhalb des Messintervalls. In dem gewonnnen optimierten Datensatz stecken aufgrund der berücksichtigten Randbedingungen Informationen, die es erlauben, morphologische und/oder physiologische Strukturen innerhalb des Datensatzes gezielt hervorzuheben oder zu unterdrücken.
Diese Möglichkeit ist in der folgenden Ausfuhrungsform in optimaler Art und Weise gegeben, wobei entsprechende Operationen auch in der vorliegenden Ausfuhrungsform möglich sind.
2. Anpassen an ein mathematisches Modell
In einer zweiten bevorzugten Ausführungsform des erfindungsgemäßen Verfahrens wird ein mathematisches Modell verwendet, um die optimierten Strukturweite in Schritt c) zu berechnen.
Diese Ausführungsform des erfindungsgemäßen Verfahrens umfässt die folgenden Schritte: cl) Bereitstellen eines mathematischen Modells, das das zeitliche Verhalten des Strukturwertes in den Bereichen des Körpers beschreibt; c2) für jeden Teilbereich: Anpassen mindestens eines Parameters des Modells an die gemessenen Strukturwerte und Ermitteln einer Modellfunktion, die den zeitlichen Verlauf der gemessenen Strukturwerte als Ergebnis eines mathematischen Optimierungsverfahrens optimal wiedergibt, wobei die Modellfunktion optimierte
Strukturwerte liefert und wobei durch das Optimierungsverfahren auch optimierte Modellparameter gewonnen werden.
Das mathematische Modell repräsentiert die Randbedingungen, die in Schritt b) des erfindungsgemäßen Verfahrens aufgestellt worden sind.
Als mathematisches Modell wird - je nach untersuchtem Körperbereich und den physikalisch-biologisch-chemischen Eigenschaften eines eventuell applizierten Hilfsmittels wie eines z.B. eines Tracers oder Kontrastmittels - vorzugsweise ein Ein- oder Mehr- Kompartimenten-Modell verwendet.
Solche Modelle sind dem Fachmann der Pharmakokinetik hinreichend bekannt (siehe z.B. Molecuiar Imaging: Computer Reconstruction and Practice, Proceedings of the NATO Advanced Study Institute on Molecuiar Imaging from Physical Principles to Computer Reconstruction and Practice, Springer-Verlag 2006 / Physiologically based pharmacokinetic modelling; ed. by M. B. Reddy et al. \ Wiley-Interscience 2005 / Peter L. Bonate: Pharmacokinetic-Pharmacodynamic Modeling and Simulation; 2nd ed., Springer- Verlag 2011).
Bei solchen Modellen wird der betrachtete Kötperbereich als ein aus einem oder mehreren Kompartimenten aufgebauter Körper betrachtet. Dabei wird für jede zeitliche Änderung des Strukturwertes ein Kompartiment im Modell verwendet. So verteilt sich beispielsweise ein Tracer nach einer Bolus-Applikation im Blut eines Patienten in einer für den Patienten und den Tracer charakteristischen Art und Weise und Geschwindigkeit und wird allmählich eliminiert und ggf. metabolisiert.
Ein weiteres Kompartiment wird beispielsweise für das Modell benötigt, wenn der Tracer das Gefäßsystem aufgrund seiner physiologisch-chemischen Eigenschaften verlassen und extravasieren kann. Für alle Effekte oder physiologischen Funktionen, die zu einer zeitlichen Änderung des Strukturwertes in dem betrachteten Datensatz fuhren, ist ein Kompartiment in der Modellfunktion vorzusehen. Um das zeitliche Verhalten der Strukturwerte mit Hilfe des Modells bestmöglich nachzubilden, können verschiedene mathematische Verfahren angewendet werden. So kann eine Modellfunktion beispielsweise durch Lösung der für das Modell aufstellbaren Differentialgleichungen erhalten werden, wie es bei pharmakokinetischen Modellierungen durchgeführt wird.
Die Modellfunktion kann aber auch durch Simulation der zeitlichen Entwicklung der betrachteten Strukturwerte über die Messzeit erhalten werden. Durch Variation der ModeUfunktionsparameter ist hierbei eine mathematische Anpassung der Modellfunktion an das zeitliche Verhalten der Strukturwerte möglich.
Die Ermittlung einer Modellfunktion durch Anpassung an ein mathematisches Modell wird im erfmdungsgemfißen Verfahren vorzugsweise mit dem Simulationsansatz durchgeführt Das Ergebnis ist eine Modellfunktion, die das zeitliche Verhalten der Strukturwerte im mathematischen Sinne optimal wiedergibt Die Modellfunktion stellt optimierte Strukturwerte zu beliebigen Zeitpunkten innerhalb des Messintervalls bereit, da die Modellfunktion eine kontinuierliche Zeitkurve darstellt und nicht aus diskreten Werten besteht Weiterhin resultiert aus der genannten Verfahrensvariante für jeden Teilbereich des gescannten Körpers ein Datensatz von optimierten Parametern, der den Einftuss eines jeden Kompartiments auf den zeitlichen Verlauf des Strukturwertes angibt
Damit ist es möglich, die Beitrage der einzelnen Kompartimente hervorzuheben, zu vennindern oder ganz wegzulassen.
Dies kann dadurch erfolgen, dass bei der Berechnung des Datensatzes für einen beliebigen Zeitpunkt innerhalb der Messzeit nicht alle durch die Anpassungsrechung ermittelten optimierten Werte der Modellparameter verwendet werden. Durch Eingrenzung des Wertebereiches eines oder mehrerer Parameter kann gezielt der Beitrag eines oder mehrerer Kompartimente beeinflusst werden.
So kann beispielsweise bei einer kontrastmittelgestützten MR-Tomografie am Patienten die Kontrastierung des Gefäßsystems in dem ausgegebenen Datensatz nach Bedarf unterdrückt oder hervorgehoben werden. Das Ergebnis der Modellanpassung ist somit ein Datensatz mit optimierten Strukturwerten und ein Datensatz mit zugehörigen Modeüparametern, mit dem sich der optimierte Datensatz in verschiedenen für das Verständnis der Untersuchungsdaten nützlichen Varianten ausgeben lässt.
Oben wurde der Einfachheit halber angenommen, dass sich der Körperbereich bei der Erzeugung des ersten, auf Messwerten beruhenden Datensatzes in Bezug zur Messvorrichtung nicht bewegt hat. Hat er sich hingegen bewegt, so sind zeitliche Änderungen der Strukturwerte nicht allein auf Änderungen des strukturellen oder funktionellen Zustands des betrachteten Körperbereichs zurückzufuhren sondern auch darauf, dass sich die betrachteten Teilbereiche im Lauf der Zeit gegenüber der Messvorrichtung verschieben. Wenn diese zeitlichen Änderungen des Strukturwertes nicht mit den Randbedingungen vereinbar sind, werden sie durch das beschriebene Verfahren reduziert oder eliminiert. Dies gilt insbesondere für Strukturwcrtänderungen, die durch Bewegungen verursacht werden, welche schneller sind, als die betrachteten zeitlichen Änderungen des Strukturwertes oder welche oszillatorischen Charakter haben, wie beispielsweise die Bewegung des Herzmuskels.
Da ungewollte Bewegungen des Körpers während des Scanvorgangs immer zu einer Verfälschung der Darstellung des gescannten Körpers führen können, ist es grundsatzlich von Vorteil, diese bereits im ersten, auf Messwerten beruhenden Datensatz erkennen zu können und zu reduzieren oder zu eliminieren. Weist der erste Datensatz jedoch einen zu großen räumlichen Rauschanteil auf, so kann eine Bewegungskorrektur auch auf Basis des optimierten Datensatzes durchgeführt werden, d.h. nach der Durchführung des erfindungsgemäßen Verfahrens, sofern die Bewegung nicht schon durch das erfindungs gemäße Verfahren genügend reduziert ist.
In Schritt d) des erfindungsgemäßen Verfahrens erfolgt die Ausgabe eines optimierten Datensatzes. Der optimierte Datensatz repräsentiert einen Bereich im untersuchten Körper. Üblicherweise stimmt der Bereich in Schritt d) mit dem Bereich in Schritt a) überein. Es ist jedoch auch denkbar, dass der Bereich in Schritt d) nur einen Teilbereich des Bereichs aus Schritt a) darstellt Es ist denkbar, dass im Rahmen oder im Anschluss der Berechnung der optimierten Strukturwerte in Schritt c) oder durch eine Bewegungskorrektur Teilbereiche verworfen worden sind. Dies gilt insbesondere für Randbereiche des Datensatzes, die aufgrund von Bewegung eventuell nicht in allen Messzeitintervallen räumlich übereinstimmen.
Der optimierte Datensatz basiert auf den optimierten Strukturwerten aus Schritt c). Daher kann Schritt d) nur im Anschluss an Schritt c) erfolgen.
Der optimierte Datersatz kann in Form einer oder mehrerer zwei- oder dreidimensionaler Darstellungen des Bereichs des Körpers auf einem Bildschirm oder als Ausdruck ausgegeben werden. Ebenso ist es denkbar, dass die Ausgabe auf einem Datenträger in Form von maschinenlesbaren Daten erfolgt.
Der optimierte Datensatz, der mittels des erfindungs gemäßen Verfahrens erzeugt worden ist, ist ebenfalls ein Gegenstand der vorliegenden Erfindung. Ein weiterer Gegenstand der vorliegenden Erfindung ist ein Computerprogrammprodukt mit auf einem maschinenlesbaren Träger gespeichertem Programmcode zur Durchführung des erfindungs gemäßen Verfahrens auf einem Computer.
Das erfindungsgemäße Verfahren eignet sich zur Optimierung aller bekannten 3D- Aufnahmen oder Tomografieaufnahmen, wie beispielsweise zur Optimierung von SPECT-, PET-, CT- oder MRT- Aufnahmen, oder Messdaten aus einem 3D- oder 4D-Ultraschall- Verfahren oder der optischen Tomografie (siehe einschlägige Literatur wie z.B.: Ashok Kliurana, Nirvikar Dahiya: 3D & 4D Vitrasound - A Text and Atlas, Jaypee Brothers Medical Pubiishers (P) Ltd., 2004; R. Weissleder et alz. Molecular Imaging - Principles and Practice, Peopie's Medical Publishing House, USA, 2010; G. B. Saha; Basics of PET Imaging, 2nd edition, Springer 2010; S. A. Jackson, R. M. Thomas; CT, MRT, Ultraschall auf einen Blick, Elsevier 2009; Olaf Dössel: Bildgebende Verfahren in der Medizin, Springer- Verlag Berlin Heidelberg New York, 2000). Mit Hilfe des erfindungsgemäßen Verfahrens lassen sich überraschenderweise aus einer Sequenz von gemessenen Tomografie-Aufnahmen in der Regel deutlich rauschreduzierte Tomografie-Aufnahmen erzeugen, ohne dass die Kinetik der Messdaten verloren geht, wie beispielsweise bei der Erstellung der sogenannten MIP (Maximum Intensity Projection) oder der Mittelung über alle Einzelscans. Bewegungen, die während der Messzeit im gescannten Körper oder in Teilbereichen des gescannten Körpers vorkommen, werden durch das erfindungsgemäße Verfahren in vielen Fallen vermindert, was insbesondere bei stark verrauschten Datensätzen vorteilhaft ist Bildunschärfen, wie sie bei statischen Aufnahmen mit nur einem Datensatz pro Gesamtmesszeit unumgänglich sind, werden mit dem erfindungsgemäßen Verfahren reduziert und die räumliche Auflösung liegt näher an der physikalisch möglichen Auflösung des Scangerätes.
Es lassen sich je nach Bedarf Darstellungen eines Körperbereichs erzeugen, in denen morphologische und/oder physiologische Strukturen gezielt hervorgehoben oder unterdrückt sind. Dies erlaubt beispielsweise die Erstellung besserer Diagnosen.
Die Erfindung wird in der Figurenbeschreibung (Fig. 1 bis 4) und anhand eines Beispiels näher erläutert ohne es hierauf zu beschränken.
Beispiel
Die nachfolgende Erläuterung des erfindungsgemäßen Verfahrens wird für den Fall der abschnittsweisen Glättung geführt Es sei ein zeitlicher Verlauf eines Strukturwertes für einen diskreten räumlichen Teilbereich aus einem tomografi sehen PET-Datensatz gegeben, wie er in Figur la dargestellt ist
Zu Beginn dieses Zeitverlaufs ist ein Signalabfall zu erkennen, wie er nach Applikation und Anflutung des Tracers in vivo zu erwarten ist Danach durchläuft die Kurve anscheinend noch ein Maximum, bevor sie zum Ende der Scanzeit auf einen niedrigen Wert abfällt. Allem Überlagert ist das für PET-Daten nicht untypische Rauschen aufgrund der Statistik der Zerfallsereignisse.
Einen solchen Verlauf würde man für einen Thrombus-Tracer erwarten, der ein Hauptmaximum in der Datenkurve aufgrund der Anflutung und Auswaschung des Tracers nach Applikation und ein weiteres Maximum aufgrund einer möglichen Anreicherung des Tracers in oder an eventuell vorhandene Thromben im Gefäßraum haben könnte. Dementsprechend werden die Randbedingungen für diesen Fall mit einem Haupt- und einem Nebenmaximum in der Strukturwert-Zeit-Kurve gewählt Die Längen der Abschnitte, welche für die abschnittsweise Glättung benötigt werden, sind in Figur lb eingetragen. Sie lassen sich grob aus der Messkurve ablesen. Schnelle Änderung des Strukturwertes am Anfang der Kurve erfordern kurze Abschnitte, wohingegen lange Abschnitte für das sich über einen längeren Zeitraum erstreckende Nebenmaximum zu wählen sind. Bei Messungen, die in der Kombination Tracer oder Kontrastmittel und untersuchte Spezies nicht zum ersten Mal durchgeführt werden, sind die möglichen Änderungen des Strukturwertes und damit auch die Abschnittslängen bekannt und können dementsprechend ausgewählt werden. Für den Fall der Anpassung der Messdaten an ein pharmakologisches Modell gilt Analoges.
Als nächstes werden die in den verschiedenen zeitlichen Abschnitten liegenden Strukturwerte pro Abschnitt gemittelt und gemäß der gewählten Randbedingungen für ein Hauptmaximum und maximal ein Nebenmaximum, wenn nötig, in der Höhe des Wertes korrigiert In der vorliegenden Stnikturwertkurve ist aus diesem Grund der etwas höhere Mittelwert des vorletzten Abschnitts (Minute 44-52) auf den Mittelwert des drittletzten Abschnitts (Minute 36-44) herunterzukorrigieren, da es aufgrund der Randbedingungen außer dem deutlich größeren Nebenmaximum bei kleiner 20 Minuten kein weiteres Maximum in der Kurve geben darf.
Zum Schluss wurde mathematisch eine Ausgleichskurve durch die errechneten Mittelwerte der Abschnitte gelegt (siehe Figur 1 c) und damit ein optimierter Datensatz erstellt
In den Figuren 2 bis 4 ist beispielhaft ein Ausschnitt aus einem Messdatensatz in den anatomisch üblichen Ebenen dargestellt Figur 2 zeigt den Datensatz ohne Bearbeitung mit dem erfindungsgemäßen Verfahren. Im Vergleich hierzu ist in Figur 3 die erfolgte Rauschminderung mit dem erfindungsgemäßen Verfahren anhand leichter erkennbarer Strukturen und wesentlich weniger Einzelspots ersichtlich. In Figur 4 wird die in Figur 3 erkennbare Struktur bestätigt Durch die Mittelung über alle Messzeitintervalle lässt der in Figur 4 dargestellte Datensatz jedoch keine Rückschlüsse mehr auf die Kinetik der Tracerverteilung im Scankörper zu, im Gegensatz zum Datensatz aus Figur 3. Figurenbeschreibung
Figur 1: Darstellung eines beispielhaften zeitlichen Verlaufs der Tracerkonzentration wahrend eines in vivo PET-Scans in einem diskreten Teilbereich eines PET- Datensatzes
a) ohne Rauschmindcrung durch das erfindungsgemäße Verfahren, b) ohne Rauschminderung durch das erfindungsgemäß e Verfahren und mit zusätzlich eingezeichneten geeigneten Abschnitten für die Abschnittsmittelung nach Schritt c2) der abschnittsweisen Glättung (horizontale Balken) und c) nach Anwendung des erfindungsgemäßen Verfahrens.
Die Abschnittsbalken in Figur lb sind jeweils in der Höhe des aus der Abschnittsmittelung erhalten Wertes eingetragen. Der Start des PET-Scans sei direkt nach Applikation des Tracers erfolgt
Figur 2: Darstellung der anatomischen Ansichten
(a) transversal,
(b) coronal und
(c) sagittal
aus einem in vivo 3 D -PET- Scan.
Der Scan wurde an einem C nomolgus Affen nach Applikation eines Thrombus-Tracers aus der PET-Tracer-Forschung mit einem Kleintier-PET- Scanner aufgenommen. Dargestellt ist der Messdatensatz Nummer 28 von 60 nacheinander durchgeführten Scans ohne Rauschminderung durch das erfindungsgemäße Verfahren. Die Messdauer jedes Messdatensatzes beträgt 1 Minute. Die Messung aller Datensätze erfolgte nacheinander ohne Pause. Die Ebenen für die dargestellten Ansichten sind identisch mit denen aus Figur 3a-c und Figur 4a-c. Die in den Figuren erkennbaren Kreuze stellen die Cursorposition in dem erfindungsgemäßen Computerprogrammprodukt dar, mit welchem die Figuren erstellt wurden.
Figur 3: Darstellung der anatomischen Ansichten
(a) transversal,
(b) coronal und
(c) sagittal aus einem in vivo 3D-PET-Scan.
Der Scan wurde an einem Cynomolgus Affen nach Applikation eines Thrombus-Tracers aus der PET-Tracer-Forschung mit einem Kleintier-PET- Scanner aufgenommen. Dargestellt ist der Messdatensatz Nummer 28 von 60 nacheinander durchgeführten Scans nach Anwendung des erfindungs gemäßen Verfahrens. Die Messdauer jedes Messdatensatzes beträgt 1 Minute. Die Messung aller Datensätze erfolgte nacheinander ohne Pause. Die Ebenen für die dargestellten Ansichten sind identisch mit denen aus Figur 2a-c und Figur 4a-c. Die in den Figuren erkennbaren Kreuze stellen die Cursorposition in dem erfindungsgemäßen Computerprograrnmprodukt dar, mit welchem die Figuren erstellt wurden.
Darstellung der anatomischen Ansichten
(a) transversal,
(b) coronal und
(c) sagittal
aus einem in vivo 3D-PET-Scan.
Der Scan wurde an einem Cynomolgus Affen nach Applikation eines Thrombus-Tracers aus der PET-Tracer-Forschung mit einem Kleintier-PET- Scanner aufgenommen. Dargestellt ist die Mittelung aller 60 Einzeldatensätze, die während der Gesamtmesszeit gescannt wurden. Die Messdauer jedes Messdatensatzes beträgt 1 Minute. Die Messung aller Datensätze erfolgte nacheinander ohne Pause. Die Einzeldatensätze sind nicht mit dem erfindungsgemäßen Verfahren bearbeitet worden. Die Ebenen für die dargestellten Ansichten sind identisch mit denen aus Figur 2a-c und Figur 3a-c. Die in den Figuren erkennbaren Kreuze stellen die Cursorposition in dem erfindungsgemäßen Computeiprograinmprodukt dar, mit welchem die Figuren erstellt wurden.

Claims

Patentansprüche
1. Verfahren zur Erzeugung optimierter Tomografie-Aufhahmen, mindestens umfassend die Schritte:
a) Bereitstellen eines Datensatzes, der einen Bereich im Korper eines Patienten wahrend einer Messzeit repräsentiert, wobei die Repräsentation des Körperbereichs im Datensatz in eine Mehrzahl von diskreten Teilbereichen unterteilt ist, wobei die Messzeit im Datensatz in eine Mehrzahl von diskreten Messintervallen unterteilt ist, wobei jedem Teilbereich zu jedem Messintervall ein diskreter Strukturwert zugeordnet ist; b) Aufstellen von Randbedingungen über den zu erwartenden zeitlichen Verlauf einer Strukturgröße in dem Bereich des Körpers wahrend der Messzeit; c) Berechnen von optimierten Strukturwerten für jeden einzelnen Teilbereich auf Basis von Strukturwerten des einzelnen Teilbereichs zu zeitlich aufeinander folgenden Messintervallen unter Berücksichtigung der Randbedingungen; d) Ausgabe eines optimierten Datensatzes, der einen Bereich im Körper zu beliebig wählbaren Zeitpunkten wahrend der Messzeit repräsentiert und der auf den optimierten Strukturwerten basiert
2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass für jeden Teilbereich in Schritt c) die folgenden Operationen durchgeführt werden: cl) Unterteilen der Messzeit in eine Mehrzahl von Abschnitten, wobei die einzelnen Abschnitte umso kürzer sind, je größer die Änderung der Strukturwerte in einem Bereich der Messzeit ist; c2) Mitteln der Strukturwerte für jeden Teilbereich innerhalb eines jeden Abschnitts; c3) Einpassen einer Ausgleichskurve in die gemittelten Strukturwerte, wobei die Ausgleichskurve optimierte Strukturwerte liefert Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass die Grüße eines jeden Abschnitts in Schritt cl) umgekehrt proportional zum Betrag der ersten Ableitung der
Strukturwerte nach der Zeit ist Verfahren nach Anspruch 2 oder 3, dadurch gekennzeichnet, dass die Abschnitte in Schritt cl) so gestaltet sind, dass sich jeweils zwei zeitlich aufeinander folgende Abschnitte in ihren Randbereichen Überlappen. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass in Schritt c) die folgenden Operationen durchgeftlhrt werden: cl) Bereitstellen eines mathematischen Modells, das das zeitliche Verhalten des
Strukturwertes in den Bereichen des Körpers beschreibt; c2) für jeden Teilbereich: Anpassen mindestens eines Parameters des Modells an die gemessenen Strukturwerte und Ermitteln einer ModeUfunktion, die den zeitlichen Verlauf der gemessenen Strukturwerte als Ergebnis eines mathematischen Optimierungsverfahrens optimal wiedergibt, wobei die Modellfunktion optimierte Strukturwerte liefert und wobei durch das Optimierungsverfahren auch optimierte Modellparameter gewonnen werden; Verfahren nach Anspruch S, dadurch gekennzeichnet, dass das mathematische Modell ein pharmakokinetisches Ein- oder Mehr-Kompartimenten-Modell ist Verfahren nach einem der Ansprache 1 bis 6, dadurch gekennzeichnet, dass der erste Datensatz aus Messungen resultiert, die an einem lebenden Organismus vorgenommen wurden. Verfahren nach einem der Ansprüche 1 bis 6, dadurch gekennzeichnet, dass der erste Datensatz aus Messungen resultiert, die an einem nicht lebenden Objekt vorgenommen wurden. Verfahren nach einem der Ansprüche 1 bis 8, dadurch gekennzeichnet, dass es sich bei dem ersten Datensatz um SPECT-, PET-, CT-, MRT-Aumahmen, oder einen Messdatensatz aus einem 3D- oder 4D-Ultraschall- Verfahren oder der optischen Tomografie handelt
Verfahren nach einem der Ansprüche 1 bis 9, dadurch gekennzeichnet, dass in dem optimierten Datensatz Strukturwerte auf Basis der Randbedingungen gezielt verändert werden, um morphologische und oder physiologische Strukturen hervorzuheben oder zu unterdrücken.
Optimierter Datensatz, erzeugt durch ein Verfahren nach einem der Ansprüche 1 bis 10.
Computerprogrammprodukt mit Programmcodemitteln zur Ausführung Verfahrens nach einem der Ansprüche 1 bis 10 auf einem Computersystem.
EP12784529.5A 2011-10-25 2012-10-24 Verfahren zur erzeugung optimierter tomografie-aufnahmen Withdrawn EP2770911A1 (de)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102011085180A DE102011085180A1 (de) 2011-10-25 2011-10-25 Verfahren zur Erzeugung optimierter Tomografie-Aufnahmen
PCT/EP2012/071035 WO2013060716A1 (de) 2011-10-25 2012-10-24 Verfahren zur erzeugung optimierter tomografie-aufnahmen

Publications (1)

Publication Number Publication Date
EP2770911A1 true EP2770911A1 (de) 2014-09-03

Family

ID=47177956

Family Applications (1)

Application Number Title Priority Date Filing Date
EP12784529.5A Withdrawn EP2770911A1 (de) 2011-10-25 2012-10-24 Verfahren zur erzeugung optimierter tomografie-aufnahmen

Country Status (16)

Country Link
US (1) US20140294275A1 (de)
EP (1) EP2770911A1 (de)
JP (1) JP2014535048A (de)
KR (1) KR20140131500A (de)
CN (1) CN104144649A (de)
AU (1) AU2012330480A1 (de)
BR (1) BR112014009244A8 (de)
CA (1) CA2853188A1 (de)
DE (1) DE102011085180A1 (de)
IL (1) IL231618A0 (de)
IN (1) IN2014MN00918A (de)
MX (1) MX2014004276A (de)
RU (1) RU2014120903A (de)
SG (1) SG11201401780UA (de)
WO (1) WO2013060716A1 (de)
ZA (1) ZA201402962B (de)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105717087B (zh) * 2016-03-10 2019-05-14 天津大学 螺旋离散扫描式荧光剂药代动力学参数直接成像方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020010550A1 (en) * 1998-09-14 2002-01-24 George M. Grass Pharmacokinetic-based drug design tool and method
DE102007046579B3 (de) * 2007-09-27 2009-01-29 Siemens Ag Verfahren zur Detektion von Bewegungen und Korrektur von Bewegungen in tomographischen und projektiven Aufnahmeserien und Tomographie- beziehungsweise Projektionssystem zur Durchführung dieses Verfahrens
WO2011070465A2 (en) * 2009-12-10 2011-06-16 Koninklijke Philips Electronics, N.V. Method and apparatus for using time of flight information to detect and correct for motion in imaging scans
US8761478B2 (en) * 2009-12-15 2014-06-24 General Electric Company System and method for tomographic data acquisition and image reconstruction
US20110148928A1 (en) * 2009-12-17 2011-06-23 General Electric Company System and method to correct motion in gated-pet images using non-rigid registration
FR2957441B1 (fr) * 2010-03-10 2016-01-01 Commissariat Energie Atomique Procede d'extraction simultanee de la fonction d'entree et des parametres pharmacocinetiques d'un principe actif.
JP5339562B2 (ja) * 2010-03-30 2013-11-13 独立行政法人放射線医学総合研究所 核医学イメージング装置の画像化方法、システム、核医学イメージグシステム及び放射線治療制御システム
CN102151142B (zh) * 2011-04-14 2012-08-15 华中科技大学 一种正电子发射断层成像中的运动门控方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2013060716A1 *

Also Published As

Publication number Publication date
WO2013060716A1 (de) 2013-05-02
ZA201402962B (en) 2015-04-29
BR112014009244A2 (pt) 2017-06-13
CN104144649A (zh) 2014-11-12
US20140294275A1 (en) 2014-10-02
CA2853188A1 (en) 2013-05-02
IL231618A0 (en) 2014-05-28
RU2014120903A (ru) 2015-12-10
JP2014535048A (ja) 2014-12-25
DE102011085180A1 (de) 2013-04-25
AU2012330480A1 (en) 2014-04-17
BR112014009244A8 (pt) 2017-06-20
MX2014004276A (es) 2014-07-09
KR20140131500A (ko) 2014-11-13
IN2014MN00918A (de) 2015-05-01
SG11201401780UA (en) 2014-09-26

Similar Documents

Publication Publication Date Title
DE102007061935A1 (de) Verfahren zur Qualitätssteigerung von computertomographischen Aufnahmeserien durch Bildverarbeitung und CT-System mit Recheneinheit
EP3111417B1 (de) Rauschreduktion in tomogrammen
DE102013103832A1 (de) Dämpfungskorrektur in der Positronen-Emissions-Tomographie unter Verwendung von Magnetresonanztomographie
DE102005023907A1 (de) Verfahren zur Ermittlung von Positronen-Emissions-Messinformationen im Rahmen der Positronen-Emissions-Tomographie
DE102012205935B4 (de) Verfahren zur Aufnahme eines vierdimensionalen Angiographie-Datensatzes
DE102012204980A1 (de) Verfahren zur Rekonstruktion von CT-Bildern mit Streustrahlenkorrektur, insbesondere für Dual-Source CT-Geräte
DE102011083647A1 (de) Verfahren, Rechensystem und CT-System zur Erzeugung eines bewegungskompensierten CT-Bilddatensatzes eines sich teilweise und zyklisch bewegenden Untersuchungsobjektes
DE102011086771A1 (de) Computertomographieanlage und Verfahren zum Ermitteln von Volumeninformationen zu einem Körper
DE102017217599A1 (de) Medizinische Informationsverarbeitungsvorrichtung, Röntgen-CT-Vorrichtung und medizinisches Informationsverarbeitungsverfahren
DE102015205959A1 (de) Ermitteln der Geschwindigkeit eines Fluids mit Hilfe eines bildgebenden Verfahrens
DE102012218289A1 (de) Bewegungskorrigierte PET-Bilder
DE102012212774A1 (de) Verfahren zur Korrektur von Metallartefakten und Röntgeneinrichtung
DE102009030110A1 (de) Verfahren zur Ermittlung der Ventilation einer Lunge
DE102008047825A1 (de) Verfahren, Vorrichtung und Computerprogrammprodukt zum Darstellen verschiedener Bilder einer Kavität
CH702961A2 (de) Verfahren zum Abgleich von mit unterschiedlichen Systemen, z.B. Tomografen, aufgenommenen Bilddaten.
EP2242023A1 (de) Verfahren zur bewegungskompensierten Rekonstruktion eines dreidimensionalen endgültigen Bilddatensatzes
DE102007045313B4 (de) Verfahren zur getrennten dreidimensionalen Darstellung von Arterien und Venen in einem Untersuchungsobjekt
DE102006015749B4 (de) Verfahren zur Ermittlung von Positronen-Emissions-Messinformationen eines Körperbereichs eines Untersuchungsobjekts sowie zugehörige Vorrichtung
DE102009007236A1 (de) CT-Bildrekonstruktion eines sich bewegenden Untersuchungsobjektes
DE102010040041B3 (de) Verfahren zur Korrektur von durch zeitliche Veränderungen von Schwächungswerten auftretenden Artefakten
DE102016215976A1 (de) Ermittelung einer klinischen Kenngröße mit einer Kombination unterschiedlicher Aufnahmemodalitäten
DE102004006548B4 (de) Verfahren zur Planung der Strahlentherapie eines Patienten und CT-System hierzu und zur Erstellung von CT-Aufnahmen
DE102015207894A1 (de) Ermitteln der Geschwindigkeit eines Fluids mit Hilfe eines bildgebenden Verfahrens
WO2013060716A1 (de) Verfahren zur erzeugung optimierter tomografie-aufnahmen
DE102008020814A1 (de) Verfahren zur Ermittlung und Darstellung von Perfusionsparametern in der Tomographie

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20140502

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

DAX Request for extension of the european patent (deleted)
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20141216