RTK  1.4.0
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 
51 
53  itkNewMacro(Self)
54 
55 
57 
59  itkSetMacro(Radius, typename TImage::RegionType::SizeType)
60  itkGetMacro(Radius, typename TImage::RegionType::SizeType)
61 
63  itkSetMacro(RegularizationWeights, typename TImage::PixelType)
64  itkGetMacro(RegularizationWeights, typename TImage::PixelType)
65 
66 protected:
67  SeparableQuadraticSurrogateRegularizationImageFilter();
68  virtual ~SeparableQuadraticSurrogateRegularizationImageFilter() ITK_OVERRIDE {}
69 
71  itk::DataObject::Pointer MakeOutput(itk::ProcessObject::DataObjectPointerArraySizeType idx) override;
72 
74 #if ITK_VERSION_MAJOR<5
75  void ThreadedGenerateData(const typename TImage::RegionType& outputRegionForThread, itk::ThreadIdType itkNotUsed(threadId)) ITK_OVERRIDE;
76 #else
77  void DynamicThreadedGenerateData(const typename TImage::RegionType& outputRegionForThread) ITK_OVERRIDE;
78 #endif
79  void GenerateInputRequestedRegion() ITK_OVERRIDE;
81 
83  typename TImage::PixelType GreenPriorFirstDerivative(typename TImage::PixelType pix);
84  typename TImage::PixelType GreenPriorSecondDerivative(typename TImage::PixelType pix);
86 
88  typename TImage::RegionType::SizeType m_Radius;
91  typename TImage::PixelType m_RegularizationWeights;
92 
93 private:
94  SeparableQuadraticSurrogateRegularizationImageFilter(const Self &); //purposely not implemented
95  void operator=(const Self &); //purposely not implemented
96 
97 };
98 } //namespace RTK
99 
100 
101 #ifndef ITK_MANUAL_INSTANTIATION
102 #include "rtkSeparableQuadraticSurrogateRegularizationImageFilter.hxx"
103 #endif
104 
105 #endif
TImage::PixelType GreenPriorSecondDerivative(typename TImage::PixelType pix)
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 ...
itk::DataObject::Pointer MakeOutput(itk::ProcessObject::DataObjectPointerArraySizeType idx) override
#define itkSetMacro(name, type)