1 /*=========================================================================
2 Program: vv http://www.creatis.insa-lyon.fr/rio/vv
5 - University of LYON http://www.universite-lyon.fr/
6 - Léon Bérard cancer center http://www.centreleonberard.fr
7 - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr
9 This software is distributed WITHOUT ANY WARRANTY; without even
10 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
11 PURPOSE. See the copyright notices for more information.
13 It is distributed under dual licence
15 - BSD See included LICENSE.txt file
16 - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ===========================================================================**/
25 #include "clitkCommon.h"
28 #include <vtkImageData.h>
29 #include <vtkTransform.h>
30 #include <vtkDataObject.h>
31 #include <vtkStreamingDemandDrivenPipeline.h>
32 #include <vtkInformation.h>
33 #include <vtkVersion.h>
38 //--------------------------------------------------------------------
43 //--------------------------------------------------------------------
46 //--------------------------------------------------------------------
53 //--------------------------------------------------------------------
55 //--------------------------------------------------------------------
60 //--------------------------------------------------------------------
62 //--------------------------------------------------------------------
68 //--------------------------------------------------------------------
70 //--------------------------------------------------------------------
71 void vvImage::AddVtkImage(vtkImageData* input)
74 // Note that we're simply adding a new image to the vector.
75 // mItkToVtkConverters is therefore not being updated, but
76 // up to here it's not being used anyway...
78 #if VTK_MAJOR_VERSION <= 5
79 int* extent = input->GetWholeExtent();
81 int* extent = input->GetInformation()->Get(vtkDataObject::DATA_EXTENT());
84 if (extent[4] != extent[5])
86 else if (extent[3] != extent[4])
88 else if (extent[0] != extent[1])
91 mVtkImages.push_back(input);
94 //--------------------------------------------------------------------
96 //--------------------------------------------------------------------
97 int vvImage::GetNumberOfSpatialDimensions()
99 return mImageDimension;
101 //--------------------------------------------------------------------
103 //--------------------------------------------------------------------
104 int vvImage::GetNumberOfDimensions() const
106 if (mVtkImages.size()) {
107 if (IsTimeSequence())
108 return mImageDimension+1;
110 return mImageDimension;
114 //--------------------------------------------------------------------
116 //--------------------------------------------------------------------
117 void vvImage::GetScalarRange(double* range)
119 assert(mVtkImages.size());
120 double * temp = mVtkImages[0]->GetScalarRange();
123 for (unsigned int i=1; i<mVtkImages.size(); i++) {
124 temp = mVtkImages[i]->GetScalarRange();
125 if (temp[0] < range[0]) range[0]=temp[0];
126 if (temp[1] > range[1]) range[1]=temp[1];
129 //--------------------------------------------------------------------
131 //--------------------------------------------------------------------
132 std::string vvImage::GetScalarTypeAsITKString()
134 // WARNING VTK pixel type different from ITK Pixel type
135 std::string vtktype = mVtkImages[0]->GetScalarTypeAsString();
136 if (vtktype == "unsigned char") return "unsigned_char";
137 if (vtktype == "unsigned short") return "unsigned_short";
138 if (vtktype == "unsigned int") return "unsigned_int";
141 //--------------------------------------------------------------------
143 //--------------------------------------------------------------------
144 int vvImage::GetNumberOfScalarComponents()
146 return mVtkImages[0]->GetNumberOfScalarComponents();
148 //--------------------------------------------------------------------
150 //--------------------------------------------------------------------
151 int vvImage::GetScalarSize()
153 return mVtkImages[0]->GetScalarSize();
155 //--------------------------------------------------------------------
157 //--------------------------------------------------------------------
158 std::vector<double> vvImage::GetSpacing()
160 std::vector<double> spacing;
161 int dim = this->GetNumberOfDimensions();
162 for (int i = 0; i < dim; i++) {
164 spacing.push_back(mTimeSpacing);
166 spacing.push_back(mVtkImages[0]->GetSpacing()[i]);
170 //--------------------------------------------------------------------
171 std::vector<double> vvImage::GetOrigin() const
173 std::vector<double> origin;
174 int dim = this->GetNumberOfDimensions();
175 for (int i = 0; i < dim; i++) {
177 origin.push_back(mTimeOrigin);
179 origin.push_back(mVtkImages[0]->GetOrigin()[i]);
183 //--------------------------------------------------------------------
185 //--------------------------------------------------------------------
186 std::vector<int> vvImage::GetSize()
188 std::vector<int> size0;
189 int dim = this->GetNumberOfDimensions();
190 for (int i = 0; i < dim; i++) {
192 size0.push_back(mVtkImages.size());
194 size0.push_back(mVtkImages[0]->GetDimensions()[i]);
198 //--------------------------------------------------------------------
200 //--------------------------------------------------------------------
201 unsigned long vvImage::GetActualMemorySize()
203 unsigned long size = 0;
204 for (unsigned int i = 0; i < mVtkImages.size(); i++) {
205 size += mVtkImages[i]->GetActualMemorySize();
209 //--------------------------------------------------------------------
212 //--------------------------------------------------------------------
213 bool vvImage::IsTimeSequence() const
215 return mVtkImages.size()>1;
217 //--------------------------------------------------------------------
220 //--------------------------------------------------------------------
221 const std::vector<vtkImageData*>& vvImage::GetVTKImages()
225 //--------------------------------------------------------------------
227 //--------------------------------------------------------------------
228 vtkImageData* vvImage::GetFirstVTKImageData()
230 return mVtkImages[0];
232 //--------------------------------------------------------------------
234 //--------------------------------------------------------------------
235 bool vvImage::IsScalarTypeInteger()
237 assert(mVtkImages.size()> 0);
238 int t = mVtkImages[0]->GetScalarType();
239 return IsScalarTypeInteger(t);
241 //--------------------------------------------------------------------
243 //--------------------------------------------------------------------
244 bool vvImage::IsScalarTypeInteger(int t)
246 if ((t == VTK_BIT) ||
248 (t == VTK_UNSIGNED_CHAR) ||
250 (t == VTK_UNSIGNED_SHORT) ||
252 (t == VTK_UNSIGNED_INT) ||
254 (t == VTK_UNSIGNED_LONG)) {
260 //--------------------------------------------------------------------
262 //--------------------------------------------------------------------
263 const std::vector< vtkSmartPointer<vtkTransform> >& vvImage::GetTransform()
265 return this->mTransform;
267 //--------------------------------------------------------------------
270 //--------------------------------------------------------------------
271 bool vvImage::HaveSameSizeAndSpacingThan(vvImage * other)
274 for(int i=0; i<GetNumberOfDimensions(); i++) {
275 if ( GetSize()[i] != other->GetSize()[i]) same = false;
276 if ( GetSpacing()[i] != other->GetSpacing()[i]) same = false;
280 //--------------------------------------------------------------------
281 itk::MetaDataDictionary* vvImage::GetFirstMetaDataDictionary()
283 return mDictionary[0];
285 //--------------------------------------------------------------------
288 #endif // VVIMAGE_CXX