RTK  2.0.1
Reconstruction Toolkit
rtkFourDToProjectionStackImageFilter.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 rtkFourDToProjectionStackImageFilter_h
19 #define rtkFourDToProjectionStackImageFilter_h
20 
21 #include <itkExtractImageFilter.h>
22 #include <itkPasteImageFilter.h>
23 #include <itkMultiplyImageFilter.h>
24 
25 #include "rtkConstantImageSource.h"
28 
29 namespace rtk
30 {
98 template< typename ProjectionStackType, typename VolumeSeriesType>
99 class FourDToProjectionStackImageFilter : public itk::ImageToImageFilter<ProjectionStackType, ProjectionStackType>
100 {
101 public:
102  ITK_DISALLOW_COPY_AND_ASSIGN(FourDToProjectionStackImageFilter);
103 
108 
110  using VolumeType = ProjectionStackType;
111 
113  itkNewMacro(Self)
114 
115 
117 
119  void SetInputVolumeSeries(const VolumeSeriesType* VolumeSeries);
120 
123  void SetInputProjectionStack(const ProjectionStackType* Projection);
124 
126  using ForwardProjectionFilterType = rtk::ForwardProjectionImageFilter< ProjectionStackType, ProjectionStackType >;
127  using PasteFilterType = itk::PasteImageFilter<ProjectionStackType, ProjectionStackType>;
132 
135 
137  virtual void SetGeometry(GeometryType::Pointer _arg);
138 
140  void SetWeights(const itk::Array2D<float> _arg);
141 
144 
146  virtual void SetSignal(const std::vector<double> signal);
147 
148 protected:
149  FourDToProjectionStackImageFilter();
150  ~FourDToProjectionStackImageFilter() override = default;
151 
152  typename VolumeSeriesType::ConstPointer GetInputVolumeSeries();
153  typename ProjectionStackType::Pointer GetInputProjectionStack();
154 
156  void GenerateData() override;
157 
158  void GenerateOutputInformation() override;
159 
160  void GenerateInputRequestedRegion() override;
161 
168 
172  typename ConstantProjectionStackSourceType::OutputImageRegionType m_PasteRegion;
173  std::vector<double> m_Signal;
174 
175 };
176 } //namespace ITK
177 
178 
179 #ifndef ITK_MANUAL_INSTANTIATION
180 #include "rtkFourDToProjectionStackImageFilter.hxx"
181 #endif
182 
183 #endif
Base class for forward projection, i.e. accumulation along x-ray lines.
ConstantProjectionStackSourceType::Pointer m_ConstantProjectionStackSource
Generate an n-dimensional image with constant pixel values.
virtual void SetSignal(const std::vector< double > signal)
STL namespace.
VolumeSeriesType::ConstPointer GetInputVolumeSeries()
virtual void SetGeometry(GeometryType::Pointer _arg)
Projection geometry for a source and a 2-D flat panel.
ProjectionStackType::Pointer GetInputProjectionStack()
void SetForwardProjectionFilter(const typename ForwardProjectionFilterType::Pointer _arg)
void SetWeights(const itk::Array2D< float > _arg)
void SetInputProjectionStack(const ProjectionStackType *Projection)
Interpolates (linearly) in a 3D+t sequence of volumes to get a 3D volume.
ForwardProjectionFilterType::Pointer m_ForwardProjectionFilter
void SetInputVolumeSeries(const VolumeSeriesType *VolumeSeries)
ConstantProjectionStackSourceType::OutputImageRegionType m_PasteRegion
Implements part of the 4D reconstruction by conjugate gradient.