I tried what you suggested, but it didn&#39;t work. I got the same result though.<div><br></div><div>Ming<br><br><div class="gmail_quote">On Wed, Nov 4, 2009 at 9:56 AM, Bill Lorensen <span dir="ltr">&lt;<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Try making the InputPIxelType short<br>
<br>
 typedef short InputPixelType;<br>
<div><div></div><div class="h5"><br>
<br>
On Wed, Nov 4, 2009 at 10:38 AM, Ming Chao &lt;<a href="mailto:mingchao2005@gmail.com">mingchao2005@gmail.com</a>&gt; wrote:<br>
&gt; Hi ITKers.<br>
&gt; Recently I was asking questions here about how to write out an multiframe<br>
&gt; dicom file. Thanks for all the help.<br>
&gt; As a matter of fact, it is not so hard to save a 3D volumetric image into a<br>
&gt; multiframe dicom file. Please see the attached code. But I do encounter a<br>
&gt; problem in the output dicom file. Basically the range of the voxel values of<br>
&gt; the output file  much narrower than expected ( I attached the images from<br>
&gt; the output dicom and original dicom). I double checked the data types and<br>
&gt; they seemed fine. Could anyone help? thanks.<br>
&gt; Ming<br>
&gt;<br>
&gt; int main( int argc, char* argv[] ) {<br>
&gt; // Verify the number of parameters in the command line<br>
&gt; if( argc &lt; 4 ) {<br>
&gt; std::cerr &lt;&lt; &quot;Usage: &quot; &lt;&lt; std::endl;<br>
&gt; std::cerr &lt;&lt; argv[0] &lt;&lt; &quot; vtkImage InputDicomImage &quot;;<br>
&gt; std::cerr &lt;&lt; &quot; OutputImage \n&quot;;<br>
&gt; return EXIT_FAILURE;<br>
&gt;     }<br>
&gt; // the following is for an input 2D dicom file<br>
&gt; typedef unsigned short InputPixelType;<br>
&gt; const unsigned int   InputDimension = 3;<br>
&gt; typedef itk::Image&lt; InputPixelType, InputDimension &gt; InputImageType;<br>
&gt; typedef itk::ImageFileReader&lt; InputImageType &gt; InputReaderType;<br>
&gt; InputReaderType::Pointer reader1 = InputReaderType::New();<br>
&gt; reader1-&gt;SetFileName( argv[1] );<br>
&gt; typedef itk::GDCMImageIO          InputImageIOType;<br>
&gt; InputImageIOType::Pointer gdcmImageIO = InputImageIOType::New();<br>
&gt; reader1-&gt;SetImageIO( gdcmImageIO );<br>
&gt; reader1-&gt;Update();<br>
&gt; std::cout &lt;&lt; &quot; (1) first dicom file is read in successfully ....&quot; &lt;&lt;<br>
&gt; std::endl;<br>
&gt; // The following is a print out of the tag info<br>
&gt; typedef itk::MetaDataDictionary   DictionaryType;<br>
&gt; const  DictionaryType &amp; dictionary = gdcmImageIO-&gt;GetMetaDataDictionary();<br>
&gt; typedef itk::MetaDataObject&lt; std::string &gt; MetaDataStringType;<br>
&gt; DictionaryType::ConstIterator itr = dictionary.Begin();<br>
&gt; DictionaryType::ConstIterator end = dictionary.End();<br>
&gt; while( itr != end )<br>
&gt;     {<br>
&gt; itk::MetaDataObjectBase::Pointer  entry = itr-&gt;second;<br>
&gt; MetaDataStringType::Pointer entryvalue =<br>
&gt; dynamic_cast&lt;MetaDataStringType *&gt;( entry.GetPointer() ) ;<br>
&gt; if( entryvalue )<br>
&gt; {<br>
&gt; std::string tagkey   = itr-&gt;first;<br>
&gt; std::string tagvalue = entryvalue-&gt;GetMetaDataObjectValue();<br>
&gt; std::cout &lt;&lt; tagkey &lt;&lt;  &quot; = &quot; &lt;&lt; tagvalue &lt;&lt; std::endl;<br>
&gt; }<br>
&gt; ++itr;<br>
&gt;     }<br>
&gt; // the following is for a 3D volume file<br>
&gt; typedef float PixelType;<br>
&gt; const unsigned int   Dimension = 3;<br>
&gt; typedef itk::Image&lt; PixelType, Dimension &gt; ImageType;<br>
&gt; typedef itk::ImageFileReader&lt; ImageType  &gt;  ReaderType;<br>
&gt;     ReaderType::Pointer reader2 = ReaderType::New();<br>
&gt;     reader2-&gt;SetFileName( argv[2] );<br>
&gt; reader2-&gt;Update();<br>
&gt; std::cout &lt;&lt; &quot; (2) the 3D volume file is read in successfully ....&quot; &lt;&lt;<br>
&gt; std::endl;<br>
&gt;     // the following is for an output 3D dicom file<br>
&gt; typedef float OutputPixelType;<br>
&gt; const unsigned int   OutputDimension = 3;<br>
&gt; typedef itk::Image&lt; OutputPixelType, OutputDimension &gt; OutputImageType;<br>
&gt; typedef itk::ImageSeriesWriter&lt; ImageType, OutputImageType &gt;<br>
&gt; OutputWriterType;<br>
&gt; OutputWriterType::Pointer writer = OutputWriterType::New();<br>
&gt; writer-&gt;SetFileName( argv[3] );<br>
&gt; writer-&gt;SetInput( reader2-&gt;GetOutput() );<br>
&gt; writer-&gt;SetMetaDataDictionary( reader1-&gt;GetMetaDataDictionary() );<br>
&gt; gdcmImageIO-&gt;KeepOriginalUIDOn();<br>
&gt; writer-&gt;SetImageIO( gdcmImageIO );<br>
&gt; writer-&gt;Update();<br>
&gt; std::cout &lt;&lt; &quot; (3) the output dicom file is saved successfully and we are<br>
&gt; done....&quot; &lt;&lt; std::endl;<br>
&gt;<br>
&gt; return 0;<br>
&gt;<br>
</div></div>&gt; _____________________________________<br>
&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;<br>
&gt; Visit other Kitware open-source projects at<br>
&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;<br>
&gt; Kitware offers ITK Training Courses, for more information visit:<br>
&gt; <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
&gt;<br>
&gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;<br>
&gt; Follow this link to subscribe/unsubscribe:<br>
&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt;<br>
&gt;<br>
</blockquote></div><br></div>