Tomographic images by proton Computed Tomography ... - IEEE Xplore

2 downloads 118 Views 824KB Size Report
Tomographic images by proton Computed. Tomography system for proton therapy applications. Valeria Sipala, Mara Bruzzi, Mirko Brianzi, Marta Bucciolini, ...
MICI5.S-305

2011 IEEE Nuclear Science Symposium Conference Record

Tomographic images by proton Computed Tomography system for proton therapy applications Valeria Sipala, Mara Bruzzi, Mirko Brianzi, Marta Bucciolini, Giuseppe Antonio Pablo Cirrone, Carlo Civinini, Giacomo Cuttone, Domenico Lo Presti, Stefania Pallotta, Nunzio Randazzo, Francesco Romano, Monica Scaringella, Concetta Stancampiano, Cinzia Talamonti, Mauro Tesi.

Abstract-Proton therapy is a highly precise form of cancer treatment,

which

requires

accurate

knowledge

of

the

dose

silicon

tracker

and

a

Y AG:Ce

calorimeter.

The

image

reconstruction technique uses data of the single particle

delivered to the patient and verification of the correct patient

trajectory in order to overcome the limits due to the Multiple

position

Coulomb Scattering (MCS) of the proton in the matter. The

to

avoid

damage

to

critical

normal

tissues.

The

development of pCT (proton Computed Tomography) system represents an important feature for precise proton radiation treatment

planning

measurement

of

because

the

it

proton

could

stopping

permit power

the

direct

distribution,

improving the accuracy in dose calculus, and the patient's

apparatus was tested at Laboratori Nazionali del Sud (LNS) using a 60MeV proton beam. During this test radiographic images of PMMA phantoms were taken. Phantoms with different density zones have been used to estimate the system

to

performances (spatial and energy resolution). Moreover, a

demonstrate the capability to acquire, during treatments in

tomographic image was acquired mounting a cylindrical

position.

A

pCT

prototype

was

manufactured

in

order

proton therapy centers, radiographic and tomographic images

phantom, with embedded in-homogeneities, on a rotating

according to clinical demands.

system to provide multiple projection angles. Preliminary results

I.

INTRODUCTION

Computed Tomography

(pCT)

is

an

imaging

Ptechnique able to reconstruct 3D images containing proton ROTON

obtained using the Filter Back Projection

(FBP)

reconstruction algorithm. Different reconstruction algorithms, based on the tracking information, are presently under test.

stopping power information that can be used to simulate and II.

plan the proton treatment. Nowadays X-ray tomography (X­ CT) is used for proton treatment planning but there are fundamental

differences

in

the

interaction

mechanisms

The

apparatus,

ApPARATUS ARCHITECTURE developed by

Italian

PRIMA

(PRoton

IMAging) collaboration, is able to acquire tomography images

between X-rays and protons. The conversion of the Hounsfield

by proton beam. The goal is to demonstrate that it meets

coefficients from the X-ray CT into proton stopping power

clinical demands for the use of protons in radiotherapy

introduces an intrinsic error (1]. A pCT system, using the same therapeutic proton beam, could produce images avoiding these conversion uncertainties resulting in a better dose accuracy. Starting from these potential advantages of pCT and from the growing interest about proton therapy treatments, a 2 pCT prototype with 5x5cm active area has been designed, manufactured and tested [2-7]. This prototype consists of a Manuscript received November 14, 2011. This work was supported by INFN (CSNS) and MIUR]RIN 2006. V.Sipala is with the Istituto Nazionale di Fisica Nucleare (INFN) of Catania and University of Sassari, 4, Via Piandanna, 1-07100 Sassari, Italy (telephone: +39079229488, e-mail: [email protected]). M. Bruzzi, M. Scaringella and M.Tesi are with INFN of Florence and Energetic Department, University of Florence, Italy. M. Bucciolini, S. Pallotta and C. Talamonti are with INFN of Florence and Clinical Physiopathology Department, University of Florence, Italy. C. Civinini and M. Brianzi are with INFN of Florence, Italy. N.Randazzo is with INFN of Catania, Italy. D. Lo Presti and C. Stancampiano are with INFN of Catania Department of Physics - University of Catania Italy. G.A.P. Cirrone and G. Cuttone are with INFN- Laboratori Nazionali del Sud, Catania, Italy. F.Romano is with INFN- Laboratori Nazionali del Sud and Centro Studi e Ricerche e Museo Storico della Fisica "Enrico Fermi", Roma, Italy

978-1-4673-0120-6/11/$26.00 ©2011 IEEE

treatments, spatial resolution lower than Imm and density energy resolution lower than 1% [8]. The principal problem in the image quality is due to the Multiple Coulomb Scattering (MCS) of the proton in matter. The pCT prototype manufactured by our collaboration is based on the single proton tracking with the aim of using this information in the reconstruction technique. The apparatus consists of a silicon tracker and a Y AG:Ce calorimeter. Silicon detectors in front of and behind the object to be imaged provide spatial coordinates and direction of single protons and Y AG:Ce scintillating crystals measure its residual energy. From these measurements,

by a semi analytical

algorithm, the Most Likely Path (MLP) [9] is calculated and applied in the image reconstruction.

2 The system has a field of view of about 5x5cm , acquisition 5 times in the order of lOs (lOkHz, 10 events), it has been optimized to detect protons with energy in the range of 200250MeV (the energy needed to traverse a human head and trunk respectively) and to sustain potentially IMHz data acquisition rate. Apparatus details were described in previous

3505

works [2-7] where first results obtained with a 60MeV proton beam at LNS were discussed. Fig. 1 shows the picture of the LNS, more quantitative measurements were carried out to

Dhole Dsmm

determine the spatial resolution and the electronic density

DlOmm

apparatus located along the beam line. During the test beam at

resolution.

Radiographic projections of phantoms were 7 performed acquiring 10 events. Details of the results will be presented in the next section.

Fig. 3. Draft of the PMMA phantom irradiated with 60MeV proton beam which of the radiographic image has been analyzed. Different density zone are obtained with holes of different depth. The dimensions of holes and the reciprocal distances are shown.

Fig. 1. Picture of the pCT prototype located in the beam line at LNS. Four x-y planes of the tracker are visible. Tomographic images were acquired using the setup shown in Fig. 2. A cylindrical phantom with two empty zones of known sizes was located between the tracker planes P2 and P3. To provide multiple projection angles, the phantom was mounted on a rotating support. The object was irradiated at angles between 0 and 360 degrees in steps of 10 degrees and 6 for each angle 2x10 events were acquired.

Data have been acquired and analyzed off-line. In Fig. 4 the phantom radiographic image obtained is shown. The spatial coordinates of the single proton path, acquired in the plane nearest to the phantom, have been correlated with residual energy value measured by the calorimeter. The maximum energy

value

is

about

55MeV,

whereas

25MeV

is

the

minimum value: SRIM code [10] confirms that the residual energy values measured are comparable with energy loss values obtained in simulation phase. Four different density zones can be distinguished and the spatial resolution is able to reconstruct all holes. In order to understand the spatial resolution of pCT apparatus, an image profile is obtained

Rotating

starting of this image and shown in Fig. 5. The good spatial resolution is supported by these data: between 55MeV and 20MeV the sigma is equal to 3 strips (600flm).

Beam

P1

P2

Phantom P3

P4

Calorimeter

Fig. 2. Schematic display of the setup used for tomographic images acquisition. The rotating phantom is positioned between four planes (PI-P4) of the tracker, while the calorimeter is at the end of the particle path.

Ill. A.

RESULTS

Radiography images

PMMA

phantoms

have

been

placed

inside

the

pCT

apparatus along the proton beam line at LNS. The main aim was to acquire and to reconstruct radiographic images. In particular a PMMA phantom with different density zone has 7 been irradiated and 10 events have been acquired. The phantom is a cylinder with 2 cm of thickness and holes of different depth. Four different zones were manufactured: empty hole and holes 5, 10 and 15 mm thick. The Fig. 3 shows the draft of the phantom just described.

Fig. 4. Color map of PMMA phantom irradiated at LNS with 60MeV proton beam (see Fig. 3). With different colors are indicated the residual proton energy value (the values are expressed in MeV.

3506

energia

vs x,

IV.

135