[Gate-users] 3-D image reconstruction using STIR from GATE Results

[Kris Thielemans]

Hi Guifen

I do not see your ecat7 sinogram on our ftp. Can you put it there again?
Don't forget the "cd /in/kris", and add a "bin" to be sure.

See questions below


> 
> Hi,kris
> I simulate a point source at (0,0,0) and at (10,0,0) in GE 
> Discovery LS PET and in Ecat 962 systems using Gate. I want 
> to reconstruct the image using two methods.here are my steps:
> Method one:(GE Discovery LS PET )
> <1> simulate a point source (see CylindricalPET.mac ) <2> 
> [root at localhost Gate_Reconstruction]# root
>    root [0] .x C_GE_Sinogram_full.C

I don't know this one. Where does it come from? (sorry, my GATE knowledge is
close to 0).

> <3> /usr/local/gate/STIR/bin/display_projdata proj.hs ......

how did you get proj.hs?

> WARNING: Interfile warning: inconsistent 
> effective_central_bin_size_in_cm Value in header is 0.388 
> while I expect 0.221734 from the inner ring radius etc 
> Ignoring value in header  the results iare shown in 
> CylindricalPET_0.png and CylindricalPET_10.png.
>  
> Method two:(Ecat 962)
> <1> simulate a point source (see PET_Ecat.mac)
> <2> [root at localhost my_test]# 
> /usr/local/gate/STIR/bin/ifheaders_for_ecat7     MySinogramFile.S
> Attempt all data-sets (Y) or single data-set (N) [Y/N D:Y]:
> Data number ? [0,8 D:0]:
> WARNING: ECAT7 IO: Bin size from header.x_resolution (2.2471) 
> does not agree with expected value 2.25 for scanner ECAT 962. 
> Using bin size from header...
> 
> <3>[root at localhost my_test]# 
> /usr/local/gate/STIR/bin/display_projdata 
> MySinogramFile_S_f1g1d0b0.hs Which segment number do you want 
> to display[-2,2 D:0]:
> Display as SegmentByView (0) or BySinogram (1)?[0,1 D:0]: 1 
> Maximum in color scale (default is actual max)[0,10 D:5]:
> Displaying
> 
> The results are shown in Ecat_0.png and Ecat_10.png.
>  

I'm guessing the CylindricalPET_0 and 10.png files are the output for the
Discovery?

If so, they look just as wrong, and they have nothing to do with ECAT7
file-format. This starts to look like you have a local GATE problem.

> I'm very looking forwand to your help!
> Thanks in advance!
> guifen
> 
> 
> 
> 
> "Thielemans, Kris" <kris.thielemans at csc.mrc.ac.uk> 写道：
> 
> 	Hi Claude
> 	 
> 	I didn't know that iv_volumetool could display 
> sinograms. Thanks for the info!
> 	 
> 	So, I guess we either see a problem with 
> ifheaders_for_ecat7 or with the Guifen's local GATE 
> installation. Guifen, maybe you could drop your ECAT7 file at 
> ftp.hammersmithimanet.com <ftp://ftp.hammersmithimanet.com/>  
> (cd /in/kris) and let me know. I can have a quick look at it here.
> 	 
> 	Kris
> 
> ________________________________
> 
> 	From: gate-users-bounces at lists.healthgrid.org on behalf 
> of Claude Comtat
> 	Sent: Thu 12/04/2007 06:08
> 	To: Guifen Lin; gate-users at lists.healthgrid.org
> 	Subject: Re: 回复： Re: [Gate-users] 3-D image 
> reconstruction using STIR from GATE Results
> 	
> 	
> 	No, I did not use ifheaders_for_ecat7. I used an ecat7 dedicated
> 	software developped by the University of 
> Louvain-la-Neuve, Belgium,
> 	(iv_volumetool) to display directly the ecat7 sinogram 
> file generated by
> 	GATE. As far as GATE is concerned, there is no problem 
> with your macro
> 	and the ecat7 sinogram file is correct.
> 	
> 	Claude
> 	
> 	Guifen Lin wrote:
> 	
> 	> Hi,Claude
> 	>
> 	> Does you you do as my steps?I means using 
> "ifheaders_for_ecat7"
> 	> and than display the sinogram.
> 	> guifen
> 	>
> 	> Claude Comtat <claude.comtat at cea.fr> 写道：
> 	>
> 	>     Hi Guifen,
> 	>
> 	>     I just run your macro. First, your point source 
> is not at the center
> 	>     (x=0,y=0):
> 	>
> 	>     /gate/source/F18LineSource/gps/centre 10.0 0.0 0.0 cm
> 	>
> 	>     Anyhow, the resulting sinogram (MySinogramFile.S) 
> just looks fine.
> 	>     I do not see any particular problem. Maybe your 
> problem is related to
> 	>     the way you display the sinogram.
> 	>
> 	>     All the best,
> 	>
> 	>     Claude
> 	>
> 	
> 	
> 	_______________________________________________
> 	Gate-users mailing list
> 	Gate-users at lists.healthgrid.org
> 	http://lists.healthgrid.org/mailman/listinfo/gate-users
> 	
> 
> 
> ________________________________
> 
> 抢注雅虎免费邮箱-3.5G容量，20M附件！ <http://cn.mail.yahoo.com>  
>