[Castor-users] TR: SPECT reading Interfile

[COMTAT Claude]

Dear Michael,

 

radial : I

axial : J

azimuthal : K

 

You should not store your projection data as it is usually done :

 

  for (k=0 ;k<K ;k++) for (j=0 ;j<J ;j++) for (i=0 ;i<I ;i++) fwrite(Mat(i,j,k));

 

and then perform an OSEM reconstruction. In CASToR, there is no assumption at all about the organization of the data. The code will read sequentially the data : if there are N subsets, first subset will contain Mat(i=0,j=0,k=0), Mat(i=N,j=0,k=0), Mat (i=2*N,j=0,k=0), … and the resulting subsets will not be well balanced.

 

You can rather write your matrix as : 

 

  for (j=0 ;j<J ;j++) for (i=0 ;i<I ;i++) for (k=0 ;k<K ;k++)  fwrite(Mat(i,j,k)) ;

 

Then, first subset will contain Mat(i=0,j=0,k=0), Mat(i=0,j=0,k=N), Mat(i=0,j=0,k=2*N), … and the resulting subsets will be better balanced.

 

I did not get your figure 9.

 

Best,

 

Claude

De : Michael Ljungberg [mailto:michael.ljungberg at med.lu.se] 
Envoyé : mardi 25 avril 2017 08:18
À : COMTAT Claude
Objet : Re: [Castor-users] SPECT reading Interfile

 

Hi Claude

 

Sorry but I am still puzzle about this. You mentioning I should write the histogram events with the inner loop over the radial or axial coordinate. Radial – what does this mean? Next sentences mention the azimuthal coordinate. What does this mean relative to the previous question? 

 

I have a float matrix Mat(I,J,K) where I,J is the projections bins for a specific SPECT camera position (projection ID ) K. The index I goes from left to right and J goes axially from the patient head and down.  

 

>From Figure 9 it looks like I need to flip the axial orientation because you start from the legs and goes to head (lower left corner to upper left corner for the axial rows. 

Any clarification about how to write the dafa file based on the above would be very helpful because the osem if of course of interest to be able to use.

Hälsningar / Best regards

Michael

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20170425/a9b2e7cd/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 5587 bytes
Desc: not available
URL: <http://lists.castor-project.org/pipermail/castor-users/attachments/20170425/a9b2e7cd/attachment.bin>