[Insight-users] DICOMWriter for wide dynamic intensity range

[Li, Wen]

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.

Now the problem is gdcm::Rescaler may not work correctly. It stopped at gdcmRescaler::InverseRescaleFunction without finishing it.

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.

Wen
________________________________
From: insight-users-bounces at itk.org [insight-users-bounces at itk.org] on behalf of Li, Wen [WLi at childrensnational.org]
Sent: Monday, August 19, 2013 11:24 PM
To: insight-users at itk.org
Subject: [Insight-users] DICOMWriter for wide dynamic intensity range

Hi Insight users,

I am talking about PET images. The data I am dealing with has a wide dynamic intensity range.
Now I am trying to write it as a DICOM image.

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.

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.

BTW, I define my writer as
typedef itk::ImageSeriesWriter< OutputImageType, OutputImageType2D > SeriesWriterType;
OutputImageType::PixelType is float and OutputImageType2D::PixelType is signed short.

Thanks,
Wen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20130821/6ee5d721/attachment.htm>