<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=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:Verdana;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
        {mso-style-priority:99;
        mso-style-link:"Balloon Text Char";
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:8.0pt;
        font-family:"Tahoma","sans-serif";}
span.hoenzb
        {mso-style-name:hoenzb;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        font-family:"Calibri","sans-serif";
        color:#1F497D;}
span.BalloonTextChar
        {mso-style-name:"Balloon Text Char";
        mso-style-priority:99;
        mso-style-link:"Balloon Text";
        font-family:"Tahoma","sans-serif";
        mso-fareast-language:EN-GB;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri","sans-serif";
        mso-fareast-language:EN-US;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:70.85pt 70.85pt 2.0cm 70.85pt;}
div.WordSection1
        {page:WordSection1;}
--></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-GB" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">It doesn’t work.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Filter.GetOuput() already has the new origin (I can print it and it’s right). Something wrong is in the writing process.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">By the way, I noticed this problem occurs for change of the origin and of rescaling intensities. It does not occur for changes in size and spacing.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Fabio</span><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span lang="EN-US" style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> Dženan Zukić [mailto:dzenanz@gmail.com]
<br>
<b>Sent:</b> Dienstag, 20. September 2016 14:46<br>
<b>To:</b> D'Isidoro Fabio<br>
<b>Cc:</b> Timothee Evain; insight-users@itk.org<br>
<b>Subject:</b> Re: [ITK-users] [ITK] itk Python: origin of DICOM slice does not change<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif"">Hi Fabio,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif""><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif"">did you try setting the origin directly on the result image:<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">filter->Update();</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">outImage=filter->GetOutput();</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">outImage->SetOrigin(new_origin);</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">writer_output->SetInput(outImage);</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif""><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif"">Regards,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif"">Dženan<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana","sans-serif""><o:p> </o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">On Tue, Sep 20, 2016 at 8:34 AM, D'Isidoro Fabio <<a href="mailto:fisidoro@ethz.ch" target="_blank">fisidoro@ethz.ch</a>> wrote:<o:p></o:p></p>
<p class="MsoNormal">Hi,<br>
<br>
I am currently using the command:<br>
<br>
writer_output.UseInputMetaDataDictionaryOff()<br>
<br>
If I don't use it, the image is saved with default origin (0,0) mm and default spacing (1., 1.).<br>
<br>
According to the manual for UseInputMetaDataDictionaryOff(), "this flag defines whether the MetaDataDictionary to use will be the one from the input image or the one already set in the ImageIO object." My understanding is that if I set it off, the DICOM slice
 will be written with current values (and not the input values), which is what I want. But this seems not to happen.<br>
<br>
Relevant parts of my code:<br>
<br>
filter.SetOutputOrigin( new_origin )<br>
writer_output.SetInput(filter.GetOutput())<br>
writer_output.UseInputMetaDataDictionaryOff()<br>
writer_output.SetImageIO(gdcmImageIO)<br>
writer_output.Update()<br>
<br>
Maybe the function UseInputMetaDataDictionaryOff() does not work properly for the itk with Python wrapping?<br>
<br>
Thanks for your answer.<br>
<span style="color:#888888"><br>
<span class="hoenzb">Fabio.</span></span><o:p></o:p></p>
<div>
<div>
<p class="MsoNormal"><br>
-----Original Message-----<br>
From: Timothee Evain [mailto:<a href="mailto:tevain@telecom-paristech.fr">tevain@telecom-paristech.fr</a>]<br>
Sent: Dienstag, 20. September 2016 10:35<br>
To: D'Isidoro Fabio<br>
Cc: <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
Subject: Re: [ITK] [ITK-users] itk Python: origin of DICOM slice does not change<br>
<br>
Hi,<br>
<br>
You refer to the DICOMImageReadWrite example. If you don't change the writing part, I think your problem comes from the fact that you assign the read metadata directly to the output ( ->UseInputMetaDataDictionaryOff() ) so the original origin will be saved.<br>
<br>
HTH,<br>
<br>
Tim<br>
<br>
----- Mail original -----<br>
De: "D'Isidoro Fabio" <<a href="mailto:fisidoro@ethz.ch">fisidoro@ethz.ch</a>><br>
À: <a href="mailto:insight-users@itk.org">insight-users@itk.org</a><br>
Envoyé: Lundi 19 Septembre 2016 22:09:39<br>
Objet: [ITK] [ITK-users] itk Python: origin of DICOM slice does not change<br>
<br>
<br>
<br>
Hi,<br>
<br>
<br>
<br>
Following the code ResampleImageFilter2.cxx, I am applying an Identity Transform to a single DICOM slice, I change the origin to<br>
<br>
<br>
<br>
output_origin = reader.GetOutput().GetOrigin() + offset<br>
<br>
<br>
<br>
and I write the transformed DICOM slice into a new DICOM slice (following the code DICOMImageReadAndWrite.cxx).<br>
<br>
<br>
<br>
I visualize the written DICOM slice with a viewer (it’s an itk-vtk segmentation software).<br>
<br>
<br>
<br>
The transformed DICOM slice looks as it should. Surprisingly however the viewer shows the original input values as origin. This is confirmed by the fact that if I read with an itk script the written DICOM slice again, and I print the information of the DICOM
 slice, the origin is still the original one and not the modified one.<br>
<br>
<br>
<br>
How could this be?<br>
<br>
<br>
<br>
Original DICOM slice:<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Transformed DICOM slice:<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Thank you for help,<br>
<br>
<br>
<br>
Fabio.<br>
<br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">
http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.php" target="_blank">http://www.kitware.com/products/protraining.php</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/insight-users" target="_blank">http://public.kitware.com/mailman/listinfo/insight-users</a><br>
<br>
_______________________________________________<br>
Community mailing list<br>
<a href="mailto:Community@itk.org">Community@itk.org</a><br>
<a href="http://public.kitware.com/mailman/listinfo/community" target="_blank">http://public.kitware.com/mailman/listinfo/community</a><br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.php" target="_blank">http://www.kitware.com/products/protraining.php</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/insight-users" target="_blank">http://public.kitware.com/mailman/listinfo/insight-users</a><o:p></o:p></p>
</div>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</body>
</html>