RTK  2.0.1
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 
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 SetForwardProjectionFilter (ForwardProjectionType _arg) override;
165 
167  void SetBackProjectionFilter (BackProjectionType itkNotUsed(_arg)) override {itkExceptionMacro(<< "Backprojection cannot be changed");}
168 
170  itkGetMacro(TruncationCorrection, double);
171  itkSetMacro(TruncationCorrection, double);
173 
175  itkGetMacro(HannCutFrequency, double);
176  itkSetMacro(HannCutFrequency, double);
178 
180  itkGetMacro(HannCutFrequencyY, double);
181  itkSetMacro(HannCutFrequencyY, double);
183 
185  itkGetMacro(ProjectionSubsetSize, unsigned int);
186  itkSetMacro(ProjectionSubsetSize, unsigned int);
188 
190  itkSetMacro(DisableDisplacedDetectorFilter, bool)
191  itkGetMacro(DisableDisplacedDetectorFilter, bool)
192 
193 protected:
195  ~IterativeFDKConeBeamReconstructionFilter() override = default;
196 
197  void GenerateInputRequestedRegion() override;
198 
199  void GenerateOutputInformation() override;
200 
201  void GenerateData() override;
202 
205 #if ITK_VERSION_MAJOR<5
206  void VerifyInputInformation() override {}
207 #else
208  void VerifyInputInformation() const override {}
209 #endif
210 
211 
223 
230 
231 private:
234 
236  unsigned int m_NumberOfIterations;
237 
239  double m_Lambda;
240 }; // end of class
241 
242 } // end namespace rtk
243 
244 #ifndef ITK_MANUAL_INSTANTIATION
245 #include "rtkIterativeFDKConeBeamReconstructionFilter.hxx"
246 #endif
247 
248 #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.
Implements Feldkamp, David and Kress cone-beam reconstruction.
Mother class for cone beam reconstruction filters which need runtime selection of their forward and b...
#define itkSetMacro(name, type)