Reconstruction of pet data with attenuation correction factor

Dear CASToR users and developers,

I’m trying to reconstruct GATE simulation data with attenuation correction factors in data file, but it ended with an error as shown below:

by the way, i’m trying to reconstruct emission image with attenuation corrected by transmission PET, so the attenuation correction factors is the ratio of the blank transmission sinogram to transmission with object one’s.

and here is my code for castor-recon:

castor-recon -df ./xxx.Cdh -opti MLEM -it 5:4 -proj joseph -fout ./xxx/xxx -th 0 -dim 170,170,71 -vox 1,1,1.98

Thank you very much in advance for your help.
seanlam

Dear Seanlam,

According to my understanding, the answer is in the first line of CASToR output : « Attenuation correction is included in the data file while it is not in the sensitivity computation ! ». Apparently, you are reconstructing list-mode data. To compute the sensitivity image, you need to provide CASToR either a normalization datafile with attenuation correction factor (see section 6.1 and Table 4 in the general documentation, option -norm) or an attenuation image (option -atn).

Regards,

Claude COMTAT
BioMaps
Service Hospitalier Frédéric Joliot
91401 Orsay, France
www.biomaps.universite-paris-saclay.fr

Claude COMTAT
BioMaps
Service Hospitalier Frédéric Joliot
91401 Orsay, France
www.biomaps.universite-paris-saclay.fr

Hi Claude,
First of all, thank you very much for your reply.
I’ve try your advice by providing an attenuation image and it did solved the error. But the reconstruction image seems to be the same img of no acf datafile reconstructed with the attenuation image. The acf contained in the datafile doesn’t seem to have any effect.
Could you further explain the role of attenuation correction factor in datafile? and I will try providing a normalization datafile with acf also, to see any difference.
thank you for your kindly help.
Best wishes
sean

Hi Claude
it worked, by reconstructing with the normalization datafile with attenuation factor.
thanks for you help!
Regards
sean

Great !

Claude COMTAT
BioMaps
Service Hospitalier Frédéric Joliot
91401 Orsay, France
www.biomaps.universite-paris-saclay.fr