1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkMaskNegatedImageFilter.h,v $
Language: C++
Date: $Date: 2007-09-27 11:36:40 $
Version: $Revision: 1.7 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef __itkMaskNegatedImageFilter_h
#define __itkMaskNegatedImageFilter_h
#include "itkBinaryFunctorImageFilter.h"
#include "itkNumericTraits.h"
namespace itk
{
/** \class MaskNegatedImageFilter
* \brief Implements an operator for pixel-wise masking of the input
* image with the negative of a mask.
*
* This class is parametrized over the types of the
* input image type, the mask image type and the type of the output image.
* Numeric conversions (castings) are done by the C++ defaults.
*
* The pixel type of the input 2 image must have a valid defintion of the
* operator != with zero . This condition is required because internally this
* filter will perform the operation
*
* if pixel_from_mask_image != 0
* pixel_output_image = 0
* else
* pixel_output_image = pixel_input_image
*
* The pixel from the input 1 is cast to the pixel type of the output image.
*
* Note that the input and the mask images must be of the same size.
*
* \warning Any pixel value other than 0 will not be masked out.
*
* \sa MaskImageFilter
* \ingroup IntensityImageFilters Multithreaded
*/
namespace Functor {
template< class TInput, class TMask, class TOutput=TInput >
class MaskNegatedInput
{
public:
typedef typename NumericTraits< TInput >::AccumulateType AccumulatorType;
MaskNegatedInput() {};
~MaskNegatedInput() {};
bool operator!=( const MaskNegatedInput & ) const
{
return false;
}
bool operator==( const MaskNegatedInput & other ) const
{
return !(*this != other);
}
inline TOutput operator()( const TInput & A, const TMask & B)
{
if (B != NumericTraits< TMask >::Zero )
{
return NumericTraits< TOutput >::Zero;
}
else
{
return static_cast<TOutput>( A );
}
}
};
}
template <class TInputImage, class TMaskImage, class TOutputImage=TInputImage>
class ITK_EXPORT MaskNegatedImageFilter :
public
BinaryFunctorImageFilter<TInputImage,TMaskImage,TOutputImage,
Functor::MaskNegatedInput<
typename TInputImage::PixelType,
typename TMaskImage::PixelType,
typename TOutputImage::PixelType> >
{
public:
/** Standard class typedefs. */
typedef MaskNegatedImageFilter Self;
typedef BinaryFunctorImageFilter<TInputImage,TMaskImage,TOutputImage,
Functor::MaskNegatedInput<
typename TInputImage::PixelType,
typename TMaskImage::PixelType,
typename TOutputImage::PixelType>
> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Runtime information support. */
itkTypeMacro(MaskNegatedImageFilter,
BinaryFunctorImageFilter);
#ifdef ITK_USE_CONCEPT_CHECKING
/** Begin concept checking */
itkConceptMacro(MaskEqualityComparableCheck,
(Concept::EqualityComparable<typename TMaskImage::PixelType>));
itkConceptMacro(InputConvertibleToOutputCheck,
(Concept::Convertible<typename TInputImage::PixelType,
typename TOutputImage::PixelType>));
/** End concept checking */
#endif
protected:
MaskNegatedImageFilter() {}
virtual ~MaskNegatedImageFilter() {}
private:
MaskNegatedImageFilter(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
};
} // end namespace itk
#endif
|