RTK  2.0.1
Reconstruction Toolkit
rtkGetNewtonUpdateImageFilter.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 #ifndef rtkGetNewtonUpdateImageFilter_h
19 #define rtkGetNewtonUpdateImageFilter_h
20 
21 #include "itkImageToImageFilter.h"
22 #include "rtkMacro.h"
23 
24 namespace rtk
25 {
41 template< class TGradient,
43 class GetNewtonUpdateImageFilter : public itk::ImageToImageFilter<TGradient, TGradient>
44 {
45 public:
46  ITK_DISALLOW_COPY_AND_ASSIGN(GetNewtonUpdateImageFilter);
47 
52 
54  itkNewMacro(Self)
55 
56 
58 
60  static constexpr unsigned int nChannels = TGradient::PixelType::Dimension;
61 
63  using dataType = typename TGradient::PixelType::ValueType;
64 
66  void SetInputGradient(const TGradient* gradient);
67  void SetInputHessian(const THessian* hessian);
69 
70 protected:
71  GetNewtonUpdateImageFilter();
72  ~GetNewtonUpdateImageFilter() override = default;
73 
74  void GenerateInputRequestedRegion() override;
75 
77 #if ITK_VERSION_MAJOR<5
78  void ThreadedGenerateData(const typename TGradient::RegionType& outputRegionForThread, itk::ThreadIdType itkNotUsed(threadId)) override;
79 #else
80  void DynamicThreadedGenerateData(const typename TGradient::RegionType& outputRegionForThread) override;
81 #endif
82 
83 
85  typename TGradient::ConstPointer GetInputGradient();
86  typename THessian::ConstPointer GetInputHessian();
88 
89 };
90 } //namespace RTK
91 
92 
93 #ifndef ITK_MANUAL_INSTANTIATION
94 #include "rtkGetNewtonUpdateImageFilter.hxx"
95 #endif
96 
97 #endif
typename TGradient::PixelType::ValueType dataType
void SetInputGradient(const TGradient *gradient)
THessian::ConstPointer GetInputHessian()
void GenerateInputRequestedRegion() override
static constexpr unsigned int nChannels
Computes update from gradient and Hessian in Newton&#39;s method.
void ThreadedGenerateData(const typename TGradient::RegionType &outputRegionForThread, itk::ThreadIdType) override
TGradient::ConstPointer GetInputGradient()
unsigned int ThreadIdType
void SetInputHessian(const THessian *hessian)