Jerome Piovano

Research in Medical Imaging and Computer Vision

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

2D_Anisotropic_Blur.C

anisotropic Blur of a 2D Image

#include <Image.H>
#include <DericheFilter.h>
#include <Globals.h>

#include <image/RecFilters.H>
#include <image/IIR.H>


int main(int argc, char* argv[])
{
#if 0
        using namespace Images;

        // image to blur
        BaseImage<2, float> image;
        //std::ifstream f("image/tiger.inr");
        std::ifstream f("image/synth3/bruit/000000.inr");
        f >> format("inrimage") >> image;

        // anisotropic variance
        //float sigma[2] = {20, 2};
        float sigma = 400000;

        DericheFilter<2,float> deriche(sigma);

        timeval _time;
        levelset::Globals::initChrono(_time); 

        // filter the image
        deriche.filter(image);
        
        levelset::Globals::printChrono("\n      => Anisotropic Blur : ", _time);
        levelset::Globals::initChrono(_time);  

        std::ofstream f1("2D_Anisotropic_Blur.inr"); 
        f1 << format("inrimage") << image; 
        f1.close();

        return 0;
#endif
#if 1 //comparaison blue theo, et cimg

        typedef           Images::Filters::IIR::Separable::GaussianFilter<4,double> DericheFilter_2;
        
        // image to blur
        BaseImage<2, float> image;
        std::ifstream f("image/tiger.inr");
        //std::ifstream f("image/synth3/bruit/000000.inr");
        f >> format("inrimage") >> image;

        // anisotropic variance
        float sigma = 20;

        DericheFilter_2 deriche(sigma, 0);
        DericheFilter<2,float> deriche_moi(sigma);

        timeval _time;

        // filter the image theo
        for (int i = 0 ; i < 20 ; i++){
                BaseImage<2, float> image2 = image.copy();
        
                levelset::Globals::initChrono(_time);  
                image2 = Images::Filter<BaseImage<2, float> >(image2, deriche);
                levelset::Globals::printChrono("\n      => Anisotropic Blur : ", _time);
        
                std::ostringstream name;
                name.str("");
                name <<  "2D_Anisotropic_Blur_" << std::setw(6) << std::setfill('0') << i << ".inr";
                std::ofstream f1(name.str().c_str());
                f1 << format("inrimage") << image2;
                f1.close();
        }

        //filter the image CImg
        for (int i = 0 ; i < 20 ; i++){
                BaseImage<2, float> image2 = image.copy();
        
                levelset::Globals::initChrono(_time);  
                deriche_moi.filter(image2);
                levelset::Globals::printChrono("\n      => Anisotropic Blur : ", _time);
        
                std::ostringstream name;
                name.str("");
                name <<  "2D_Anisotropic_Blur_" << std::setw(6) << std::setfill('0') << i << ".inr";
                std::ofstream f1(name.str().c_str());
                f1 << format("inrimage") << image2;
                f1.close();
        }

        

#endif
}

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