Hi Mr Malekzadeh

You can use MATLAB or other tools to read these files. As I remember, GATE output for Dosimetry was float32 and exactly the same dimension as the original digital phantom you used for simulation(it depends on wich module you used for Dose report in GATE). Also, GATE output are based on voxel dose, so if you want to extract the organ dose you should consider this.

Hi Dr. Parach,

I really appreciate your help. It has a procedure as you mentioned. Fortunately, I could develope a MATLAB code to read the file.

The whole MATLAB code is here:

#For reading a dose file:

File=fopen('filepath', 'rb');

Data=fread(File, 'float');

A=reshape(Data, 128,128,243);

#For reading a phantom:

File2= fopen ( 'phanompath', 'rb');

Data2= fread(File2, 'int8');

A2= reshape (Data2, 128, 128, 243);

imlook3d( A2); #for phantom view

#For dose extraction: [cGy unit]

d= sum (Data(A2==2)); #num 2 is for example the ID of brain.

#For isodose

A3= 100* A2./max(A2(:));# normalization

a= A3(:,:,87); #cross section selection

contour (a); # isodose curve plotting

Also, I need to appreciate Dr. Asgaree.

Regards,