1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
/*=========================================================================
Program: Visualization Toolkit
Module: vtkImageIterator.txx
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notice for more information.
=========================================================================*/
// Include blockers needed since vtkImageIterator.h includes this file
// when VTK_NO_EXPLICIT_TEMPLATE_INSTANTIATION is defined.
#ifndef vtkImageIterator_txx
#define vtkImageIterator_txx
#include "vtkImageIterator.h"
#include "vtkImageData.h"
//----------------------------------------------------------------------------
template <class DType>
vtkImageIterator<DType>::vtkImageIterator()
{
this->Pointer = 0;
this->EndPointer = 0;
this->SpanEndPointer = 0;
this->SliceEndPointer = 0;
}
//----------------------------------------------------------------------------
template <class DType>
void vtkImageIterator<DType>::Initialize(vtkImageData *id, int *ext)
{
this->Pointer = static_cast<DType *>(id->GetScalarPointerForExtent(ext));
id->GetIncrements(this->Increments[0], this->Increments[1],
this->Increments[2]);
id->GetContinuousIncrements(ext,this->ContinuousIncrements[0],
this->ContinuousIncrements[1],
this->ContinuousIncrements[2]);
this->EndPointer =
static_cast<DType *>(id->GetScalarPointer(ext[1],ext[3],ext[5]))
+this->Increments[0];
// if the extent is empty then the end pointer should equal the beg pointer
if (ext[1] < ext[0] || ext[3] < ext[2] || ext[5] < ext[4])
{
this->EndPointer = this->Pointer;
}
this->SpanEndPointer =
this->Pointer + this->Increments[0]*(ext[1] - ext[0] + 1);
this->SliceEndPointer =
this->Pointer + this->Increments[1]*(ext[3] - ext[2] + 1);
}
//----------------------------------------------------------------------------
template <class DType>
vtkImageIterator<DType>::vtkImageIterator(vtkImageData *id, int *ext)
{
this->Initialize(id, ext);
}
//----------------------------------------------------------------------------
template <class DType>
void vtkImageIterator<DType>::NextSpan()
{
this->Pointer += this->Increments[1];
this->SpanEndPointer += this->Increments[1];
if (this->Pointer >= this->SliceEndPointer)
{
this->Pointer += this->ContinuousIncrements[2];
this->SpanEndPointer += this->ContinuousIncrements[2];
this->SliceEndPointer += this->Increments[2];
}
}
#endif
|
