RTK  2.6.0
Reconstruction Toolkit
rtkTotalVariationDenoiseSequenceImageFilter.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright RTK Consortium
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * https://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 
19 #ifndef rtkTotalVariationDenoiseSequenceImageFilter_h
20 #define rtkTotalVariationDenoiseSequenceImageFilter_h
21 
22 #include "rtkConstantImageSource.h"
23 
24 #include <itkExtractImageFilter.h>
25 #include <itkPasteImageFilter.h>
26 #include <itkCastImageFilter.h>
27 
28 #ifdef RTK_USE_CUDA
30 #else
32 #endif
33 
34 namespace rtk
35 {
78 template <typename TImageSequence>
80  : public itk::ImageToImageFilter<TImageSequence, TImageSequence>
81 {
82 public:
83  ITK_DISALLOW_COPY_AND_MOVE(TotalVariationDenoiseSequenceImageFilter);
84 
89 
91  itkNewMacro(Self);
92 
94 #ifdef itkOverrideGetNameOfClassMacro
95  itkOverrideGetNameOfClassMacro(TotalVariationDenoiseSequenceImageFilter);
96 #else
98 #endif
99 
100 
102  itkGetMacro(Gamma, double);
103  itkSetMacro(Gamma, double);
105 
106  itkGetMacro(NumberOfIterations, int);
107  itkSetMacro(NumberOfIterations, int);
108 
109  void
110  SetDimensionsProcessed(bool * arg);
111 
114 #ifdef RTK_USE_CUDA
115  typedef typename std::conditional<
116  std::is_same<TImageSequence, CPUImageSequenceType>::value,
117  itk::Image<typename TImageSequence::PixelType, TImageSequence::ImageDimension - 1>,
118  itk::CudaImage<typename TImageSequence::PixelType, TImageSequence::ImageDimension - 1>>::type ImageType;
119  typedef typename std::conditional<std::is_same<TImageSequence, CPUImageSequenceType>::value,
121  CudaTotalVariationDenoisingBPDQImageFilter>::type TVDenoisingFilterType;
122 #else
123  using ImageType = itk::Image<typename TImageSequence::PixelType, TImageSequence::ImageDimension - 1>;
125 #endif
130 
131 protected:
133  ~TotalVariationDenoiseSequenceImageFilter() override = default;
134 
136  void
137  GenerateData() override;
138 
139  void
140  GenerateOutputInformation() override;
141  void
142  GenerateInputRequestedRegion() override;
143 
150 
152  typename TImageSequence::RegionType m_ExtractAndPasteRegion;
153 
155  double m_Gamma{ 1. };
156  int m_NumberOfIterations{ 1 };
157  bool m_DimensionsProcessed[ImageType::ImageDimension];
158 };
159 } // namespace rtk
160 
161 
162 #ifndef ITK_MANUAL_INSTANTIATION
163 # include "rtkTotalVariationDenoiseSequenceImageFilter.hxx"
164 #endif
165 
166 #endif
Generate an n-dimensional image with constant pixel values.
Applies 3D total variation denoising to a 3D + time sequence of images.
Applies a total variation denoising, only alm_SingularValueThresholdFilterong the dimensions specifie...
#define itkSetMacro(name, type)
typename itk::Image< typename TImageSequence::PixelType, TImageSequence::ImageDimension > CPUImageSequenceType