ITK  4.6.0
Insight Segmentation and Registration Toolkit
itkImageRegion.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright Insight Software 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 /*=========================================================================
19  *
20  * Portions of this file are subject to the VTK Toolkit Version 3 copyright.
21  *
22  * Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
23  *
24  * For complete copyright, license and disclaimer of warranty information
25  * please refer to the NOTICE file at the top of the ITK source tree.
26  *
27  *=========================================================================*/
28 #ifndef __itkImageRegion_h
29 #define __itkImageRegion_h
30 
31 #include "itkRegion.h"
32 
33 #include "itkSize.h"
34 #include "itkContinuousIndex.h"
35 #include "vnl/vnl_math.h"
36 
37 namespace itk
38 {
39 // Forward declaration of ImageBase so it can be declared a friend
40 // (needed for PrintSelf mechanism)
41 template< unsigned int VImageDimension >
42 class ImageBase;
43 
68 template< unsigned int VImageDimension >
69 class ImageRegion:public Region
70 {
71 public:
73  typedef ImageRegion Self;
74  typedef Region Superclass;
75 
77  itkTypeMacro(ImageRegion, Region);
78 
80  itkStaticConstMacro(ImageDimension, unsigned int, VImageDimension);
81 
84  itkStaticConstMacro( SliceDimension, unsigned int,
85  ( ImageDimension - ( ImageDimension > 1 ) ) );
86 
88  static unsigned int GetImageDimension()
89  { return ImageDimension; }
90 
98 
102 
105 
107  virtual typename Superclass::RegionType GetRegionType() const ITK_OVERRIDE
109 
112  ImageRegion();
113 
116  virtual ~ImageRegion();
117 
120  ImageRegion(const Self & region):Region(region), m_Index(region.m_Index), m_Size(region.m_Size) {}
121 
124  ImageRegion(const IndexType & index, const SizeType & size)
125  { m_Index = index; m_Size = size; }
126 
130  ImageRegion(const SizeType & size)
131  { m_Size = size; m_Index.Fill(0); }
132 
135  void operator=(const Self & region)
136  { m_Index = region.m_Index;m_Size = region.m_Size; }
137 
139  void SetIndex(const IndexType & index)
140  { m_Index = index; }
141 
143  const IndexType & GetIndex() const { return m_Index; }
146 
149  void SetSize(const SizeType & size)
150  { m_Size = size; }
151 
153  const SizeType & GetSize() const { return m_Size; }
156 
159  void SetSize(unsigned int i, SizeValueType sze)
160  { m_Size[i] = sze; }
161  SizeValueType GetSize(unsigned int i) const
162  { return m_Size[i]; }
164 
167  void SetIndex(unsigned int i, IndexValueType sze)
168  { m_Index[i] = sze; }
169  IndexValueType GetIndex(unsigned int i) const
170  { return m_Index[i]; }
172 
174  IndexType GetUpperIndex() const;
175 
177  void SetUpperIndex( const IndexType & idx );
178 
180  void ComputeOffsetTable(OffsetTableType offsetTable) const;
181 
183  bool
184  operator==(const Self & region) const
185  {
186  bool same = ( m_Index == region.m_Index );
187  same = same && ( m_Size == region.m_Size );
188  return same;
189  }
191 
193  bool
194  operator!=(const Self & region) const
195  {
196  bool same = ( m_Index == region.m_Index );
197  same = same && ( m_Size == region.m_Size );
198  return !same;
199  }
201 
203  bool
204  IsInside(const IndexType & index) const
205  {
206  for ( unsigned int i = 0; i < ImageDimension; i++ )
207  {
208  if ( index[i] < m_Index[i] )
209  {
210  return false;
211  }
212  if ( index[i] >= ( m_Index[i] + static_cast< IndexValueType >( m_Size[i] ) ) )
213  {
214  return false;
215  }
216  }
217  return true;
218  }
220 
225  template< typename TCoordRepType >
226  bool
228  {
229  for ( unsigned int i = 0; i < ImageDimension; i++ )
230  {
231  if ( Math::RoundHalfIntegerUp< IndexValueType >(index[i]) < static_cast< IndexValueType >( m_Index[i] ) )
232  {
233  return false;
234  }
235  // bound is the last valid pixel location
236  const TCoordRepType bound = static_cast< TCoordRepType >(
237  m_Index[i] + m_Size[i] - 0.5 );
239 
240  /* Note for NaN: test using negation of a positive test in order
241  * to always evaluate to true (and thus return false) when index[i]
242  * is NaN. The cast above to integer via RoundHalfIntegerUp will cast
243  * NaN into a platform-dependent value (large negative, -1 or large
244  * positive, empirically). Thus this test here is relied on
245  * to 'catch' NaN's. */
246  if ( ! (index[i] <= bound) )
247  {
248  return false;
249  }
250  }
251  return true;
252  }
253 
258  bool
259  IsInside(const Self & region) const
260  {
261  IndexType beginCorner = region.GetIndex();
262 
263  if ( !this->IsInside(beginCorner) )
264  {
265  return false;
266  }
267  IndexType endCorner;
268  SizeType size = region.GetSize();
269  for ( unsigned int i = 0; i < ImageDimension; i++ )
270  {
271  endCorner[i] = beginCorner[i] + static_cast< OffsetValueType >( size[i] ) - 1;
272  }
273  if ( !this->IsInside(endCorner) )
274  {
275  return false;
276  }
277  return true;
278  }
279 
283 
287  void PadByRadius(OffsetValueType radius);
288 
289  void PadByRadius(const IndexValueArrayType radius);
290 
291  void PadByRadius(const SizeType & radius);
292 
297  bool ShrinkByRadius(OffsetValueType radius);
298 
299  bool ShrinkByRadius(const IndexValueArrayType radius);
300 
301  bool ShrinkByRadius(const SizeType & radius);
302 
307  bool Crop(const Self & region);
308 
312  SliceRegion Slice(const unsigned int dim) const;
313 
314 protected:
319  virtual void PrintSelf(std::ostream & os, Indent indent) const ITK_OVERRIDE;
320 
321 private:
324 
326  friend class ImageBase< VImageDimension >;
327 };
328 
329 template< unsigned int VImageDimension >
330 std::ostream & operator<<(std::ostream & os, const ImageRegion< VImageDimension > & region);
331 } // end namespace itk
332 
333 #ifndef ITK_MANUAL_INSTANTIATION
334 #include "itkImageRegion.hxx"
335 #endif
336 
337 #endif
void SetSize(const SizeType &size)
void SetUpperIndex(const IndexType &idx)
void ComputeOffsetTable(OffsetTableType offsetTable) const
IndexValueType GetIndex(unsigned int i) const
const IndexType & GetIndex() const
Represent the offset between two n-dimensional indexes in a n-dimensional image.
Definition: itkOffset.h:55
Size< itkGetStaticConstMacro(ImageDimension) > SizeType
bool ShrinkByRadius(OffsetValueType radius)
signed long OffsetValueType
Definition: itkIntTypes.h:154
bool IsInside(const IndexType &index) const
Index< itkGetStaticConstMacro(ImageDimension) > IndexType
void PadByRadius(OffsetValueType radius)
An image region represents a structured region of data.
bool operator==(const Self &region) const
static unsigned int GetImageDimension()
OffsetType::OffsetValueType OffsetValueType
static const unsigned int ImageDimension
unsigned long SizeValueType
Definition: itkIntTypes.h:143
SizeValueType GetSize(unsigned int i) const
IndexValueType IndexValueArrayType[ImageDimension]
ImageRegion< itkGetStaticConstMacro(SliceDimension) > SliceRegion
bool Crop(const Self &region)
ImageRegion(const SizeType &size)
IndexType & GetModifiableIndex()
SizeType & GetModifiableSize()
void Fill(IndexValueType value)
Definition: itkIndex.h:292
virtual ~ImageRegion()
IndexType::IndexValueType IndexValueType
void SetSize(unsigned int i, SizeValueType sze)
const SizeType & GetSize() const
A region represents some portion or piece of data.
Definition: itkRegion.h:64
SizeType::SizeValueType SizeValueType
IndexType::OffsetType OffsetType
void operator=(const Self &region)
bool operator!=(const Self &region) const
bool IsInside(const Self &region) const
Base class for templated image classes.
Definition: itkImageBase.h:112
ImageRegion(const Self &region)
A templated class holding a point in n-Dimensional image space.
SliceRegion Slice(const unsigned int dim) const
Control indentation during Print() invocation.
Definition: itkIndent.h:49
virtual Superclass::RegionType GetRegionType() const ITK_OVERRIDE
ImageRegion Self
IndexType GetUpperIndex() const
static const unsigned int SliceDimension
void SetIndex(unsigned int i, IndexValueType sze)
OffsetValueType OffsetTableType[ImageDimension+1]
SizeValueType GetNumberOfPixels() const
itk::OffsetValueType OffsetValueType
Definition: itkOffset.h:69
void SetIndex(const IndexType &index)
virtual void PrintSelf(std::ostream &os, Indent indent) const ITK_OVERRIDE
ImageRegion(const IndexType &index, const SizeType &size)
bool IsInside(const ContinuousIndex< TCoordRepType, VImageDimension > &index) const