<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"> <head> <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> <meta name="Generator" content="Microsoft Word 14 (filtered medium)"> <style></style> </head> <body lang="FR" link="blue" vlink="purple"> <div class="WordSection1"> Hi Johannes,<o:p></o:p> <o:p> </o:p> The problem is that some data that are required to reconstruct an image are missing in the ECAT7 header. An example is the ring diameter. It is not possible to store the ring diameter in a ECAT7 file (neither the number of detectors per ring).<o:p></o:p> As STIR requires this information, it has a database for each known scanner model. What is printed at the end of your email is not the content of the ECAT7 header, but the content of STIR database:<o:p></o:p> <o:p> </o:p> WARNING: ECAT7 IO: Bin size from header.x_resolution (2.00714) does not agree with expected value 2.25<o:p></o:p> for scanner ECAT 962. Using bin size from header...<o:p></o:p> Info for file yourSinogram.S<o:p></o:p> ProjDataInfoCylindricalNoArcCorr := <o:p></o:p> Scanner parameters:= <o:p></o:p> Scanner type := ECAT 962<o:p></o:p> Number of rings := 32<o:p></o:p> Number of detectors per ring := 576<o:p></o:p> Inner ring diameter (cm) := 82.4<o:p></o:p> Average depth of interaction (cm) := 0.7<o:p></o:p> Distance between rings (cm) := 0.485<o:p></o:p> Default bin size (cm) := 0.200714<o:p></o:p> View offset (degrees) := 0<o:p></o:p> Maximum number of non-arc-corrected bins := 288<o:p></o:p> Default number of arc-corrected bins := 288<o:p></o:p> Number of blocks per bucket in transaxial direction := 3<o:p></o:p> Number of blocks per bucket in axial direction := 4<o:p></o:p> Number of crystals per block in axial direction := 8<o:p></o:p> Number of crystals per block in transaxial direction := 8<o:p></o:p> Number of detector layers := 1<o:p></o:p> Number of crystals per singles unit in axial direction := 8<o:p></o:p> Number of crystals per singles unit in transaxial direction := 24<o:p></o:p> end scanner parameters:=<o:p></o:p> <o:p> </o:p> As you set the scanner model to zero in your ECAT7 file (#/gate/output/ecat7/system 962 ), STIR does not recognize it. As a consequence, STIR assumes it is an HR+ (default scanner):<o:p></o:p> <o:p> </o:p> WARNING: ECAT7 IO: Couldn't determine the scanner <o:p></o:p> (Main_header.system_type=0), defaulting to 962.<o:p></o:p> This might give dramatic problems.<o:p></o:p> <o:p> </o:p> You should define a new scanner model ID in GATE, like /gate/output/ecat7/system 492, and create in STIR the corresponding database. I do not know how to create your own database in STIR.<o:p></o:p> <o:p> </o:p> Claude<o:p></o:p> <o:p> </o:p> <div> <div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm"> De : gate-users-bounces@lists.opengatecollaboration.org [mailto:gate-users-bounces@lists.opengatecollaboration.org] De la part de Johannes Anderl Envoy� : mercredi 2 octobre 2013 10:49 � : gate-users@lists.opengatecollaboration.org Objet : [Gate-users] ECAT7 output with user defined scanner<o:p></o:p> </div> </div> <o:p> </o:p> <div> Dear all, I have a problem concerning the ECAT7 Output. I would like to simulate a user defined scanner and not the E962 scanner type. The created yourSinogram.S file still contains the default values (ring diameter, detectors per ring,...) for the E962 scanner model. Is there a possbility to get this ECAT7 *.S output with the user defined scanner values? The code for the user defined scanner type can be seen below: /gate/world/daughters/name ecat /gate/world/daughters/insert cylinder /gate/ecat/setMaterial Air /gate/ecat/geometry/setRmax 44.1 cm /gate/ecat/geometry/setRmin 42.1 cm /gate/ecat/geometry/setHeight 21.7685 cm /gate/ecat/vis/forceWireframe # B L O C K /gate/ecat/daughters/name block /gate/ecat/daughters/insert box /gate/block/placement/setTranslation 431.0 0.0 0.0 mm /gate/block/geometry/setXLength 20.0 mm /gate/block/geometry/setYLength 52.17127 mm /gate/block/geometry/setZLength 52.17127 mm /gate/block/setMaterial Air /gate/block/vis/forceWireframe # C R Y S T A L /gate/block/daughters/name crystal /gate/block/daughters/insert box /gate/crystal/geometry/setXLength 20.0 mm /gate/crystal/geometry/setYLength 4 mm /gate/crystal/geometry/setZLength 4 mm /gate/crystal/setMaterial LSO /gate/crystal/vis/setColor yellow # R E P E A T C R Y S T A L /gate/crystal/repeaters/insert cubicArray /gate/crystal/cubicArray/setRepeatNumberX 1 /gate/crystal/cubicArray/setRepeatNumberY 13 /gate/crystal/cubicArray/setRepeatNumberZ 13 /gate/crystal/cubicArray/setRepeatVector 0. 4.0142725 4.0142725 mm # R E P E A T BLOCK (reproduziert die Box entlang eines Rings) /gate/block/repeaters/insert linear /gate/block/linear/setRepeatNumber 4 /gate/block/linear/setRepeatVector 0. 0. 55.17127 mm /gate/block/repeaters/insert ring /gate/block/ring/setRepeatNumber 48 The ecat7 code can be seen below: /gate/output/ecat7/enable /gate/output/ecat7/verbose 2 /gate/output/ecat7/setFileName yourSinogram /gate/output/ecat7/describe /gate/output/ecat7/mashing 2 /gate/output/ecat7/span 103 /gate/output/ecat7/maxringdiff 51 #/gate/output/ecat7/system 962 Looking at the values from the yoursinogram.S using list_projdata_info from STIR shows: WARNING: ECAT7 IO: Couldn't determine the scanner (Main_header.system_type=0), defaulting to 962. This might give dramatic problems. WARNING: ECAT7 IO: Bin size from header.x_resolution (2.00714) does not agree with expected value 2.25 for scanner ECAT 962. Using bin size from header... Info for file yourSinogram.S ProjDataInfoCylindricalNoArcCorr := Scanner parameters:= Scanner type := ECAT 962 Number of rings := 32 Number of detectors per ring := 576 Inner ring diameter (cm) := 82.4 Average depth of interaction (cm) := 0.7 Distance between rings (cm) := 0.485 Default bin size (cm) := 0.200714 View offset (degrees) := 0 Maximum number of non-arc-corrected bins := 288 Default number of arc-corrected bins := 288 Number of blocks per bucket in transaxial direction := 3 Number of blocks per bucket in axial direction := 4 Number of crystals per block in axial direction := 8 Number of crystals per block in transaxial direction := 8 Number of detector layers := 1 Number of crystals per singles unit in axial direction := 8 Number of crystals per singles unit in transaxial direction := 24 end scanner parameters:= Thanks in advance, Johannes <o:p></o:p> </div> </div> </body> </html>