<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<title></title>
</head>
<body>
<div name="messageBodySection">
<div dir="auto">Hi Daniela,
<div dir="auto">yes, dicomfiles contains the filenames which match the pattern '*.dcm’, so the dicom files. Up to that point, this has nothing to do with gatetools, by the way. Just a standard way to read in filenames in python. A priori, the filenames are not sorted, but that should not be a problem for you. The read_dicom function in gatetools looks for the <span style="background-color: var(--backgroundColor);">SliceLocation dicom attribute and, if present, sorts the slices accordingly. </span></div>
<div dir="auto"><span style="background-color: var(--backgroundColor);"><br /></span></div>
<div dir="auto">On to your error: The error looks like the file is not where you think it is. My code works only if you execute it in the folder where the dicom files are. If, f<span style="color: var(--textColor); background-color: var(--backgroundColor);">or example, you write</span></div>
<div dir="auto"><span style="color: var(--textColor); background-color: var(--backgroundColor);"><br /></span></div>
<div dir="auto"><span style="color: var(--textColor); background-color: var(--backgroundColor);">dicomfiles = fnmatch.filter(os.listdir(‘/Users/nkrah/dicoms‘), ‘*dcm’)</span></div>
<div dir="auto"><span style="color: var(--textColor); background-color: var(--backgroundColor);"><br /></span></div>
<div dir="auto"><span style="color: var(--textColor); background-color: var(--backgroundColor);">the list dicomfiles will contain the names in the directory </span><span style="caret-color: rgb(39, 39, 40); color: var(--textColor); background-color: var(--backgroundColor);">/Users/nkrah/dicoms but not their full path. If I run the script from /Users/nkrah/scripts and feed the dicomfiles list to gt.read_dicom(dicomfiles), python will try to find the files in the current directory, i.e., where the script is executed. </span></div>
<div dir="auto"><span style="caret-color: rgb(39, 39, 40); color: var(--textColor); background-color: var(--backgroundColor);"><br /></span></div>
<div dir="auto"><span style="caret-color: rgb(39, 39, 40); color: var(--textColor); background-color: var(--backgroundColor);">That might explain your problem. </span><span style="caret-color: rgb(39, 39, 40); color: var(--textColor); background-color: var(--backgroundColor);">Can you share the code you are using or are you using the exact same lines I have provided?</span></div>
<div dir="auto"><span style="caret-color: rgb(39, 39, 40); color: var(--textColor); background-color: var(--backgroundColor);"><br /></span></div>
<div dir="auto">Try using absolute paths: </div>
<div dir="auto"><br /></div>
<div dir="auto">dicomfiles_names = fnmatch.filter(os.listdir('./'), '*.dcm’)</div>
<div dir="auto">dicomfiles_absolute = [os.path.abspath(df) for df in dicomfiles_names] # this compact syntax is called ‘list comprehension’ in python. </div>
<div dir="auto"><br /></div>
<div dir="auto">… and feed this as input to the read_dicom function. </div>
<div dir="auto"><br /></div>
<div dir="auto">I know these things are a bit frustrating at the beginning, but once you get a hang of how to deal with files and so on in python, it’s quite straight forward. </div>
</div>
</div>
<div name="messageSignatureSection"><br />
<div class="matchFont">Cheers,<br />
Nils</div>
</div>
<div name="messageReplySection">On 30 Jan 2020 at 17:40 +0100, Botnariuc, Daniela <daniela.botnariuc.19@ucl.ac.uk>, wrote:<br />
<blockquote type="cite" class="spark_quote" style="margin: 5px 5px; padding-left: 10px; border-left: thin solid #1abc9c;">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Dear Nils,</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Thank you very much for the quick reply, I really appreciate it.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">I tried to run the code you suggested but I still get an error.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">According to the code, the ‘dicomfiles’ variable is an array which lists all the files in the folder: [‘dicom1.dcm’, ‘dicom2.dcm’, ‘dicom3.dcm’, ….] - it actually lists the files in a random order.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">When I run the code, I get the following error:</span></p>
<p class="MsoNormal" style="text-indent:36.0pt"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">no such file or directory: ‘dicom 17.dcm' (being ‘dicom17.dcm’ the first file appearing in my list)</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">I appreciate any help with this issue.</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> </span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Thank you</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US">Daniela</span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D;mso-fareast-language:EN-US"> </span></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif" xml:lang="EN-US">From:</span></b> <span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri",sans-serif" xml:lang="EN-US">Nils Krah <nils.krah@creatis.insa-lyon.fr><br />
<b>Sent:</b> 30 January 2020 14:10<br />
<b>To:</b> gate-users@lists.opengatecollaboration.org; Botnariuc, Daniela <daniela.botnariuc.19@ucl.ac.uk><br />
<b>Subject:</b> Re: [Gate-users] Run GATETOOLS from Python script</span></p>
</div>
</div>
<p class="MsoNormal"> </p>
<div name="messageBodySection">
<div>
<p class="MsoNormal">Hi Daniela, </p>
<div>
<p class="MsoNormal">within a python script, you need to call the python function image_convert. It takes as input an ITK image. </p>
</div>
<div>
<p class="MsoNormal">So, say you have a series of dicom files representing a 3D image and you want to convert it to float, you could do: </p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">import os</p>
</div>
<div>
<p class="MsoNormal">import fnmatch</p>
</div>
<div>
<p class="MsoNormal">import itk</p>
</div>
<div>
<p class="MsoNormal">import gatetools as gt</p>
</div>
<div>
<p class="MsoNormal"><br />
<br /></p>
</div>
<div>
<p class="MsoNormal"># get file in current folder and filter out those ending on .dcm</p>
</div>
<div>
<p class="MsoNormal">dicomfiles = fnmatch.filter(os.listdir('./'), '*.dcm’)</p>
</div>
<div>
<p class="MsoNormal"># read them in and get them as ITK image</p>
</div>
<div>
<p class="MsoNormal">itk_image = gt.read_dicom(dicomfiles)</p>
</div>
<div>
<p class="MsoNormal"># pass the ITK image to the convert function and specify to which type you want to convert</p>
</div>
<div>
<p class="MsoNormal">converted_image = gt.image_convert(itk_image, pixeltype=‘float’)</p>
</div>
<div>
<p class="MsoNormal"># you get again an ITK image and yu can continue manipulating it as you wish.</p>
</div>
<div>
<p class="MsoNormal"># You may dump it with:</p>
</div>
<div>
<p class="MsoNormal">itk.imwrite(converted_image, ‘converted_image.mhd’)</p>
</div>
<div>
<p class="MsoNormal"><br />
<br /></p>
</div>
<div>
<p class="MsoNormal">(I haven’t tested this bit of code, but it should work.)</p>
</div>
<div>
<p class="MsoNormal"><br />
<br /></p>
</div>
<div>
<p class="MsoNormal">Technical remark for python newcomers: </p>
</div>
<div>
<p class="MsoNormal">When writing gt. preceding the function name, you are referring to the namespace of the imported module. The module’s name gatetools has been aliased to gt for brevity. That is not mandatory, but often practical. </p>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<p class="MsoNormal">Hope that helps. </p>
</div>
</div>
</div>
<div name="messageSignatureSection">
<p class="MsoNormal"> </p>
<div>
<p class="MsoNormal">Best,<br />
Nils</p>
</div>
</div>
<div name="messageReplySection">
<p class="MsoNormal">On 30 Jan 2020 at 12:56 +0100, Botnariuc, Daniela <<a href="mailto:daniela.botnariuc.19@ucl.ac.uk">daniela.botnariuc.19@ucl.ac.uk</a>>, wrote:<br />
<br /></p>
<blockquote style="margin: 5px 5px; padding-left: 10px; border-left: thin solid #e67e22;" class="spark_quote">
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Dear GATE users,</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I have been using the gatetools (thank you for this great contribution!) from the command line but I would like to run them from a python script. I am not as familiar with python as I should be to understand the code in detail. I tried to transform the bash scripts that run the tools from the command line to python files but I am having problems with the type of input/output the functions require.</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">For example, for the gt_image_convert, how should I write my input (path to the dicom folder) and output in python to run the function?</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">I thought it would be something like: (but I am missing something)</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">gt_image_convert('/home/gate/Documents/dicom_folder','output.mhd','float')</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Thank you. I look forward to hearing from you.</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Kind regards,</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Daniela</p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> </p>
</div>
<p class="MsoNormal">_______________________________________________<br />
Gate-users mailing list<br />
<a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.opengatecollaboration.org</a><br />
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users">http://lists.opengatecollaboration.org/mailman/listinfo/gate-users</a></p>
</blockquote>
</div>
</div>
</blockquote>
</div>
</body>
</html>