mitkSubtractImageFilter.h

/*=========================================================================

  Program:   3DMed
  Date:      $Date: 2014-02-25 18:30:00 +0800 $
  Version:   $Version: 4.6.0 $
  Copyright: MIPG, Institute of Automation, Chinese Academy of Sciences

=========================================================================*/


#ifndef __mitkSubtractImageFilter_h
#define __mitkSubtractImageFilter_h

#include "mitkVolumeToVolumeFilter.h"
#include "mitkInterpolateFilter.h"

#define MITK_SUBSTRACT_MODE_USERDEFINED     0
#define MITK_SUBSTRACT_MODE_SUB             1
#define MITK_SUBSTRACT_MODE_ADD             2
#define MITK_SUBSTRACT_MODE_SUBCOLOR        3
#define MITK_SUBSTRACT_MODE_ADDCOLOR        4
#define MITK_SUBSTRACT_MODE_CHECKERBOARD    5
#define MITK_SUBSTRACT_MODE_MIN             6
#define MITK_SUBSTRACT_MODE_MAX             7
#define MITK_SUBSTRACT_MODE_AVERAGE         8
#define MITK_SUBSTRACT_MODE_MULTIPLY        9
#define MITK_SUBSTRACT_MODE_EXPLOGMULTI     10  
#define MITK_SUBSTRACT_MODE_DIVIDE          11
#define MITK_SUBSTRACT_MODE_AND             12
#define MITK_SUBSTRACT_MODE_OR              13
#define MITK_SUBSTRACT_MODE_XOR             14

//for single input
#define MITK_SUBSTRACT_MODE_LOG_10          101
#define MITK_SUBSTRACT_MODE_LOG_E           102
#define MITK_SUBSTRACT_MODE_SQUAREROOT      103
#define MITK_SUBSTRACT_MODE_SQUARED         104
#define MITK_SUBSTRACT_MODE_ABS             105
#define MITK_SUBSTRACT_MODE_COMPLEMENT      106
#define MITK_SUBSTRACT_MODE_NOT             107


class MITK_REGISTRATION_API mitkSubtractImageFilter : public mitkVolumeToVolumeFilter
{
public:
    MITK_TYPE(mitkSubtractImageFilter, mitkVolumeToVolumeFilter)    
    virtual void PrintSelf(ostream &os);

    mitkSubtractImageFilter(int mode);

    mitkSubtractImageFilter();

    void SetInput1(mitkVolume *inData) { m_InData = inData; }

    void SetInput2(mitkVolume *inData) { m_InData2 = inData; }

    void SetInterpolator(mitkInterpolateFilter* filter) { m_Interpolator = filter; }

    mitkInterpolateFilter* GetInterpolator();

    void SetSubtractMode(int mode) { m_SubtractMode = mode; }

    void SetChannelMixMode(int mode) { m_ChannelMixMode = mode; }

    void SetAdjustOutputFlag(bool f) { m_FlagAdjustOutput = f; }

    void SetOutputDatatype(int dataType) {m_OutputDataType = dataType;}

    void SetCheckImageDim(int row, int col, int dep=1);

    void SetMaxOutputValue(ScalarPixelType max) { m_MaxOutputValue = max; }

    void SetMinOutputValue(ScalarPixelType min) { m_MinOutputValue = min; }

    void SetSubtractKernel(void (*kernel)(VectorPixelType& ip,VectorPixelType& op))
    {
        SubtractKernel1 = kernel;
    }

    void SetSubtractKernel(void (*kernel)(VectorPixelType& ip1,VectorPixelType& ip2,VectorPixelType& op))
    {
        SubtractKernel2 = kernel;
    }

protected:
    virtual ~mitkSubtractImageFilter();
    virtual bool Execute();
    bool Update();
    void GenerateCheckImage();
    bool BindSubtractKernel();
    void AdjustOutput();
    void _init();

    void (*SubtractKernel1)(VectorPixelType&,VectorPixelType&);
    void (*SubtractKernel2)(VectorPixelType&,VectorPixelType&,VectorPixelType&);

    mitkRCPtr<mitkVolume>               m_MaskImage;
    mitkRCPtr<mitkVolume>               m_InData2;
    mitkRCPtr<mitkInterpolateFilter>    m_Interpolator;

    int                                 m_SubtractMode;
    //    -1: fill n ch with ch[0] 
    // 0~n-1: swap ch[0] to ch[i]
    int                                 m_ChannelMixMode;
    int                                 m_OutputDataType;

    float                               m_Spacings[3];
    int                                 m_Dimensions[3];
    unsigned int                        m_CheckImageDim[3];
    unsigned int                        m_NumberOfChannels;
    unsigned int                        m_ImageDimension;

    ScalarPixelType                     m_MinOutputValue;
    ScalarPixelType                     m_MaxOutputValue;

    bool                                m_FlagNeedInterpolation;
    bool                                m_FlagAdjustOutput;

private:
    mitkSubtractImageFilter(const mitkSubtractImageFilter&);
    void operator = (const mitkSubtractImageFilter&);

};


//#define DEFINED_mitkSubtractImageFilter



#endif

---