<html dir="ltr"><head><meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1" /><style id="owaParaStyle" type="text/css">
<!--
p
 {margin-top:0;
 margin-bottom:0}
-->
P {margin-top:0;margin-bottom:0;}</style></head><body ocsi="0" fpstyle="1">Hi David,<br>
<br>
have you tried the metafile (.mhd, .hdr) file format? It usually works better.<br>
<br>
Seasons greetings,<br>
Hermann<br><br><div class="gmail_quote">Am 23. Dezember 2014 17:59:13 MEZ, schrieb David Larsson <davlars@kth.se>:<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">






<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Dear GATE-users,<br />
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt"><br />
I am trying to implement a voxelized phantom into my GATE-simulations for the first time, but seem to encounter problems in the Interfile-read of my script.<br />
<br />
Using a simple imageconverter ((x)medcon) I've created a test-interfile (.h33 + .i33) with a 10x10x10 matrix with voxel values ranging from 0 to 10. The conversion seems to work fine, and when I open to generated interfile-set in a random image viewer (ImageJ/MATLAB)
 I can see that the voxel values are preserved (ranging from 0 - 10).<br />
<br />
However, when I run my script I notice that the loader does not import the voxels as ranging from 0-10, instead I can see (through the error messages) that all voxels are given one value (in this specific case 9280). The original values and the corresponding
 input structure is therefore lost.<br />
<br />
I've tried to compare my generated code to the (working) example with the hof_brain but cannot seem to understand the problem. I've provided the essential parts of the mac-files and the interfile-header of my testfile below. Can anyone identify what the problem
 might be?<br />
<br />
----------------------------main.mac----------------------------------------<br />
<br />
# WORLD <br />
/gate/world/geometry/setXLength       400. cm<br />
/gate/world/geometry/setYLength       400. cm<br />
/gate/world/geometry/setZLength       400. cm<br />
<br />
# VOXEL PHANTOM<br />
<br />
/gate/world/daughters/name                                                       teststack<br />
/gate/world/daughters/insert                                                      fictitiousVoxelMap<br />
/gate/teststack/geometry/insertReader                                      interfile<br />
/gate/teststack/interfileReader/insertTranslator                       range<br />
/gate/teststack/interfileReader/rangeTranslator/readTable     attenuation_range.dat   
<br />
/gate/teststack/interfileReader/rangeTranslator/describe        1<br />
/gate/teststack/interfileReader/readFile                                     testintf.h33<br />
/gate/teststack/placement/setTranslation                                  0. 0. 0. mm<br />
/gate/teststack/placement/setRotationAxis                                1 0 0<br />
/gate/teststack/placement/setRotationAngle                             0 deg<br />
/gate/teststack/attachVoxelPhantomSD<br />
<br />
-------------------attenuation_range.dat--------------------------<br />
<br />
5<br />
0 1     Air        false    0.0    0.0    0.0    0.0<br />
1 2     Water        true    1.0    0.0    0.0    0.0<br />
2 3     Water        true    1.0    0.0    0.0    0.0<br />
4 9280     Water        true    1.0    0.0    0.0    0.2<br />
9280 9280     Air    true    1.0    0.0    0.0    0.2     #This changed for visualisation of incorrect values<br />
<br />
------------------------- testintf.h33--------------------------<br />
<br />
!INTERFILE := <br />
!imaging modality := nucmed <br />
!originating system := (X)MedCon <br />
!version of keys := 3.3 <br />
date of keys := 1996:09:24 <br />
conversion program := (X)MedCon <br />
program author := Erik Nolf <br />
program version := 0.13.0 <br />
program date := 2013:06:23 <br />
; <br />
!GENERAL DATA := <br />
original institution := NucMed <br />
!data offset in bytes := 0 <br />
!name of data file := m000-testnew.i33 <br />
patient name := Unknown <br />
!patient ID := Unknown <br />
patient dob := 0000:00:00 <br />
patient sex := Unknown <br />
!study ID := Unknown <br />
exam type := Unknown <br />
data compression := none <br />
data encode := none <br />
organ := Unknown <br />
isotope := Unknown <br />
dose := 0 <br />
NUD/Patient Weight [kg] := 0.00 <br />
NUD/imaging modality := NM <br />
NUD/activity := 0 <br />
NUD/activity start time := 00:00:00 <br />
NUD/isotope half life [hours] := 0.000000 <br />
; <br />
!GENERAL IMAGE DATA := <br />
!type of data := Tomographic <br />
!total number of images := 10 <br />
study date := 0000:00:00 <br />
study time := 00:00:00 <br />
imagedata byte order := LITTLEENDIAN <br />
process label :=  <br />
; <br />
number of energy windows := 1 <br />
; <br />
energy window [1] := <br />
energy window lower level [1] := <br />
energy window upper level [1] := <br />
flood corrected := N <br />
decay corrected := N <br />
; <br />
!SPECT STUDY (general) := <br />
number of detector heads := 1 <br />
; <br />
!number of images/energy window := 10 <br />
!process status := Reconstructed <br />
!matrix size [1] := 10 <br />
!matrix size [2] := 10 <br />
!number format := long float <br />
!number of bytes per pixel := 8 <br />
scaling factor (mm/pixel) [1] := +10.000000e+00 <br />
scaling factor (mm/pixel) [2] := +10.000000e+00 <br />
!number of projections := 10 <br />
!extent of rotation :=  <br />
!time per projection (sec) := 0 <br />
study duration (sec) := 0 <br />
!maximum pixel count := +1.000000e+01 <br />
patient orientation := Unknown <br />
patient rotation := Unknown <br />
; <br />
!SPECT STUDY (reconstructed data) := <br />
method of reconstruction := Unknown <br />
!number of slices := 10 <br />
number of reference frame := 0 <br />
slice orientation := unknown <br />
slice thickness (pixels) := +10.000000e+00 <br />
centre-centre slice separation (pixels) := +1.000000e+00 <br />
filter name := Unknown <br />
filter parameters := Cutoff <br />
method of attenuation correction := measured <br />
scatter corrected := N <br />
oblique reconstruction := N <br />
!END OF INTERFILE :=<br />
<br />
----------------------------------------------------------------------------<br />
<br />
Best regards,<br />
<br />
David  <br />
</div>
</div>
</div>
</div>


<p style="margin-top: 2.5em; margin-bottom: 1em; border-bottom: 1px solid #000"></p><pre class="k9mail"><hr /><br />Gate-users mailing list<br />Gate-users@lists.opengatecollaboration.org<br /><a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></pre></blockquote></div></body></html>