[Castor-users] few questions about castor

[tmerlin]

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/60fccdd6/attachment.html>