Inclusion of organ motion in IMRT optimization using probabilistic treatment planning [Elektronische Ressource] / presented by Jan Unkelbach

Sujets

• Gesundheit

Informations

• Publié par: ruprecht-karls-universitat_heidelberg
• Publié le: 01 janvier 2006
• Nombre de lectures: 35
• Langue: Deutsch
• Poids de l'ouvrage: 3 Mo

Extrait

Dissertation
submitted to the
Combined Faculties for the Natural Sciences and for Mathematics
of the Ruperto-Carola University of Heidelberg, Germany
for the degree of
Doctor of Natural Sciences
presented by
Diplom-Physiker: Jan Unkelbach
born in: Darmstadt, Germany
Oral examination: February, 15th 2006Inclusion of organ motion
in IMRT optimization
using probabilistic treatment planning
Referees: Prof. Dr. Uwe Oelfke
Prof. Dr. Josef BilleZusammenfassung
Beruc ksichtigung von Organbewegungen in der IMRT Optimierung
durch probabilistische Bestrahlungsplanung
Die vorliegende Arbeit beschreibt ein Verfahren zur Beruc ksichtigung von Organbewegun-
geninderBestrahlungsplanungfur diefraktionierteintensit atsmodulierteStrahlentherapie
(IMRT). Organbewegungen werden durch ein mathematisches Modell beschrieben welches
dieGrundlagederBestrahlungsplan-Optimierungdarstellt. DasModellenth altZufallsvari-
ablenumdiestochastischeNaturvonOrganbewegungenzubeschreiben. Dievorausgesagte
DosisverteilungimPatientenmussdaherebenfallsalsZufallsvariableaufgefasstwerdenund
wird durch einen Erwartungswert der Dosis und dessen Varianz charakterisiert. Zur Op-
timierung des Bestrahlungsplans unter Beruc ksichtigung des Bewegungsmodells wird der
Erwartungswert einer quadratischen Kostenfunktion minimiert, der sich als Summe der
Dosisvarianz und derhen Di erenz von Solldosis und Erwartungswert der Dosis
darstellen l asst. Die daraus resultierenden Bestrahlungspl ane haben die Eigenschaft, dass
Bereiche in denen Tumorgewebe nur relativ selten lokalisiert ist mit einer geringeren Do-
sis bestrahlt werden. Dies wird ausgeglichen durch eine Dosisub erh ohung in benachbarten
Bereichen, sodass sichdurch denEin uss derBewegung imVerlaufdergesamten Therapie
eine n aherungsweise homogene Gesamtdosisverteilung im Tumor ergibt. DasVerfahren er-
laubt eine potentiell bessere Schonung von angrenzenden gesunden Geweben im Vergleich
zur Sicherheitsrand-Methode.
Abstract
Inclusion of organ motion in IMRT optimization
using probabilistic treatment planning
The presented thesis describes an o -line approach to incorporate organ motion into the
treatment plan optimization for fractionated intensity modulated radiotherapy (IMRT).
Organ movement is described in terms of a mathematical model that represents the basis
of the treatment plan optimization process. The motion model contains random variables
in order to describe the stochastic nature of organ movements. As a consequence, the
predicted dose distribution in the patient must be considered as a random variable as well.
It is characterized by the expectation value and the variance of the dose. For treatment
plan optimization incorporating the motion model, the expectation value of a quadratic
cost function is minimized, which can be expressed as the sum of the variance of the dose
and the quadratic di erence of expected and prescribed dose. The resulting treatment
plans show a reduction of the dose in regions where tumor tissue is only rarely present.
This is compensated for by delivering a higher dose to neighboring regions that are mostly
occupied by tumor tissue. Due to organ movement during the course of treatment, a
widely homogeneous cumulative dose distribution is delivered to the tumor. This method,
compared to the standard safety margin approach, potentially allows for a better sparing
of healthy tissues from dose burden.ii
Acknowledgments
I would like to thank my supervisor Prof. Uwe Oelfke for providing me the opportunity
to work on this subject. I always enjoyed working with him and his group. I appreciate
hisfriendlypersonality,whichisoneofthereasonsforthepleasantatmosphereinhisgroup.
I am grateful to all colleagues in the department who made my stay at the DKFZ a very
enjoyable period of my life, both at work and at various social events. In particular, I
would like to thank Simeon Nill and Jan Wilkens for their exceptional helpfulness and
their valuable support in countless situations.
IwouldliketothankProf. WolfgangSchlegelforhissupportandforprovidingaproductive
and friendly work environment in his department. I also thank Prof. Josef Bille as the
second referee of this thesis. I am grateful to Annette Seeber for her careful proof reading
and to Daniel Maleike for his work as a diploma student.iii
Preface
This PhD thesis was carried out at the German Cancer Research Center (DKFZ) between
November2002andDecember2005. Thethesisdescribesanapproachtoincorporateorgan
movements thatoccurduringthe course offractionatedradiotherapy intotheoptimization
of a treatment plan for intensity modulated radiotherapy. This approach is referred to as
probabilistic treatment planning.
The manuscript does not contain a general introduction to radiotherapy and assumes the
reader to be familiar with the concept of intensity modulated radiation therapy (IMRT).
A brief introduction to di eren t methods to deal with organ motion in radiotherapy is
provided in order to integrate this work into the eld of Adaptive Radiotherapy (section
1.1).
Parts of this work have already been published as a series of three peer-reviewed papers.
Chapter 2 is mainly covered by the papers [1] and [2], and section 3.2 is widely covered
by the paper [3]. Parts of the work were also presented at v e international conferences
[4, 5, 6, 7, 8]. For the contribution at the ICCR Congress 2004 in South Korea [5], a
young investigator’s prize was awarded. Under the supervision of the author and Prof.
Uwe Oelfke, one diploma thesis was prepared which is closely related to Chapter 3 of this
work [9].Contents
1 Introduction 1
1.1 Organ motion in radiotherapy . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.1.1 Classi cation of organ movements . . . . . . . . . . . . . . . . . . . 1
1.1.2 Dealing with organ motion . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Philosophy of probabilistic treatment planning . . . . . . . . . . . . . . . . 3
1.3 A mathematical paradigm . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4 A probabilistic optimization postulate . . . . . . . . . . . . . . . . . . . . . 8
2 Probabilistic optimization for interfractional motion 11
2.1 The motion model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.1.1 Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.1.2 Maximum likelihood estimate of the most probable data model . . . 13
2.1.3 The concept of Bayesian inference . . . . . . . . . . . . . . . . . . . 14
2.2 An idealized geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.2.1 The static solution of the inverse problem . . . . . . . . . . . . . . 19
2.3 The optimization problem . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.4 Probabilistic treatment planning incorporating random errors . . . . . . . 21
2.4.1 Optimization of the expectation value . . . . . . . . . . . . . . . . . 23
2.5 Incorporating systematic errors into the optimization . . . . . . . . . . . . 24
2.6 Accounting for an uncertain magnitude of motion . . . . . . . . . . . . . . 27
2.6.1 Inclusion of population based knowledge . . . . . . . . . . . . . . . 30
2.7 Adapting uence pro les . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.7.1 Comparison to section 2.5 . . . . . . . . . . . . . . . . . . . . . . . 35
2.7.2 Further comments . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3 Application to prostate cancer 39
3.1 Preliminaries . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
3.2 Accounting for random errors . . . . . . . . . . . . . . . . . . . . . . . . . 41
3.2.1 Implementation of probabilistic optimization for random errors . . . 41
3.2.2 Treatment planning for prostate patients . . . . . . . . . . . . . . . 47
3.2.3 Results: N =30. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
3.2.4 N =1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48
3.2.5 Results: N !1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51vi CONTENTS
3.2.6 Results: small and large organ movements . . . . . . . . . . . . . . 52
3.2.7 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.3 Dose calculation in the presence of systematic errors . . . . . . . . . . . . . 56
3.3.1 A simulation and visualization tool for assessing dose uncertainties. 56
3.3.2 Implementation into the optimization tool KonRad . . . . . . . . . 59
3.4 Di culties in treatment plan evaluation and visualization . . . . . . . . . . 60
3.5 The impact of systematic errors on the optimized treatment plan . . . . . 62
3.6 The of di eren t sources of uncertainty . . . . . . . . . . . . . . . . 65
4 Dealing with respiratory motion 69
4.1 The idealized geometry . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
4.2 The respiratory motion model . . . . . . . . . . . . . . . . . . . . . . . . . 70
4.3 Uncertainties due to a nite irradiation time . . . . . . . . . . . . . . . . . 71
4.3.1 The expectation value of the dose . . . . . . . . . . . . . . . . . . . 72
4.3.2 The variance due to a nite irradiation time . . . . . . . . . . . . . 73
4.3.3 Approximation of the standard deviation . . . . . . . . . . . . . . . 75
4.4 Uncertainties due to variations in amplitude and exhale position . . . . . . 77
4.4.1 The expectation value of the dose . . . . . . . . . . . . . . . . . . . 77