RTK  2.1.0
Reconstruction Toolkit
rtkIterativeFDKConeBeamReconstructionFilter.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  * http://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 rtkIterativeFDKConeBeamReconstructionFilter_h
20 #define rtkIterativeFDKConeBeamReconstructionFilter_h
21 
22 #include <itkMultiplyImageFilter.h>
23 #include <itkSubtractImageFilter.h>
26 
27 #include "rtkConstantImageSource.h"
33 
34 namespace rtk
35 {
36 
106 template <class TInputImage, class TOutputImage = TInputImage, class TFFTPrecision = double>
108  : public rtk::IterativeConeBeamReconstructionFilter<TInputImage, TOutputImage>
109 {
110 public:
111  ITK_DISALLOW_COPY_AND_ASSIGN(IterativeFDKConeBeamReconstructionFilter);
112 
118 
120  using InputImageType = TInputImage;
121  using OutputImageType = TOutputImage;
122 
123  using ForwardProjectionType = typename Superclass::ForwardProjectionType;
124  using BackProjectionType = typename Superclass::BackProjectionType;
125 
136 
138  itkNewMacro(Self);
139 
142 
144  itkGetModifiableObjectMacro(Geometry, ThreeDCircularProjectionGeometry);
145  itkSetObjectMacro(Geometry, ThreeDCircularProjectionGeometry);
147 
149  itkGetMacro(NumberOfIterations, unsigned int);
150  itkSetMacro(NumberOfIterations, unsigned int);
152 
154  itkGetMacro(Lambda, double);
155  itkSetMacro(Lambda, double);
157 
159  itkGetMacro(EnforcePositivity, bool);
160  itkSetMacro(EnforcePositivity, bool);
162 
164  void
165  SetBackProjectionFilter(BackProjectionType itkNotUsed(_arg)) override
166  {
167  itkExceptionMacro(<< "Backprojection cannot be changed");
168  }
169 
171  itkGetMacro(TruncationCorrection, double);
172  itkSetMacro(TruncationCorrection, double);
174 
176  itkGetMacro(HannCutFrequency, double);
177  itkSetMacro(HannCutFrequency, double);
179 
181  itkGetMacro(HannCutFrequencyY, double);
182  itkSetMacro(HannCutFrequencyY, double);
184 
186  itkGetMacro(ProjectionSubsetSize, unsigned int);
187  itkSetMacro(ProjectionSubsetSize, unsigned int);
189 
191  itkSetMacro(DisableDisplacedDetectorFilter, bool);
192  itkGetMacro(DisableDisplacedDetectorFilter, bool);
194 
195 protected:
197  ~IterativeFDKConeBeamReconstructionFilter() override = default;
198 
200  void
201  VerifyPreconditions() ITKv5_CONST override;
202 
203  void
204  GenerateInputRequestedRegion() override;
205 
206  void
207  GenerateOutputInformation() override;
208 
209  void
210  GenerateData() override;
211 
214  void
215  VerifyInputInformation() const override
216  {}
217 
219  typename Superclass::ForwardProjectionPointerType m_ForwardProjectionFilter;
229 
236 
237 private:
240 
242  unsigned int m_NumberOfIterations;
243 
245  double m_Lambda;
246 }; // end of class
247 
248 } // end namespace rtk
249 
250 #ifndef ITK_MANUAL_INSTANTIATION
251 # include "rtkIterativeFDKConeBeamReconstructionFilter.hxx"
252 #endif
253 
254 #endif
Generate an n-dimensional image with constant pixel values.
Weigting for displaced detectors.
Analytical projection of a BoxShape.
Projection geometry for a source and a 2-D flat panel.
#define itkSetMacro(name, type)
TOutputImage OutputImageType
Implements Feldkamp, David and Kress cone-beam reconstruction.
Mother class for cone beam reconstruction filters which need runtime selection of their forward and b...