<div dir="ltr">It is different but related. The handling of rescale/slope in the current implementation is tricky.<br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 21, 2013 at 10:39 AM, Li, Wen <span dir="ltr"><<a href="mailto:WLi@childrensnational.org" target="_blank">WLi@childrensnational.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Thank you for your responding.
<div>The bug you reported is related with my issue but slightly different.</div>
<div><br>
</div>
<div>I can't even confirm whether the output dicom retains the specific RescaleIntercept and RescaleSlope I set. The writer simple stopped working at the first slice when gdcm is doing inverse rescaling. Of course I use a double RescaleSlope like normally PET
images need. I suspect it is the memory (pixel type) issue. I will investigate it.</div>
<div><br>
</div>
<div>Wen<br>
<div>
<div style="font-size:16px;font-family:Times New Roman">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> Bill Lorensen [<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a>]<br>
<b>Sent:</b> Wednesday, August 21, 2013 10:24 AM<br>
<b>To:</b> Li, Wen<br>
<b>Cc:</b> <a href="mailto:insight-users@itk.org" target="_blank">insight-users@itk.org</a><br>
<b>Subject:</b> Re: [Insight-users] DICOMWriter for wide dynamic intensity range<br>
</font><br>
</div><div><div class="h5">
<div></div>
<div>
<div dir="ltr">I have an outstanding bug request that deals with similar issues. I hope to get to it soon...<br>
<a href="https://issues.itk.org/jira/browse/ITK-3194" target="_blank">https://issues.itk.org/jira/browse/ITK-3194</a><br>
<br>
</div>
<div class="gmail_extra"><br>
<br>
<div class="gmail_quote">On Wed, Aug 21, 2013 at 9:38 AM, Li, Wen <span dir="ltr">
<<a href="mailto:WLi@childrensnational.org" target="_blank">WLi@childrensnational.org</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Only ITK in Debug mode allows manipulating RescaleIntercept and RescaleSlope AND if output pixel type of the dicom writer is float or double. That's fine.
<div><br>
</div>
<div>Now the problem is gdcm::Rescaler may not work correctly. It stopped at gdcmRescaler::InverseRescaleFunction without finishing it. </div>
<div><br>
</div>
<div>I suppose that function handles inverse rescaling for float input pixel type. I don't expect anyone to look into that problem for me but just a little frustrated.</div>
<div><br>
</div>
<div>Wen<br>
<div>
<div style="font-size:16px;font-family:Times New Roman">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> <a href="mailto:insight-users-bounces@itk.org" target="_blank">
insight-users-bounces@itk.org</a> [<a href="mailto:insight-users-bounces@itk.org" target="_blank">insight-users-bounces@itk.org</a>] on behalf of Li, Wen [<a href="mailto:WLi@childrensnational.org" target="_blank">WLi@childrensnational.org</a>]<br>
<b>Sent:</b> Monday, August 19, 2013 11:24 PM<br>
<b>To:</b> <a href="mailto:insight-users@itk.org" target="_blank">insight-users@itk.org</a><br>
<b>Subject:</b> [Insight-users] DICOMWriter for wide dynamic intensity range<br>
</font><br>
</div>
<div>
<div>
<div></div>
<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">Hi Insight users,
<div><br>
</div>
<div>I am talking about PET images. The data I am dealing with has a wide dynamic intensity range.</div>
<div>Now I am trying to write it as a DICOM image. </div>
<div><br>
</div>
<div>As we all know, DICOM only stores 16bit integer as pixel value. I tried to manipulate DICOM tag 0028|1052 and 0028|1053 for the RescaleIntercept and RescaleSlope. Failed. They were always reset to be 0 and 1.</div>
<div><br>
</div>
<div>Is there any suggestion about how to store the original intensity values instead of having them chopped off to be 16bit integer. In my case, I probably need float non-1 RescaleSlope. </div>
<div><br>
</div>
<div>BTW, I define my writer as </div>
<div>typedef itk::ImageSeriesWriter< OutputImageType, OutputImageType2D > SeriesWriterType;</div>
<div>OutputImageType::PixelType is float and OutputImageType2D::PixelType is signed short.</div>
<div><br>
</div>
<div>Thanks,</div>
<div>Wen</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<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://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
<br>
</blockquote>
</div>
<br>
<br clear="all">
<br>
-- <br>
Unpaid intern in BillsBasement at noware dot com<br>
</div>
</div>
</div></div></div>
</div>
</div>
</div>
</div>
</blockquote></div><br><br clear="all"><br>-- <br>Unpaid intern in BillsBasement at noware dot com<br>
</div>