/home/neoflo/smb4k/SERVEUR/Millie/trunk/src/PredefinedKernel.cpp

Aller à la documentation de ce fichier.

/******************************************************************************
 *       __    _    _   _       _       _   _____                             *
 *       | \  / |  | | | |     | |     | |  | ___|                            *
 *       |  \/  |  | | | |     | |     | |  | |_                              *
 *       |      |  | | | |     | |     | |  |  _|                             *
 *       | |\/| |  | | | |__   | |__   | |  | |__                             *
 *       |_|  |_|  |_| |____|  |____|  |_|  |____|                            *
 * __________________________________________________________________________ *
 *                 Multifunctional Library For Image Processing               *
 *                                                                            *
 *                                                                            *
 *                                                                            *
 *      (c) Copyright 2007 by Humbert Florent                                 *
 *                                                                            *
 *      This program is free software; you can redistribute it and/or modify  *
 *      it under the terms of the GNU General Public License as published by  *
 *      the Free Software Foundation; only version 2 of the License.          *
 *                                                                            *
 *      This program is distributed in the hope that it will be useful,       *
 *      but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *      MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *      GNU General Public License for more details.                          *
 *                                                                            *
 *      You should have received a copy of the GNU General Public License     *
 *      along with this program; if not, write to the Free Software           *
 *      Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA             *
 *      02111-1307, USA.                                                      *
 ******************************************************************************/

#include <cmath>

#include "PredefinedKernel.hpp"

using namespace Millie;

static const float floatSobelVerticalKernel [] =
  {
    -1, -2 , -1,
    0,  0,   0,
    1,  2,   1
  };

static const float floatSobelHorizontalKernel []=
  {
    -1, 0, 1,
    -2, 0, 2,
    -1, 0, 1
  };

static const float floatPrewittVerticalKernel [] =
  {
    -1, -1 , -1,
    0,  0,   0,
    1,  1,   1
  };

static const float floatPrewittHorizontalKernel []=
  {
    -1, 0, 1,
    -1, 0, 1,
    -1, 0, 1
  };


static const float floatLaplaceKernel [] =
  {
    0, 1, 0,
    1,  -4, 1,
    0, 1, 0
  };

static const float floatMDIFKernel []  =
  {
    0, 1, 0, -1, 0,
    1, 2, 0, -2, -1,
    1, 3, 0, -3, -1,
    1, 2, 0, -2, -1,
    0, 1, 0, -1, 0
  };

static const float floatMediumSmoothKernel [] =
  {
    1.0f / 9.0f, 1.0f / 9.0f, 1.0f / 9.0f,
    1.0f / 9.0f,1.0f / 9.0f,1.0f / 9.0f,
    1.0f / 9.0f,1.0f / 9.0f,1.0f / 9.0f
  };

KernelSobelVertical::KernelSobelVertical()
    :
    Kernel(3,3,1,1, floatSobelVerticalKernel)
{}

KernelSobelHorizontal::KernelSobelHorizontal()
    : Kernel(3,3 ,1,1, floatSobelHorizontalKernel)
{}

KernelPrewittVertical::KernelPrewittVertical()
    : Kernel(3,3,1,1, floatPrewittVerticalKernel)
{}

KernelPrewittHorizontal::KernelPrewittHorizontal()
    : Kernel(3,3 ,1,1, floatPrewittHorizontalKernel)
{}

KernelLaplace::KernelLaplace()
    : Kernel(3,3,1,1, floatLaplaceKernel)
{}

KernelMDIF::KernelMDIF()
    : Kernel(5,5,2,2, floatMDIFKernel)
{}


KernelMediumSmooth::KernelMediumSmooth()
    : Kernel(3,3,1,1, floatMediumSmoothKernel)
{}



KernelGaussian::KernelGaussian(int rayon, double sigma)
    : Kernel(2*rayon+1,2*rayon+1,rayon,rayon)
{
  if(sigma <0.0f)
    throw IllegalArgument("KernelGaussian : sigma negative");

  float gaussianKernelFactor=0;

  for(int ky=-rayon;ky<=rayon;ky++)
  {
    for(int kx=-rayon;kx<=rayon;kx++)
    {
      float e = exp( - (kx*kx+ky*ky) / (2*sigma*sigma) );
      gaussianKernelFactor += e;
      setValue(kx+rayon, ky+rayon, e);
    }
  }

  for(int ky = 0; ky< getHeight(); ky++)
   for(int kx = 0; kx < getWidth(); kx++) {
      float valeur = (*this)(kx,ky);
      setValue(kx, ky, valeur/ gaussianKernelFactor);
   }

}

---