vlvoxelop_sobel.h

Go to the documentation of this file.

// $Id: vlvoxelop_sobel.h,v 1.1 2003/08/12 04:18:55 sarang Exp $
//
// OpenVL - A Library for working with volumetric datasets.
//          http://openvl.sf.net
//
// Copyright (C) 2000-2003  Sarang Lakare <sarang#users.sf.net>
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Library General Public
// License as published by the Free Software Foundation; either
// version 2 of the License, or (at your option) any later version.
//
// This library 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
// Library General Public License for more details.
//
// You should have received a copy of the GNU Library General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA.
//
// Please report all bugs and problems to openvl-devel@lists.sf.net or
// at the OpenVL website.
//

#ifndef _vlVoxelOpSobel_h
#define _vlVoxelOpSobel_h

#include "vlvoxeloperator.h"
#include "vlvarlistadmin.h"
#include "vloffset.h"

template <typename DataType, vlLayoutType Layout>
class vlVoxelOpSobel : public vlVoxelOperator<DataType, Layout>
{
public:
  vlVoxelOpSobel()
  {
    // for X direction
    m_offsetsX[0] = vlOffset(1,1,0);
    m_offsetsX[1] = vlOffset(1,-1,0);
    m_offsetsX[2] = vlOffset(1,0,1);
    m_offsetsX[3] = vlOffset(1,0,-1);

    m_offsetsX[4] = vlOffset(-1,1,0);
    m_offsetsX[5] = vlOffset(-1,-1,0);
    m_offsetsX[6] = vlOffset(-1,0,1);
    m_offsetsX[7] = vlOffset(-1,0,-1);

    // for Y direction
    m_offsetsY[0] = vlOffset(1,1,0);
    m_offsetsY[1] = vlOffset(-1,1,0);
    m_offsetsY[2] = vlOffset(0,1,1);
    m_offsetsY[3] = vlOffset(0,1,-1);

    m_offsetsY[4] = vlOffset(-1,-1,0);
    m_offsetsY[5] = vlOffset(1,-1,0);
    m_offsetsY[6] = vlOffset(0,-1,1);
    m_offsetsY[7] = vlOffset(0,-1,-1);

    // for Z direction
    m_offsetsZ[0] = vlOffset(1,0,1);
    m_offsetsZ[1] = vlOffset(-1,0,1);
    m_offsetsZ[2] = vlOffset(0,1,1);
    m_offsetsZ[3] = vlOffset(0,-1,1);

    m_offsetsZ[4] = vlOffset(-1,0,-1);
    m_offsetsZ[5] = vlOffset(1,0,-1);
    m_offsetsZ[6] = vlOffset(0,1,-1);
    m_offsetsZ[7] = vlOffset(0,-1,-1);
  }

  vlVoxelOpType type() { return vlVoxelOp::Sobel; };

  std::string name() { return ("Sobel"); };

  vlLayoutType layout() { return Layout; };

  vlVarList * config() { return (&m_varList); };

  bool getValue(vlVolIterConst<DataType, Layout> & iter, DataType & value)
  { return (false); };

  bool getValue(vlVolIterConst<DataType, Layout> & iter, vlVoxelOpValue & value)
  {
    double dx = iter.getRelativeX(1);
    dx += iter.getRelative(m_offsetsX[0]);
    dx += iter.getRelative(m_offsetsX[1]);
    dx += iter.getRelative(m_offsetsX[2]);
    dx += iter.getRelative(m_offsetsX[3]);

    dx -= iter.getRelativeX(-1);
    dx -= iter.getRelative(m_offsetsX[4]);
    dx -= iter.getRelative(m_offsetsX[5]);
    dx -= iter.getRelative(m_offsetsX[6]);
    dx -= iter.getRelative(m_offsetsX[7]);

    double dy = iter.getRelativeY(1);
    dy += iter.getRelative(m_offsetsY[0]);
    dy += iter.getRelative(m_offsetsY[1]);
    dy += iter.getRelative(m_offsetsY[2]);
    dy += iter.getRelative(m_offsetsY[3]);

    dy -= iter.getRelativeY(-1);
    dy -= iter.getRelative(m_offsetsY[4]);
    dy -= iter.getRelative(m_offsetsY[5]);
    dy -= iter.getRelative(m_offsetsY[6]);
    dy -= iter.getRelative(m_offsetsY[7]);

    double dz = iter.getRelativeZ(1);
    dz += iter.getRelative(m_offsetsZ[0]);
    dz += iter.getRelative(m_offsetsZ[1]);
    dz += iter.getRelative(m_offsetsZ[2]);
    dz += iter.getRelative(m_offsetsZ[3]);

    dz -= iter.getRelativeZ(-1);
    dz -= iter.getRelative(m_offsetsZ[4]);
    dz -= iter.getRelative(m_offsetsZ[5]);
    dz -= iter.getRelative(m_offsetsZ[6]);
    dz -= iter.getRelative(m_offsetsZ[7]);


    value.setNormal3f(vlNormal3f(dx, dy, dz));

    value.setDl(sqrt(dx*dx+dy*dy+dz*dz));

    return (true);
  };

protected:
  vlVarListAdmin m_varList;

  vlOffset m_offsetsX[8], m_offsetsY[8], m_offsetsZ[8];

};

#endif // _vlVoxelOpSobel_h

---