RTK  2.6.0
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  * 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 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>
98 class ITK_TEMPLATE_EXPORT ADMMWaveletsConjugateGradientOperator : public ConjugateGradientOperator<TOutputImage>
99 {
100 public:
101  ITK_DISALLOW_COPY_AND_MOVE(ADMMWaveletsConjugateGradientOperator);
102 
107 
109  itkNewMacro(Self);
110 
112 #ifdef itkOverrideGetNameOfClassMacro
113  itkOverrideGetNameOfClassMacro(ADMMWaveletsConjugateGradientOperator);
114 #else
116 #endif
117 
118 
120  using BackProjectionFilterPointer = typename BackProjectionFilterType::Pointer;
121 
123  using ForwardProjectionFilterPointer = typename ForwardProjectionFilterType::Pointer;
124 
126  using AddFilterType = itk::AddImageFilter<TOutputImage>;
127 
129 
131  void
132  SetBackProjectionFilter(const BackProjectionFilterPointer _arg);
133 
135  void
136  SetForwardProjectionFilter(const ForwardProjectionFilterPointer _arg);
137 
139  void
140  SetGeometry(const ThreeDCircularProjectionGeometry::Pointer _arg);
141 
143  itkSetMacro(Beta, float);
144 
146  itkSetMacro(DisableDisplacedDetectorFilter, bool);
147  itkGetMacro(DisableDisplacedDetectorFilter, bool);
149 
150 protected:
152  ~ADMMWaveletsConjugateGradientOperator() override = default;
153 
155  void
156  GenerateData() override;
157 
161 
162  typename AddFilterType::Pointer m_AddFilter;
167 
168  float m_Beta;
170 
174  void
175  VerifyInputInformation() const override
176  {}
177 
180  void
181  GenerateInputRequestedRegion() override;
182  void
183  GenerateOutputInformation() override;
184 };
185 } // namespace rtk
187 
188 
189 #ifndef ITK_MANUAL_INSTANTIATION
190 # include "rtkADMMWaveletsConjugateGradientOperator.hxx"
191 #endif
192 
193 #endif
Weigting for displaced detectors.
typename ForwardProjectionFilterType::Pointer ForwardProjectionFilterPointer
#define itkSetMacro(name, type)
Implements the operator A used in the conjugate gradient step of ADMM reconstruction with wavelets re...
typename BackProjectionFilterType::Pointer BackProjectionFilterPointer