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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkImage.txx,v $
Language: C++
Date: $Date: 2007-10-14 22:15:56 $
Version: $Revision: 1.100 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
Portions of this code are covered under the VTK copyright.
See VTKCopyright.txt or http://www.kitware.com/VTKCopyright.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 notices for more information.
=========================================================================*/
#ifndef __itkImage_txx
#define __itkImage_txx
#include "itkImage.h"
#include "itkProcessObject.h"
namespace itk
{
/**
*
*/
template<class TPixel, unsigned int VImageDimension>
Image<TPixel, VImageDimension>
::Image()
{
m_Buffer = PixelContainer::New();
}
//----------------------------------------------------------------------------
template<class TPixel, unsigned int VImageDimension>
void
Image<TPixel, VImageDimension>
::Allocate()
{
unsigned long num;
this->ComputeOffsetTable();
num = this->GetOffsetTable()[VImageDimension];
m_Buffer->Reserve(num);
}
template<class TPixel, unsigned int VImageDimension>
void
Image<TPixel, VImageDimension>
::Initialize()
{
//
// We don't modify ourselves because the "ReleaseData" methods depend upon
// no modification when initialized.
//
// Call the superclass which should initialize the BufferedRegion ivar.
Superclass::Initialize();
// Replace the handle to the buffer. This is the safest thing to do,
// since the same container can be shared by multiple images (e.g.
// Grafted outputs and in place filters).
m_Buffer = PixelContainer::New();
}
template<class TPixel, unsigned int VImageDimension>
void
Image<TPixel, VImageDimension>
::FillBuffer (const TPixel& value)
{
const unsigned long numberOfPixels =
this->GetBufferedRegion().GetNumberOfPixels();
for(unsigned int i=0; i<numberOfPixels; i++)
{
(*m_Buffer)[i] = value;
}
}
template<class TPixel, unsigned int VImageDimension>
void
Image<TPixel, VImageDimension>
::SetPixelContainer(PixelContainer *container)
{
if (m_Buffer != container)
{
m_Buffer = container;
this->Modified();
}
}
//----------------------------------------------------------------------------
template<class TPixel, unsigned int VImageDimension>
void
Image<TPixel, VImageDimension>
::Graft(const DataObject *data)
{
// call the superclass' implementation
Superclass::Graft( data );
if ( data )
{
// Attempt to cast data to an Image
const Self * imgData;
try
{
imgData = dynamic_cast<const Self *>( data );
}
catch( ... )
{
return;
}
if ( imgData )
{
// Now copy anything remaining that is needed
this->SetPixelContainer( const_cast< PixelContainer * >
(imgData->GetPixelContainer()) );
}
else
{
// pointer could not be cast back down
itkExceptionMacro( << "itk::Image::Graft() cannot cast "
<< typeid(data).name() << " to "
<< typeid(const Self *).name() );
}
}
}
/**
*
*/
template<class TPixel, unsigned int VImageDimension>
void
Image<TPixel, VImageDimension>
::PrintSelf(std::ostream& os, Indent indent) const
{
Superclass::PrintSelf(os,indent);
os << indent << "PixelContainer: " << std::endl;
m_Buffer->Print(os, indent.GetNextIndent());
// m_Origin and m_Spacing are printed in the Superclass
}
} // end namespace itk
#endif
|
