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
139
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkShiftScaleImageFilter.h,v $
Language: C++
Date: $Date: 2006-03-31 14:31:04 $
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 __itkShiftScaleImageFilter_h
#define __itkShiftScaleImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkArray.h"
namespace itk
{
/** \class ShiftScaleImageFilter
* \brief Shift and scale the pixels in an image.
*
* ShiftScaleImageFilter shifts the input pixel by Shift (default 0.0)
* and then scales the pixel by Scale (default 1.0). All computattions
* are performed in the precison of the input pixel's RealType. Before
* assigning the computed value to the output pixel, the value is clamped
* at the NonpositiveMin and max of the pixel type.
* \ingroup IntensityImageFilters
*
*/
template <class TInputImage, class TOutputImage>
class ITK_EXPORT ShiftScaleImageFilter:
public ImageToImageFilter<TInputImage,TOutputImage>
{
public:
/** Standard class typedefs. */
typedef ShiftScaleImageFilter Self;
typedef ImageToImageFilter<TInputImage,TOutputImage> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Typedef to describe the output and input image region types. */
typedef typename TInputImage::RegionType InputImageRegionType;
typedef typename TOutputImage::RegionType OutputImageRegionType;
/** Typedef to describe the pointer to the input/output. */
typedef typename TInputImage::Pointer InputImagePointer;
typedef typename TOutputImage::Pointer OutputImagePointer;
/** Typedef to describe the type of pixel. */
typedef typename TInputImage::PixelType InputImagePixelType;
typedef typename TOutputImage::PixelType OutputImagePixelType;
/** Typedef to describe the output and input image index and size types. */
typedef typename TInputImage::IndexType InputImageIndexType;
typedef typename TInputImage::SizeType InputImageSizeType;
typedef typename TInputImage::OffsetType InputImageOffsetType;
typedef typename TOutputImage::IndexType OutputImageIndexType;
typedef typename TOutputImage::SizeType OutputImageSizeType;
typedef typename TOutputImage::OffsetType OutputImageOffsetType;
/** Type to use form computations. */
typedef typename NumericTraits<OutputImagePixelType>::RealType RealType;
/** Image related typedefs. */
itkStaticConstMacro(ImageDimension, unsigned int,
TInputImage::ImageDimension ) ;
/** Run-time type information (and related methods). */
itkTypeMacro(ShiftScaleImageFilter, ImageToImageFilter);
/** Set/Get the amount to Shift each Pixel. The shift is followed by a Scale. */
itkSetMacro(Shift,RealType);
itkGetMacro(Shift,RealType);
/** Set/Get the amount to Scale each Pixel. The Scale is applied after the Shift. */
itkSetMacro(Scale,RealType);
itkGetMacro(Scale,RealType);
/** Get the number of pixels that underflowed and overflowed. */
itkGetMacro(UnderflowCount,long);
itkGetMacro(OverflowCount,long);
#ifdef ITK_USE_CONCEPT_CHECKING
/** Begin concept checking */
itkConceptMacro(OutputHasNumericTraitsCheck,
(Concept::HasNumericTraits<OutputImagePixelType>));
itkConceptMacro(InputPlusRealTypeCheck,
(Concept::AdditiveOperators<InputImagePixelType, RealType, RealType>));
itkConceptMacro(RealTypeMultiplyOperatorCheck,
(Concept::MultiplyOperator<RealType>));
/** End concept checking */
#endif
protected:
ShiftScaleImageFilter();
~ShiftScaleImageFilter();
void PrintSelf(std::ostream& os, Indent indent) const;
/** Initialize some accumulators before the threads run. */
void BeforeThreadedGenerateData ();
/** Tally accumulated in threads. */
void AfterThreadedGenerateData ();
/** Multi-thread version GenerateData. */
void ThreadedGenerateData (const OutputImageRegionType&
outputRegionForThread,
int threadId) ;
private:
ShiftScaleImageFilter(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
RealType m_Shift;
RealType m_Scale;
long m_UnderflowCount;
long m_OverflowCount;
Array<long> m_ThreadUnderflow;
Array<long> m_ThreadOverflow;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkShiftScaleImageFilter.txx"
#endif
#endif
|
