Jerome Piovano

Research in Medical Imaging and Computer Vision

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

2D_Parzen.C

Parzen Region Segmentation of a 2D Image

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

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

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

        // initialization of the segmentations objects
        FastMarchingInitializer<2,float> init;                               // Fast Marching Initializer
        LevelSet<2,float>                phi(image.shape(), (double) 2.0);   // Level Set of band thickness = 2
        Segmentation<2,float>            segm(&phi, &init, &image);          // Segmentation

        // Initalization of the level set
        //phi.init_small_spheres(10, 30, init);    // parameters : radius, windows
         phi.init_small_spheres(10, 25, init);

        // initialization of the modules
        RegularizationModule<2,float>      regul(1.f);                // parameters : weight
        ParzenRegionModule<2,float,float>  parzen(&image, 1.f, 20);   // parameters : image, weight, variance of the kernel

        segm.add_module(regul);
        segm.add_module(parzen);

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

        timeval _time;
        Globals::initChrono(_time); 

        segm.segment(1000, "parzen_", neg, pos);

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

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

        return 0;
}

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