<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>