[Castor-users] few questions about castor

[Emanuele Ripiccini]

Dear Thibaud,

Many thanks. Our scanner has about 1.5M detectors, so it is normal that it takes all that time.

Cheers,

Emanuele

****************************************************************
Dr Emanuele Ripiccini
Assistant Post-doctorant

Université de Genève
Département de physique nucléaire et corpusculaire
Quai Ernest-Ansrmet 24, 1205 Genève<x-apple-data-detectors://1/1>

Phone nr : +41 223796151<tel:+41%20223796151>
Pavillon des isotopes
Office nr : 124
****************************************************************


Le 20 juin 2018 à 17:37, tmerlin <Thibaut.Merlin at univ-brest.fr<mailto:Thibaut.Merlin at univ-brest.fr>> a écrit :


Dear Emanuele,

You can enable the verbose to see the progression of the sensitivity image generation (add the following option : -vb 2)

This computation performs a dual loop on all the detector elements to compute all possible lines of response, depending on the potential angle restrictions in the geometry file and datafile. So the duration will mostly depend on the total number of detectors and the image dimensions.

Once the sensitivity image is computed, it could be re-used for other reconstructions (provided these reconstructions use the same dimensions, projectors and image-based operations such as psf, etc..) using the option -sens sensitivity_file.hdr.

I will forward this answer to our user mailing-list as this is a question that other users could potentially have.

Best,
Thibaut


On 20/06/2018 14:03, Emanuele Ripiccini wrote:
Dear Thibaut,

many thanks again for your precious answers. I successfully created the geometry LUT, the geometry header, the binary data file and the header data file. I executed castor using my geometry and my data with this command

./castor-recon -df /Users/emanuele/Desktop/example_ttpet/TTPETdata.cdh -opti MLEM -it 5:1 -proj joseph -conv gaussian,0.7,0.7,1.::psf -dim 200,200,200 -vox 0.25,0.25,0.25 -dout image_test


I observed this output

sOutputManager::InitOutputDirectory() -> Output path is 'image_test'
castor-recon() -> Launch reconstruction from CASToR version 2.0.1.
sScannerManager::Initialize() -> From scanner TTPET
oProjectorManager::Initialize() -> Initialize projectors and projection lines
oOptimizerManager::Initialize() -> Initialize optimizer and penalty
oImageConvolverManager::Initialize() -> Initialize image convolvers
oSensitivityGenerator::Launch() -> Start the sensitivity computation

but after this it has been 3 hours and half more or less. You told me that it would take time given our geometry, but I don’t know if it is normal or not and I should see some output to follow the status.

best regards,

Emanuele


****************************************************************
Dr Emanuele Ripiccini
Assistant Post-doctorant

Université de Genève
Département de physique nucléaire et corpusculaire
Quai Ernest-Ansrmet 24, 1205 Genève

Phone nr : +41 223796151
Pavillon des isotopes
Office nr : 124
****************************************************************



--
Thibaut MERLIN -- PhD

Docteur en Imagerie Médicale au Laboratoire de Traitement de l'Information Medicale (LaTIM - INSERM UMR 1101)
Bâtiment 1
CHRU Morvan - 2, Av. Foch
29609 Brest CEDEX - FRANCE
Tel: 06.75.12.24.90
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20180620/aafc1f79/attachment-0001.html>