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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkTwoOutputExampleImageFilter.h,v $
Language: C++
Date: $Date: 2007-08-20 07:50:04 $
Version: $Revision: 1.11 $
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 __itkTwoOutputExampleImageFilter_h
#define __itkTwoOutputExampleImageFilter_h
#include "itkImageToImageFilter.h"
namespace itk
{
/** \class TwoOutputExampleImageFilter
* \brief Example of a filter that produce two outputs.
*
* TwoOutputExampleImageFilter sets image values to a user-specified "outside"
* value (by default, "black") if the image values are below, above, or
* between simple threshold values. The filter can produce two outputs,
* one the inverse of the other. (GetOutput() returns an image whose
* pixels satisfy the threshold values and are passed to the output
* unchanged (and those that don't are marked with the outside user-value);
* GetInverseOutput() returns an image in which pixels satisfying the
* threshold are marked "outside", and the other pixel values are passed
* through.)
*
* The pixels must support the operators >= and <=.
*
* \ingroup IntensityImageFilters Multithreaded
*/
template <class TImage>
class ITK_EXPORT TwoOutputExampleImageFilter:
public ImageToImageFilter<TImage, TImage>
{
public:
/** Standard class typedefs. */
typedef TwoOutputExampleImageFilter Self;
typedef ImageToImageFilter<TImage,TImage> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Typedef to describe the type of pixel. */
typedef typename TImage::PixelType PixelType;
/** Run-time type information (and related methods). */
itkTypeMacro(TwoOutputExampleImageFilter, ImageToImageFilter);
/** Set the "outside" pixel value. The default value
* NumericTraits<PixelType>::Zero. */
itkSetMacro(OutsideValue,PixelType);
/** Get the "outside" pixel value. */
itkGetMacro(OutsideValue,PixelType);
/** The values greater than or equal to the value are set to OutsideValue. */
void ThresholdAbove(PixelType &thresh);
/** The values less than or equal to the value are set to OutsideValue. */
void ThresholdBelow(PixelType &thresh);
/** The values outside the range are set to OutsideValue. */
void ThresholdOutside(PixelType &lower, PixelType &upper);
typedef typename Superclass::InputImageConstPointer InputImageConstPointer;
/** Some typedefs to handle the second output. */
typedef TImage OutputImageType;
typedef typename OutputImageType::Pointer OutputImagePointer;
typedef typename OutputImageType::RegionType OutputImageRegionType;
typedef typename OutputImageType::PixelType OutputImagePixelType;
/** Get the image output of this process object. */
OutputImagePointer GetInverseOutput()
{ return static_cast<TImage *>(this->ProcessObject::GetOutput(1)); }
/** Set the image output of this process object. */
void SetInverseOutput(OutputImageType *output)
{ this->SetNthOutput(1, output); };
#ifdef ITK_USE_CONCEPT_CHECKING
/** Begin concept checking */
itkConceptMacro(ComparableCheck,
(Concept::Comparable<PixelType>));
itkConceptMacro(OStreamWritableCheck,
(Concept::OStreamWritable<PixelType>));
/** Begin concept checking */
#endif
protected:
TwoOutputExampleImageFilter();
~TwoOutputExampleImageFilter() {};
void PrintSelf(std::ostream& os, Indent indent) const;
/** TwoOutputExampleImageFilter can be implemented as a multithreaded filter.
* Therefore, this implementation provides a ThreadedGenerateData() routine
* which is called for each processing thread. The output image data is
* allocated automatically by the superclass prior to calling
* ThreadedGenerateData(). ThreadedGenerateData can only write to the
* portion of the output image specified by the parameter
* "outputRegionForThread"
*
* \sa ImageToImageFilter::ThreadedGenerateData(),
* ImageToImageFilter::GenerateData() */
void ThreadedGenerateData(const OutputImageRegionType& outputRegionForThread,
int threadId );
private:
TwoOutputExampleImageFilter(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
PixelType m_OutsideValue;
PixelType m_Lower;
PixelType m_Upper;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkTwoOutputExampleImageFilter.txx"
#endif
#endif
|
