<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div>the code c++ is :<br>#if defined(_MSC_VER)<br>#pragma warning ( disable : 4786 )<br>#endif<br><br>#ifdef __BORLANDC__<br>#define ITK_LEAN_AND_MEAN<br>#endif<br>#include "itkHessian3DToVesselnessMeasureImageFilter.h"<br>#include "itkHessianRecursiveGaussianImageFilter.h"<br>#include "itkSymmetricSecondRankTensor.h"<br>#include "itkImage.h"<br>#include "itkImageFileReader.h"<br>#include "itkImageFileWriter.h"<br><br>int main( int argc, char *argv[] )<br>{<br><br> if( argc < 3 )<br> {<br> std::cerr << "Usage: inputImage outputImage [sigma] [alpha_1] [alpha_2]" << std::endl;<br> }<br> <br> const unsigned int
Dimension = 3;<br> typedef double InputPixelType;<br> typedef float OutputPixelType;<br> <br> typedef itk::Image< InputPixelType, Dimension > InputImageType;<br> typedef itk::Image< OutputPixelType, Dimension > OutputImageType;<br> <br> typedef itk::HessianRecursiveGaussianImageFilter< <br> InputImageType > HessianFilterType;<br> typedef
itk::Hessian3DToVesselnessMeasureImageFilter<<br> OutputPixelType > VesselnessMeasureFilterType;<br> typedef itk::ImageFileReader< InputImageType > ReaderType;<br> typedef itk::ImageFileWriter< OutputImageType > WriterType;<br> <br> HessianFilterType::Pointer hessianFilter = HessianFilterType::New();<br> VesselnessMeasureFilterType::Pointer vesselnessFilter = <br> VesselnessMeasureFilterType::New();<br><br> ReaderType::Pointer reader = ReaderType::New();<br> WriterType::Pointer writer = WriterType::New();<br> <br> reader->SetFileName( argv[1] );<br> hessianFilter->SetInput(
reader->GetOutput() );<br> if( argc >= 4 )<br> {<br> hessianFilter->SetSigma( static_cast< double >(atof(argv[3])) );<br> }<br> vesselnessFilter->SetInput( hessianFilter->GetOutput() );<br> writer->SetInput( vesselnessFilter->GetOutput() );<br> writer->SetFileName( argv[2] );<br> if( argc >= 5 )<br> {<br> vesselnessFilter->SetAlpha1( static_cast< double >(atof(argv[4])));<br> }<br> if( argc >= 6 )<br> {<br> vesselnessFilter->SetAlpha2( static_cast< double >(atof(argv[5])));<br> }<br><br> writer->Update();<br> return EXIT_SUCCESS;<br>}<br><br> <br>I'm testing the vesselness filter on an MRA image .<br></div><div style="font-family: times new roman,new york,times,serif;
font-size: 12pt;"><br><div style="font-family: arial,helvetica,sans-serif; font-size: 13px;"><font size="2" face="Tahoma"><hr size="1"><b><span style="font-weight: bold;">De :</span></b> John Drescher <drescherjm@gmail.com><br><b><span style="font-weight: bold;">À :</span></b> Syrine Sahmim <syrine.sahmim@yahoo.fr><br><b><span style="font-weight: bold;">Envoyé le :</span></b> Jeudi, 20 Août 2009, 13h48mn 45s<br><b><span style="font-weight: bold;">Objet :</span></b> Re: Re : needs a lot of recommandations<br></font><br>On Thu, Aug 20, 2009 at 5:58 AM, Syrine Sahmim<<a ymailto="mailto:syrine.sahmim@yahoo.fr" href="mailto:syrine.sahmim@yahoo.fr">syrine.sahmim@yahoo.fr</a>> wrote:<br>> the result that i had is on the joined file resultat<br>><br><br>The last picture was what I expected except for the crash. If you post<br>you code I will take a look at this.<br><br>John<br></div></div></div><br>
</body></html>