DE102009048070B4 - Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens - Google Patents

Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens Download PDF

Info

Publication number
DE102009048070B4
DE102009048070B4 DE102009048070.6A DE102009048070A DE102009048070B4 DE 102009048070 B4 DE102009048070 B4 DE 102009048070B4 DE 102009048070 A DE102009048070 A DE 102009048070A DE 102009048070 B4 DE102009048070 B4 DE 102009048070B4
Authority
DE
Germany
Prior art keywords
image
target object
control parameters
images
movement
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.)
Active
Application number
DE102009048070.6A
Other languages
English (en)
Other versions
DE102009048070A1 (de
Inventor
Dr. Bert Christoph
Dr. Rietzel Eike
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.)
GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
Original Assignee
GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
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 GSI Helmholtzzentrum fuer Schwerionenforschung GmbH filed Critical GSI Helmholtzzentrum fuer Schwerionenforschung GmbH
Priority to DE102009048070.6A priority Critical patent/DE102009048070B4/de
Priority to EP10750085A priority patent/EP2482728A1/de
Priority to PCT/EP2010/062292 priority patent/WO2011038999A1/de
Publication of DE102009048070A1 publication Critical patent/DE102009048070A1/de
Application granted granted Critical
Publication of DE102009048070B4 publication Critical patent/DE102009048070B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1049Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • G06T7/0014Biomedical image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/32Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4064Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
    • A61B6/4092Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam for producing synchrotron radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1048Monitoring, verifying, controlling systems and methods
    • A61N5/1049Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam
    • A61N2005/1061Monitoring, verifying, controlling systems and methods for verifying the position of the patient with respect to the radiation beam using an x-ray imaging system having a separate imaging source
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1085X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
    • A61N2005/1087Ions; Protons
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/103Treatment planning systems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing

Abstract

Verfahren zur Bildbearbeitung, wobei- ein erstes von einem medizinischen Bildaufnahmegerät (CT1) erfasstes Bild (PL-CT) eines Objektes, welches ein Zielobjekt enthält, vorliegt, umfassend Bildpunkte und ihnen zugeordnete Bildpunktwerte, und ein zu einem anderen Zeitpunkt erfasstes zweites Bild (K-CT) des Objektes vorliegt,- das erste Bild (PL-CT) und das zweite Bild (K-CT) hinsichtlich des Zielobjektes miteinander verglichen werden,- unter Berücksichtigung des Vergleichs eine Bewegung des Zielobjektes innerhalb des ersten Bildes (PL-CT) durchgeführt wird, indem eine Neuzuordnung von Bildpunktwerten zu Bildpunkten vorgenommen wird,- als erstes Bild (PL-CT) eine zeitliche Abfolge von Bildern vorliegt, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen,- das zweite Bild (K-CT) nur einen Zeitpunkt abbildet, und die Neuzuordnung für alle Bilder (PL-CT) der zeitlichen Abfolge vorgenommen wird, wobei von einer bei Aufnahme des zweiten Bildes (K-CT) gegenüber der Aufnahme der Abfolge der ersten Bilder (PL-CT) unveränderten Form der Bewegungstrajektorie ausgegangen wird.

Description

  • Die Erfindung betrifft ein Verfahren zur Bildbearbeitung, welches insbesondere in einem Verfahren zur Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage einsetzbar ist. Ferner betrifft die Erfindung eine Vorrichtung zur Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage. Eine derartige Vorrichtung bzw. ein derartiges Verfahren kann insbesondere im Rahmen der Partikeltherapie Einsatz finden, beispielsweise im Rahmen der Therapieplanung, bei der im Vorfeld einer Bestrahlung Steuerparameter ermittelt werden, die es erlauben, anschließend während der Bestrahlung ein Zielobjekt gemäß bestimmten Vorgaben zu bestrahlen.
  • Die Partikeltherapie ist ein etabliertes Verfahren zur Behandlung von Gewebe, insbesondere von Tumorerkrankungen. Bestrahlungsverfahren, wie sie in der Partikeltherapie für Patienten eingesetzt werden, finden auch in nichttherapeutischen Gebieten Anwendung. Hierzu gehören beispielsweise Forschungsarbeiten, etwa zur Produktentwicklung, im Rahmen der Partikeltherapie an nicht-lebenden Phantomen oder Körpern, Bestrahlungen von Materialien, etc. Bei der Partikeltherapie werden geladene Partikel wie z.B. Protonen oder Kohlenstoff- oder andere Ionen auf hohe Energien beschleunigt, zu einem Partikelstrahl geformt und über ein Hochenergiestrahltransportsystem zu einem oder mehreren Bestrahlungsräumen geführt. In einem Bestrahlungsraum wird das zu bestrahlende Objekt in einem Zielvolumen mit dem Partikelstrahl bestrahlt. Unter einem Zielvolumen wird ein Raumbereich des zu bestrahlenden Objektes verstanden, innerhalb dessen eine bestimmte Dosis deponiert werden soll.
  • Die Druckschrift US 2009/0116616 A1 zeigt ein System und ein Verfahren aus dem Feld der Röntgentherapie zum optimierten Bereitstellen einer Strahlentherapiebehandlung. Die Bereitstellung der Behandlung wird in Echtzeit optimiert um verschieden Faktoren wie anatomische und physiologische Änderungen des Patienten und Änderungen der Maschinenkonfiguration zu berücksichtigen.
  • Aus der Druckschrift WO 00/ 67 202 A1 geht ein Verfahren zur Kombination eines ersten digitalen radiographischen Bildes mit einem zweiten digitalen radiographischen Bild hervor Das finale Bild entspricht dem ersten Bild, in welchem Pixel entsprechend eines Grauanteils innerhalb eines „CT-Intervalls“ einer digitalen Bearbeitung unterzogen werden.
  • Die Druckschrift US 2009/0087124 A1 zeigt schließlich ein Verfahren aus dem Feld der Röntgenstrahltherapie, bei welchem ein erstes und ein zweites Bild erhalten wird und ein „Deformationsfeld“ anhand der beiden Bilder bestimmt wird. Aus dem Deformationsfeld wird weiter eine Transformationsmatrix bestimmt.
  • Abhängig von der Energie des Partikelstrahls dringt dieser in das Objekt ein und wechselwirkt mit diesem. In einem relativ eng umschriebenen Bereich (Bragg Peak Region) weist das Tiefendosisprofil ein Maximum auf, da es hier besonders viele Wechselwirkungsprozesse gibt. Partikelstrahlen sind aus diesem Grund für die Bestrahlung eines Zielvolumens besonders vorteilhaft, da sie am Ende ihrer Reichweite ein Maximum der Energiedeposition aufweisen, den so genannten Bragg-Peak. So können auch eingebettete Strukturen effektiv bestrahlt werden, ohne die einbettende Umgebung allzu stark zu schädigen. Dies hat also den Vorteil, dass das das Zielvolumen umgebende Gebiet verglichen mit anderen Bestrahlungsarten effektiv geschont werden kann.
  • Insbesondere mit einem gescannten Partikelstrahl, d.h. einem fein kollimierten Partikelstrahl mit einer Halbwertsbreite von 3-10 mm, der durch laterale magnetische Ablenkung und Veränderung der Strahlenergie sukzessive auf unterschiedliche Bereiche des Zielvolumens gerichtet wird und somit das Zielobjekt scannt, kann eine genaue Dosisapplikation erreicht werden.
  • Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren und eine entsprechende Vorrichtung aufzuzeigen, welche im Rahmen der Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage einsetzbar sind. Ferner sollen ein Verfahren zur Bestrahlung eines Zielobjektes mit Partikeln vorgestellt werden, sowie ein Computerprogramm und ein Computerprogrammprodukt.
  • Diese Aufgabe wird durch ein Verfahren mit den Merkmalen des Anspruchs 1, sowie durch eine Vorrichtung, durch ein Computerprogramm und ein Computerprogrammprodukt mit Merkmalen von nebengeordneten Ansprüchen gelöst. Vorteilhafte Ausgestaltungen und Weiterbildungen sind Gegenstand von Unteransprüchen.
  • Das erfindungsgemäße Verfahren dient der Bildbearbeitung. Hierbei liegt ein erstes von einem medizinischen Bildaufnah megerät erfasstes Bild eines Objektes, welches ein Zielobjekt enthält, vor. Dieses erste Bild umfasst Bildpunkte und ihnen zugeordnete Bildpunktwerte. Ferner liegt ein zu einem anderen Zeitpunkt erfasstes zweites Bild des Objektes vor. Das erste Bild und das zweite Bild werden hinsichtlich des Zielobjektes miteinander verglichen. Unter Berücksichtigung des Vergleichs wird eine Bewegung des Zielobjektes innerhalb des ersten Bildes durchgeführt, indem eine Neuzuordnung von Bildpunktwerten zu Bildpunkten vorgenommen wird. Das erste Bild liegt dabei als eine zeitliche Abfolge von Bildern vor, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen. Das zweite Bild bildet dabei nur einen Zeitpunkt ab. Die Neuzuordnung wird für alle Bilder der zeitlichen Abfolge vorgenommen, wobei von einer bei Aufnahme des zweiten Bildes gegenüber der Aufnahme der Abfolge der ersten Bilder unveränderten Form der Bewegungstrajektorie ausgegangen wird.
  • Die Bearbeitung betrifft also ein Bild, welches von einem medizinischen Bildaufnahmegerät erfasst wurde. Als Bildaufnahmegerät kommt insbesondere in Computertomographie-Gerät in Betracht. In diesem Fall handelt es sich bei dem ersten Bild um ein CT-Bild. Alternativ ist auch der Einsatz andersartiger Bildaufnahmegeräte möglich, insbesondere der Einsatz anderer Tomographie-Geräte wie z.B. eines NMR-Gerätes. Das zweite Bild kann mit einem gegenüber dem ersten Bild gleich- oder andersartigen medizinischen Bildaufnahmegerät erfasst werden.
  • Bei dem ersten Bild kann es sich um ein zweidimensionales Bild oder ein dreidimensionales Volumenbild handeln. Ein solches Bild ist gegliedert in einzelne Bildpunkte, bei einem zweidimensionalen Bild in so genannte Pixel bzw. bei einem dreidimensionalen Bild in so genannte Voxel. Zu jedem Bildpunkt gehört ein Bildpunktwert; diese werden in CT-Bildern üblicherweise in HU (Hounsfield Units) angegeben. Im Falle eines zweidimensionalen Bildes kann ein Schnittbild vorliegen, oder eine Projektion, wie sie z.B. durch planares Röntgen oder Röntgenfluorosokopie entsteht.
  • Das erste Bild zeigt ein Objekt oder einen Ausschnitt aus diesem Objekt. Insbesondere zeigt das erste Bild ein Zielobjekt, welches Bestandteil des Objektes ist. Bei diesem Zielobjekt kann es sich z.B. um einen Tumor eines Patienten handeln. Es ist jedoch auch möglich, dass das Objekt nicht lebendig ist.
  • Der Bearbeitung des ersten Bildes geht ein Vergleich des ersten Bildes mit dem zweiten Bild voraus. Dieser Vergleich betrifft insbesondere das Zielobjekt, welches Bestandteil des von den Bildern abgebildeten Objektes ist. Da die beiden Datensätze, d.h. das erste Bild und das zweite Bild, zu verschiedenen Zeitpunkten erfasst wurden, ist es möglich, dass das Zielobjekt zwischen diesen beiden Aufnahmen eine Veränderung erfahren hat. Informationen über diese Veränderung konnen durch den Vergleich ermittelt werden.
  • Nach dem Vergleich wird das Zielobjekt innerhalb des ersten Bildes bewegt. Diese Bewegung kann aus einer Verschiebung und/oder Drehung und/oder Formveränderung bestehen. Nach der erfolgten Bewegung existiert somit eine überarbeitete Version des ersten Bildes. Dieses überarbeitete erste Bild sieht so aus, als hätte das Zielobjekt nicht mehr die Lage und/oder Orientierung und/oder Form gemäß dem ursprünglichen ersten Bild, sondern eine hiervon abweichende Lage und/oder Orientierung und/oder Form.
  • In Weiterbildung der Erfindung wird durch die Bewegung des Zielobjektes innerhalb des ersten Bildes dieses an das zweite Bild angenähert. Dies ist insbesondere vorteilhaft, wenn das zweite Bild einen aktuelleren Zustand des Zielobjektes abbildet. Auf diese Weise kann das erste Bild also aktualisiert werden. Die Annäherung des Zielobjektes an das zweite Bild kann vollständig erfolgen oder auch nur teilweise.
  • Besonders vorteilhaft ist es, wenn das zweite Bild zu einem späteren Zeitpunkt erfasst wurde als das erste Bild und das Zielobjekt zwischenzeitlich eine Bewegung innerhalb des Objektes erfahren hat. Diese tatsächliche Bewegung des Zielobjektes kann durch die erfindungsgemaße Bewegung des Zielobjektes innerhalb des ersten Bildes teilweise oder vollständi nachgeführt werden.
  • In Ausgestaltung der Erfindung werden bei der Neuzuordnung dem Zielobjekt entsprechende Bildpunktwerte zu gemäß dem ursprünglichen ersten Bild nicht das Zielobjekt darstellenden Bildpunkten zugeordnet. Es existieren Bildpunktwerte, welche dem Zielobjekt entsprechen, und solche, welche nicht dem Zielobjekt entsprechen. Beispielsweise kann das Zielobjekt besonders große Bildpunktwerte aufweisen, d.h. Werte, welche deutlich größer als diejenigen der Umgebung des Zielobjektes sind. In der vorliegenden Ausgestaltung findet ein Belegen von Bildpunkten, welche ursprünglich außerhalb des Zielobjektes lagen, durch Bildpunktwerte des Zielobjektes statt. Es erfolgt also ein Überschreiben von nicht-Zielobjekt-Bildpunktwerten durch Zielobjekt-Bildpunktwerte. Zusätzlich oder alternativ ist auch der umgekehrte Fall möglich, d.h. bei der Neuzuordnung werden nicht dem Zielobjekt entsprechende Bildpunktwerte zu gemäß dem ursprünglichen ersten Bild das Zielobjekt darstellenden Bildpunkten zugeordnet.
  • Vorteilhaft ist es, wenn bei der Neuzuordnung ein oberer und/oder ein unterer Bildpunktwertgrenzwert für nicht das Zielobjekt darstellende Bildpunkte eingesetzt werden. Es werden also Bildpunkte, welche bei dem überarbeiteten ersten Bild außerhalb des Zielobjektes liegen, mit bestimmten Bildpunktwerten belegt, welche einem oberen und/oder unteren Grenzwert darstellen.
  • Einer Ausgestaltung der Erfindung gemäß wird bei dem Vergleich des ersten Bildes mit dem zweiten Bild ein Algorithmus zur Optimierung eines Ähnlichkeitsmaßes eingesetzt. Hierbei können virtuell verschiedene Bewegungen des Zielobjektes innerhalb des ersten Bildes durchgeführt werden und jeweils das Ähnlichkeitsmaß zwischen diesem durch die virtuelle Bewegung erhaltenen ersten Bild und dem zweiten Bild bestimmt werden.
  • Es ist vorteilhaft, wenn der Vergleich des ersten Bildes mit dem zweiten Bild auf das Zielobjekt beschrankt ist. In diesem Fall können andere Bestandteile des Objektes außer Betracht bleiben. Dies ermöglicht es einerseits, den Vergleich zwischen den beiden Bildern rasch durchzuführen. Andererseits muss das zweite Bild nur hinsichtlich des Zielobjektes bestimmte Qualitätsanforderungen erfüllen, und nicht hinsichtlich der restlichen von dem zweiten Bild abgebildeten Bestandteile des Objektes.
  • Einer Weiterbildung der Erfindung gemäß handelt es sich bei dem ersten Bild um ein das Zielobjekt enthaltendes Volumenbild und bei dem zweiten Bild um ein zweidimensionales Bild des Objektes, z.B. um ein einen Schnitt durch das Zielobjekt enthaltendes Bild des Objektes. Im Falle eines Schnitts kann der direkte Vergleich zwischen diesen beiden Bildern also nur auf diesen Schnitt durch das Zielobjekt beschränkt sein. Hieraus können jedoch Schlussfolgerungen für die Lage und/oder Orientierung und/oder Form des Zielobjektes außerhalb des Schnitts gezogen werden.
  • Insbesondere kann es sich sowohl bei dem ersten als auch bei dem zweiten Bild jeweils um ein CT-Bild handeln, wobei das zweite Bild gegenüber dem ersten CT-Bild von schlechterer Qualität ist. Dies ermöglicht es, das Objekt bei Aufnahme des zweiten Bildes einer niedrigeren Röntgenstrahlungsdosis auszusetzen oder z.B. die Cone-Beam-Computertomographie (kV oder MV) einzusetzen, welche derzeit noch im Vergleich zum konventionellen CT von minderer Qualität ist.
  • Bei dem zweiten Bild kann es sich um ein Fluoroskopie-Bild oder ein Röntgenbild oder ein Ultraschallbild oder ein NMR-Bild handeln. Es ist also möglich, dass das erste Bild hinsichtlich der Aufnahmemodalität grundsätzlich unterschiedlich gegenüber dem zweiten Bild ist. Dies kann voraussetzen, dass vor dem Vergleich das erste und/oder das zweite Bild bearbeitet werden oder für den Vergleich geeignete Größen hieraus berechnet werden müssen, z.B. ein DRR (Digitally Reconstructed Radiograph).
  • In Weiterbildung einer für das Verständnis der Erfindung hilfreichen Ausführung liegt als erstes Bild eine zeitliche Abfolge von Bildern vor, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen, und als zweites Bild liegt eine zeitliche Abfolge von Bildern vor, welche auch eine Bewegungstrajektorie des Zielobjektes aufzeigen, und die Neuzuordnung wird für alle Bilder der zeitlichen Abfolge vorgenommen. Es wird also jedes der ersten Bilder der zeitlichen Abfolge verglichen mit einem entsprechenden zweiten Bild, woraufhin eine Bewegung des Zielobjektes innerhalb des jeweiligen ersten Bildes erfolgt. Auf diese Weise kann die Bewegungstrajektorie der ersten Bilder derjenigen der zeitlichen Abfolge von zweiten Bildern angeglichen werden. Das bislang für den Fall eines einzigen ersten Bildes und eines einzigen zweiten Bildes beschriebene Verfahren wird also auf jedes vorliegende Paar von erstem Bild und zweitem Bild angewandt, wobei jedes Paar jeweils einen bestimmten Zeitpunkt innerhalb der Bewegungstrajektorie abbildet.
  • Im Rahmen der Erfindung liegt hierbei als erstes Bild eine zeitliche Abfolge von Bildern vor, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen, und das zweite Bild bildet nur einen Zeitpunkt ab. Die Neuzuordnung wird für alle Bilder der zeitlichen Abfolge vorgenommen, wobei von einer bei Aufnahme des zweiten Bildes gegenüber der Aufnahme der ersten Bilder unveränderten Form der Bewegungstrajektorie ausgegangen wird. Es findet ein Vergleich zwischen dem ersten Bild und dem zweiten Bild also nur für einen einzigen Zeitpunkt innerhalb der zeitlichen Abfolge statt. Die Resultate dieses Vergleichs werden jedoch auch für die Überarbeitung der anderen ersten Bilder der zeitlichen Abfolge herangezogen. Hierbei wird ausgenutzt, dass sich das Zielobjekt zwar zwischen der Erfassung der ersten Bilder und des zweiten Bildes bewegt haben mag, dass die Bewegungstrajektorie hinsichtlich ihrer Form jedoch unverändert ist.
  • Bei der erfindungsgemäßen Verwendung des beschriebenen Verfahrens zur Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage erfolgt die Bestimmung der Steuerparameter auf Grundlage des durch die Neuzuordnung von Bildpunktwerten veränderten ersten Bildes.
  • Besonders vorteilhaft ist es hierbei, wenn zuerst eine Bestimmung der Steuerparameter auf Grundlage des ersten Bildes erfolgt, und im Anschluss eine Neubestimmung oder eine Verifikation oder eine Veränderung der bestimmten Steuerparameter auf Grundlage des durch die Neuzuordnung von Bildpunktwerten veränderten ersten Bildes.
  • Die erfindungsgemäße Vorrichtung zur Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage umfasst einen Eingang zum Empfangen eines ersten Bildes eines Objektes, welches das Zielobjekt enthält, umfassend Bildpunkte und ihnen zugeordnete Bildpunktwerte, sowie einen Eingang zum Empfangen eines zu einem anderen Zeitpunkt erfassten zweiten Bildes des Objektes, sowie einen Bildbearbeitungsbestandteil zum Vergleichen des ersten Bildes mit dem zweiten Bild hinsichtlich des Zielobjektes und zum Durchführen einer Bewegung des Zielobjektes innerhalb des ersten Bildes, indem unter Berücksichtigung des Vergleichs eine Neuzuordnung von Bildpunktwerten zu Bildpunkten vorgenommen wird, sowie schließlich einen Ermittlungsbestandteil zum Bestimmen der Steuerparameter für die Partikelbestrahlungsanlage auf Grundlage des durch die Neuzuordnung von Bildpunktwerten veränderten ersten Bildes. Das erste Bild liegt dabei als eine zeitliche Abfolge von Bildern vor, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen. Das zweite Bild bildet dabei nur einen Zeitpunkt ab. Die Neuzuordnung wird für alle Bilder der zeitlichen Abfolge vorgenommen, wobei von einer bei Aufnahme des zweiten Bildes gegenüber der Aufnahme der Abfolge der ersten Bilder unveränderten Form der Bewegungstrajektorie ausgegangen wird.
  • Die erfindungsgemäße Vorrichtung kann wie obenstehend hinsichtlich des erfindungsgemäßen Verfahrens erläutert weitergebildet und ausgestaltet werden.
  • Das erfindungsgemäße Computerprogramm verfügt über Programmcode-Mittel, die geeignet sind, ein Verfahren der oben beschriebenen Art durchzuführen, wenn das Computerprogramm auf einem Computer ausgeführt wird.
  • Das erfindungsgemäße Computerprogrammprodukt umfasst auf einem computerlesbaren Datenträger gespeicherte Programmcode-Mittel, die geeignet sind, ein Verfahren der oben beschriebenen Art durchzuführen, wenn das Computerprogramm auf einem Computer ausgeführt wird.
  • Im Folgenden wird die Erfindung anhand eines Ausführungsbeispiels näher erläutert. Dabei zeigen:
    • 1: einen schematischen Überblick über den Aufbau einer Partikeltherapieanlage,
    • 2: in Blockdiagrammform eine Bestrahlungsanlage.
  • 1 zeigt in schematischer Darstellung einen Aufbau einer Partikeltherapieanlage. Diese wird zur Bestrahlung eines auf einer Positioniervorrichtung 12 angeordneten Patienten 14 mit einem Strahl aus Partikeln 16 eingesetzt, der im Folgenden als Partikelstrahl 16 bezeichnet ist. Der Partikelstrahl 16 ist ein Strahl aus Partikeln mit einem definierten Querschnitt und einem definierten, in der Regel schmalen, Spektrum der Partikelenergie. Die Partikelenergie ist die Energie eines einzelnen Partikels vor bzw. beim Eintritt in das zu bestrahlende Objekt.
  • Insbesondere kann ein tumorerkranktes Gewebe des Patienten 14 mit dem Partikelstrahl 16 bestrahlt werden. Es ist ebenfalls möglich, die Partikelbestrahlungsanlage zur Bestrahlung eines nicht-lebenden Objektes 18, insbesondere eines Wasserphantoms, einzusetzen. Die Bestrahlung des Wasserphantoms 18 erfolgt beispielsweise zu Zwecken der Überprüfung und Verifizierung von Bestrahlungsparametern vor und/oder nach einer erfolgten Bestrahlung eines Patienten 14. Es kann ferner vorgesehen werden, andere zu bestrahlende Objekte, insbesondere Versuchsaufbauten wie beispielsweise Zellkulturen oder Bakterienkulturen zu Forschungszwecken mit dem Partikelstrahl 16 zu bestrahlen.
  • Als Partikel werden vornehmlich Teilchen wie beispielsweise Protonen, Pionen, Heliumionen, Kohlenstoffionen oder Ionen anderer Elemente eingesetzt. Üblicherweise werden derartige Partikel in einer Partikelquelle bzw. Ionenquelle 20 erzeugt. Der von der Ionenquelle 20 erzeugte Ionenstrahl oder Partikelstrahl wird in dem Vorbeschleuniger 22 auf ein erstes Energieniveau beschleunigt. Der Vorbeschleuniger 22 ist beispielsweise ein Linearbeschleuniger. Anschließend werden die Partikel in einen weiteren Beschleuniger 26, beispielsweise einen Kreisbeschleuniger, insbesondere ein Synchrotron oder Zyklotron, eingespeist. In dem Beschleuniger 26 wird der Partikelstrahl auf eine zur Bestrahlung benötigte Energie beschleunigt. Nachdem der Partikelstrahl den Beschleuniger 26 verlassen hat, transportiert ein Hochenergiestrahl-Transportsystem 28 den Partikelstrahl in einen oder mehrere Bestrahlungsräumen 30, 30', 30'', wobei dort beispielsweise die Positioniervorrichtung 12 - etwa eine Patientenliege - mit dem Patienten 14 oder dem nicht-lebenden Objekt 18 angeordnet ist.
  • In dem Bestrahlungsraum 30 oder 30' erfolgt die Bestrahlung des zu bestrahlenden Objektes 14 oder 18 von einer festen Richtung aus, wobei das zu bestrahlende Objekt 14 oder 18 raumfest angeordnet ist. Diese Bestrahlungsräume 30, 30' werden als „fixed beam“-Räume bezeichnet. Im Behandlungsraum 30'' ist eine um eine Achse 32 beweglich angeordnete, vorzugsweise drehbar angeordnete Gantry 34 vorgesehen. Mittels der Gantry 34 kann das zu bestrahlende Objekt 14 oder 18 von verschiedenen Richtungen aus bestrahlt werden. Hierbei wird der Partikelstrahl 16 mittels der in der Gantry 34 angeordneten Gantrystrahlführung 36 um das zu bestrahlende Objekt 14 oder 18 gedreht. Es sind in 1 stellvertretend für die unterschiedlichen Positionen der Gantrystrahlführung 36 der Gantry 34 eine erste Position 38 und eine zweite Position 38' gezeigt. Selbstverständlich sind auch Zwischenpositionen für die Gantrystrahlführung 36, die aus Gründen der Übersichtlichkeit nicht gezeigt sind, auf zumindest einem Halbkreis oberhalb des zu bestrahlenden Objektes 14 oder 18 in einer gedachten Kugel möglich. Somit kann das zu bestrahlende Zielvolumen von mehreren Richtungen aus senkrecht zur Achse 32 bestrahlt werden. Dies ist aus geometrischen Gründen vorteilhaft.
  • Im Bestrahlungsraum 30, 30' tritt der Partikelstrahl aus einem als Strahlauslass 40, 40' bezeichneten Ende eines Vakuumsystems der Hochenergiestrahlführung 28 aus und wird auf das zu bestrahlende Zielvolumen im zu bestrahlenden Objekt 14 oder 18 gelenkt. Das Zielvolumen ist hierbei üblicherweise in einem Isozentrum 42, 42' des jeweiligen Bestrahlungsraums 30, 30' angeordnet.
  • Der anhand der 1 dargestellte Grundaufbau einer Partikeltherapieanlage ist beispielhaft für Partikeltherapieanlagen, kann aber auch hiervon abweichen. Die nachfolgend beschriebenen Ausführungen sind sowohl im Zusammenhang mit der anhand von 1 dargestellten als auch mit anderen Partikeltherapieanlagen einsetzbar.
  • Üblicherweise soll, z.B. bei der Bestrahlung eines Tumors, eine bestimmte Verteilung der Dosis, also eine Zieldosisverteilung, insbesondere eine biologisch effektive Zieldosisverteilung für bestimmte Ionen, über das Zielvolumen erreicht werden. Die Zieldosisverteilung wird z.B. als deponierte Energie pro Volumeneinheit quantifiziert. Verbreitet ist die Angabe der Dosis als Joule pro Kilogramm, dies entspricht der Einheit Gray.
  • Die Steuerung der Bestrahlung, d.h. insbesondere die Lenkung von Partikeln bestimmter Energie in eine bestimmte Richtung, so dass eine bestimmte Dosis an einem bestimmen Zielpunkt innerhalb des Zielvolumens deponiert wird, erfolgt durch die Bestrahlungsvorrichtung BV, welche schematisch in 2 dargestellt ist. Hierzu kann die Bestrahlungsvorrichtung BV z.B. Partikelstrahlablenkeinrichtungen, welche den Partikelstrahl horizontal und/oder vertikal abzulenken vermögen, steuern. Im in der 1 dargestellten System befindet sich die Bestrahlungsvorrichtung BV vor den Strahlauslassen 40, 40'.
  • Vor der Bestrahlung wird eine Bestrahlungsplanung durchgeführt, um dann die Bestrahlung, d.h. das Scannen bzw. Abtasten des Zielvolumens mit dem Partikelstrahl, entsprechend des erstellten Bestrahlungsplans mittels der Bestrahlungsvorrichtung BV durchzuführen und zu steuern. Die Bestrahlungsplanung stellt also die Bestimmung von Steuerparametern zur späteren Steuerung des Bestrahlungsvorgangs dar. Die Bestrahlungsplanung wird mittels einer Bestrahlungsplanungsvorrichtung BPV durchgeführt. Die Bestrahlungsplanungsvorrichtung BPV ist z.B. ein Arbeitsplatz-Computer, eine Workstation oder ein anderer Rechner. Die Bestrahlungsplanungsvorrichtung BPV bestimmt die Steuerparameter, die von der Bestrahlungsvorrichtung BV zur späteren Steuerung und Durchführung der Bestrahlung eingesetzt werden. Die Bestrahlungsvorrichtung BV steuert also den Bestrahlungsverlauf entsprechend der von der Bestrahlungsplanungsvorrichtung BPV ermittelten Steuerparameter.
  • Die Steuerparameter werden zur Durchführung der Bestrahlung an die Bestrahlungsvorrichtung BV weitergeleitet. Selbstverständlich muss die Bestrahlungsplanungsvorrichtung BPV nicht physikalisch an die Bestrahlungsvorrichtung BV angebunden sein. Vielmehr ist es beispielsweise auch möglich, dass die Berechnungsergebnisse der Bestrahlungsplanungsvorrichtung BPV über einen Datenträger zur Bestrahlungsvorrichtung BV transferiert werden. Auch kann zwischen Durchführung der Bestrahlungsplanung und der Bestrahlung ein gewisser Zeitraum, z.B. von mehreren Tagen, liegen.
  • Der Ablauf der Bestrahlungsplanung ist der folgende: zunächst wird mittels eines Computertomographen CT1 ein Volumenbild des zu bestrahlenden Objektes erfasst. Dieses Volumenbild wird im Folgenden als Planungs-CT PL-CT bezeichnet. Das Planungs-CT PL-CT kann hierbei das gesamte zu bestrahlende Objekt abbilden, oder auch nur einen Bereich des Objektes, in welchem sich das Zielvolumen befindet. Zusätzlich zum Planungs-CT PL-CT ist auch der Einsatz eines Kernspin-Tomographen oder anderer diagnostischer Bildgebungsvorrichtungen möglich.
  • Es wird beispielhaft und ohne Einschränkung der Erfindung davon ausgegangen, dass es sich bei dem Zielvolumen um einen Tumor eines Patienten handelt. Aus dem Planungs-CT PL-CT kann die Lage und Ausdehnung des zu bestrahlenden Tumors ermittelt werden. Die Bestrahlungsplanungsvorrichtung BPV erhält Daten des Computertomographen CT1. Hierbei ist es möglich, dass der Bestrahlungsplanungsvorrichtung BPV die erfassten Rohdaten, die Daten in bereits teilweise verarbeiteter Form, oder auch die bereits rekonstruierten Bilder PL-CT des Patienten zur Verfügung gestellt werden.
  • Die Bestrahlungsplanungsvorrichtung BPV ist ferner durch ihre Benutzerschnittstelle I dazu ausgebildet, dass ein Anwender ihr die zu applizierende Dosisverteilung vorgibt. Der Anwender kann also für die Zielpunkte des Zielvolumens eine biologisch wirksame Dosis angeben, die so genannte Zieldosisverteilung. Die Bestrahlungsplanungsvorrichtung BPV ermittelt anhand des Planungs-CTs PL-CT und der vorgegebenen Zieldosisverteilung die Steuerparameter für die Bestrahlung. Insbesondere ist zu ermitteln, wie viele Partikel einer bestimmten Energie in eine bestimmte Richtung auszustrahlen sind. Es ist möglich, dass der Bestrahlungsplanungsvorrichtung BPV weitere Informationen zur Verfügung gestellt werden, welche in die Bestimmung der Steuerparameter für die Bestrahlung eingehen; dies ist für das Verständnis der Erfindung nicht relevant.
  • Um die Reichweite von Partikeln und somit den Ort innerhalb des Patienten, an welchem die jeweilige Dosis deponiert wird, genau zu bestimmen, muss die Zusammensetzung der Materie bzw. des Gewebes, welches die Partikel durchdringen, genau bekannt sein. Hierzu besteht das Planungs-CT PL-CT aus kalibrierten HU (Hounsfield Units), welche den linearen Absorptionskoeffizienten des Materials für Röntgenstrahlung in Bezug auf denjenigen von Wasser angeben. Das Planungs-CT PL-CT besteht also aus Bildpunkten, wobei jedem Bildpunkt ein Bildpunktwert in Form eines HU-Wertes zugeordnet ist. Der HU-Wert ist korreliert mit der Energieabgabe der Partikel bei der Bestrahlung.
  • Nach der Erfassung des Planungs-CTs PL-CT und vor der Durchführung der Bestrahlung kann eine Bewegung des Tumors stattfinden. Eine solche Bewegung kann eine Verschiebung, Drehung oder Formveränderung, sowie eine Kombination hieraus umfassen. Die zwischen verschiedenen Behandlungen bzw. zwischen der Aufnahme des Planungs-CTs PL-CT und einer Behandlung auftretenden Bewegungen werden als interfraktionelle Bewegungen bezeichnet. Sie kommen z.B. zustande durch eine Veränderung des Füllstandes des Magens oder eine Bewegung des Darms.
  • Um derartige interfraktionelle Bewegungen bei der Bestrahlung zu berücksichtigen, ist es vorteilhaft, die Lage des Tumors vor der Bestrahlung erneut zu überprüfen. Denn würde diese Überprüfung unterbleiben und der auf einer nicht mehr aktuellen Tumorposition erstellte Bestrahlungsplan eingesetzt werden, könnte es zu einer falschen Bestrahlung des Patienten kommen, d.h. einer Bestrahlung, die nicht der vorgegebenen Zieldosisverteilung entspricht, wodurch insbesondere der Tumor nicht wie geplant sterilisiert wird und als ungewünschter Nebeneffekt unnötig viel gesundes Gewebe beschädigt werden kann.
  • Es wird daher direkt vor der Bestrahlung erneut ein Bild des Patienten aufgenommen. Hierzu wird die Bildgebungseinrichtung CT2 verwendet; die Ausgestaltung dieser Bildgebungseinrichtung CT2 und/oder die Art der rekonstruierten Bilder kann sich von dem Planungs-CT PL-CT unterscheiden. So sind z.B. zwei- oder dreidimensionale CT-Aufnahmen oder Cone-Beam CT-Aufnahmen möglich, Fluoroskopie-Aufnahmen, Röntgenaufnahmen, . Ultraschallaufnahmen, oder NMR-Aufnahmen, jeweils z.B. aus mehreren Richtungen und/oder in Kombination miteinander. Es können sowohl volumetrische als auch 2D-Datensätze zum Einsatz kommen. Das von der Bildgebungseinrichtung CT2 erhaltene Bild wird im Folgenden - unabhängig davon, ob es ein zwei- oder dreidimensionales Bild ist und auf welche Weise die Daten hierfür erfasst werden - als Kontrollbild K-CT bezeichnet.
  • Von Vorteil bei der im Folgenden erläuterten Vorgehensweise ist, dass das Kontrollbild K-CT der Bildgebungseinrichtung CT2 im Gegensatz zum Planungs-CT PL-CT nicht in Form von kalibrierten HU vorliegen muss. Dies ermöglicht es, Bildgebungseinrichtungen CT2 einzusetzen, welche auf anderen physikalischen Vorgängen als die Computertomographie beruhen und somit keine HU-Werte liefern, wie z.B. NMR. Ferner ermöglicht dies auch, Kontrollbilder K-CT zu verwenden, welche aufgrund von Artefakten oder physikalischen Eigenschaften keine für eine Bestrahlungsplanung verwendbaren HU-Werte enthalten. Ein Beispiel hierfür sind CT-Bilder mit Bewegungsartefakten, welche durch Atem- und/oder Herzbewegung zustande kommen. Diese wären von verminderter oder gegebenenfalls sogar nicht von ausreichender Qualität, um sie als Planungs-CT PL-CT einzusetzen, für die Kontrolle bzw. Überarbeitung des Bestrahlungsplans reicht die Qualität jedoch aus, da nur die Verschiebung des Zielvolumens detektiert und auf das Planungs-CT angewendet werden muss, die verwendeten HUs jedoch vom Planungs-CT PL-CT selbst stammen. Ferner kann der durch das Kontrollbild K-CT abgebildete Ausschnitt des Patienten kleiner als derjenige des Planungs-CTs PL-CT sein. Letzteres ist insbesondere der Fall, wenn als Kontrollbild K-CT eine zweidimensionale Abbildung verwendet wird; denn bei dem Planungs-CT PL-CT handelt es sich um eine Volumenaufnahme.
  • Bevor das Kontrollbild K-CT aufgenommen wird, wird überprüft, dass die Patientenposition derjenigen bei Aufnahme des Planungs-CTs PL-CT entspricht. Besteht das Kontrollbild K-CT beispielsweise aus einem CT-Schnittbild des Patienten, so muss genau bekannt sein, in welcher Höhe innerhalb des dreidimensionalen Planungs-CTs PL-CT dieses Schnittbild anzuordnen ist. Diese korrekte äußere Patientenposition kann mit an sich bekannten Maßnahmen sichergestellt werden, wie z.B. der Verwendung von Markern auf der Patientenoberfläche, der Registrierung der gesamten Patientenanatomie in Bezug auf den Übergang Körper-Luft oder der knöchernen Anatomie.
  • Das Kontrollbild K-CT wird nun verwendet, um den von der Bestrahlungsplanungsvorrichtung BPV berechneten Bestrahlungsplan zu kontrollieren. Diese Kontrolle kann zu einer Verifikation, d.h. zu einer unveränderten Übernahme des Bestrahlungsplans für die Bestrahlung, oder auch zu einer Modifikation, d.h. zu einer Änderung des Bestrahlungsplans vor der Bestrahlung, führen. Die beiden Datensätze, d.h. das Planungs-CT PL-CT und das Kontrollbild K-CT werden miteinander verglichen. Hierbei wird ermittelt, wie der Tumor im Planungs-CT PL-CT verschoben und/oder rotiert und/oder hinsichtlich seiner Form verändert werden müsste, um mit demjenigen des Kontrollbildes K-CT möglichst gut übereinzustimmen. Hierzu können verschiedene mathematische Algorithmen zum Einsatz kommen. Vorteilhafterweise wird im Rahmen eines Optimierungsverfahrens ein Ähnlichkeitsmaß definiert und nach einem Maximum für die Ähnlichkeit zwischen dem Kontrollbild K-CT und dem Planungs-CT PL-CT, welches den bereits verschobenen und/oder rotierten und/oder hinsichtlich seiner Form veränderten Tumor enthält, gesucht. Die besten Ergebnisse werden selbstverständlich erreicht, wenn es sich bei dem Kontrollbild K-CT wie auch bei dem Planungs-CT PL-CT um eine dreidimensionale Aufnahme handelt. Das Vorgehen ist jedoch hierauf nicht beschränkt; auch bei einem zweidimensionalen Kontrollbild K-CT oder insbesondere der Kombination von verschiedenen Aufnahmemodalitaten für das Kontrollbild K-CT kann ein Abgleich mit der Lage, Orientierung und gegebenenfalls auch der Form des Tumors innerhalb des Planungs-CTs PL-CT stattfinden.
  • Es ist möglich, dass sich der Vergleich von Kontrollbild K-CT und Planungs-CT PL-CT auf den Tumor, d.h. dessen Position, Orientierung und Form, beschränkt. Denn das den Tumor umgebende Gewebe und die restliche Anatomie des Patienten verandern sich i.d.R. nur wenig. Betrachtet man beispielsweise einen Lungentumor, so gilt, dass sich die Struktur des den Tumor umgebenden gesunden Weichteilgewebes im Brustbereich, sowie die Form der Rippen üblicherweise nicht oder nur über einen längeren Zeitraum ändern. Diese Bestandteile können bei einem Abgleich von Kontrollbild K-CT und Planungs-CT PL-CT also außer Acht gelassen werden. Der Vorteil hiervon liegt darin, dass für die von der Lage des Tumors nicht betroffenen Teile des Patienten die genauen HU-Werte des Planungs-CTs PL-CT unverändert beibehalten werden können.
  • Bei dem Vergleich zwischen den beiden Datensätzen des Planungs-CTs PL-CT und des Kontrollbildes K-CT wird der Tumor innerhalb des Planungs-CTs PL-CT virtuell bewegt. Hinsichtlich dieser Bewegung des Tumors innerhalb des Planungs-CTs PL-CT ist es möglich, verschiedene Arten von Bewegungen zuzulassen: zum einen existiert die Möglichkeit der rigiden Bewegung. -In diesem Fall wird die Form des Tumors festgehalten und es wird eine Rotation und/oder eine Translation durchgeführt. Alternativ ist auch eine affine Bewegung möglich. In diesem Fall ist neben der Rotation und Translation auch eine Scherung und Skalierung des Tumors möglich. Eine dritte Alternative ist die elastische Bewegung. In diesem Fall sind neben den oben genannten Veränderungen des Tumors beliebige Arten von Deformationen zugelassen.
  • Nachdem durch die virtuelle Bewegung des Tumors innerhalb des Planungs-CTs PL-CT festgestellt wurde, welche Lage, Orientierung und gegebenenfalls auch welche Form der Tumor innerhalb des Planungs-CTs PL-CT haben müsste, um möglichst gut mit dem Kontrollbild K-CT übereinzustimmen, wird das Planungs-CT PL-CT dementsprechend angepasst. Auf diese Weise wird das Planungs-CT CT-PL so geändert, dass es der aktuellen Anatomie entspricht, wie sie durch das Kontrollbild K-CT abgebildet wird. Dies erfolgt, indem Bildpunktwerte, d.h. HU-Werte, des Planungs-CTs PL-CT geändert werden: diejenigen Bildpunkte des Planungs-CTs PL-CT, an welchen sich gemäß dem Vergleich mit dem Kontrollbild K-CT Tumorgewebe befindet, werden mit dem HU-Wert des Tumorgewebes belegt. Dieser HU-Wert ist aus dem ursprünglichen Planungs-CT PL-CT bekannt; gegebenenfalls kann aus mehreren Bildpunkten des Tumorvolumens ein Mittelwert gebildet und dieser verwendet werden. Für manche Bildpunkte bedeutet dies ein Überschreiben der HU-Werte von gesundem Gewe4 be durch diejenigen von Tumorgewebe. Für andere Bildpunkte hingegen mag keine Veränderung der HU-Werte eintreten; hierbei handelt es sich um diejenigen Bildpunkte, bei welchen sowohl gemäß dem ursprünglichen Planungs-CT PL-CT als auch gemäß der aktuellen Anatomie sich (kein) Tumorgewebe befindet.
  • Im Falle eines Lungentumors bedeutet die Änderung des Planungs-CTs PL-CT, dass weniger dichtes Lungengewebe durch die wasserähnliche Dichte von Tumorgewebe (dies entspricht ca. 40 HU) und umgekehrt dichtes Tumorgewebe durch weniger dichtes gesundes Lungengewebe (dies entspricht ca. -800 HU) überschrieben wird. Die „leeren“ Stellen, von denen der Tumor wegverschoben wurde, wird mit der mittleren Dichte von gesundem Lungengewebe aufgefüllt. Der aktuelle Wert des gesunden Lungengewebes ist wenig kritisch, da sich aufgrund seiner geringen Dichte nur kleine Einflüsse auf die Partikelreichweite ergeben.
  • Es werden also die HU-Werte innerhalb des Planungs-CTs PL-CT geändert, so dass der Tumor „künstlich“ an die der aktuellen Anatomie entsprechende Position verschoben wird. Entsprechendes gilt auch für eine Rotation und Formveränderung des Tumors.
  • Das Kontrollbild K-CT wird also nicht direkt zur Berechnung des Bestrahlungsplans verwendet. Vielmehr dient es dazu, das Planungs-CT PL-CT zu verändern, um im Anschluss dieses erneut für eine Berechnung eines neuen Bestrahlungsplans oder die Adaptierung des ursprünglichen Bestrahlungsplans einsetzen zu können. Diese indirekte Verwendung des Kontrollbildes K-CT erlaubt es, die oben beispielhaft aufgezählten verschiedenen Arten von Kontrollbildern K-CT bzw. Bildgebungseinrichtungen CT2 einzusetzen.
  • Nachdem das Planungs-CT PL-CT wie beschrieben überarbeitet wurde, wird erneut ein Bestrahlungsplan errechnet. Dieser basiert auf der aktuellen Anatomie, so dass er für die Bestrahlung verwendet werden kann. Anstelle einer vollständigen Neuberechnung kann auch das Verifizieren, dass der bisherige Bestrahlungsplan verwendbar ist, treten, oder eine Adaption des berechneten Bestrahlungsplans.
  • Um robuster gegenüber Unwägbarkeiten zu werden, ist es sinnvoll, bei Uberschreiben der HU-Werte obere und/oder untere Grenzwerte zu verwenden. Dadurch kann die Reichweite der Partikel bei einer auf dem überarbeiteten Planungs-CT PL-CT basierenden Bestimmung des Bestrahlungsplans beeinflusst werden. Wenn man beispielsweise sicherstellen möchte, dass das distale Ende des Tumors vollständig von der Bestrahlung erfasst wird, ist es vorteilhaft, für das in Partikelstrahlrichtung vor dem Tumor liegende Gewebe einen oberen Grenzwert, d.h. große HU-Werte, zu verwenden. Auf diese Weise wird gewährleistet, dass die Eindringtiefe der Partikel bis zum distalen Ende des Tumors reicht. Soll andererseits sichergestellt werden, dass bestimmtes Gewebe keine Bestrahlungsdosis erhält, kann für das vor diesem sensiblen Gewebe von dem Partikelstrahl passierte Gewebe ein unterer Grenzwert, d.h. niedrige HU-Werte, verwendet werden.
  • Außer der bisher beschriebenen interfraktionellen Bewegung eines Tumors stellt auch die intrafraktionelle Bewegung, d.h. die Bewegung des Tumors während der Bestrahlung, ein Problem dar. Diese Bewegung kommt hauptsächlich durch die Atmung des Patienten zustande. Bei der Bestrahlung von bewegten Zielvolumina werden üblicherweise Methoden eingesetzt, um diese Bewegung zu kompensieren. Denn ohne eine Berücksichtigung der Bewegung kann die gewünschte Dosisverteilung im zu bestrahlenden Objekt in der Regel nicht erreicht werden. Bekannte Verfahren zur Berücksichtigung der Bewegung sind die zur Bewegung synchronisierte Strahlapplikation des Gating, bei welchem der Strahl nur eingeschaltet wird, wenn ein bestimmter Bewegungszustand und somit eine bestimmte Position des Tumore vorliegt, und des Tracking, bei welchem der Strahl der Bewegung des Tumors nachgeführt wird, sowie das Rescanning.
  • Die vorgestellte Vorgehensweise kann auch auf die Berücksichtigung der intrafraktionellen Bewegung angewandt werden. Das Planungs-CT PL-CT enthält, wenn man die intrafraktionelle Bewegung während der Bestrahlung berücksichtigen mochte, eine zeitliche Abfolge von Bildern, so dass die Bewegungstrajektorie des Tumors bekannt ist. Es handelt sich in diesem Fall also um eine vierdimensionale CT-Aufnahme. Das Kontrollbild K-CT kann auch eine zeitliche Abfolge von Kontrollbildern enthalten, so dass die Bewegungstrajektorie des Planungs-CTs PL-CT vollständig mit derjenigen der Kontrollbilder K-CT verglichen werden kann. Durch Vergleich der beiden Datensätze kann auf die oben für den statischen Fall beschriebene Weise die Bewegungstrajektorie innerhalb des Planungs-CTs PL-CT den aktuellen Gegebenheiten angepasst werden.
  • Oftmals ist es jedoch so, dass sich die Lage eines intrafraktionell bewegten Tumors zwar auch interfraktionell ändert, z.B. durch eine Änderung des Füllstandes des Magens, die Form der Bewegungstrajektorie jedoch annähernd gleich bleibt. Dies wurde z.B. in der Veröffentlichung
    Jan-Jakob Sonke, Ph. D. Joos Lebesque, Ph. D. Marcel van Herk: „Variability of Four-Dimensional Computed Tomography Patient Models", Int. J. Radiation Oncology Biol. Phys., Vol. 70, No. 2, pp.590-598, 2008
    festgestellt. Es genügt daher, wenn das Kontrollbild K-CT die Position des Tumors zu einem definierten Zeitpunkt innerhalb des Bewegungszyklus abbildet. Dieser definierte Zeitpunkt kann dann für einen Abgleich mit dem Planungs-CT PL-CT verwendet werden. Die Lage, Orientierung und gegebenenfalls auch die Form des Tumors wird im Planungs-CT PL-CT der aktuellen Anatomie angepasst, und von dieser angepassten Anordnung des Tumors aus verläuft die aus dem ursprünglichen Planungs-CT PL-CT bekannte Bewegungstrajektorie.
  • Die Erfindung wurde voranstehend an einem Ausführungsbeispiel beschrieben. Es versteht sich, dass zahlreiche Änderungen und Modifikationen möglich sind, ohne dass der Rahmen der Erfindung verlassen wird.

Claims (21)

  1. Verfahren zur Bildbearbeitung, wobei - ein erstes von einem medizinischen Bildaufnahmegerät (CT1) erfasstes Bild (PL-CT) eines Objektes, welches ein Zielobjekt enthält, vorliegt, umfassend Bildpunkte und ihnen zugeordnete Bildpunktwerte, und ein zu einem anderen Zeitpunkt erfasstes zweites Bild (K-CT) des Objektes vorliegt, - das erste Bild (PL-CT) und das zweite Bild (K-CT) hinsichtlich des Zielobjektes miteinander verglichen werden, - unter Berücksichtigung des Vergleichs eine Bewegung des Zielobjektes innerhalb des ersten Bildes (PL-CT) durchgeführt wird, indem eine Neuzuordnung von Bildpunktwerten zu Bildpunkten vorgenommen wird, - als erstes Bild (PL-CT) eine zeitliche Abfolge von Bildern vorliegt, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen, - das zweite Bild (K-CT) nur einen Zeitpunkt abbildet, und die Neuzuordnung für alle Bilder (PL-CT) der zeitlichen Abfolge vorgenommen wird, wobei von einer bei Aufnahme des zweiten Bildes (K-CT) gegenüber der Aufnahme der Abfolge der ersten Bilder (PL-CT) unveränderten Form der Bewegungstrajektorie ausgegangen wird.
  2. Verfahren nach Anspruch 1, wobei es sich bei dem ersten Bild (PL-CT) um ein CT-Bild handelt.
  3. Verfahren nach Anspruch 1 oder 2, wobei die Bewegung des Zielobjektes innerhalb des ersten Bildes (PL-CT) aus einer Verschiebung und/oder Drehung und/oder Formveränderung besteht.
  4. Verfahren nach einem der Ansprüche 1 bis 3, wobei durch die Bewegung des Zielobjektes innerhalb des ersten Bildes (PL-CT) dieses an das zweite Bild (K-CT) angenähert wird.
  5. Verfahren nach einem der Ansprüche 1 bis 4, wobei das zweite Bild (K-CT) zu einem späteren Zeitpunkt erfasst wurde als das erste Bild (PL-CT) und das Zielobjekt zwischenzeitlich eine Bewegung innerhalb des Objektes erfahren hat.
  6. Verfahren nach einem der Ansprüche 1 bis 5, wobei bei der Neuzuordnung dem Zielobjekt entsprechende Bildpunktwerte zu gemäß dem ursprünglichen ersten Bild (PL-CT) nicht das Zielobjekt darstellenden Bildpunkten zugeordnet werden.
  7. Verfahren nach einem der Ansprüche 1 bis 6, wobei bei der Neuzuordnung nicht dem Zielobjekt entsprechende Bildpunktwerte zu gemäß dem ursprünglichen ersten Bild (PL-CT) das Zielobjekt darstellenden Bildpunkten zugeordnet werden.
  8. Verfahren nach einem der Ansprüche 1 bis 7, wobei bei der Neuzuordnung ein oberer und/oder ein unterer Bildpunktwertgrenzwert für nicht das Zielobjekt darstellende Bildpunkte eingesetzt werden.
  9. Verfahren nach einem der Ansprüche 1 bis 8, wobei bei dem Vergleich des ersten Bildes (PL-CT) mit dem zweiten Bild (K-CT) ein Algorithmus zur Optimierung eines Ähnlichkeitsmaßes eingesetzt wird.
  10. Verfahren nach einem der Ansprüche 1 bis 9, wobei der Vergleich des ersten Bildes (PL-CT) mit dem zweiten Bild (K-CT) auf das Zielobjekt beschränkt ist.
  11. Verfahren nach einem der Ansprüche 1 bis 10, wobei es sich bei dem ersten Bild (PL-CT) um ein das Zielobjekt enthaltendes Volumenbild und bei dem zweiten Bild (K-CT) um ein zweidimensionales Bild des Objektes handelt.
  12. Verfahren nach einem der Ansprüche 1 bis 11, wobei es sich bei dem ersten Bild (PL-CT) und dem zweiten Bild (K-CT) jeweils um ein CT-Bild handelt, und das zweite gegenüber dem ersten von schlechterer Qualität ist.
  13. Verfahren nach einem der Ansprüche 1 bis 12, wobei es sich bei dem zweiten Bild (K-CT) um ein Fluoroskopie-Bild oder ein Röntgenbild oder ein Ultraschallbild oder ein NMR-Bild handelt.
  14. Verwendung des Verfahren nach einem der Ansprüche 1 bis 13 zur Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage, wobei die Bestimmung der Steuerparameter auf Grundlage des durch die Neuzuordnung von Bildpunktwerten veränderten ersten Bildes (PL-CT) erfolgt.
  15. Verwendung nach Anspruch 14, wobei zuerst eine Bestimmung der Steuerparameter auf Grundlage des ersten Bildes (PL-CT) erfolgt, und im Anschluss eine Neubestimmung oder eine Verifikation oder eine Veränderung der bestimmten Steuerparameter auf Grundlage des durch die Neuzuordnung von Bildpunktwerten veränderten ersten Bildes (PL-CT) erfolgt.
  16. Vorrichtung (BPV) zur Bestimmung von Steuerparametern für eine Partikelbestrahlungsanlage (BV), wobei mit der Partikelbestrahlungsanlage (BV) verschiedene Dosiswerte an unterschiedlichen Zielpunkten in einem Zielobjekt deponierbar sind, umfassend: einen Eingang zum Empfangen eines ersten von einem medizinischen Bildaufnahmegerät (CT1) erfassten Bildes (PL-CT) eines Objektes, welches das Zielobjekt enthält, umfassend Bildpunkte und ihnen zugeordnete Bildpunktwerte, einen Eingang zum Empfangen eines zu einem anderen Zeitpunkt erfassten zweiten Bildes (K-CT) des Objektes, einen Bildbearbeitungsbestandteil zum Vergleichen des ersten Bildes (PL-CT) mit dem zweiten Bild (K-CT) hinsichtlich des Zielobjektes und zum Durchführen einer Bewegung des Zielobjektes innerhalb des ersten Bildes (PL-CT), indem unter Berücksichtigung des Vergleichs eine Neuzuordnung von Bildpunktwerten zu Bildpunkten vorgenommen wird, und einen Ermittlungsbestandteil zum Bestimmen der Steuerparameter für die Partikelbestrahlungsanlage auf Grundlage des durch die Neuzuordnung von Bildpunktwerten veränderten ersten Bildes (PL-CT), wobei als erstes Bild (PL-CT) eine zeitliche Abfolge von Bildern vorliegt, welche eine Bewegungstrajektorie des Zielobjektes aufzeigen, das zweite Bild (K-CT) nur einen Zeitpunkt abbildet, und die Neuzuordnung für alle Bilder (PL-CT) der zeitlichen Abfolge vorgenommen wird, wobei von einer bei Aufnahme des zweiten Bildes (K-CT) gegenüber der Aufnahme der Abfolge der ersten Bilder (PL-CT) unveränderten Form der Bewegungstrajektorie ausgegangen wird.
  17. Partikelbestrahlungsanlage mit einer Vorrichtung (BPV) zur Bestimmung von Steuerparametern nach Anspruch 16.
  18. Verfahren zur Bestrahlung eines Zielobjektes mit Partikeln unter Verwendung von Steuerparametern zur Steuerung einer Partikelbestrahlungsanlage, wobei die Steuerparameter gemäß Anspruch 14 oder 15 ermittelt sind.
  19. Verfahren nach Anspruch 18, wobei das Zielobjekt Bestandteil eines nicht-lebenden Objektes ist.
  20. Computerprogramm mit Programmcode-Mitteln, um das Verfahren nach einem der Ansprüche 1 bis 13 durchzuführen, wenn das Computerprogramm auf einem Computer ausgeführt wird.
  21. Computerprogrammprodukt, umfassend auf einem computerlesbaren Datenträger gespeicherte Programmcode-Mittel eines Computerprogramms, um das Verfahren nach einem der Ansprüche 1 bis 13 durchzuführen, wenn das Computerprogramm auf einem Computer ausgeführt wird.
DE102009048070.6A 2009-10-01 2009-10-01 Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens Active DE102009048070B4 (de)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE102009048070.6A DE102009048070B4 (de) 2009-10-01 2009-10-01 Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens
EP10750085A EP2482728A1 (de) 2009-10-01 2010-08-24 Optimierung von steuerparametern für eine partikelbestrahlungsanlage unter berücksichtigung von interfraktionellen bewegungen eines zielvolumens
PCT/EP2010/062292 WO2011038999A1 (de) 2009-10-01 2010-08-24 Optimierung von steuerparametern für eine partikelbestrahlungsanlage unter berücksichtigung von interfraktionellen bewegungen eines zielvolumens

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102009048070.6A DE102009048070B4 (de) 2009-10-01 2009-10-01 Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens

Publications (2)

Publication Number Publication Date
DE102009048070A1 DE102009048070A1 (de) 2011-04-07
DE102009048070B4 true DE102009048070B4 (de) 2020-08-20

Family

ID=43244782

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102009048070.6A Active DE102009048070B4 (de) 2009-10-01 2009-10-01 Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens

Country Status (3)

Country Link
EP (1) EP2482728A1 (de)
DE (1) DE102009048070B4 (de)
WO (1) WO2011038999A1 (de)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000067202A1 (en) * 1999-04-29 2000-11-09 Ge Medical Systems Sa Ct/mr fused image presentation algorithm
US20090087124A1 (en) * 2007-09-28 2009-04-02 Varian Medical Systems Finland Radiation systems and methods using deformable image registration
US20090116616A1 (en) * 2007-10-25 2009-05-07 Tomotherapy Incorporated System and method for motion adaptive optimization for radiation therapy delivery

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006018761A1 (en) * 2004-08-13 2006-02-23 Koninklijke Philips Electronics N.V. Radiotherapeutic treatment plan adaptation
US8989349B2 (en) * 2004-09-30 2015-03-24 Accuray, Inc. Dynamic tracking of moving targets
WO2007014106A2 (en) * 2005-07-22 2007-02-01 Tomotherapy Incorporated System and method of delivering radiation therapy to a moving region of interest

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000067202A1 (en) * 1999-04-29 2000-11-09 Ge Medical Systems Sa Ct/mr fused image presentation algorithm
US20090087124A1 (en) * 2007-09-28 2009-04-02 Varian Medical Systems Finland Radiation systems and methods using deformable image registration
US20090116616A1 (en) * 2007-10-25 2009-05-07 Tomotherapy Incorporated System and method for motion adaptive optimization for radiation therapy delivery

Also Published As

Publication number Publication date
DE102009048070A1 (de) 2011-04-07
WO2011038999A1 (de) 2011-04-07
EP2482728A1 (de) 2012-08-08

Similar Documents

Publication Publication Date Title
DE102006044139B4 (de) Strahlentherapieanlage und Verfahren zur Anpassung eines Bestrahlungsfeldes für einen Bestrahlungsvorgang eines zu bestrahlenden Zielvolumens eines Patienten
EP2158939B1 (de) Verfahren und Vorrichtung zur Auswahl eines Bestrahlungsplans sowie Bestrahlungsanlage
EP2352556B1 (de) Bestrahlung von zumindest zwei zielvolumen
DE102010061121B4 (de) Bestrahlungsphantom mit zumindest einer Bewegungsvorrichtung zur Bewegung eines ersten Teilbereichs
EP2344250B1 (de) Bestrahlung eines zielvolumens unter berücksichtigung eines zu schützenden volumens
DE102011000204B4 (de) Erstellung einer Bestrahlungsplanung unter Berücksichtigung der Auswirkungen zumindest einer Unsicherheit
DE102008036478A1 (de) Vorrichtung und Verfahren zur Auswertung einer Aktivitätsverteilung sowie Bestrahlungsanlage
DE102008019128A1 (de) Vorrichtung zur Durchführung einer Bestrahlung und Verfahren zur Überwachung einer solchen
WO2010043340A1 (de) Vorrichtung und verfahren zur bestimmung von steuerparametern für eine bestrahlungsanlage, bestrahlungsanlage und bestrahlungsverfahren
DE10250655A1 (de) Patientenpositionierungssystem, das eine Oberflächenfotogrammetrie verwendet
DE102014210458B4 (de) Ermittlung einer Position einer zu bestrahlenden Zielregion eines Patienten in einer Bestrahlungseinrichtung
DE102011083854B4 (de) Zeitaufgelöste Tomosynthesebildgebung
EP2790784B1 (de) Verfahren und vorrichtung zur erstellung eines bestrahlungsplans bei bewegtem zielvolumen ohne bewegungskompensation
EP1348394A1 (de) Planungs- bzw. Navigationsunterstützung durch generische und erfasste Patientdaten mit zweidimensionaler Anpassung
DE102011085773B4 (de) Dosisrekonstruktion bei der Strahlentherapie
EP2366428B1 (de) Verfahren zum Registrieren eines ersten Abbildungsdatensatzes zu einem zweiten Abbildungsdatensatz
DE10139934B4 (de) Integrierte 3D Online Bildgebung für die Positionierung bei der Strahlentherapie
DE102009048070B4 (de) Optimierung von Steuerparametern für eine Partikelbestrahlungsanlage unter Berücksichtigung von interfraktionellen Bewegungen eines Zielvolumens
EP2246096B1 (de) Verfahren und vorrichtung zur bestimmung bevorzugter ausrichtungen eines behandlungsstrahlerzeugers
DE102011081422B4 (de) Bildgestützte Strahlentherapie
DE102010034917B4 (de) Verfahren und Vorrichtung zum Ermitteln des Einflusses von Abschattungen auf eine Röntgenbildgebung, sowie entsprechender Datenträger und Computertomograph
DE102011080364B4 (de) Kontrastmittel-verstärkte Bildgebung bei der Strahlentherapie
EP2482927B1 (de) Verfahren zur erstellung eines bestrahlungsplans, verfahren zur anpassung und optimierung einer bestrahlungsplanung, bestrahlungsplanungseinrichtung sowie bestrahlungsanlage
DE102011085946B4 (de) Bildgebungsverfahren bei einem Strahlentherapiegerät und Strahlentherapiegerät
DE102006012970A1 (de) Verfahren zur Ermittlung des Übereinstimmungsgrades zwischen medizinischen Bildern und Vorrichtungen zur Durchführung des Verfahrens

Legal Events

Date Code Title Description
OP8 Request for examination as to paragraph 44 patent law
R082 Change of representative

Representative=s name: BLUMBACH ZINNGREBE, DE

R081 Change of applicant/patentee

Owner name: GSI HELMHOLTZZENTRUM FUER SCHWERIONENFORSCHUNG, DE

Free format text: FORMER OWNERS: GSI HELMHOLTZZENTRUM FUER SCHWERIONENFORSCHUNG GMBH, 64291 DARMSTADT, DE; SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE

Effective date: 20130215

Owner name: GSI HELMHOLTZZENTRUM FUER SCHWERIONENFORSCHUNG, DE

Free format text: FORMER OWNER: GSI HELMHOLTZZENTRUM FUER SCHWER, SIEMENS AKTIENGESELLSCHAFT, , DE

Effective date: 20130215

R082 Change of representative

Representative=s name: BLUMBACH ZINNGREBE PATENT- UND RECHTSANWAELTE, DE

Effective date: 20130215

Representative=s name: BLUMBACH ZINNGREBE, DE

Effective date: 20130215

Representative=s name: BLUMBACH ZINNGREBE PATENT- UND RECHTSANWAELTE , DE

Effective date: 20130215

R082 Change of representative

Representative=s name: BLUMBACH ZINNGREBE PATENT- UND RECHTSANWAELTE, DE

Representative=s name: BLUMBACH ZINNGREBE PATENT- UND RECHTSANWAELTE , DE

R016 Response to examination communication
R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final