Jerome Piovano

Research in Medical Imaging and Computer Vision

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

2D_LocalStats.C

Local Stats Region Segmentation of a 2D Image

#include <Segmentation.h>
#include <SegmentationModule/RegularizationModule.h>
#include <SegmentationModule/LocalStatsRegionModule.h>
#include <SegmentationModule/GaussianRegionModule.h>

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

        BaseImage<2, float> image;
        std::ifstream f("image/synthetic.inr");
        f >> format("inrimage") >> image;

        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 object

        // Initalization of the level set
        //LevelSet<2,float>::Index center(230,170);
        phi.init_border(100);

        // phi.init_small_spheres(20, 50);
        // initialization of the modules

        //note variance = 100 quasi equivalent to chan vese.

        RegularizationModule<2,float>          regul(1.f);                       // parameters : weight
        LocalStatsRegionModule<2,float,float>  local_stats( &image, 0.1f, 10);   // parameters : image, weight, smooth
        //GaussianRegionModule<2,float,float>    gauss(&image, 1.f);               // parameters : image, weight

        local_stats.use_local_variance() = false; 
        //local_stats.sigma() = 20;

        segm.add_module(regul);
        //segm.add_module(gauss);
        segm.add_module(local_stats);
 
        segm.segment(1000, "localstats_");

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


        return 0;
}

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