From 54de2903ffa1d7f444d2a38609de58f8f2ab0d2d Mon Sep 17 00:00:00 2001 From: David Sarrut Date: Fri, 4 Nov 2011 10:22:45 +0100 Subject: [PATCH] First version of the builder of relpos database --- ...itkRelativePositionDataBaseBuilderFilter.h | 153 ++++++++++++++++++ 1 file changed, 153 insertions(+) create mode 100644 itk/clitkRelativePositionDataBaseBuilderFilter.h diff --git a/itk/clitkRelativePositionDataBaseBuilderFilter.h b/itk/clitkRelativePositionDataBaseBuilderFilter.h new file mode 100644 index 0000000..bad3d83 --- /dev/null +++ b/itk/clitkRelativePositionDataBaseBuilderFilter.h @@ -0,0 +1,153 @@ +/*========================================================================= + Program: vv http://www.creatis.insa-lyon.fr/rio/vv + + Authors belong to: + - University of LYON http://www.universite-lyon.fr/ + - Léon Bérard cancer center http://www.centreleonberard.fr + - CREATIS CNRS laboratory http://www.creatis.insa-lyon.fr + + This software is distributed WITHOUT ANY WARRANTY; without even + the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR + PURPOSE. See the copyright notices for more information. + + It is distributed under dual licence + + - BSD See included LICENSE.txt file + - CeCILL-B http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html + ===========================================================================**/ + +#ifndef CLITKRelativePositionDataBaseBuilderFILTER_H +#define CLITKRelativePositionDataBaseBuilderFILTER_H + +// clitk +#include "clitkFilterWithAnatomicalFeatureDatabaseManagement.h" +#include "clitkFilterBase.h" +#include "clitkRelativePositionAnalyzerFilter.h" +#include "clitkRelativePositionDataBase.h" + +// itk +#include + +namespace clitk { + + //-------------------------------------------------------------------- + /* + Analyze the relative position of a Target (mask image) according + to some Object, in a given Support. Indicate the main important + position of this Target according the Object. + */ + //-------------------------------------------------------------------- + + template + class RelativePositionDataBaseBuilderFilter: + public virtual FilterBase, + public clitk::FilterWithAnatomicalFeatureDatabaseManagement, + public itk::ImageToImageFilter + { + + public: + /** Standard class typedefs. */ + typedef itk::ImageToImageFilter Superclass; + typedef RelativePositionDataBaseBuilderFilter Self; + typedef itk::SmartPointer Pointer; + typedef itk::SmartPointer ConstPointer; + + /** Method for creation through the object factory. */ + itkNewMacro(Self); + + /** Run-time type information (and related methods). */ + itkTypeMacro(RelativePositionDataBaseBuilderFilter, ImageToImageFilter); + + /** Some convenient typedefs. */ + typedef typename ImageType::ConstPointer ImageConstPointer; + typedef typename ImageType::Pointer ImagePointer; + typedef typename ImageType::RegionType RegionType; + typedef typename ImageType::PixelType PixelType; + typedef typename ImageType::SpacingType SpacingType; + typedef typename ImageType::SizeType SizeType; + typedef typename ImageType::IndexType IndexType; + typedef typename ImageType::PointType PointType; + + /** ImageDimension constants */ + itkStaticConstMacro(ImageDimension, unsigned int, ImageType::ImageDimension); + FILTERBASE_INIT; + typedef itk::Image FloatImageType; + + // Inputs + itkGetConstMacro(SupportName, std::string); + itkSetMacro(SupportName, std::string); + + itkGetConstMacro(TargetName, std::string); + itkSetMacro(TargetName, std::string); + + void AddObjectName(std::string s) { m_ObjectNames.push_back(s); } + std::string & GetObjectName(int i) { return m_ObjectNames[i]; } + int GetNumberOfObjects() { return m_ObjectNames.size(); } + + // Options + itkGetConstMacro(BackgroundValue, PixelType); + itkSetMacro(BackgroundValue, PixelType); + + itkGetConstMacro(ForegroundValue, PixelType); + itkSetMacro(ForegroundValue, PixelType); + + itkGetConstMacro(NumberOfBins, int); + itkSetMacro(NumberOfBins, int); + + itkGetConstMacro(NumberOfAngles, int); + itkSetMacro(NumberOfAngles, int); + + itkGetConstMacro(AreaLossTolerance, double); + itkSetMacro(AreaLossTolerance, double); + + // For debug + void PrintOptions(); + + // I dont want to verify inputs information + virtual void VerifyInputInformation() { } + + protected: + RelativePositionDataBaseBuilderFilter(); + virtual ~RelativePositionDataBaseBuilderFilter() {} + + PixelType m_BackgroundValue; + PixelType m_ForegroundValue; + + ImagePointer m_Support; + ImagePointer m_Object; + ImagePointer m_Target; + + std::string m_SupportName; + std::string m_TargetName; + std::vector m_ObjectNames; + + int m_NumberOfAngles; + std::vector m_ListOfAngles; + std::vector m_ListOfDirections; + + int m_NumberOfBins; + double m_AreaLossTolerance; + + virtual void GenerateOutputInformation(); + virtual void GenerateInputRequestedRegion(); + virtual void GenerateData(); + + typename FloatImageType::Pointer + ComputeFuzzyMap(ImageType * object, ImageType * target, double angle); + + void + ComputeOptimalThresholds(FloatImageType * map, ImageType * target, int bins, double tolerance, + double & threshold, double & reverseThreshold); + private: + RelativePositionDataBaseBuilderFilter(const Self&); //purposely not implemented + void operator=(const Self&); //purposely not implemented + + }; // end class + //-------------------------------------------------------------------- + +} // end namespace clitk +//-------------------------------------------------------------------- + +#include "clitkRelativePositionDataBaseBuilderFilter.txx" + +#endif -- 2.47.1