[Insight-users] How to save the result of itkHessianRecursiveGaussianImageFilter

[Luis Ibanez]

Hi zlf,

You may want to try saving that file using the MetaImage or the VTK
file formats.

Both of them would be able to manage pixels with multiple components.

In the case of the Hessian, you are actually generating an image with
six components per pixel. (there is only six, because the Hessian matrix
is symmetric).

You could load those files into Paraview (www.paraview.org) if you
need to visualize them. (probably want to extract their individual
components inside paraview as well)


      Regards,


            Luis


----------------------------------------------------------------------------
On Sun, Oct 18, 2009 at 10:14 PM, 庄凌峰 <jxdw_zlf at yahoo.com.cn> wrote:

> Hello
> I am new to itk.
> I want to employ HessianRecursiveGaussianImageFilter in my program.
> In order to see its result,I modified DicomSeriesReadImageWrite2.cxx. Wanna
> to result its output to a dicom file.
> But it failed with
>   itk::ExceptionObject (013FDEAC)
>   Location: "void __thiscall itk::GDCMImageIO::Write(const void *)"
>   File: ..\..\..\Insight\Code\IO\itkGDCMImageIO.cxx
>   Line: 1659
>   Description: itk::ERROR: GDCMImageIO(02E72FE0): DICOM does not support
> RGBPixels with components != 3
>
> The changed made to DicomSeriesReadImageWrite2.cxx is listed below.
>  typedef itk::HessianRecursiveGaussianImageFilter< ImageType >
> myFilterType;
>  typedef myFilterType::OutputImageType myHessianImageType;
>  // Create a  Filter
>  myFilterType::Pointer filter = myFilterType::New();
>
>  // Connect the input images
>  filter->SetInput( reader->GetOutput() );
>  // Select the value of Sigma
>  filter->SetSigma( 2.5 );
>
>  // Execute the filter
>  filter->Update();
>
>  writer->SetInput( filter->GetOutput() );
>
> Thanks!
> zlf
>
> ------------------------------
> 好玩贺卡等你发，邮箱贺卡全新上线！<http://cn.rd.yahoo.com/mail_cn/tagline/card/*http://card.mail.cn.yahoo.com/>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20091021/d863e145/attachment.htm>