[Gate-users] Voxelized source not filled

[J.J.M. Kouwenberg]

Dear all,

Im having trouble with getting GATE to recognize my voxelized source properly. I want to add a time activity curve to a simple cylindrical source. Since this is only possible for voxelized source, I generated a cylinder using STIR and replaced 'scaling factor (mm/pixel) [3]' by 'slice thickness (pixels)' and '!matrix size [2]' by 'number of slices' so GATE can read the file. The voxilized source is attached to a simple cylinder geometry (non-voxelized). I define the source via:
/gate/source/addSource Art voxel
/gate/source/Art/gps/particle e+
/gate/source/Art/gps/energytype Arb
/gate/source/Art/attachTo fem_Art
/gate/source/Art/reader/insert interfile
/gate/source/Art/interfileReader/translator/insert range
/gate/source/Art/interfileReader/SetTimeActivityTablesFrom ../../Phantoms/AIF/range_art.dat
/gate/source/Art/interfileReader/SetTimeSampling 1 s
/gate/source/Art/interfileReader/readFile ../../Phantoms/AIF/voxel_source.h33 
/gate/source/Art/setPosition -5. -5. -10. mm

GATE then outputs the following
activity range  [0 - 1e+07] is associated  to Time Activity Curve read from file ../../Phantoms/AIF/art.dat
GateSourceVoxelImageReader::ReadFile : fileName: ../../Phantoms/AIF/voxel_source.h33
 Header read from       '../../Phantoms/AIF/voxel_source.h33'
 Data file name         '../../Phantoms/AIF/voxel_source.i33'
 Nb of planes:           200
 Nb of pixels per plane: 100 100
 Pixel size:             0.1 0.1
 Slice thickness:        0.1
 Matrix size:            10 10
 Data type:              UNSIGNED INTEGER
 Data byte order:	 LITTLEENDIAN

Source ---------------> Art
  ID                  : 0
  type                : 
  startTime (s)       : 0
  time (s)            : 0
  forcedUnstable      : 1
  forcedLifetime (s)  : 5886.2
  verboseLevel        : 0
----------------------- 
  Voxel reader ---------->   number of voxels       : 0
  total activity (Bq)    : 0
  position  (mm)         : 0 0 0
  voxel size  (mm)    : 0.1 0.1 0.1
Source --------------> Art
  ID                 : 0
  type               : 
  activity (Bq)      : 0
  startTime (s)      : 0
  time (s)           : 0
  forcedUnstable     : 1
  forcedHalfLife (s) : 4080
  verboseLevel       : 0
  relative to vol    : fem_Art

Here you can see the problem. GATE's voxel reader somehow reads zero voxels, while the header file is correctly read and xmedcon shows that the voxel source is correctly saved.

Can someone with experience in voxelized source please help me out? I've been trying for quite a while but nothing seems to solve this issue.

With kind regards,

Jasper J.M. Kouwenberg
Medical Physics Graduate student
Nuclear Medicine



E-mail j.kouwenberg at erasmusmc.nl | Phone +316 2532 9809