Slicer  4.10
Slicer is a multi-platform, free and open source software package for visualization and medical image computing
vtkPolyDataToFractionalLabelmapFilter.h
Go to the documentation of this file.
1 /*==============================================================================
2 
3  Copyright (c) Laboratory for Percutaneous Surgery (PerkLab)
4  Queen's University, Kingston, ON, Canada. All Rights Reserved.
5 
6  See COPYRIGHT.txt
7  or http://www.slicer.org/copyright/copyright.txt for details.
8 
9  Unless required by applicable law or agreed to in writing, software
10  distributed under the License is distributed on an "AS IS" BASIS,
11  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
12  See the License for the specific language governing permissions and
13  limitations under the License.
14 
15  This file was originally developed by Kyle Sunderland, PerkLab, Queen's University
16  and was supported through the Applied Cancer Research Unit program of Cancer Care
17  Ontario with funds provided by the Ontario Ministry of Health and Long-Term Care
18 
19  This file is a modified version of vtkPolyDataToImageStencil.h
20 
21 ==============================================================================*/
22 
23 #ifndef vtkPolyDataToFractionalLabelmapFilter_h
24 #define vtkPolyDataToFractionalLabelmapFilter_h
25 
26 
27 // VTK includes
28 #include <vtkPolyDataToImageStencil.h>
29 #include <vtkSmartPointer.h>
30 #include <vtkCellArray.h>
31 #include <vtkSetGet.h>
32 #include <vtkMatrix4x4.h>
33 #include <vtkCellLocator.h>
34 
35 // Segmentations includes
36 #include <vtkOrientedImageData.h>
37 
38 // std includes
39 #include <map>
40 
41 #include "vtkSegmentationCoreConfigure.h"
42 
43 // Define the datatype and fractional constants for fractional labelmap conversion based on the value of VTK_FRACTIONAL_DATA_TYPE
44 #define VTK_FRACTIONAL_DATA_TYPE VTK_CHAR
45 
46 #if VTK_FRACTIONAL_DATA_TYPE == VTK_UNSIGNED_CHAR
47  #define FRACTIONAL_DATA_TYPE VTK_TYPE_NAME_UNSIGNED_CHAR
48  #define FRACTIONAL_MIN 0
49  #define FRACTIONAL_MAX 216
50  #define FRACTIONAL_STEP_SIZE 1
51 #elif VTK_FRACTIONAL_DATA_TYPE == VTK_CHAR
52  #define FRACTIONAL_DATA_TYPE VTK_TYPE_NAME_CHAR
53  #define FRACTIONAL_MIN -108
54  #define FRACTIONAL_MAX 108
55  #define FRACTIONAL_STEP_SIZE 1
56 #elif VTK_FRACTIONAL_DATA_TYPE == VTK_FLOAT
57  #define FRACTIONAL_DATA_TYPE VTK_TYPE_NAME_FLOAT
58  #define FRACTIONAL_MIN 0.0
59  #define FRACTIONAL_MAX 1.0
60  #define FRACTIONAL_STEP_SIZE (1.0/216.0)
61 #endif
62 
63 class vtkSegmentationCore_EXPORT vtkPolyDataToFractionalLabelmapFilter :
64  public vtkPolyDataToImageStencil
65 {
66 private:
67  std::map<double, vtkSmartPointer<vtkCellArray> > LinesCache;
68  std::map<double, vtkSmartPointer<vtkPolyData> > SliceCache;
69  std::map<double, vtkIdType*> PointIdsCache;
70  std::map<double, vtkIdType> NptsCache;
71  std::map<double, vtkSmartPointer<vtkIdTypeArray> > PointNeighborCountsCache;
72 
73  vtkCellLocator* CellLocator;
74 
75  vtkOrientedImageData* OutputImageTransformData;
76  int NumberOfOffsets;
77 
78 public:
80  vtkTypeMacro(vtkPolyDataToFractionalLabelmapFilter, vtkPolyDataToImageStencil);
81 
82  virtual vtkOrientedImageData* GetOutput();
83  virtual void SetOutput(vtkOrientedImageData* output);
84 
85  void SetOutputImageToWorldMatrix(vtkMatrix4x4* imageToWorldMatrix);
86  void GetOutputImageToWorldMatrix(vtkMatrix4x4* imageToWorldMatrix);
87 
88  using Superclass::GetOutputOrigin;
89  double* GetOutputOrigin() VTK_OVERRIDE;
90  void GetOutputOrigin(double origin[3]) VTK_OVERRIDE;
91 
92  void SetOutputOrigin(double origin[3]) VTK_OVERRIDE;
93  void SetOutputOrigin(double x, double y, double z) VTK_OVERRIDE;
94 
95  using Superclass::GetOutputSpacing;
96  double* GetOutputSpacing() VTK_OVERRIDE;
97  void GetOutputSpacing(double spacing[3]) VTK_OVERRIDE;
98 
99  void SetOutputSpacing(double spacing[3]) VTK_OVERRIDE;
100  void SetOutputSpacing(double x, double y, double z) VTK_OVERRIDE;
101 
102 
104  void DeleteCache();
105 
106  vtkSetMacro(NumberOfOffsets, int);
107  vtkGetMacro(NumberOfOffsets, int);
108 
109 protected:
112 
113  virtual int RequestData(vtkInformation *, vtkInformationVector **,
114  vtkInformationVector *) VTK_OVERRIDE;
115  vtkOrientedImageData *AllocateOutputData(vtkDataObject *out, int* updateExt);
116  virtual int FillOutputPortInformation(int, vtkInformation*) VTK_OVERRIDE;
117 
123  void FillImageStencilData(vtkImageStencilData *output, vtkPolyData* closedSurface, int extent[6]);
124 
128  void AddBinaryLabelMapToFractionalLabelMap(vtkImageData* binaryLabelMap, vtkImageData* fractionalLabelMap);
129 
135  void PolyDataCutter(vtkPolyData *input, vtkPolyData *output,
136  double z);
137 
138 private:
140  void operator=(const vtkPolyDataToFractionalLabelmapFilter&); // Not implemented.
141 };
142 
143 #endif
Image data containing orientation information.