[Insight-users] Hessian Matrix

[Chris McIntosh]

Hello,

Try:

try
{
writer->Update();
}
catch( itk::ExceptionObject & err )
{
std::cerr << "ExceptionObject caught !" << std::endl;
std::cerr << err << std::endl;
return EXIT_FAILURE;
}

instead of just:
writer->Update();

This will catch the exception (hopefully), and give some more
information about whats happened.

-Chris

예동해 wrote:

> Thank you for your help.
>
> I have tried what you told.
>
> I complied VesselnessMeasureImageFilter.cxx
> However, while debugging, I'm in trouble to use it.
>
> The exception message is "Unhandled exception in xxx.exe
> (KERNEL32.DLL): 0xE06D7363 MICROSOFT C++ Exception."
>
> My inputimage and output image's file format is *.raw.
>
> This exception occurs in "writer->update();"
>
> What can I do in order to solve the pb?
>
> Please help me out.
>
>     --- Original Message ---
>     *From : *"Chris McIntosh"<cmcintos.sfu at gmail.com>
>     *To : *"예동해"<eastsun3 at snu.ac.kr>
>     *Date : *2007/01/18 목요일 오후 4:49:24
>     *Subject : *Re: [Insight-users] Hessian Matrix
>
>     Hello,
>
>     To my knowledge the itk::Hessian3DToVesselnessMeasureImageFilter<
>     TPixel >
>     class will compute the Hessian, which can then be input to the
>     Hessian3DToVesselnessMeasureImageFilter in the ./Code/Basic filters.
>
>     It appears from its header file that it defaultly sets the output
>     image type
>     to: "Image< SymmetricSecondRankTensor< ITK_TYPENAME
>     <itkWin32Header_8h.html#a0> NumericTraits< ITK_TYPENAME
>     TInputImage::PixelType>::RealType,
>     ::itk::GetImageDimension<TInputImage>::ImageDimension >,
>     ::itk::GetImageDimension<TInputImage>::ImageDimension >"
>
>     which is what's needed for the
>     Hessian3DToVesselenssMeasureImageFilter
>
>     http://www.itk.org/Doxygen/html/itkHessianRecursiveGaussianImageFilter_8h-source.html
>
>     http://www.itk.org/Doxygen/html/classitk_1_1Hessian3DToVesselnessMeasureImageFilter.html
>
>
>
>     So in theory you need only do something *like* this: (note I have not
>     compiled or tested this)
>
>     typedef itk::HessianRecursiveGaussianImageFilter< TInputImage >
>     hrgType;
>     typedef itk::Hessian3DToVesselnessMeasureImageFilter<
>     TypeOutputPixel>
>     hvmType;
>
>     hrgType::Pointer hessian = hrgType::New();
>     hessian->SetInput(reader->GetOutput()); //Pipe the reader into the
>     hessian filter
>     hessian->SetSigma(5.0);// Will find vessels of radius 5
>
>     hvmType::Poiner vesselFilter = hvmType::New();
>     vesselFilter->SetInput(hessian->GetOutput()); //Set the hessian as
>     the
>     input to the vessel filter
>     vesselFilter->SetAlpha1(0.5);//
>     vesselFilter->SetAlpha2(2.0);//
>     http://splweb.bwh.harvard.edu:8000/pages/papers/yoshi/node3.html#SECTION00021000000000000000
>
>
>     writer->SetInput(vesselFilter->GetOutput());
>     writer->Update();
>
>     For more details check the included ITK example in:
>     InsightToolkit-2.6.0-src/Examples/Filtering/VesselnessMeasureImageFilter.cxx
>
>
>     HTH,
>
>     Chris
>
>
>     예동해 wrote:
>
>     > Hi all.
>     >
>     > I am trying to compute 3D hessian matrix for vessel enhancement.
>     >
>     > I found some functions that are supposed to be useful for my goal.
>     >
>     > This is Hessian3DToVesselnessMeasureImageFilter in the ./Code/Basic
>     > filters.
>     >
>     > But, I don't know how to use it.
>     >
>     > Please show the way and Help me out.
>     >
>     > I'm looking forward to your reply.
>     >
>     > Thank you.
>     >
>     >------------------------------------------------------------------------
>
>     >
>     >_______________________________________________
>     >Insight-users mailing list
>     >Insight-users at itk.org
>     >http://www.itk.org/mailman/listinfo/insight-users
>     >
>     >
>