Filter Library

Camera

Interface Physics

ScanFilterMath.cpp

/*
 * ScanFilterMath.h - math filter.
 *
 * This file is part of the Camera Filter Library.
 * Computer Aided Measurement Environment for Realtime Atomic imaging (Camera)
 *
 * Copyright (C) 2004-2005, Leiden Probe Microscopy.
 * Copyright (C) 2004-2005, Universiteit Leiden.
 *
 * Authors: M. Seynen (original), Martin J. Moene
 *
 * $Id: ScanFilterMath.cpp 400 2006-09-28 13:34:02Z moene $
 */

/*
 * configuration defines:
 */

// currently nothing

/*
 * end of configuration defines.
 */

#include "stdafx.h"                     // for common (pre-compiled) headers
//#include <cfl/stdafx.h>                 // for common (pre-compiled) headers
#include <cfl/resource.h>               // for dialog resources

#include <cfl/FilterDlg_Math.h>         // for class CFilterDlg_Math
#include <cfl/ScanFilterMath.h>         // for class CScanFilterMath

#pragma warning( disable: 4018 4100 4116 4146 4244 4663)
#include <cil/Math.h>                   // for lpm::min()
#pragma warning( default: 4018 4100 4116 4146 4244 4663)

#include <algorithm>                    // for std::transform()

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

#define CFL_SCHEMAVERSION_FILTERMATH 0          // permanent storage schema version
                                                // this filter's name
LPCTSTR CScanFilterMath::m_lpcsFilterName = _T("Math");
                                        // this filter's short name
LPCTSTR CScanFilterMath::m_lpcsShortFilterName = CScanFilterMath::m_lpcsFilterName;
                                                // serialization
IMPLEMENT_SERIAL( CScanFilterMath, CScanFilter, CFL_SCHEMAVERSION_FILTERMATH )

/*
 * default math operation: a - b.
 */
const CScanFilterMath::MathOperationType CScanFilterMath::def_MathOperation = MATH_SUBAB;

/*
 * constructor.
 */
CScanFilterMath::CScanFilterMath() :
   m_lpsbIn2( NULL ),
   m_MathOperation( def_MathOperation )
{
   m_csFilterName = m_lpcsFilterName;

   ReadFilterSettings();
}

/*
 * destructor.
 */
CScanFilterMath::~CScanFilterMath()
{
//   Q_LOG( _T("CScanFilterMath::~CScanFilterMath()") );

   ; // do nothing
}

/*
 * store or retrieve the object's settings;
 * see also CScanFilterNull::Serialize().
 */
void CScanFilterMath::Serialize(CArchive& ar)
{
   CScanFilter::Serialize( ar );

   if ( ar.IsStoring() )
   {
      ar << static_cast< int >( m_MathOperation );
   }
   else
   {
      int op;
      ar >> op; m_MathOperation = static_cast< MathOperationType >( op );
   }
};

/*
 * \brief configure filter with settings as provided by the application on the
 * filterlist window (e.g. via registry); see also CScanFilterNull::ReadFilterSettings().
 */
void CScanFilterMath::ReadFilterSettings()
{
//   Q_LOG( _T("CScanFilterMath::ReadFilterSettings()") );

   CWinApp* pApp = AfxGetApp();

   if ( Q_INVALID( NULL == pApp ) )
      return ;

   /*
    * parameter string:
    */
   SetParameters( pApp->GetProfileString( gCflRegistrySubkey, m_lpcsFilterName, _T( "0" ) ) );
}

/*
 * \brief save filter settings e.g. for use with the filterlist window (e.g. via registry);
 * see also ReadFilterSettings().
 */
void CScanFilterMath::WriteFilterSettings() const
{
//   Q_LOG( _T("CScanFilterMath::WriteFilterSettings()") );

   CWinApp* pApp = AfxGetApp();

   if ( Q_INVALID( NULL == pApp ) )
      return ;

   /*
    * parameter string:
    */
   LPCTSTR pStr = GetParameters();

   if ( pStr )
   {
      pApp->WriteProfileString( gCflRegistrySubkey, m_lpcsFilterName, pStr );
      delete const_cast<LPTSTR>( pStr );
   }
}

/*
 * provide current filter parameter(s) for filter-scripting capture
 * in the main application;
 * see also CScanFilterNull::GetParameters().
 */
LPCTSTR CScanFilterMath::GetParameters() const
{
//   Q_LOG( _T("CScanFilterMath::GetParameters()") );

   /*
    * format parameters and return a c-string on the heap:
    */
   CString csParameters;

   csParameters.Format( _T( "%d" ), m_MathOperation );

   return strcpy( new TCHAR[ csParameters.GetLength() + 1 ], csParameters );
}

/*
 * set parameters for filter-script execution of filter;
 * see also CScanFilterNull::SetParameters().
 */
BOOL CScanFilterMath::SetParameters( LPCTSTR lpParameters )
{
//   Q_LOG( _T("CScanFilterMath::SetParameters()") );

   if ( Q_INVALID( NULL == lpParameters && "CScanFilterMath (SetParameters): parameters expected, none provided (NULL).") )
   {
      return FALSE;
   }

   /*
    * set defaults, read parameter string with factor and optional offset:
    */
   m_MathOperation = def_MathOperation;

   int nfields = _stscanf( lpParameters, _T( "%d" ), &m_MathOperation );

   if ( Q_INVALID( 1 > nfields && "CScanFilterMath (SetParameters): math operation expected, got none." ) )
   {
      return FALSE;
   }

   return TRUE;
}

/**
 * the number of input buffers.
 */
DWORD CScanFilterMath::GetNrOfInBuffers() const
{
//   Q_LOG( _T("CScanFilterMath::SetParameters()") );

   return 2;
}

/**
 * return the requested input buffer.
 */
CScanBaseBuffer *CScanFilterMath::GetInputBuffer( UINT nNr )
{
//   Q_LOG( _T("CScanFilterMath::GetInputBuffer()") );

   Q_ASSERT( 0 <= nNr && nNr < GetNrOfInBuffers() && "ensure valid input buffer nr" );

   if ( 1 == nNr )
   {
      return m_lpsbIn2;
   }

   return m_lpsbIn;
}

/**
 * set the given input buffer.
 */
void CScanFilterMath::SetInputBuffer( UINT nNr, CScanBaseBuffer *lpsbIn )
{
//   Q_LOG( _T("CScanFilterMath::SetInputBuffer()") );

   Q_ASSERT( 0 <= nNr && nNr < GetNrOfInBuffers() && "ensure valid input buffer nr" );
   Q_ASSERT( NULL != lpsbIn && "ensure valid input buffer" );

   if ( 0 == nNr )
   {
      m_lpsbIn = lpsbIn;
   }
   else if ( 1 == nNr )
   {
      m_lpsbIn2 = lpsbIn;
   }
}

/*
 * create filter dialog and pre Apply() filter to view result;
 * see also CScanFilterNull::RunModeless().
 */
BOOL CScanFilterMath::RunModeless( CWnd* pParentWnd, CDocument* pDoc )
{
//   Q_LOG( _T("CScanFilterMath::RunModeless()") );

   CFilterDlg_MathPtr pDlg = new CFilterDlg_Math();

   if ( Q_INVALID( NULL == pDlg ) )
      return FALSE;

   pDlg->SetParameters( this, pDoc );   // (void) method

   /*
    * Create returns non-zero if dialog was created and initialized succesfully.
    * Note: no pDlg->DestroyWindow() must be done.
    */
   pDlg->Create( IDD_FILTERDLG_MATH, pParentWnd );

   m_pDlg = pDlg;

   BOOL bRet = Apply();

   if ( bRet )
   {
      pDlg->UpdateView();
   }

   return bRet;
}

/*
 * check parameters, take care of a properly sized output buffer,
 * set its name and copy filter parameters and process;
 * see also CScanFilterNull::ApplyCore().
 */
BOOL CScanFilterMath::Apply()
{
//   Q_LOG( _T("CScanFilterMath::Apply()") );

   if ( Q_INVALID( NULL == m_lpsbIn && "ensure first input buffer") )
      return FALSE;

   if ( Q_INVALID( NULL == m_lpsbIn2 && "ensure second input buffer" ) )
      return FALSE;

   if ( Q_INVALID( NULL == m_lpsbOut && "ensure output buffer" ) )
      return FALSE;

   if ( Q_INVALID( !m_lpsbIn->IsCompleteFrame() && "ensure complete frame" ) )
      return FALSE;

   /*
    * copy the filter settings (not the data):
    */
   Q_RETURN( m_lpsbOut->CreateOutputBufferFor( *m_lpsbIn ) );

   m_lpsbOut->m_csBufferName.Format( _T("%s - %s"), m_lpsbIn->m_csBufferName, m_lpcsShortFilterName );

   if ( IsInteractive() )
   {
      m_pDlg->SetDlgItemText( IDC_BUFFERINNAME2, m_lpsbIn2->m_csBufferName );
      m_pDlg->SetDlgItemText( IDC_BUFFEROUTNAME, m_lpsbOut->m_csBufferName );
   }

   /*
    * limit buffer access to the smallest one: this also handles backtrace removal:
    */
   SizeElementType size = lpm::min( m_lpsbIn->GetFrameSize(), m_lpsbIn2->GetFrameSize(), m_lpsbOut->GetFrameSize() );

   /*
    * process first part (L2R):
    */
   Process(
      m_lpsbIn->Data(),
      m_lpsbIn->Data() + size,
      m_lpsbIn2->Data(),
      m_lpsbOut->Data()
   );

   /*
    * if we have a double buffer, L2R and R2L, process that part too.
    */
   if ( m_lpsbIn->IsBidirectionalScan() && m_lpsbIn2->IsBidirectionalScan() && m_lpsbOut->IsBidirectionalScan() )
   {
      Process(
         m_lpsbIn->Data()  + m_lpsbIn->GetFrameSize(),
         m_lpsbIn->Data()  + m_lpsbIn->GetFrameSize() + size,
         m_lpsbIn2->Data() + m_lpsbIn2->GetFrameSize(),
         m_lpsbOut->Data() + m_lpsbOut->GetFrameSize()
      );
   }

   /*
    * recalculate all temporary information:
    */
   m_lpsbOut->m_dwParameterMask = 0;

   return TRUE;
}

/**
 * anonymous namespace to hide class AddOp, SubOp, InvSubOp from other modules:
 */
namespace
{
   /**
    * binary functor to perform calculation on a + b in std::transform.
    */
   template <typename T>
   class AddOp : public std::binary_function<T, T, T>
   {
   public:
      /**
       * visit elements; return a + b.
       */
      inline T operator()( const T a, const T b ) const { return static_cast<T>( a + b ); }
   };

   /**
    * unary functor to perform calculation a - b in std::transform.
    */
   template <typename T>
   class SubOp : public std::binary_function<T, T, T>
   {
   public:
      /**
       * visit elements; return a - b.
       */
      inline T operator()( const T a, const T b ) const { return static_cast<T>( a - b ); }
   };

   /**
    * unary functor to perform calculation b - a in std::transform.
    */
   template <typename T>
   class InvSubOp : public std::binary_function<T, T, T>
   {
   public:
      /**
       * visit elements; return b - a.
       */
      inline T operator()( const T a, const T b ) const { return static_cast<T>( b - a ); }
   };

   /**
    * binary functor to perform calculation on a * b in std::transform.
    */
   template <typename T>
   class MulOp : public std::binary_function<T, T, T>
   {
   public:
      /**
       * visit elements; return a + b.
       */
      inline T operator()( const T a, const T b ) const { return static_cast<T>( a * b ); }
   };

}  // end anonymous namespace

/*
 * apply the operation to the buffers.
 */
void CScanFilterMath::Process( ConstPointer in1, ConstPointer end1, ConstPointer in2, Pointer out )
{
   switch( m_MathOperation )
   {
      case MATH_ADD:
         std::transform( in1, end1, in2, out, AddOp< ValueType >() );
         break;

      case MATH_SUBAB:
         std::transform( in1, end1, in2, out, SubOp< ValueType >() );
         break;

      case MATH_SUBBA:
         std::transform( in1, end1, in2, out, InvSubOp< ValueType >() );
         break;

      case MATH_MUL:
         std::transform( in1, end1, in2, out, MulOp< ValueType >() );
         break;

   };
}

/*
 * end of file
 */

---