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
|
/*=========================================================================
*
* Copyright NumFOCUS
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* https://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#ifndef itkMeanProjectionImageFilter_h
#define itkMeanProjectionImageFilter_h
#include "itkProjectionImageFilter.h"
#include "itkConceptChecking.h"
namespace itk
{
/** \class MeanProjectionImageFilter
* \brief Mean projection
*
* This class was contributed to the Insight Journal by Gaetan Lehmann.
* The original paper can be found at
* https://www.insight-journal.org/browse/publication/71
*
* \author Gaetan Lehmann. Biologie du Developpement et de la Reproduction,
* INRA de Jouy-en-Josas, France.
*
* \sa ProjectionImageFilter
* \sa MedianProjectionImageFilter
* \sa MinimumProjectionImageFilter
* \sa StandardDeviationProjectionImageFilter
* \sa SumProjectionImageFilter
* \sa BinaryProjectionImageFilter
* \sa MaximumProjectionImageFilter
* \ingroup ITKImageStatistics
*/
namespace Functor
{
template <typename TInputPixel, typename TAccumulate>
class MeanAccumulator
{
public:
using RealType = typename NumericTraits<TInputPixel>::RealType;
MeanAccumulator(SizeValueType size) { m_Size = size; }
~MeanAccumulator() { m_Size = SizeValueType{}; }
inline void
Initialize()
{
m_Sum = TAccumulate{};
}
inline void
operator()(const TInputPixel & input)
{
m_Sum = m_Sum + input;
}
inline RealType
GetValue()
{
return ((RealType)m_Sum) / m_Size;
}
TAccumulate m_Sum;
SizeValueType m_Size;
};
} // namespace Functor
template <typename TInputImage,
typename TOutputImage,
typename TAccumulate = typename NumericTraits<typename TOutputImage::PixelType>::AccumulateType>
class MeanProjectionImageFilter
: public ProjectionImageFilter<TInputImage,
TOutputImage,
Functor::MeanAccumulator<typename TInputImage::PixelType, TAccumulate>>
{
public:
ITK_DISALLOW_COPY_AND_MOVE(MeanProjectionImageFilter);
using Self = MeanProjectionImageFilter;
using Superclass = ProjectionImageFilter<TInputImage,
TOutputImage,
Functor::MeanAccumulator<typename TInputImage::PixelType, TAccumulate>>;
using Pointer = SmartPointer<Self>;
using ConstPointer = SmartPointer<const Self>;
using InputImageType = TInputImage;
using InputPixelType = typename InputImageType::PixelType;
using OutputImageType = TOutputImage;
using OutputPixelType = typename OutputImageType::PixelType;
/** \see LightObject::GetNameOfClass() */
itkOverrideGetNameOfClassMacro(MeanProjectionImageFilter);
/** Method for creation through the object factory. */
itkNewMacro(Self);
#ifdef ITK_USE_CONCEPT_CHECKING
// Begin concept checking
itkConceptMacro(InputPixelToOutputPixelTypeGreaterAdditiveOperatorCheck,
(Concept::AdditiveOperators<OutputPixelType, InputPixelType, OutputPixelType>));
itkConceptMacro(InputHasNumericTraitsCheck, (Concept::HasNumericTraits<InputPixelType>));
// End concept checking
#endif
protected:
MeanProjectionImageFilter() = default;
~MeanProjectionImageFilter() override = default;
}; // end MeanProjectionImageFilter
} // end namespace itk
#endif
|