300 likes | 433 Views
Benchmark of the 4D Treatment Planning at PSI. Danilo Andrea Zanin. 1. Danilo Andrea Zanin, June 2010. y z x. y z x. motion file time displacement 1s Δ x Δ y Δ z 2s Δ x Δ y Δ z 3s Δ x Δ y Δ z 4s Δ x Δ y Δ z.
E N D
Benchmark of the 4D Treatment Planning at PSI Danilo Andrea Zanin 1 Danilo Andrea Zanin, June 2010
y z x y z x motion file timedisplacement 1s ΔxΔyΔz 2s ΔxΔyΔz 3s ΔxΔyΔz 4s ΔxΔyΔz ... motion file timedisplacement 1s ΔxΔyΔz 2s ΔxΔyΔz 3s ΔxΔyΔz 4s ΔxΔyΔz ... ... 4s ... 4s 4D Treatment Planning: Version 1 Antje-Christin Knopf, May 2010 Body moves rigidly 2 Antje-Christin Knopf, May 2010
z z y y x x 1s motion map 1s motion map ... 4s ... 4s 2s motion map 2s motion map 3s motion map 3s motion map 4s motion map 4s motion map 1s motion map 1s motion map ... 4s ... 4s 2s motion map 2s motion map 3s motion map 3s motion map 4s motion map 4s motion map • 4D Treatment Planning: Version 2 • Antje-Christin Knopf, May 2010 • Body described with voxels • It will take into account • rigid/non-rigid • amplitude • period • Phase Motion divided in steps Motion matrices Third version Varying densities z z y y x x 3 Antje-Christin Knopf, May 2010
Motivation • Control the rigid motion • Benchmark the implementations for the non rigid motion 4
Motion Matrices • Defined in the gantry coords system: s,t,u • Divide the motion in txt-file steps Generated with MatLab 5
Motion types • Fixed body • Instant Motion The body is misplaced • Continuous Motion Natural motions 6
Instant motion in the beam direction (s) and in orthogonal to it (t, u) • Motion along “s” will not influence the calculated dose only air change! • Motion orthogonal to “s” will move the calculated dose in the opposite direction new tissue will be exposed to the pencil beam. 7
Instant motion in u-direction of 3 cm instant motion of 3 cm uChange_150610_3 Reference 9
Summary of the instant motion in u-direction 3 cm as upper limit 10
Constant motion • Unidirectional Motion • Sinusoidal Motion Constant motion in one direction will spread the dose but less than the instant one 11
Constant motion in u-dir. of 0.2 cm 2 cm Constant motion of 0.2 cm uChange_170610_m1 Instant motion of 2cm uChange_110610_10 12
Sinusoidal motion in u-dir. of 1 cm +3/-3cm Sinusoidal motion of 1 cm uSway_170610_5 Reference 13
Instant union/division • AIR gap VS Tissue superposition • AIR gap: less tissue exposed dose calculation reduced • Tissue superposition: more tissue exposeddose calculation enhanced 15
Instant union/division in u-direction of 3 cm Frontal view Instant division of 3 cm uDevide_160610_3 Instant union of 3 cm uUnion_160610_3 Sagittal view 16
Instant union/division in s-direction of 3 cm Instant union of 3 cm sDevide_170610_1 Instant division of 3 cm sDevide_170610_2 No changes observed! 18
Other Special Motions • Instant slide on an plane • Instant explosion / implosion • Pulse 19
Instant slide in t-direction of 3 cm Transaxial view Instant slide of +3/-3 cm sSlide_170610_1 Reference Sagittal view 20
4D treatment planning Instant implosion in all-direction of 3 cm Instant explosion of 3 cm explode_210610_2 Reference 21
4D treatment planning – Instant explosion in all-direction of 3 cm Instant explosion in all-direction of 3 cm Instant explosion of 3 cm explode_210610_1 Reference 22
Conclusions • Good response to rigid body motions • Good adaptability to non rigid body motions • Displacement threshold (3 cm) • Do not care about motions in depth • Do not care about density changes 23
Relation between coordinates systems Gantry coords syst. Patient coords syst. 26
References for all the results • All motion matrices are stored in: “/asm/pdata/EXP/AKT030/Motionvectors_DZ/matlab” • All the dose are calculated with: “/asm/apps/dev/ptplan04/WIDGETS/DOSE_FIELD/4d/DZ/script_M1.sh” and are saved in:“/asm/pdata/EXP/AKT030/Zanin in the same directory can be found also a file: “info.txt” containing a short description of the different calculated dose 27