// @author Leonardo Florez Valencia
// @email florez-l@javeriana.edu.co
// =========================================================================
-
#ifndef __CTBronchi__MoriLabelling__h__
#define __CTBronchi__MoriLabelling__h__
-#include <itkImageToImageFilter.h>
-#include <itkHessianRecursiveGaussianImageFilter.h>
-#include <itkHessian3DToVesselnessMeasureImageFilter.h>
-
-#include <itkImageFileWriter.h>
-
+#include <fpa/Filters/BaseMarksInterface.h>
+#include <fpa/Filters/Image/SeedsFromLabelsInterface.h>
+#include <fpa/Filters/Image/DefaultTraits.h>
+#include <fpa/Filters/Image/RegionGrow.h>
+#include <fpa/Functors/RegionGrow/BinaryThreshold.h>
namespace CTBronchi
{
/**
*/
- template< class _TInputImage, class _TLabelImage >
- class MoriLabelling
- : public itk::ImageToImageFilter< _TLabelImage, _TLabelImage >
+ template< class _TImage, class _TLabels >
+ class MoriLabellingTraits
+ : public fpa::Filters::Image::DefaultTraits< _TImage, _TLabels, typename _TLabels::PixelType >
{
public:
- typedef MoriLabelling Self;
- typedef itk::ImageToImageFilter< _TLabelImage, _TLabelImage > Superclass;
- typedef itk::SmartPointer< Self > Pointer;
- typedef itk::SmartPointer< const Self > ConstPointer;
+ typedef fpa::Filters::Image::DefaultTraits< _TImage, _TLabels, typename _TLabels::PixelType > Superclass;
+ typedef typename Superclass::TInternalTraits TInternalTraits;
+ typedef typename Superclass::TMarksImage TMarksImage;
+ typedef typename Superclass::TFilterInterface TFilterInterface;
- typedef _TInputImage TInputImage;
- typedef _TLabelImage TLabelImage;
- typedef typename TInputImage::PixelType TPixel;
- typedef typename TLabelImage::PixelType TLabel;
- typedef typename TLabelImage::RegionType TRegion;
+ typedef fpa::Filters::BaseMarksInterface< TInternalTraits > TMarksInterface;
+ typedef fpa::Filters::Image::SeedsFromLabelsInterface< TInternalTraits > TSeedsInterface;
+ };
+ /**
+ */
+ template< class _TImage, class _TLabels, class _TScalarImage >
+ class MoriLabelling
+ : public fpa::Filters::Image::RegionGrow< _TImage, _TLabels, typename _TLabels::PixelType, CTBronchi::MoriLabellingTraits< _TImage, _TLabels > >
+ {
public:
- itkNewMacro( Self );
- itkTypeMacro( MoriLabelling, itk::ImageToImageFilter );
+ typedef CTBronchi::MoriLabellingTraits< _TImage, _TLabels > TTraits;
+ fpaTraitsMacro( typename TTraits );
- itkGetConstMacro( UpperThreshold, TPixel );
- itkSetMacro( UpperThreshold, TPixel );
+ typedef fpa::Filters::Image::RegionGrow< _TImage, _TLabels, TMark, TTraits > Superclass;
+ typedef MoriLabelling Self;
+ typedef itk::SmartPointer< Self > Pointer;
+ typedef itk::SmartPointer< const Self > ConstPointer;
- itkGetConstMacro( InsideValue, TLabel );
- itkSetMacro( InsideValue, TLabel );
+ typedef fpa::Functors::RegionGrow::BinaryThreshold< TInputValue > TFunctor;
- itkGetConstMacro( InsideLabel, TLabel );
- itkSetMacro( InsideLabel, TLabel );
+ public:
+ itkNewMacro( Self );
+ itkTypeMacro( MoriLabelling, fpa::Filters::Image::RegionGrow );
- itkGetConstMacro( OutsideLabel, TLabel );
- itkSetMacro( OutsideLabel, TLabel );
+ itkGetConstMacro( VesselnessThreshold, double );
+ itkSetMacro( VesselnessThreshold, double );
- public:
- const TLabelImage* GetInputLabelImage( ) const;
- void SetInputLabelImage( TLabelImage* image );
+ itkGetConstMacro( UpperThreshold, TInputValue );
+ itkSetMacro( UpperThreshold, TInputValue );
- const TInputImage* GetInputRawImage( ) const;
- void SetInputRawImage( TInputImage* image );
+ ivqITKInputMacro( InputLabels, _TLabels );
+ ivqITKInputMacro( InputVesselness, _TScalarImage );
protected:
- MoriLabelling( );
- virtual ~MoriLabelling( );
-
- virtual void BeforeThreadedGenerateData( ) override
+ MoriLabelling( )
+ : Superclass( ),
+ m_VesselnessThreshold( 0.05 ),
+ m_UpperThreshold( -650 )
{
- this->Superclass::BeforeThreadedGenerateData( );
-
- const TInputImage* raw = this->GetInputRawImage( );
- typename TInputImage::SpacingType spac = raw->GetSpacing( );
- double sigma = spac[ 0 ];
- for( unsigned int d = 1; d < TInputImage::ImageDimension; ++d )
- sigma = ( spac[ d ] < sigma )? spac[ d ]: sigma;
- sigma *= 1.5;
-
- typedef itk::HessianRecursiveGaussianImageFilter< TInputImage > _THessian;
- typename _THessian::Pointer hessian = _THessian::New( );
- hessian->SetInput( raw );
- hessian->SetSigma( sigma );
-
- typedef itk::Hessian3DToVesselnessMeasureImageFilter< double > _TVesselness;
- typename _TVesselness::Pointer vesselness = _TVesselness::New( );
- vesselness->SetInput( hessian->GetOutput( ) );
- vesselness->SetAlpha1( 0.5 );
- vesselness->SetAlpha2( 2.0 );
-
- typename itk::ImageFileWriter< typename _TVesselness::OutputImageType >::Pointer w =
- itk::ImageFileWriter< typename _TVesselness::OutputImageType >::New( );
- w->SetInput( vesselness->GetOutput( ) );
- w->SetFileName( "vessel.mhd" );
- w->Update( );
+ ivqITKInputConfigureMacro( InputLabels, _TLabels );
+ ivqITKInputConfigureMacro( InputVesselness, _TScalarImage );
+ this->m_Functor = TFunctor::New( );
+ this->SetPredicate( this->m_Functor );
}
- virtual void ThreadedGenerateData( const TRegion& region, itk::ThreadIdType threadId ) override;
+ virtual ~MoriLabelling( )
+ {
+ }
private:
// Purposely not implemented.
Self& operator=( const Self& other );
protected:
- TPixel m_UpperThreshold;
- TLabel m_InsideValue;
- TLabel m_InsideLabel;
- TLabel m_OutsideLabel;
+ typename TFunctor::Pointer m_Functor;
+ double m_VesselnessThreshold;
+ TInputValue m_UpperThreshold;
};
} // ecapseman
-#ifndef ITK_MANUAL_INSTANTIATION
-# include <CTBronchi/MoriLabelling.hxx>
-#endif // ITK_MANUAL_INSTANTIATION
-
-#endif // __CTBronchi__MoriLabelling__h__
+#endif // __CTBronchi__Functions__h__
// eof - $RCSfile$