RTK  1.4.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  * 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:
105 
107  itkNewMacro(Self)
108 
109 
110  itkTypeMacro(rtkADMMWaveletsConjugateGradientOperator, ConjugateGradientOperator)
111 
112  typedef rtk::BackProjectionImageFilter< TOutputImage, TOutputImage > BackProjectionFilterType;
113  typedef typename BackProjectionFilterType::Pointer BackProjectionFilterPointer;
114 
115  typedef rtk::ForwardProjectionImageFilter< TOutputImage, TOutputImage > ForwardProjectionFilterType;
116  typedef typename ForwardProjectionFilterType::Pointer ForwardProjectionFilterPointer;
117 
119  typedef itk::AddImageFilter<TOutputImage> AddFilterType;
120 
122 
124  void SetBackProjectionFilter (const BackProjectionFilterPointer _arg);
125 
127  void SetForwardProjectionFilter (const ForwardProjectionFilterPointer _arg);
128 
130  void SetGeometry(const ThreeDCircularProjectionGeometry::Pointer _arg);
131 
133  itkSetMacro(Beta, float)
134 
136  itkSetMacro(DisableDisplacedDetectorFilter, bool)
137  itkGetMacro(DisableDisplacedDetectorFilter, bool)
138 
139 protected:
141  virtual ~ADMMWaveletsConjugateGradientOperator() ITK_OVERRIDE {}
142 
144  void GenerateData() ITK_OVERRIDE;
145 
149 
155 
156  float m_Beta;
158 
162  void VerifyInputInformation() ITK_OVERRIDE {}
163 
166  void GenerateInputRequestedRegion() ITK_OVERRIDE;
167  void GenerateOutputInformation() ITK_OVERRIDE;
169 
170 private:
171  ADMMWaveletsConjugateGradientOperator(const Self &); //purposely not implemented
172  void operator=(const Self &); //purposely not implemented
173 
174 };
175 } //namespace RTK
176 
177 
178 #ifndef ITK_MANUAL_INSTANTIATION
179 #include "rtkADMMWaveletsConjugateGradientOperator.hxx"
180 #endif
181 
182 #endif
Base class for forward projection, i.e. accumulation along x-ray lines.
Weigting for displaced detectors.
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)
#define itkSetMacro(name, type)
void SetBackProjectionFilter(const BackProjectionFilterPointer _arg)