There is no error. The problem is I cant see the final image.:-)<div>The output image does not have to be a dicom image.</div><div><br></div><div><br><div class="gmail_quote">2009/5/4 Mathieu Malaterre <span dir="ltr">&lt;<a href="mailto:mathieu.malaterre@gmail.com">mathieu.malaterre@gmail.com</a>&gt;</span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">What is the error message ? I am guessing something like:<br>
<br>
    &quot;A Floating point buffer was passed but the stored pixel type was<br>
not specified.&quot; &quot;This is currently not supported&quot;<br>
<br>
I do not know of any DICOM Media Storage class that can store floating<br>
point data. If you are dealing with CT Image Storage Class, you can<br>
define the Rescale Slope / Intercept to use so that your 32bits<br>
floating point is stored as a 16bits integer type with an appropriate<br>
linear function (floating point parameters).<br>
<br>
2cts<br>
<br>
2009/5/3 Selim Hasan &lt;<a href="mailto:selimhasan@gmail.com">selimhasan@gmail.com</a>&gt;:<br>
<div><div></div><div class="h5">&gt; Hi again,<br>
&gt;<br>
&gt; I added the function that i use to make float number a volume. I tested<br>
&gt;  this code with adding only one pixel value, but it didnt work.<br>
&gt; I cant also write as dicom image, this functions works without error if i<br>
&gt; change the final image to .vtk but I can not see anything in the image.<br>
&gt; Can someone help me to solve this problem.<br>
&gt;                 typedef itk::GDCMImageIO ImageIOType;<br>
&gt; ImageIOType::Pointer gdcmImageIO = ImageIOType::New();<br>
&gt; typedef itk::Image&lt; float, 3 &gt; ImageType;<br>
&gt; typedef itk::ImageFileWriter&lt; ImageType &gt; WriterType;<br>
&gt; ImageType::Pointer image = ImageType::New();<br>
&gt; WriterType::Pointer writer = WriterType::New();<br>
&gt; cout&lt;&lt;&quot;MakeImage&quot;&lt;&lt;endl;<br>
&gt; ImageType::IndexType start;<br>
&gt; ImageType::SizeType  size;<br>
&gt; size[0]  = 144;  // size along X<br>
&gt; size[1]  = 144;  // size along Y<br>
&gt; size[2]  = 60;  // size along Z<br>
&gt; start[0] =   0;  // first index on X<br>
&gt; start[1] =   0;  // first index on Y<br>
&gt; start[2] =   0;  // first index on Z<br>
&gt; ImageType::RegionType region;<br>
&gt; region.SetSize( size );<br>
&gt; region.SetIndex( start );<br>
&gt; image-&gt;SetRegions( region );<br>
&gt; image-&gt;Allocate();<br>
&gt; ImageType::PixelType  initialValue = 0;<br>
&gt; image-&gt;FillBuffer( initialValue );<br>
&gt; cout&lt;&lt;&quot;MakeImage&quot;&lt;&lt;endl;<br>
&gt; ImageType::PixelType   pixelValue;<br>
&gt; ImageType::IndexType pixelIndex;<br>
&gt;<br>
&gt; pixelIndex[0] = 27;   // x position<br>
&gt; pixelIndex[1] = 29;   // y position<br>
&gt; pixelIndex[2] = 37;   // z position<br>
&gt; pixelValue=0.52356; // TEST<br>
&gt; image-&gt;SetPixel(  pixelIndex,   pixelValue );<br>
&gt; cout&lt;&lt;&quot;MakeImage&quot;&lt;&lt;endl;<br>
&gt; writer-&gt;SetInput(image);<br>
&gt;                 writer-&gt;SetImageIO(gdcmImageIO);<br>
&gt; writer-&gt;SetFileName(&quot;C:\\selimmingwinput\\test.dcm&quot;);<br>
&gt; try<br>
&gt; {<br>
&gt; writer-&gt;Update();<br>
&gt; }<br>
&gt;  catch( itk::ExceptionObject &amp; e )<br>
&gt; {<br>
&gt; std::cerr &lt;&lt; &quot;WRITER UPDATE ERROR&quot; &lt;&lt; std::endl;<br>
&gt; std::cerr &lt;&lt; e &lt;&lt; std::endl;<br>
&gt;<br>
&gt; }<br>
&gt; On Thu, Apr 30, 2009 at 11:19 PM, C.Cagatay Bilgin &lt;<a href="mailto:bilgic@cs.rpi.edu">bilgic@cs.rpi.edu</a>&gt;<br>
&gt; wrote:<br>
&gt;&gt;<br>
&gt;&gt; Merhaba Selim,<br>
&gt;&gt; follow &quot; 4.1.1 Creating an Image &quot; Software manual to create an image.<br>
&gt;&gt; Then iterate on the image you just created and follow 4.1.3 Accessing<br>
&gt;&gt; Pixel Data<br>
&gt;&gt; to write the pixel values you read from the text file. Make sure the<br>
&gt;&gt; spacing<br>
&gt;&gt; values and the order you iterate (row,col,depth vs col,row,depth...)are<br>
&gt;&gt; consistent<br>
&gt;&gt; with the given input.<br>
&gt;&gt; Regards,<br>
&gt;&gt; C. Cagatay Bilgin<br>
&gt;&gt; On Apr 30, 2009, at 3:04 AM, Selim Hasan wrote:<br>
&gt;&gt;<br>
&gt;&gt; Hello all,<br>
&gt;&gt;<br>
&gt;&gt; I have a problem for writing an image(3D). I have an 144*144*60 values<br>
&gt;&gt; inside a 3-D image, those values are writen in a txt file.<br>
&gt;&gt;<br>
&gt;&gt; I need to write a program that takes this .txt fýle as an input and writes<br>
&gt;&gt; the volume basically having the same image again. This input comes from an<br>
&gt;&gt; another software.<br>
&gt;&gt;<br>
&gt;&gt; Does anyone know how to write this to an image?<br>
&gt;&gt; Thank You<br>
&gt;&gt; --<br>
&gt;&gt; Selim Hasan<br>
&gt;&gt;<br>
&gt;&gt; Yeditepe University<br>
&gt;&gt; Computer Science Engineering<br>
&gt;&gt;<br>
&gt;&gt; <a href="mailto:selimhasan@gmail.com">selimhasan@gmail.com</a><br>
&gt;&gt; _____________________________________<br>
&gt;&gt; Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
&gt;&gt;<br>
&gt;&gt; Visit other Kitware open-source projects at<br>
&gt;&gt; <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
&gt;&gt;<br>
&gt;&gt; Please keep messages on-topic and check the ITK FAQ at:<br>
&gt;&gt; <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
&gt;&gt;<br>
&gt;&gt; Follow this link to subscribe/unsubscribe:<br>
&gt;&gt; <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Selim Hasan<br>
&gt;<br>
&gt; Yeditepe University<br>
&gt; Computer Science Engineering<br>
&gt;<br>
&gt; <a href="mailto:selimhasan@gmail.com">selimhasan@gmail.com</a><br>
&gt;<br>
&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; 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>
<br>
<br>
<br>
--<br>
</div></div><font color="#888888">Mathieu<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Selim Hasan<br><br>Yeditepe University<br>Computer Science Engineering<br><br><a href="mailto:selimhasan@gmail.com">selimhasan@gmail.com</a><br>
</div>