RTK  2.0.1
Reconstruction Toolkit
rtkSeparableQuadraticSurrogateRegularizationImageFilter.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 rtkSeparableQuadraticSurrogateRegularizationImageFilter_h
20 #define rtkSeparableQuadraticSurrogateRegularizationImageFilter_h
21 
23 #include <itkImageRegionIterator.h>
24 #include <itkImageToImageFilter.h>
25 
26 namespace rtk
27 {
28 
42 template<typename TImage>
44 {
45 public:
46  ITK_DISALLOW_COPY_AND_ASSIGN(SeparableQuadraticSurrogateRegularizationImageFilter);
47 
52 
54  itkNewMacro(Self)
55 
56 
58 
60  itkSetMacro(Radius, typename TImage::RegionType::SizeType)
61  itkGetMacro(Radius, typename TImage::RegionType::SizeType)
62 
64  itkSetMacro(RegularizationWeights, typename TImage::PixelType)
65  itkGetMacro(RegularizationWeights, typename TImage::PixelType)
66 
67 protected:
68  SeparableQuadraticSurrogateRegularizationImageFilter();
69  ~SeparableQuadraticSurrogateRegularizationImageFilter() override = default;
70 
72  itk::ProcessObject::DataObjectPointer MakeOutput(itk::ProcessObject::DataObjectPointerArraySizeType idx) override;
73  itk::ProcessObject::DataObjectPointer MakeOutput(const itk::ProcessObject::DataObjectIdentifierType &) override;
75 
77 #if ITK_VERSION_MAJOR<5
78  void ThreadedGenerateData(const typename TImage::RegionType& outputRegionForThread, itk::ThreadIdType itkNotUsed(threadId)) override;
79 #else
80  void DynamicThreadedGenerateData(const typename TImage::RegionType& outputRegionForThread) override;
81 #endif
82  void GenerateInputRequestedRegion() override;
84 
86  typename TImage::PixelType GreenPriorFirstDerivative(typename TImage::PixelType pix);
87  typename TImage::PixelType GreenPriorSecondDerivative(typename TImage::PixelType pix);
89 
91  typename TImage::RegionType::SizeType m_Radius;
94  typename TImage::PixelType m_RegularizationWeights;
95 
96 };
97 } //namespace RTK
98 
99 
100 #ifndef ITK_MANUAL_INSTANTIATION
101 #include "rtkSeparableQuadraticSurrogateRegularizationImageFilter.hxx"
102 #endif
103 
104 #endif
TImage::PixelType GreenPriorSecondDerivative(typename TImage::PixelType pix)
itk::ProcessObject::DataObjectPointer MakeOutput(itk::ProcessObject::DataObjectPointerArraySizeType idx) override
typename TPixelType::ValueType ValueType
TImage::PixelType GreenPriorFirstDerivative(typename TImage::PixelType pix)
ImageBaseType::SizeType SizeType
unsigned int ThreadIdType
void ThreadedGenerateData(const typename TImage::RegionType &outputRegionForThread, itk::ThreadIdType) override
For one-step inversion of spectral CT data by the method Mechlem2017, computes regularization term&#39;s ...
ImageBaseType::RegionType RegionType
#define itkSetMacro(name, type)