RTK  2.7.0
Reconstruction Toolkit
rtkZengBackProjectionImageFilter.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 rtkZengBackProjectionImageFilter_h
20 #define rtkZengBackProjectionImageFilter_h
21 
22 #include "rtkConfiguration.h"
24 #include "rtkMacro.h"
25 #include <itkPixelTraits.h>
26 
27 #include <itkMultiplyImageFilter.h>
28 #include <itkAddImageFilter.h>
30 #include <itkPasteImageFilter.h>
31 #include <itkResampleImageFilter.h>
32 #include <itkVector.h>
35 #include <itkExtractImageFilter.h>
39 
40 #include "rtkConstantImageSource.h"
41 
42 #include <itkVectorImage.h>
43 namespace rtk
44 {
45 
111 template <class TInputImage, class TOutputImage>
112 class ITK_TEMPLATE_EXPORT ZengBackProjectionImageFilter : public BackProjectionImageFilter<TInputImage, TOutputImage>
113 {
114 public:
117  using InputPixelType = typename TInputImage::PixelType;
118  using OutputPixelType = typename TOutputImage::PixelType;
119  using OutputImageRegionType = typename TOutputImage::RegionType;
123 
150 
152  static constexpr unsigned int InputImageDimension = TInputImage::ImageDimension;
153 
156 
158  itkNewMacro(Self);
159 
161  itkOverrideGetNameOfClassMacro(ZengBackProjectionImageFilter);
162 
164  itkGetMacro(SigmaZero, double);
165  itkSetMacro(SigmaZero, double);
167 
169  itkGetMacro(Alpha, double);
170  itkSetMacro(Alpha, double);
172 
173 
174 protected:
176  ~ZengBackProjectionImageFilter() override = default;
177 
179  void
180  GenerateInputRequestedRegion() override;
181 
182  void
183  GenerateOutputInformation() override;
184 
185  void
186  GenerateData() override;
187 
190  void
191  VerifyInputInformation() const override;
192 
208 
209 private:
210  ZengBackProjectionImageFilter(const Self &) = delete; // purposely not implemented
211  void
212  operator=(const Self &) = delete; // purposely not implemented
213 
214  double m_SigmaZero{ 1.5417233052142099 };
215  double m_Alpha{ 0.016241189545787734 };
216  VectorType m_VectorOrthogonalDetector{ 0. };
217  PointType m_centerVolume{ 0 };
218 };
219 
220 } // end namespace rtk
221 
222 #ifndef ITK_MANUAL_INSTANTIATION
223 # include "rtkZengBackProjectionImageFilter.hxx"
224 #endif
225 
226 #endif
typename ExtractImageFilterType::Pointer ExtractImageFilterPointerType
ConstantVolumeSourcePointerType m_ConstantVolumeSource
typename OuputCPUImageType::PointType PointType
typename PasteImageFilterType::Pointer PasteImageFilterPointerType
Generate an n-dimensional image with constant pixel values.
typename CustomUnaryFilterType::Pointer CustomUnaryFilterPointerType
typename TInputImage::PixelType InputPixelType
MultiplyImageFilterPointerType m_AttenuationMapMultiplyImageFilter
typename TOutputImage::PixelType OutputPixelType
#define itkSetMacro(name, type)
typename ConstantVolumeSourceType::Pointer ConstantVolumeSourcePointerType
typename RegionOfInterestFilterType::Pointer RegionOfInterestPointerType
DiscreteGaussianFilterPointeurType m_DiscreteGaussianFilter
typename MultiplyImageFilterType::Pointer MultiplyImageFilterPointerType
ChangeInformationPointerType m_AttenuationMapChangeInformation
typename ResampleImageFilterType::Pointer ResampleImageFilterPointerType
RegionOfInterestPointerType m_AttenuationMapRegionOfInterest
typename TOutputImage::RegionType OutputImageRegionType
typename DiscreteGaussianFilterType::Pointer DiscreteGaussianFilterPointeurType
MultiplyImageFilterPointerType m_MultiplyImageFilter
ResampleImageFilterPointerType m_ResampleImageFilter
typename ChangeInformationFilterType::Pointer ChangeInformationPointerType
ResampleImageFilterPointerType m_AttenuationMapResampleImageFilter
typename TransformType::Pointer TransformPointerType
typename AddImageFilterType::Pointer AddImageFilterPointerType