Jerome Piovano

Research in Medical Imaging and Computer Vision

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

3D_ChanVese.C

Chan-Vese Region Segmentation of a 3D Image

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

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

        BaseImage<3, unsigned> image;
        std::ifstream f("image/brain.inr");
        f >> image; //format("inrimage") >> image;
#if 0
        FastMarchingInitializer<3,float> init;                               // Fast Marching Initializer
        LevelSet<3,float>                phi(image.shape(), (double) 2.0);   // Level Set of band thickness = 2
        Segmentation<3,float>            segm(&phi, &init, &image);          // Segmentation

        // Initalization of the level set
        phi.init_border(60);    

        // initialization of the modules
        RegularizationModule<3,float>        regul(1.f);           // parameters : weight
        GaussianRegionModule<3,float,float>  gauss(&image, 1.f);   // parameters : image, weight

        gauss.use_variance(false);    // => false = Chan-Vese region segmentation
        gauss.set_weight(0.01f);
        gauss.set_data_term(image);

        segm.add_module(regul);
        segm.add_module(gauss);

        // launch the segmentation
        segm.segment(1000, "gaussian_");

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

#endif
        return 0;
}

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