From: Leonardo Flórez-Valencia Date: Tue, 6 Dec 2016 00:26:02 +0000 (-0500) Subject: ... X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=commitdiff_plain;ds=sidebyside;h=4f7a9ad44485de053935926e227c76980152bd97;p=FrontAlgorithms.git ... --- diff --git a/appli/experiments/CMakeLists.txt b/appli/experiments/CMakeLists.txt index 3be8e1c..65b5df1 100644 --- a/appli/experiments/CMakeLists.txt +++ b/appli/experiments/CMakeLists.txt @@ -5,4 +5,7 @@ INCLUDE_DIRECTORIES( ADD_EXECUTABLE(fpa_experiments_Skeleton_00 Skeleton_00.cxx) TARGET_LINK_LIBRARIES(fpa_experiments_Skeleton_00 ${ITK_LIBRARIES} ${VTK_LIBRARIES}) +ADD_EXECUTABLE(fpa_experiments_Hausdorff Hausdorff.cxx) +TARGET_LINK_LIBRARIES(fpa_experiments_Hausdorff ${ITK_LIBRARIES} ${VTK_LIBRARIES}) + ## eof - $RCSfile$ diff --git a/appli/experiments/Hausdorff.cxx b/appli/experiments/Hausdorff.cxx new file mode 100644 index 0000000..dce01ef --- /dev/null +++ b/appli/experiments/Hausdorff.cxx @@ -0,0 +1,48 @@ +#include +#include +#include + +typedef double TScalar; +typedef itk::Image< TScalar, 3 > TImage; +typedef itk::ImageFileReader< TImage > TReader; +typedef itk::HausdorffDistanceImageFilter< TImage, TImage > THausdorff; + +int main( int argc, char* argv[] ) +{ + if( argc < 3 ) + { + std::cerr + << "Usage: " << argv[ 0 ] + << " image1 image2" + << std::endl; + return( 1 ); + + } // fi + std::string image1_name = argv[ 1 ]; + std::string image2_name = argv[ 2 ]; + + TReader::Pointer reader1 = TReader::New( ); + reader1->SetFileName( image1_name ); + reader1->Update( ); + + TReader::Pointer reader2 = TReader::New( ); + reader2->SetFileName( image2_name ); + reader2->Update( ); + + THausdorff::Pointer hausdorff = THausdorff::New( ); + hausdorff->SetInput1( reader1->GetOutput( ) ); + hausdorff->SetInput2( reader2->GetOutput( ) ); + hausdorff->SetUseImageSpacing( true ); + hausdorff->Update( ); + + std::cout + << image1_name << " " + << image2_name << " " + << hausdorff->GetHausdorffDistance( ) << " " + << hausdorff->GetAverageHausdorffDistance( ) + << std::endl; + + return( 0 ); +} + +// eof - $RCSfile$