[Castor-users] Potential bug in castor-PetScannerLutEx

[Thibaut Merlin]

  Hi Peter,
Did you center your ring of crystal at z=0 ? Because really the axial
position just depends on where your ring is located.

Best,
Thibaut

Peter Neilson <neilson at txcorp.com> a écrit :

> Hello Thibaut,      
>    Thanks for your assistance, I’ve been able to get the script up and
> reconstructing images. However I’ve noticed something rather odd when I
> view them in amide. 
>    
>   In this case I’m still working with a single ring of axial modules
and
> rsectors to keep things simple, with point sources at z=0. I would
> expect the produced image to be centered in amide when I view it, and do
> seem to have the correct X and Y coordinates, however I find the sources
> at z = -11.63 mm. Is there some form of metadata or option to set these
> coordinates?
>    
>   Thanks,
>   -Peter Neilson
>
>
>> On May 19, 2020, at 1:38 PM, tmerlin <Thibaut.Merlin at univ-brest.fr>
>> wrote:
>>
>>
>>        Hi Peter,
>>
>>        Thank you for noticing this.
>>
>>        This code is just an example on how to make a LUT following
>> CASToR ring by ring detector indexation.
>>
>>        In order to adapt it to GATE indexing, you must change the last
>> loop which attribute the crystal index to the /cryID /variable/./ For
>> your example system, I believe it should be something like
>> /int cryID = rs*nb_axial_crystal*nb_trans_crystal + l /
>>
>>        ..to index all the crystals in the first rsector, then the next
>> one, etc. You should check that the ordering and position of the first
>> rsector/crystal are consistent between CASToR and GATE indexing.
>>
>>        Hope this helps,
>> Thibaut
>>        On 18/05/2020 22:01, Peter Neilson wrote:
>>
>>> Hello All,                   
>>>            I’m working on adapting castor-PetScannerLutEx to match
a
>>> Gate geometry that is incompatible with castor-GATEMacToGeo. However
>>> on my first run of castor-PetScannerLutEx (no changes) I get a seg
>>> fault. 
>>>          I found there was a bug at line 342, the oMatrix needs to be
>>> initialized with nb_resctors +1 oMatrix *****crystal_center = new
>>> oMatrix ****[nb_rsectors+1];
>>>           
>>>          After that seemed to get the program to compile alright, but
>>> I’m having some difficulty in getting sensible results out of it with
>>> my geometry. For now I’m using GATE with a cylindricalPET system with
>>> only a single ring of 28 rsectors, no layers and 16 x 16 crystal
>>> arrangement, nothing that castor-GateRootToCastor couldn’t handle. I
>>> end up with a non-sequential castor indexing, that I think is the
>>> cause of my problems. 
>>>           
>>>          Are there any further tutorials on
>>> using castor-PetScannerLutEx ? 
>>>           
>>>          Thanks,
>>>          -Peter Neilson
>>>
>>> _______________________________________________ Castor-users mailing
>>> list
>>>
Castor-users at lists.castor-project.orghttp://lists.castor-project.org/listinfo/castor-users
>>
>> --  Thibaut MERLIN -- PhD  Docteur en Imagerie Médicale au Laboratoire
>> de Traitement de l'Information Medicale (LaTIM - INSERM UMR 1101)
>> Institut Brestois de recherche en Bio-Santé (IBRBS) 12 Avenue Foch,
>> 29200 Brest, FRANCE Tel: 06.75.12.24.90
>> _______________________________________________
>> Castor-users mailing list
>> Castor-users at lists.castor-project.org
>> http://lists.castor-project.org/listinfo/castor-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20200528/f9e75772/attachment.html>