<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=Windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:1307855683;
mso-list-type:hybrid;
mso-list-template-ids:1376522624 67698703 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level2
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level3
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level4
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level5
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level6
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level7
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level8
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level9
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l1
{mso-list-id:1470249524;
mso-list-type:hybrid;
mso-list-template-ids:1395013084 1056983852 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l1:level1
{mso-level-tab-stop:none;
mso-level-number-position:left;
margin-left:54.0pt;
text-indent:-18.0pt;}
@list l1:level2
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
margin-left:90.0pt;
text-indent:-18.0pt;}
@list l1:level3
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
margin-left:126.0pt;
text-indent:-9.0pt;}
@list l1:level4
{mso-level-tab-stop:none;
mso-level-number-position:left;
margin-left:162.0pt;
text-indent:-18.0pt;}
@list l1:level5
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
margin-left:198.0pt;
text-indent:-18.0pt;}
@list l1:level6
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
margin-left:234.0pt;
text-indent:-9.0pt;}
@list l1:level7
{mso-level-tab-stop:none;
mso-level-number-position:left;
margin-left:270.0pt;
text-indent:-18.0pt;}
@list l1:level8
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
margin-left:306.0pt;
text-indent:-18.0pt;}
@list l1:level9
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
margin-left:342.0pt;
text-indent:-9.0pt;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-CA" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi Khalid,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I cannot speak to your specific error, but did you make sure to update the dataDir directory variable in the ICRPtoSTL.py file? Also, I’m not sure if you are referring to Andrew’s README.txt file for organ names, i.e.
‘liver’, but with the ICRP 110 data that I downloaded <a href="https://www.icrp.org/publication.asp?id=icrp%20publication%20110">
here</a>, the organ names are actually tagged with their medium number, i.e. ‘liver_30’. You can verify this by uncommenting the line “data.listOrgans()” in the ICRPtoSTL.py file.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">The scripts appear to be written for python2 and I was unable to execute the scripts at all due to an import error for VTK. I converted all scripts to python3 with “2to3 -w <name>.py” and implemented the following changes:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<ol style="margin-top:0cm" start="1" type="1">
<li class="MsoListParagraph" style="margin-left:0cm;mso-list:l0 level1 lfo2"><span lang="EN-US">In readICRP.py<br>
-On line 88, I replaced “words = re.split(“ *”, line)” with “words = re.split(r’\W+’, line)” to correct a ValueError which was not properly splitting the contents of line; and<o:p></o:p></span></li></ol>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<ol style="margin-top:0cm" start="2" type="1">
<li class="MsoListParagraph" style="margin-left:0cm;mso-list:l0 level1 lfo2"><span lang="EN-US">In dicomToSTL.py<o:p></o:p></span></li></ol>
<p class="MsoNormal" style="margin-left:36.0pt"><span lang="EN-US">-On line 71, the function expects a 3x1 vector for a kernel radius but Andrew’s implementation only passes a single integer. I admit, I do not know what this radius refers to or how it is applied,
maybe someone can shed some light on this, but replacing “organMask=sitk.BinaryClosingByReconstruction(organMask,30)” with “organMask=sitk.BinaryClosingByReconstruction(organMask)” or “organMask=sitk.BinaryClosingByReconstruction(organMask,[30,30,30])” mitigates
the issue and I notice no apparent difference in the output. <o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I have attached the modified scripts which allows me to execute a conversion, for example, with “python ICRPtoSTL.py female teeth_1” which produces STL files that can be imported into GATE.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Cheers,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Matthew<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;color:black">From:
</span></b><span style="font-size:12.0pt;color:black">Gate-users <gate-users-bounces@lists.opengatecollaboration.org> on behalf of Khalid Hussain <khalidhussain1134@gmail.com><br>
<b>Date: </b>Friday, April 30, 2021 at 14:36<br>
<b>To: </b>Said Kaddouch <saidkaddouch@gmail.com>, gate-users <Gate-users@lists.opengatecollaboration.org>, Andrew Robinson <andrew.robinson@npl.co.uk><br>
<b>Subject: </b>Re: [Gate-users] ICRP 110 AF Phantom<o:p></o:p></span></p>
</div>
<div align="center">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;border-radius:5px;overflow:hidden">
<tbody>
<tr>
<td style="padding:7.5pt 0cm 7.5pt 0cm">
<div align="center">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;background:#707372">
<tbody>
<tr>
<td valign="top" style="border:none;border-top:solid #FBE122 6.0pt;padding:3.0pt 6.0pt 3.0pt 6.0pt">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0">
<tbody>
<tr>
<td style="padding:0cm 0cm 0cm 0cm">
<p class="MsoNormal" style="line-height:12.0pt;mso-line-height-rule:exactly"><b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:white">CAUTION:</span></b><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:white"> The Sender
of this email is not from within Dalhousie.<o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:white">Can anyone advise me, how can I read ICRP110 AF ph</span>antom data to generate an image to use in GATE simulation?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">I got this error while running Andrew <a href="https://github.com/a-p-robinson/STLPhantoms">
code</a>.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><img border="0" width="497" height="66" style="width:5.177in;height:.6875in" id="_x0000_i1026" src="cid:image001.png@01D73DF4.46A545F0"><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Sincerely,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Khalid<o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Regards<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"><b><span style="color:#666666">KHALID HUSSAIN </span></b><span style="color:#666666">|</span> <b><span style="color:#BF9000">PhD Scholar</span></b><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><b><span style="font-family:"Times New Roman",serif;color:#073763">Universiti Putra Malaysia (UPM), Malaysia</span></b><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Times New Roman",serif">Email:<a href="mailto:hussain.khalid@student.upm.edu.my" target="_blank">hussain.khalid@student.upm.edu.my</a></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Times New Roman",serif">Cell: <u><span style="color:blue">+60-11-6455-2354</span></u></span><o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Fri, Apr 30, 2021 at 10:31 PM Khalid Hussain <<a href="mailto:khalidhussain1134@gmail.com" target="_blank">khalidhussain1134@gmail.com</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">Dear Said, <o:p></o:p></p>
<div>
<p class="MsoNormal">Thank you for your email.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">I think this xmedcon is used once we have voxel image and in my case, I don't have voxel phantom, I only have data downloaded from
<a href="https://www.icrp.org/publication.asp?id=ICRP%20Publication%20110" target="_blank">
icrp.org</a> which is in the form of .dat files and needs to be read by the program such as python program advised by Andrew but it gives some parsing/namespace error, which I am requesting to Andrew or anyone else who has used this data before.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Can anyone advise me, how can I read ICRP110 AF phantom data to generate an image to use in GATE simulation?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><img border="0" width="497" height="66" style="width:5.177in;height:.6875in" id="_x0000_i1025" src="cid:image001.png@01D73DF4.46A545F0"><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">Sincerely,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Khalid<o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Fri, Apr 30, 2021 at 3:26 AM Said Kaddouch <<a href="mailto:saidkaddouch@gmail.com" target="_blank">saidkaddouch@gmail.com</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">Dear khalid <o:p></o:p></p>
<div>
<p class="MsoNormal">To read the voxelised phantom you can use xmedcon to read the phantom and translate it to intefile or analyse format. After you can use Matlab to read your phantom format analyse. You can find matlab analyse reader in the internet.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Hop this help. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Best regards. <o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Mon, Apr 26, 2021, 3:27 PM Khalid Hussain <<a href="mailto:khalidhussain1134@gmail.com" target="_blank">khalidhussain1134@gmail.com</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">Dear Gate Users, <o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Is there anyone who has a python or matlab code to read the ICRP110 AF phantom?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I want to request you to please share the code to read this phantom. I am trying to translate the fortran code given in readme file but I did find the detail of some functions used in fortan code provided in readme file. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thank you for your cooperation and support.<br clear="all">
<o:p></o:p></p>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Regards<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"><b><span style="color:#666666">KHALID HUSSAIN </span></b><span style="color:#666666">|</span> <b><span style="color:#BF9000">PhD Scholar</span></b><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><b><span style="font-family:"Times New Roman",serif;color:#073763">Universiti Putra Malaysia (UPM), Malaysia</span></b><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Times New Roman",serif">Email:<a href="mailto:hussain.khalid@student.upm.edu.my" target="_blank">hussain.khalid@student.upm.edu.my</a></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Times New Roman",serif">Cell: <u><span style="color:blue">+60-11-6455-2354</span></u></span><o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal">_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org" target="_blank">Gate-users@lists.opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" target="_blank">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a><o:p></o:p></p>
</blockquote>
</div>
</blockquote>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</body>
</html>