1 #include <cpPluginsImageMeshFilters/TriangleMeshToBinaryImageFilter.h>
2 #include <cpPlugins/BoundingBox.h>
3 #include <cpPlugins/Image.h>
4 #include <cpPlugins/Mesh.h>
6 #include <cpPlugins_Instances/ImageMeshFilters.h>
8 // -------------------------------------------------------------------------
9 cpPluginsImageMeshFilters::TriangleMeshToBinaryImageFilter::
10 TriangleMeshToBinaryImageFilter( )
13 this->_AddInput( "Input" );
14 this->_AddInput( "BoundingBox", false );
15 this->_AddOutput< cpPlugins::Image >( "Output" );
17 this->m_Parameters.ConfigureAsUint( "InsideValue" );
18 this->m_Parameters.ConfigureAsUint( "OutsideValue" );
19 this->m_Parameters.ConfigureAsUint( "MinimumSize" );
21 this->m_Parameters.SetUint( "InsideValue", 1 );
22 this->m_Parameters.SetUint( "OutsideValue", 0 );
23 this->m_Parameters.SetUint( "MinimumSize", 100 );
26 // -------------------------------------------------------------------------
27 cpPluginsImageMeshFilters::TriangleMeshToBinaryImageFilter::
28 ~TriangleMeshToBinaryImageFilter( )
32 // -------------------------------------------------------------------------
33 void cpPluginsImageMeshFilters::TriangleMeshToBinaryImageFilter::
36 typedef itk::Mesh< float, 3 > _3F;
37 typedef itk::Mesh< double, 3 > _3D;
39 auto _3f = this->GetInputData( "Input" )->GetITK< _3F >( );
40 auto _3d = this->GetInputData( "Input" )->GetITK< _3D >( );
41 if ( _3f != NULL ) this->_GD0( _3f );
42 else if( _3d != NULL ) this->_GD0( _3d );
43 else this->_Error( "No valid input mesh." );
46 // -------------------------------------------------------------------------
47 template< class _TMesh >
48 void cpPluginsImageMeshFilters::TriangleMeshToBinaryImageFilter::
51 this->_GD1< _TMesh, unsigned char >( mesh );
54 // -------------------------------------------------------------------------
55 template< class _TMesh, class _TPixel >
56 void cpPluginsImageMeshFilters::TriangleMeshToBinaryImageFilter::
59 typedef cpPlugins::BoundingBox _TBB;
60 typedef itk::Image< _TPixel, _TMesh::PointDimension > _TImage;
61 typedef itk::TriangleMeshToBinaryImageFilter< _TMesh, _TImage > _TFilter;
62 typedef typename _TImage::PointType _TPoint;
64 static const unsigned int PAD = 10;
66 _TFilter* filter = this->_CreateITK< _TFilter >( );
68 auto in_bb = dynamic_cast< _TBB* >( this->GetInputData( "BoundingBox" ) );
72 auto bb = mesh->GetBoundingBox( );
73 minBB = bb->GetMinimum( );
74 maxBB = bb->GetMaximum( );
78 minBB = in_bb->GetMinimum< _TPoint >( );
79 maxBB = in_bb->GetMaximum< _TPoint >( );
83 double lx = double( maxBB[ 0 ] - minBB[ 0 ] );
84 double ly = double( maxBB[ 1 ] - minBB[ 1 ] );
85 double lz = double( maxBB[ 2 ] - minBB[ 2 ] );
86 double lm = ( lx < ly )? lx: ly;
87 lm = ( lm < lz )? lm: lz;
90 double mSpac = lm / double( this->m_Parameters.GetUint( "MinimumSize" ) );
91 typename _TImage::SpacingType spac;
93 filter->SetSpacing( spac );
96 typename _TImage::SizeType size;
97 size[ 0 ] = ( unsigned long )( std::ceil( lx / mSpac ) );
98 size[ 1 ] = ( unsigned long )( std::ceil( ly / mSpac ) );
99 size[ 2 ] = ( unsigned long )( std::ceil( lz / mSpac ) );
101 // ... add some padding pixels
105 filter->SetSize( size );
108 typename _TImage::PointType origin = minBB;
109 origin[ 0 ] -= double( PAD >> 1 ) * spac[ 0 ];
110 origin[ 1 ] -= double( PAD >> 1 ) * spac[ 1 ];
111 origin[ 2 ] -= double( PAD >> 1 ) * spac[ 2 ];
112 filter->SetOrigin( origin );
115 typename _TImage::DirectionType direction;
116 direction.SetIdentity( );
117 filter->SetDirection( direction );
119 typename _TImage::IndexType index;
121 filter->SetIndex( index );
124 filter->SetInput( mesh );
125 filter->SetInsideValue( this->m_Parameters.GetUint( "InsideValue" ) );
126 filter->SetOutsideValue( this->m_Parameters.GetUint( "OutsideValue" ) );
130 auto out = this->GetOutputData( "Output" );
131 out->SetITK( filter->GetOutput( ) );