File: itkBinaryDiamondStructuringElement.txx

package info (click to toggle)
itksnap 3.4.0-2
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 10,196 kB
  • ctags: 9,196
  • sloc: cpp: 62,895; sh: 175; makefile: 13
file content (63 lines) | stat: -rw-r--r-- 1,861 bytes parent folder | download | duplicates (4)
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
/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkBinaryDiamondStructuringElement.txx,v $
  Language:  C++
  Date:      $Date: 2008/12/10 18:42:15 $
  Version:   $Revision: 1.2 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/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 notices for more information.

=========================================================================*/
#ifndef __itkBinaryDiamondStructuringElement_txx
#define __itkBinaryDiamondStructuringElement_txx
#include "itkBinaryDiamondStructuringElement.h"

#include "itkNumericTraits.h"

namespace itk
{

// Create the structuring element
template <class TPixel, unsigned int VDimension, class TAllocator>
void
BinaryDiamondStructuringElement<TPixel, VDimension, TAllocator>
::CreateStructuringElement()
{

  unsigned int minRadius = this->GetRadius( 0 );
  for ( unsigned d = 1; d < NeighborhoodDimension; d++ )
    {
    if ( minRadius > this->GetRadius( d ) )
      {
      minRadius = this->GetRadius( d );
      }    
    }

  for ( unsigned int n = 0; n < this->Size(); n++ )
    {
    OffsetType offset = this->GetOffset( n );
    unsigned int manhattanDistance = 0;
    for ( unsigned int d = 0; d < NeighborhoodDimension; d++ )
      {
      manhattanDistance += vnl_math_abs( offset[d] );
      } 
    if ( manhattanDistance <= minRadius )
      {
      this->operator[]( n ) = NumericTraits<TPixel>::One;
      } 
    else
      {
      this->operator[]( n ) = NumericTraits<TPixel>::Zero;
      }
    }  
}

} // namespace itk

#endif