<br>Hi zlf,<br><br>You may want to try saving that file using the MetaImage or the VTK<br>file formats.<br><br>Both of them would be able to manage pixels with multiple components.<br><br>In the case of the Hessian, you are actually generating an image with<br>
six components per pixel. (there is only six, because the Hessian matrix<br>is symmetric).<br><br>You could load those files into Paraview (<a href="http://www.paraview.org">www.paraview.org</a>) if you <br>need to visualize them. (probably want to extract their individual <br>
components inside paraview as well)<br><br><br> Regards,<br><br><br> Luis<br><br> <br>----------------------------------------------------------------------------<br><div class="gmail_quote">On Sun, Oct 18, 2009 at 10:14 PM, 庄凌峰 <span dir="ltr"><<a href="mailto:jxdw_zlf@yahoo.com.cn">jxdw_zlf@yahoo.com.cn</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font-family: inherit; font-style: inherit; font-variant: inherit; font-weight: inherit; font-size: inherit; line-height: inherit; font-size-adjust: inherit; font-stretch: inherit;" valign="top">
<div>Hello<br>I am new to itk.</div>
<div>I want to employ HessianRecursiveGaussianImageFilter in my program.<br>In order to see its result,I modified DicomSeriesReadImageWrite2.cxx. Wanna to result its output to a dicom file.<br>But it failed with <br> itk::ExceptionObject (013FDEAC)<br>
Location: "void __thiscall itk::GDCMImageIO::Write(const void *)"<br> File: ..\..\..\Insight\Code\IO\itkGDCMImageIO.cxx<br> Line: 1659<br> Description: itk::ERROR: GDCMImageIO(02E72FE0): DICOM does not support RGBPixels with components != 3</div>
<div> </div>
<div>The changed made to DicomSeriesReadImageWrite2.cxx is listed below.</div>
<div> typedef itk::HessianRecursiveGaussianImageFilter< ImageType > myFilterType;<br> typedef myFilterType::OutputImageType myHessianImageType;</div>
<div> // Create a Filter <br> myFilterType::Pointer filter = myFilterType::New();</div>
<div><br> // Connect the input images<br> filter->SetInput( reader->GetOutput() ); </div>
<div> // Select the value of Sigma<br> filter->SetSigma( 2.5 ); </div>
<div><br> // Execute the filter<br> filter->Update();</div>
<div><br> writer->SetInput( filter->GetOutput() );</div>
<div><br>Thanks!</div>
<div>zlf</div></td></tr></tbody></table><br>
<hr size="1"><a href="http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/" target="_blank"> 好玩贺卡等你发，邮箱贺卡全新上线！</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>
Please keep messages on-topic and check the ITK FAQ at: <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>