Jerome Piovano

Research in Medical Imaging and Computer Vision

• Home
• Research
• Publications
• Teaching
• Biography
• Shared code and data

2D_GeodesicActiveContour.C

Geodesic Active Contour Segmentation of a 2D Image

#include <Segmentation.h>
#include <SegmentationModule/RegularizationModule.h>
#include <SegmentationModule/GeodesicActiveContourModule.h>

int main(int argc, char* argv[])
{       
        using namespace Images;
        using namespace levelset; 
        using namespace segmentation;

        BaseImage<2, float> image;
        std::ifstream f("image/head3.inr");
        f >> format("inrimage") >> image;
        Globals::normalize(image, 0.0f, 250.0f);

        FastMarchingInitializer<2,float> init;                               // Fast Marching Initializer
        LevelSet<2,float>                phi(image.shape(), (double) 2.0);   // Level Set
        Segmentation<2,float>            segm(&phi, &init, &image);          // Segmentation

        // Initalization of the level set
        phi.init_border(120, init);

        // Initialization of the module
        // IMPORTANT Do not add a regularization term using the geodesic active contour module
        GeodesicActiveContourModule<2,float,float>  gac(&image, 1.f, -1.f , 2);   // parameters : image, weight, balloon, smooth

        gac.balloon() = 2;

        segm.add_module(gac);


        // launch the segmentation
        Segmentation<2,float>::RGBPixel neg(215,152,0); 
        Segmentation<2,float>::RGBPixel pos(255,252,0); 

        timeval _time;

        std::ostringstream name;
        name.str("");  
        name << "mkdir CVPR/";
        system (name.str().c_str()); 
        name.str("");  
        name << "mkdir CVPR/gaussian_cortex";
        system (name.str().c_str());
        name.str("");  
        name << "CVPR/gaussian_cortex/geodesic";
        Globals::initChrono(_time);

        segm.segment(2000, name.str().c_str(), neg, pos);

        Globals::printChrono("\n        => Gaussian : ", _time);
        Globals::initChrono(_time); 

        std::cout << "OK" << std::endl;


        return 0;
}

For further information, please contact Jerome Piovano - Last update 2008-02-08