RTK  2.0.1
Reconstruction Toolkit
rtkADMMWaveletsConjugateGradientOperator.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 rtkADMMWaveletsConjugateGradientOperator_h
20 #define rtkADMMWaveletsConjugateGradientOperator_h
21 
22 #include <itkMultiplyImageFilter.h>
23 #include <itkAddImageFilter.h>
24 
29 
31 
32 namespace rtk
33 {
34 
97 template< typename TOutputImage >
99 {
100 public:
101  ITK_DISALLOW_COPY_AND_ASSIGN(ADMMWaveletsConjugateGradientOperator);
102 
107 
109  itkNewMacro(Self)
110 
111 
112  itkTypeMacro(rtkADMMWaveletsConjugateGradientOperator, ConjugateGradientOperator)
113 
114  using BackProjectionFilterType = rtk::BackProjectionImageFilter< TOutputImage, TOutputImage >;
116 
117  using ForwardProjectionFilterType = rtk::ForwardProjectionImageFilter< TOutputImage, TOutputImage >;
119 
121  using AddFilterType = itk::AddImageFilter<TOutputImage>;
122 
124 
127 
130 
132  void SetGeometry(const ThreeDCircularProjectionGeometry::Pointer _arg);
133 
135  itkSetMacro(Beta, float)
136 
138  itkSetMacro(DisableDisplacedDetectorFilter, bool)
139  itkGetMacro(DisableDisplacedDetectorFilter, bool)
140 
141 protected:
143  ~ADMMWaveletsConjugateGradientOperator() override = default;
144 
146  void GenerateData() override;
147 
151 
152  typename AddFilterType::Pointer m_AddFilter;
157 
158  float m_Beta;
160 
164 #if ITK_VERSION_MAJOR<5
165  void VerifyInputInformation() override {}
166 #else
167  void VerifyInputInformation() const override {}
168 #endif
169 
170 
173  void GenerateInputRequestedRegion() override;
174  void GenerateOutputInformation() override;
175 };
176 } //namespace RTK
178 
179 
180 #ifndef ITK_MANUAL_INSTANTIATION
181 #include "rtkADMMWaveletsConjugateGradientOperator.hxx"
182 #endif
183 
184 #endif
Weigting for displaced detectors.
typename ForwardProjectionFilterType::Pointer ForwardProjectionFilterPointer
Projection geometry for a source and a 2-D flat panel.
Implements the operator A used in the conjugate gradient step of ADMM reconstruction with wavelets re...
void SetForwardProjectionFilter(const ForwardProjectionFilterPointer _arg)
void SetGeometry(const ThreeDCircularProjectionGeometry::Pointer _arg)
typename BackProjectionFilterType::Pointer BackProjectionFilterPointer
#define itkSetMacro(name, type)
void SetBackProjectionFilter(const BackProjectionFilterPointer _arg)