C:/cmcintos/defOrgs/examples/Phys_LevelSetTest/Phys_LevelSetTest.cxx

#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif

#include "itkImage.h"
#include "itkGeodesicActiveContourLevelSetImageFilter.h"

#include "itkCurvatureAnisotropicDiffusionImageFilter.h"
#include "itkGradientMagnitudeRecursiveGaussianImageFilter.h"
#include "itkSigmoidImageFilter.h"
#include "itkFastMarchingImageFilter.h"
#include "itkRescaleIntensityImageFilter.h"
#include "itkBinaryThresholdImageFilter.h"
#include "itkImageFileReader.h"
#include "itkImageFileWriter.h"
#include "itkRescaleIntensityImageFilter.h"
#include "itkCastImageFilter.h"

#include "Phys_LevelSet.h"
#include "Geometric.h"
#include "Geom_vGeometry.h"

#include <string>

using namespace mial;

int main( int argc, char *argv[] )
{
        std::string inputFilename = "..\\basic\\data\\square.bmp";
        std::string outputFilename = "C:\\Phys_LevelSetTest_output.png";
        std::string nodeInputFileName = "nodesInSimple.txt";
        std::string nodeOutputFileName = "C:\\nodesOut.txt";

        // For Andy's computer:
        //std::string nodeInputFileName = "nodesInSimple.txt";
        //std::string inputFilename = "..\\build\\data\\letterU.png";
        //std::string outputFilename = "C:\\Documents and Settings\\Andy\\Desktop\\Phys_LevelSetTest_output.png";

        typedef itk::Image< double, 2 > ImageType; 

        typedef itk::Image< unsigned char, 2 > OutputImageType;

        typedef itk::ImageFileReader< ImageType >  ReaderType;
        ReaderType::Pointer reader = ReaderType::New();

        typedef itk::ImageFileWriter< OutputImageType >  WriterType;
        WriterType::Pointer writer = WriterType::New();

        typedef itk::RescaleIntensityImageFilter< 
                ImageType, 
                OutputImageType > CastFilterType;
        CastFilterType::Pointer castFilter = CastFilterType::New();

        typedef vnl_matrix< double > vType;

        reader->SetFileName(inputFilename.c_str());
        writer->SetFileName(outputFilename.c_str());

        std::cout << "This is a test of the Phys_LevelSet class." << std::endl << std::endl;
 
        Geom_vGeometry<vType, double, 2> geom;                  // geometric object
        std::cout << "Geometric object created..." << std::endl;
        //read nodes from file:
        geom.readNodesFromFile(nodeInputFileName);
        std::cout << "Nodes read from file..." << std::endl;
        std::cout << "numNodes: " << geom.numNodes() << std::endl;

        Phys_LevelSet<double, vType, 2>  physLevelSetLayer;                             // Phys_LevelSet object
        std::cout << "Phys_LevelSet object created..." << std::endl;

        try
        {
                // set the input image of the Phys_LevelSet object:
                reader->Update();
                physLevelSetLayer.setInput( reader->GetOutput() );

                // create the distance image using node information:
                std::cout << "Starting fast marching initialization of distance image..." << std::endl;
                physLevelSetLayer.initializeDistanceImg(&geom);
                std::cout << "Finished initializing distance image..." << std::endl;

                // run the level set simulation:
                std::cout << "Starting geodesic active contour level set simulate()..." << std::endl;
                physLevelSetLayer.simulate(); 
                std::cout << "Finished level set simulation..." << std::endl;
                physLevelSetLayer.printInternalValues(); // cout member variable status

                // DEBUG:
                // Output the distance image to file:
                //castFilter->SetInput( physLevelSetLayer.getDistanceImage() );
                //writer->SetInput( castFilter->GetOutput() );
                //castFilter->SetOutputMinimum( 0 );
                //castFilter->SetOutputMaximum( 255 );

                // Output the segmentation result:
                castFilter->SetInput( physLevelSetLayer.getOutputImage() );
                writer->SetInput( castFilter->GetOutput() );
                castFilter->SetOutputMinimum( 0 );
                castFilter->SetOutputMaximum( 255 );
                
                writer->Update();
                std::cout << "Output image written to " << writer->GetFileName() << std::endl;
        }
        catch(itk::ExceptionObject & err)
        {
                std::cerr << "ExceptionObject caught:" << std::endl;
                std::cerr << err << std::endl;
                return -1;
        }

        return EXIT_SUCCESS;
}

---